Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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

Phase behavior of the 38-atom Lennard-Jones cluster

We have developed a coarse-grained description of the phase behavior of the isolated 38-atom Lennard-Jones cluster (LJ{sub 38}). The model captures both the solid-solid polymorphic transitions at low temperatures and the complex cluster breakup and melting transitions at higher temperatures. For this coarse model development, we employ the manifold learning technique of diffusion mapping. The outcome of the diffusion mapping analysis over a broad temperature range indicates that two order parameters are sufficient to describe the cluster's phase behavior; we have chosen two such appropriate order parameters that are metrics of condensation and overall crystallinity. In this well-justified coarse-variable space, we calculate the cluster's free energy landscape (FEL) as a function of temperature, employing Monte Carlo umbrella sampling. These FELs are used to quantify the phase behavior and onsets of phase transitions of the LJ{sub 38} cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A second-order accurate finite volume scheme with the discrete maximum principle for solving Richards’ equation on unstructured meshes

Richards’s equation describes steady-state or transient flow in a variably saturated medium. For a medium having multiple layers of soils that are not aligned with coordinate axes, a mesh fitted to these layers is no longer orthogonal and the classical two-point flux approximation finite volume scheme is no longer accurate. Here, we propose new second-order accurate nonlinear finite volume (NFV) schemes for the head and pressure formulations of Richards’ equation. We prove that the discrete maximum principles hold for both formulations at steady-state which mimics similar properties of the continuum solution. The second-order accuracy is achieved using high-order upwind algorithms for the relative permeability. Numerical simulations of water infiltration into a dry soil show significant advantage of the second-order NFV schemes over the first-order NFV schemes even on coarse meshes. Since explicit calculation of the Jacobian matrix becomes prohibitively expensive for high-order schemes due to build-in reconstruction and slope limiting algorithms, we study numerically the preconditioning strategy introduced recently in Lipnikov et al. (2016) that uses a stable approximation of the continuum Jacobian. Lastly, numerical simulations show that the new preconditioner reduces computational cost up to 2–3 times in comparison with the conventional preconditioners.

97 MATHEMATICS AND COMPUTING↗

Rapid and multiband variability of the TeV bright active nucleus of the galaxy IC 310

Recently the radio galaxy IC 310 was identified as a γ-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E > 100 GeV) with the MAGIC telescopes. Originally classified as a head-tail radio galaxy, the nature of this object is subject of controversy since its nucleus shows blazar-like behavior. In order to understand the nature of IC 310 and the origin of the VHE emission, we studied the spectral and flux variability of IC 310 from the X-ray band to the VHE γ-ray regime. The light curve of IC 310 above 300 GeV has been measured with the MAGIC telescopes from 2009 October to 2010 February. Contemporaneous Fermi-LAT data (2008-2011) in the 10-500 GeV energy range were also analyzed. In the X-ray regime, archival observations from 2003 to 2007 with XMM-Newton, Chandra, and Swift-XRT in the 0.5-10 keV band were studied. The VHE light curve reveals several high-amplitude and short-duration flares. Day-to-day flux variability is clearly present (>5σ). The photon index between 120 GeV and 8 TeV remains at the value Γ ~ 2.0 during both low and high flux states. The VHE spectral shape does not show significant variability, whereas the flux at 1 TeV changes by a factor of ~7. Fermi-LAT detected only eight γ-ray events in the energy range 10 GeV–500 GeV in three years of observation. Moreover, the measured photon index of Γ = 1.3 ± 0.5 in the Fermi-LAT range is very hard. The X-ray measurements show strong variability in both flux and photon index. The latter varied from 1.76 ± 0.07 to 2.55 ± 0.07. The rapid variability measured in γ-rays and X-rays confirms the blazar-like behavior of IC 310. The multi-TeV γ-ray emission seems to originate from scales of less than 80 Schwarzschild radii (for a black hole mass of 2 × 10 8 M⊙) within the compact core of its FR I radio jet with orientation angle 10°-38°. The spectral energy distribution resembles that of an extreme blazar, albeit the luminosity is more than two orders of magnitude lower.

79 ASTRONOMY AND ASTROPHYSICS↗

Edge Pushing is Equivalent to Vertex Elimination for Computing Hessians

