In topology optimization of structures, materials and mechanisms, parametrization of geometry is often performed by a grey-scale density-like interpolation function. In this paper we analyze and compare the various approaches to this concept in the light of variational bounds on effective properties of composite materials. This allows us to derive simple necessary conditions for the possible realization of grey-scale via composites, leading to a physical interpretation of all feasible designs as well as the optimal design. Thus it is shown that the so-called artificial interpolation model in many circumstances actually falls within the framework of microstructurally based models. Single material and multi-material structural design in elasticity as well as in multi-physics problems is discussed.

/pdf/material-interpolation-schemes-in-topology-optimization-1ps0tr57ei.pdf

Material interpolation schemes in topology optimization

A dimension reduction method called discrete empirical interpolation is proposed and shown to dramatically reduce the computational complexity of the popular proper orthogonal decomposition (POD) method for constructing reduced-order models for time dependent and/or parametrized nonlinear partial differential equations (PDEs). In the presence of a general nonlinearity, the standard POD-Galerkin technique reduces dimension in the sense that far fewer variables are present, but the complexity of evaluating the nonlinear term remains that of the original problem. The original empirical interpolation method (EIM) is a modification of POD that reduces the complexity of evaluating the nonlinear term of the reduced model to a cost proportional to the number of reduced variables obtained by POD. We propose a discrete empirical interpolation method (DEIM), a variant that is suitable for reducing the dimension of systems of ordinary differential equations (ODEs) of a certain type. As presented here, it is applicable to ODEs arising from finite difference discretization of time dependent PDEs and/or parametrically dependent steady state problems. However, the approach extends to arbitrary systems of nonlinear ODEs with minor modification. Our contribution is a greatly simplified description of the EIM in a finite-dimensional setting that possesses an error bound on the quality of approximation. An application of DEIM to a finite difference discretization of the one-dimensional FitzHugh-Nagumo equations is shown to reduce the dimension from 1024 to order 5 variables with negligible error over a long-time integration that fully captures nonlinear limit cycle behavior. We also demonstrate applicability in higher spatial dimensions with similar state space dimension reduction and accuracy results.

/pdf/nonlinear-model-reduction-via-discrete-empirical-f3tjqvqhuf.pdf

Nonlinear Model Reduction via Discrete Empirical Interpolation

Recent advancements in computational inverse-design approaches — algorithmic techniques for discovering optical structures based on desired functional characteristics — have begun to reshape the landscape of structures available to nanophotonics. Here, we outline a cross-section of key developments in this emerging field of photonic optimization: moving from a recap of foundational results to motivation of applications in nonlinear, topological, near-field and on-chip optics. Starting with a desired optical output it is possible to use computational algorithms to inverse design devices. The approach is reviewed here with an emphasis on nanophotonics.

Inverse design in nanophotonics

Phononic band-gap materials prevent elastic waves in certain frequency ranges from propagating, and they may therefore be used to generate frequency filters, as beam splitters, as sound or vibration protection devices, or as waveguides. In this work we show how topology optimization can be used to design and optimize periodic materials and structures exhibiting phononic band gaps. Firstly, we optimize infinitely periodic band-gap materials by maximizing the relative size of the band gaps. Then, finite structures subjected to periodic loading are optimized in order to either minimize the structural response along boundaries (wave damping) or maximize the response at certain boundary locations (waveguiding).

Systematic design of phononic band-gap materials and structures by topology optimization.

Data-driven discovery is revolutionizing the modeling, prediction, and control of complex systems. This textbook brings together machine learning, engineering mathematics, and mathematical physics to integrate modeling and control of dynamical systems with modern methods in data science. It highlights many of the recent advances in scientific computing that enable data-driven methods to be applied to a diverse range of complex systems, such as turbulence, the brain, climate, epidemiology, finance, robotics, and autonomy. Aimed at advanced undergraduate and beginning graduate students in the engineering and physical sciences, the text presents a range of topics and methods from introductory to state of the art.

Data-Driven Science and Engineering: Machine Learning, Dynamical Systems, and Control

The shape of the ideal column

Extremal eigenvalue problems for two-phase conductors

Dendrites form the major components of neurons. They are complex branching structures that receive and process thousands of synaptic inputs from other neurons. It is well known that dendritic morphology plays an important role in the function of dendrites. Another important contribution to the response characteristics of a single neuron comes from the intrinsic resonant properties of dendritic membrane. In this paper we combine the effects of dendritic branching and resonant membrane dynamics by generalising the “sum-over-trips” approach (Abbott et al. in Biol Cybernetics 66, 49–60 1991). To illustrate how this formalism can shed light on the role of architecture and resonances in determining neuronal output we consider dual recording and reconstruction data from a rat CA1 hippocampal pyramidal cell. Specifically we explore the way in which an I h current contributes to a voltage overshoot at the soma.

/pdf/branching-dendrites-with-resonant-membrane-a-sum-over-trips-oj0n4khq84.pdf

Branching dendrites with resonant membrane: a “sum-over-trips” approach

The accurate simulation of a neuron’s ability to integrate distributed synaptic input typically requires the simultaneous solution of tens of thousands of ordinary differential equations. For, in order to understand how a cell distinguishes between input patterns we apparently need a model that is biophysically accurate down to the space scale of a single spine, i.e., 1 μm. We argue here that one can retain this highly detailed input structure while dramatically reducing the overall system dimension if one is content to accurately reproduce the associated membrane potential at a small number of places, e.g., at the site of action potential initiation, under subthreshold stimulation. The latter hypothesis permits us to approximate the active cell model with an associated quasi-active model, which in turn we reduce by both time-domain (Balanced Truncation) and frequency-domain (
 ${\cal H}_2$
 approximation of the transfer function) methods. We apply and contrast these methods on a suite of typical cells, achieving up to four orders of magnitude in dimension reduction and an associated speed-up in the simulation of dendritic democratization and resonance. We also append a threshold mechanism and indicate that this reduction has the potential to deliver an accurate quasi-integrate and fire model.

/pdf/low-dimensional-morphologically-accurate-models-of-39cecp0rp9.pdf

Steven J. Cox

Papers

The shape of the ideal column

Extremal eigenvalue problems for two-phase conductors

Branching dendrites with resonant membrane: a “sum-over-trips” approach

Low-dimensional, morphologically accurate models of subthreshold membrane potential.

The Generalized Gradient at a Multiple Eigenvalue