We prove the equivalence of two different Hessian evaluation algorithms in AD. The first is the Edge Pushing algorithm of Gower and Mello, which may be viewed as a second order Reverse mode algorithm for computing the Hessian. In earlier work, we have derived the Edge Pushing algorithm by exploiting a Reverse mode invariant based on the concept of live variables in compiler theory. The second algorithm is based on eliminating vertices in a computational graph of the gradient, in which intermediate variables are successively eliminated from the graph, and the weights of the edges are updated suitably. We prove that if the vertices are eliminated in a reverse topological order while preserving symmetry in the computational graph of the gradient, then the Vertex Elimination algorithm and the Edge Pushing algorithm perform identical computations. In this sense, the two algorithms are equivalent. This insight that unifies two seemingly disparate approaches to Hessian computations could lead to improved algorithms and implementations for computing Hessians. Read More: http://epubs.siam.org/doi/10.1137/1.9781611974690.ch11

Wang, Mu↗

Investigation of Dynamics in BMIM TFSA Ionic Liquid through Variable Temperature and Pressure NMR Relaxometry and Diffusometry

A comprehensive variable temperature, pressure and frequency multinuclear ( 1 H, 2 H, and 19 F) magnetic resonance study was undertaken on selectively deuterated 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide (BMIM TFSA) ionic liquid isotopologues. This study builds on our earlier investigation of the effects of increasing alkyl chain length on diffusion and dynamics in imidazolium-based TFSA ionic liquids. Fast field cycling 1 H T 1 data revealed multiple modes of motion. Through calculation of diffusion coefficient (D) values and activation energies, the low- and high-field regimes were assigned to the translational and reorientation dynamics respectively. Variable-pressure 2 H T 1 measurements reveal site-dependent interactions in the cation with strengths in the order MD 3 > CD 3 > CD 2 , indicating dissimilarities in the electric field gradients along the alkyl chain, with the CD 2 sites having the largest gradient. Additionally, the α saturation effect in T 1 vs. P was observed for all three sites, suggesting significant reduction of the short-range rapid reorientational dynamics. This reduction was also deduced from the variable pressure 1 H T 1 data, which showed an approach to saturation for both the methyl and butyl group terminal methyl sites. Pressure-dependent D measurements show independent motions for both cations and anions, with the cations having greater D values over the entire pressure range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-driven wind turbine wake modeling via probabilistic machine learning

Wind farm design primarily depends on the variability of the wind turbine wake flows to the atmospheric wind conditions and the interaction between wakes. Physics-based models that capture the wake flow field with high-fidelity are computationally very expensive to perform layout optimization of wind farms, and, thus, data-driven reduced-order models can represent an efficient alternative for simulating wind farms. In this work, we use real-world light detection and ranging (LiDAR) measurements of wind-turbine wakes to construct predictive surrogate models using machine learning. Specifically, we first demonstrate the use of deep autoencoders to find a low-dimensional latent space that gives a computationally tractable approximation of the wake LiDAR measurements. Then, we learn the mapping between the parameter space and the (latent space) wake flow fields using a deep neural network. Additionally, we also demonstrate the use of a probabilistic machine learning technique, namely, Gaussian process modeling, to learn the parameter-space-latent-space mapping in addition to the epistemic and aleatoric uncertainty in the data. Finally, to cope with training large datasets, we demonstrate the use of variational Gaussian process models that provide a tractable alternative to the conventional Gaussian process models for large datasets. Furthermore, we introduce the use of active learning to adaptively build and improve a conventional Gaussian process model predictive capability. Overall, we find that our approach provides accurate approximations of the wind-turbine wake flow field that can be queried at an orders-of-magnitude cheaper cost than those generated with high-fidelity physics-based simulations.

Deep neural networks↗

Adjusting inhomogeneous daily temperature variability using wavelet analysis

A wavelet-analysis-based homogenization method (WHM) is developed for detecting and adjusting biases of variability in a daily climate observation series, which are influential in the estimation of climate extremes and relevant trends in the time series. WHM is applied to the Central England daily mean temperature series (CET, 1772–2010) and the daily mean temperature series at Henan station in mountainous western China (HNT, 1960–2008) in order to demonstrate the usefulness of the new technique. The changes of methods for calculating daily mean temperature for CET in 1878 caused significantly reduced daily variability (DV) in the subsequent sub-series of CET. The relocation of the Chinese station in 1981 to a higher mountainous site resulted in not only enhanced DV but also in weather (weekly) time-scale variability (WV) in the subsequent sub-series of HNT. The adjustments based on WHM in the present cases improve the estimates of the long-term trends of climate extremes such as hot days, cold days, heat waves, and cold surges. A few widely applied methods of homogenization such as multiple analysis of series for homogenization (MASH), higher-order moment (HOM), RHtestsV3 with quantile-matching (QM) adjustments and two-phase regression (TPR) are also applied to HNT for comparison. TPR is not aimed at improving the homogeneity of variability in the time series; MASH does not improve daily variability either; HOM improves it in a small way; and with RHtestsV3 with QM adjustments it is considerably improved, but biases remain. In contrast, WHM improves the homogeneity of the short-term variability in the time series, resulting in reasonable assessments of climate extremes and trends in the daily observations.

54 ENVIRONMENTAL SCIENCES↗

Structural, magnetic and electrical properties of the hexagonal ferrites MFeO{sub 3} (M=Y, Yb, In)

The hexagonal ferrites MFeO{sub 3} (M=Y, Yb, In) have been studied using a combination of neutron and X-ray powder diffraction, magnetic susceptibility, dielectric measurements and {sup 57}Fe Moessbauer spectroscopy. This study confirms the previously reported crystal structure of InFeO{sub 3} (YAlO{sub 3} structure type, space group P6{sub 3}/mmc), but YFeO{sub 3} and YbFeO{sub 3} both show a lowering of symmetry to at most P6{sub 3}cm (ferrielectric YMnO{sub 3} structure type). However, Moessbauer spectroscopy shows at least two distinct Fe sites for both YFeO{sub 3} and YbFeO{sub 3} and we suggest that the best model to rationalise this involves phase separation into more than one similar hexagonal YMnO{sub 3}-like phase. Rietveld analysis of the neutron diffraction data was carried out using two hexagonal phases as a simplest case scenario. In both YFeO{sub 3} and YbFeO{sub 3}, distinct dielectric anomalies are observed near 130 K and 150 K, respectively. These are tentatively correlated with weak anomalies in magnetic susceptibility and lattice parameters, for YFeO{sub 3} and YbFeO{sub 3}, respectively, which may suggest a weak magnetoelectric effect. Comparison of neutron and X-ray powder diffraction shows evidence of long-range magnetic order in both YFeO{sub 3} and YbFeO{sub 3} at low temperatures. Due to poor sample crystallinity, the compositional and structural effects underlying the phase separation and possible magnetoelectric phenomena cannot be ascertained. - Graphical abstract: Hexagonal MFeO{sub 3} (M=Y, Yb) exhibit phase separation into two YMnO{sub 3}-like phases. Variable temperature crystallographic, electrical and magnetic studies suggest weak correlations between electrical and magnetic responses and long-range magnetic order at low temperature. Highlights: Black-Right-Pointing-Pointer Multi-technique study of multiferroic hexagonal MFeO{sub 3}. Black-Right-Pointing-Pointer Phase separation into two similar hexagonal phases. Black-Right-Pointing-Pointer Weak coupling of electrical and magnetic responses. Black-Right-Pointing-Pointer Long-range magnetic order at low T.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sr{sub 3}Bi{sub 2}(SeO{sub 3}){sub 6}·H{sub 2}O: A novel anionic layer consisting of second-order Jahn–Teller (SOJT) distortive cations

A new layered bismuth selenite hydrate material, Sr{sub 3}Bi{sub 2}(SeO{sub 3}){sub 6}·H{sub 2}O has been synthesized through a hydrothermal reaction using SrCO{sub 3}, Bi{sub 2}O{sub 3}, SeO{sub 2}, and water as reagents. The crystal structure of the reported material has been determined by single crystal X-ray diffraction analysis. The anionic layered framework of Sr{sub 3}Bi{sub 2}(SeO{sub 3}){sub 6}·H{sub 2}O consists of polyhedra of second-order Jahn–Teller (SOJT) distortive cations, Bi{sup 3+} and Se{sup 4+}. Attributable to the variable and asymmetric coordination geometry of the constituent cations, a rich structural chemistry including uni-dimensional bands and linkers is observed in the layer. The material is thermally stable up to about 390 °C and starts decomposing due to the sublimation of SeO{sub 2} above the temperature. The UV–vis diffuse reflectance spectrum suggests a band gap of 3.3 eV. Elemental analysis, infrared spectrum, local dipole moment calculations, and electronic structure calculations are also reported. - Graphical abstract: Ball-and-stick and polyhedral representation of Sr{sub 3}Bi{sub 2}(SeO{sub 3}){sub 6}·H{sub 2}O in the ab-plane. The stereoactive lone pairs are visualized around the Bi{sup 3+} and Se{sup 4+} as lobelike isosurfaces through ELF calculations with η=0.9. - Highlights: • A novel layered bismuth selenite, Sr{sub 3}Bi{sub 2}(SeO{sub 3}){sub 6}·H{sub 2}O has been hydrothermally synthesized. • The layer is only composed of second-order Jahn–Teller distortive cations in asymmetric environments. • The stereoactive lone pairs are visualized as lobelike isosurfaces through ELF calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The definitive X-ray light curve of Swift J164449.3+573451

On 2011 March 28, the Swift Burst Alert Telescope triggered on an object that had no analog in over six years of Swift operations. Follow-up observations by the Swift X-ray Telescope (XRT) found a new, bright X-ray source covering three orders of magnitude in flux over the first few days, that was much more persistent (and variable) than gamma-ray burst afterglows. Ground-based spectroscopy found a redshift of 0.35, implying extremely high luminosity, with integrated isotropic-equivalent energy output in the X-ray band alone exceeding 10{sup 53} erg in the first two weeks after discovery. Strong evidence for a collimated outflow or beamed emission was found. The observational properties of this object are unlike anything ever before observed. We interpret these unique properties as the result of emission from a relativistic jet produced in the aftermath of the tidal disruption of a main sequence star by a massive black hole (BH) in the center of the host galaxy. The source decayed slowly as the stellar remnants were accreted onto the BH, before abruptly shutting off. Here we present the definitive XRT team light curve for Swift J164449.3+573451 and discuss its implications. We show that the unabsorbed flux decayed roughly as a t{sup −1.5} power law up to 2012 August 17. The steep turnoff of an order of magnitude in 24 hr seems to be consistent with the shutdown of the jet as the accretion disk transitioned from a thick disk to a thin disk. Swift continues to monitor this source in case the jet reactivates.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection of significant cm to sub-mm band radio and γ-ray correlated variability in Fermi bright blazars

The exact location of the γ-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA programme) and γ-ray (0.1–300 GeV) ~3.5 yr light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. These results reveal: (i) the first highly significant multiband radio and γ-ray correlations (radio lagging γ rays) when averaging over the whole sample, (ii) average time delays (source frame: 76 ± 23 to 7 ± 9 d), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence τr, γ(ν) ∝ ν -1 , in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk γ-ray production region typically located within/upstream of the 3 mm core region (τ3mm, γ = 12 ± 8 d), (iv) mean distances between the region of γ-ray peak emission and the radio ‘τ = 1 photosphere’ decreasing from 9.8 ± 3.0 pc (11 cm) to 0.9 ± 1.1 pc (2 mm) and 1.4 ± 0.8 pc (0.8 mm), (v) 3 mm/γ-ray correlations in nine individual sources at a significance level where one is expected by chance (probability: 4 × 10 -6 ), (vi) opacity and ‘time lag core shift’ estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk γ-ray production region from the supermassive black hole (SMBH) of ~0.8–1.6 pc, i.e. at the outer edge of the broad-line region (BLR) or beyond. A 3 mm τ = 1 surface at ~2–3 pc from the jet base (i.e. well outside the ‘canonical BLR’) finally suggests that BLR material extends to several parsec distances from the SMBH.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectroscopic classification, variability, and SED of the Fermi -detected CSS 3C 286: the radio-loudest NLS1 galaxy?

ABSTRACT 3C 286 is a well-known calibrator source in radio astronomy. It is also one of very few compact steep-spectrum (CSS) sources detected in γ-rays. Here, we perform a detailed spectroscopic and variability analysis and present the first quasi-simultaneous optical to X-ray spectral energy distribution (SED) in order to reveal physical mechanisms that dominate its emission at different wavelengths, and arrive at a reliable optical source classification. The first main result of our study reveals several pitfalls when applying simple broad- or narrow-line Seyfert 1 (BLS1 or NLS1) classification criteria that only look at the [O iii]–Hβ complex. [O iii] and Hβ can be dominated by the same outflow components, in which case full width at half-maximum (Hβ) is no reliable classification criterion, and extinction by intrinsic or intervening material can make the highest velocity Hβ component undetectable. After careful combination of all information from ultraviolet (UV)–optical spectra along with multiwavelength data, we confirm that 3C 286 can be classified as NLS1 galaxy, with line properties and supermassive black hole mass (of order 108 M⊙ and accreting near the Eddington limit) close to the BLS1 regime, making it an important borderline object. The quasi-simultaneous SED taken with Swift shows a sharp rise in the optical–UV, implying the presence of a strong accretion-disc component with extreme UV excess, consistent with emission-line diagnostics. Finally, we report the discovery of X-ray variability of 3C 286, plausibly dominated by jet emission, and variable by at least a factor of ∼4. This result suggests to exercise caution when using 3C 286 as radio calibrator in high-resolution radio very long baseline interferometry observations.

Yao, Su↗

Simulating HFIR Core Thermal Hydraulics Using 3D-2D Model Coupling

A model utilizing interdimensional variable coupling is presented for simulating the thermal hydraulic interactions of the High Flux Isotope Reactor (HFIR) core at Oak Ridge National Laboratory (ORNL). The model s domain consists of a single, explicitly represented three-dimensional fuel plate and a simplified two-dimensional coolant channel slice. In simplifying the coolant channel, and thus the number of mesh points in which the Navier-Stokes equations must be solved, the computational cost and solution time are both greatly reduced. In order for the reduced-dimension coolant channel to interact with the explicitly represented fuel plate, however, interdimensional variable coupling must be enacted along all shared boundaries. The primary focus of this paper is in detailing the collection, storage, passage, and application of variables across this interdimensional interface. Comparisons are made showing the general speed-up associated with this simplified coupled model.

HFIR↗

Surrogate models for mixed discrete-continuous variables

Large-scale computational models have become common tools for analyzing complex man made systems. However, when coupled with optimization or uncertainty quantification methods in order to conduct extensive model exploration and analysis, the computational expense quickly becomes intractable. Furthermore, these models may have both continuous and discrete parame ters. One common approach to mitigating the computational expense is the use of response surface approximations. While well developed for models with continuous parameters, they are still new and largely untested for models with both continuous and discrete parameters. In this work, we describe and investigate the performance of three types of response surfaces developed for mixed variable models: Adaptive Component Selection and Shrinkage Operator, Treed Gaussian Process, and Gaussian Process with Special Correlation Functions. We focus our efforts on test problems with a small number of parameters of interest, a characteristic of many physics-based engineering models. We present the results of our studies and offer some insights regarding the performance of each response surface approximation method.

Swiler, Laura Painton↗

Coherent states, 6j symbols and properties of the next to leading order asymptotic expansions

We present the first complete derivation of the well-known asymptotic expansion of the SU(2) 6j symbol using a coherent state approach, in particular we succeed in computing the determinant of the Hessian matrix. To do so, we smear the coherent states and perform a partial stationary point analysis with respect to the smearing parameters. This allows us to transform the variables from group elements to dihedral angles of a tetrahedron resulting in an effective action, which coincides with the action of first order Regge calculus associated to a tetrahedron. To perform the remaining stationary point analysis, we compute its Hessian matrix and obtain the correct measure factor. Furthermore, we expand the discussion of the asymptotic formula to next to leading order terms, prove some of their properties and derive a recursion relation for the full 6j symbol.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

97 MATHEMATICS AND COMPUTING↗

Placement and Sizing of Inverter-Based Renewable Systems in Multi-Phase Distribution Networks

This study develops a tractable formulation for optimal placement and sizing of inverter-based renewable systems in multi-phase distribution networks. The goal of the formulation is to minimize the cost of inverter installation, average power import, and average distributed generation curtailment. Threephase and single-phase inverter models are presented that preserve the underlying mappings between renewable uncertainty to power injection. The uncertainty of distributed generators (DGs) and loads are characterized by a finite set of scenarios. Linear multi-phase power flow approximations are used in conjunction with scenario reduction techniques to arrive at a tractable twostage stochastic formulation for optimal DG placement and sizing. First-stage decisions are locations for DG deployment and capacity sizes, and second-stage decisions include DG real power curtailment, reactive power support, as well as feeder voltage profile. The resulting formulation is a mixed-integer second-order cone program and can be solved efficiently either by existing optimization solvers or by relaxing the binary variables to the [0,1] interval. Simulation studies on standard multi-phase IEEE test feeders promise that optimal stochastic planning of DGs reduces costs during validation, compared to a scheme where uncertainty is only represented by its average value.

24 POWER TRANSMISSION AND DISTRIBUTION↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