Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse prediction”

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 1,171 records · Page 65

Development of coated heat flux gauges for fast responding measurements

Thermal systems often exhibit transient behaviors that have important implications for the operation of the system and can be difficult to predict. For these reasons, experimental testing is often required to ensure system durability requirements are achieved. One important parameter governing the survivability of components in hot, high-stress environments is the heat flux into the part that dictates the temperature distribution for the component. However, sensors required to experimentally characterize heat fluxes in extreme environments must also be resilient. Here, this study presents the development of coated heat transfer gauges capable of robust, high-frequency measurements in turbine research facilities. The addition of a protective coating increases the durability of the gauge, but inherent of that coating is the attenuation of high-frequency temperature penetrations. As a result, this study first outlines the use of analytical solutions to define a gauge design for a specific frequency range and heat transfer, ensuring that subsurface signals can be rectified to surface conditions through inverse methods. Then, the fabrication of polyimide substrate sensors with a parylene-F coating is described. Micro surface heaters added to the custom sensors were used to determine important geometric and thermal properties necessary to calculate accurate surface heat flux. Ultimately, this work shows increased sensor robustness in a turbine test bed and experimentally validates that the frequency response of the fabricated sensors meet the design intent.

47 OTHER INSTRUMENTATION↗

Maximum respiration rates in hyporheic zone sediments are primarily constrained by organic carbon concentration and secondarily by organic matter chemistry

Abstract. River corridors are fundamental components of the Earth system, and their biogeochemistry can be heavily influenced by processes in subsurface zones immediately below the riverbed, referred to as the hyporheic zone. Within the hyporheic zone, organic matter (OM) fuels microbial respiration, and OM chemistry heavily influences aerobic and anaerobic biogeochemical processes. The link between OM chemistry and respiration has been hypothesized to be mediated by OM molecular diversity, whereby respiration is predicted to decrease with increasing diversity. Here we test the specific prediction that aerobic respiration rates will decrease with increases in the number of unique organic molecules (i.e., OM molecular richness, as a measure of diversity). We use publicly available data across the United States from crowdsourced samples taken by the Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems (WHONDRS) consortium. Our continental-scale analyses rejected the hypothesis of a direct limitation of respiration by OM molecular richness. In turn, we found that organic carbon (OC) concentration imposes a primary constraint over hyporheic zone respiration, with additional potential influences of OM richness. We specifically observed respiration rates to decrease nonlinearly with the ratio of OM richness to OC concentration. This relationship took the form of a constraint space with respiration rates in most systems falling below the constraint boundary. A similar, but slightly weaker, constraint boundary was observed when relating respiration rate to the inverse of OC concentration. These results indicate that maximum respiration rates may be governed primarily by OC concentration, with secondary influences from OM richness. Our results also show that other variables often suppress respiration rates below the maximum associated with the richness-to-concentration ratio. An important focus of future research will identify physical (e.g., sediment grain size), chemical (e.g., nutrient concentrations), and/or biological (e.g., microbial biomass) factors that suppress hyporheic zone respiration below the constraint boundaries observed here.

58 GEOSCIENCES↗

The Orbital Variation of Non-Thermal X-Ray Emission from eta Carinae

The collision of strong stellar winds in massive binary systems creates powerful shocks,which accelerates a small amount of particles to relativistic energies at the shock interface.This process is important in fundamental astrophysics in two ways;i) some particles may contribute to cosmic-rays observed around the Earth, whose origin is not well known, andii) it provides a good laboratory for particle acceleration physics as the shock occurs steadily in a predictable environment.Some relativistic particles collide with stellar photons or ambient material near the binary systemand emit non-thermal X-rays and gamma-rays, which are good probes of particle acceleration.This emission had been searched for decades but not found convincingly without sensitive high-energy telescopes.The NuSTAR telescope observed the enigmatic supermassive binary system eta Carinae multiple times after 2014 and found conclusive evidence of non-thermal emission from the star in the extremely hard X-ray band.This emission is prominent between 20-50 keV, below which thermal emission from shock colliding plasma dominates.It is relatively stable throughout the binary orbit, but it disappears near periastron when the wind colliding activity shuts off.This variation indicates that the non-thermal X-ray emission originates from the head-on wind-wind collision.The flat spectrum is consistent with inverse-Compton of stellar UV photons by accelerated electrons.The spectrum smoothly connects to a gamma-ray spectrum of a Fermi source detected around eta Carinae,suggesting that acceleration occurs up to the GeV energy.This result provide the strongest evidence so far that particle acceleration occurs at the wind colliding shock of a massive colliding wind binary system.

Hamaguchi, Kenji↗

Simultaneous Planck, Swift, and Fermi Observations of X-Ray and gamma-Ray Selected Blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and -ray bands, and we compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi Large Area Telescope (LAT), whereas 30 to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-millimetre spectral slope of blazars is quite flat, with [alpha] approximately 0 up to about 70 GHz, above which it steepens to [alpha] approximately -0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (v(sup IC) (sub (PEAK)), ranges from 10(sup 21) to 10(sup 22) HZ. The distribution of the rest-frame synchrotron peak frequency (v(sup s)(sub peak)) in the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples with (v(sup s)(sub peak) = 10(sup 13:1 plus or minus 0.1) Hz, while the mean inverse-Compton peak frequency,(v(sup IC)(sub peak) ranges from 10(sup 21) to 10(sup 22) Hz. The distributions of v(sup S)(sub peak) and of v(sup IC)(sub peak) of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at approximately 7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection effect. A comparison of our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and v(sup S)(sub peak) predicted by the blazar sequence.

Giommi, P.↗

Probabilistic Prognosis of Non-Planar Fatigue Crack Growth

Quantifying the uncertainty in model parameters for the purpose of damage prognosis can be accomplished utilizing Bayesian inference and damage diagnosis data from sources such as non-destructive evaluation or structural health monitoring. The number of samples required to solve the Bayesian inverse problem through common sampling techniques (e.g., Markov chain Monte Carlo) renders high-fidelity finite element-based damage growth models unusable due to prohibitive computation times. However, these types of models are often the only option when attempting to model complex damage growth in real-world structures. Here, a recently developed high-fidelity crack growth model is used which, when compared to finite element-based modeling, has demonstrated reductions in computation times of three orders of magnitude through the use of surrogate models and machine learning. The model is flexible in that only the expensive computation of the crack driving forces is replaced by the surrogate models, leaving the remaining parameters accessible for uncertainty quantification. A probabilistic prognosis framework incorporating this model is developed and demonstrated for non-planar crack growth in a modified, edge-notched, aluminum tensile specimen. Predictions of remaining useful life are made over time for five updates of the damage diagnosis data, and prognostic metrics are utilized to evaluate the performance of the prognostic framework. Challenges specific to the probabilistic prognosis of non-planar fatigue crack growth are highlighted and discussed in the context of the experimental results.

Leser, Patrick E.↗

Statistical tests of additional plate boundaries from plate motion inversions

The application of the F-ratio test, a standard statistical technique, to the results of relative plate motion inversions has been investigated. The method tests whether the improvement in fit of the model to the data resulting from the addition of another plate to the model is greater than that expected purely by chance. This approach appears to be useful in determining whether additional plate boundaries are justified. Previous results have been confirmed favoring separate North American and South American plates with a boundary located beween 30 N and the equator. Using Chase's global relative motion data, it is shown that in addition to separate West African and Somalian plates, separate West Indian and Australian plates, with a best-fitting boundary between 70 E and 90 E, can be resolved. These results are generally consistent with the observation that the Indian plate's internal deformation extends somewhat westward of the Ninetyeast Ridge. The relative motion pole is similar to Minster and Jordan's and predicts the NW-SE compression observed in earthquake mechanisms near the Ninetyeast Ridge.

Stein, S.↗

Tunable ferroelectricity in oxygen-deficient perovskites with Grenier structure

Using first-principles calculations, we predict that tunable ferroelectricity can be realized in oxide perovskites with the Grenier structure and ordered oxygen vacancies. Specifically, we show that R 1/3 A 2/3 FeO 2.67 solids (where R is a rare-earth ion and A an alkaline-earth cation) exhibit polar phases, with a spontaneous polarization tunable by an appropriate choice of R and A. We find that larger cations combined with small R elements lead to a maximum in the polarization and to a minimum in the energy barriers required to switch the sign of the polarization. Ferroelectricity arises from cooperative distortions of octahedral and tetrahedral units, where a combination of rotational and sliding modes controls the emergence of polarization within three-dimensional connected layers. Our results indicate that polar Grenier phases of oxide perovskites are promising materials for microelectronic applications and, in general, for the study of phenomena emerging from breaking inversion symmetry in solids.

Electronic devices↗

Quantum Dots Investigated for Solar Cells

The NASA Glenn Research Center has been investigating the synthesis of quantum dots of CdSe and CuInS2 for use in intermediate-bandgap solar cells. Using quantum dots in a solar cell to create an intermediate band will allow the harvesting of a much larger portion of the available solar spectrum. Theoretical studies predict a potential efficiency of 63.2 percent, which is approximately a factor of 2 better than any state-of-the-art devices available today. This technology is also applicable to thin-film devices--where it offers a potential four-fold increase in power-to-weight ratio over the state of the art. Intermediate-bandgap solar cells require that quantum dots be sandwiched in an intrinsic region between the photovoltaic solar cell's ordinary p- and n-type regions (see the preceding figure). The quantum dots form the intermediate band of discrete states that allow sub-bandgap energies to be absorbed. However, when the current is extracted, it is limited by the bandgap, not the individual photon energies. The energy states of the quantum dot can be controlled by controlling the size of the dot. Ironically, the ground-state energy levels are inversely proportional to the size of the quantum dots. We have prepared a variety of quantum dots using the typical organometallic synthesis routes pioneered by Ba Wendi et al., in the early 1990's. The most studied quantum dots prepared by this method have been of CdSe. To produce these dots, researchers inject a syringe of the desired organometallic precursors into heated triocytlphosphine oxide (TOPO) that has been vigorously stirred under an inert atmosphere (see the following figure). The solution immediately begins to change from colorless to yellow, then orange and red/brown, as the quantum dots increase in size. When the desired size is reached, the heat is removed from the flask. Quantum dots of different sizes can be identified by placing them under a "black light" and observing the various color differences in their fluorescence (see the photograph).

Bailey, Sheila G.↗

Baryon-Interacting Dark Matter: heating dark matter and the emergence of galaxy scaling relations

The empirical scaling relations observed in disk galaxies remain challenging for models of galaxy formation. The most striking among these is the Mass Discrepancy-Acceleration Relation (MDAR), which encodes both a tight baryonic Tully-Fisher relation (BTFR) and the observed diversity of galaxy rotation curves through the central surface density relation (CSDR). Building on our earlier work [1], we propose here that the MDAR is the result of interactions between baryons and ‘Baryon-Interacting Dark Matter’ (BIDM), which heat up the dark matter. Following a bottom-up, hydrodynamical approach, we find that the MDAR follows if: i) the BIDM equation of state approximates that of an ideal gas; ii) the BIDM relaxation time is order the Jeans time; iii) the heating rate is inversely proportional to the BIDM density. Remarkably, under these assumptions the set of hydrodynamical equations together with Poisson’s equation enjoy an anisotropic scaling symmetry. In the BIDM-dominated regime, this gives rise to an enhanced symmetry which fully captures the low-acceleration limit of the MDAR. We then show that, assuming a cored pseudo-isothermal profile at equilibrium, this set of equations gives rise to parameters reproducing the MDAR. Specifically, in the flat part of the rotation curve the asymptotic rotational velocity matches the parametric dependence of the BTFR. Moreover, in the central region of high-surface brightness galaxies, the profile reproduces the CSDR. Finally, by studying the time-dependent approach to equilibrium, we derive a global combination of the BTFR and CSDR, which matches the expectations in low surface-brightness galaxies. The form of the heating rate also makes model-independent predictions for various cosmological observables. We argue that our scenario satisfies existing observational constraints, and, intriguingly, offers a possible explanation to the EDGES anomaly.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗

Residual-based error correction for neural operator accelerated infinite-dimensional Bayesian inverse problems

We explore using neural operators, or neural network representations of nonlinear maps between function spaces, to accelerate infinite-dimensional Bayesian inverse problems (BIPs) with models governed by nonlinear parametric partial differential equations (PDEs). Neural operators have gained significant attention in recent years for their ability to approximate the parameter-to-solution maps defined by PDEs using as training data solutions of PDEs at a limited number of parameter samples. The computational cost of BIPs can be drastically reduced if the large number of PDE solves required for posterior characterization are replaced with evaluations of trained neural operators. However, reducing error in the resulting BIP solutions via reducing the approximation error of the neural operators in training can be challenging and unreliable. We provide an a priori error bound result that implies certain BIPs can be ill-conditioned to the approximation error of neural operators, thus leading to inaccessible accuracy requirements in training. To reliably deploy neural operators in BIPs, we consider a strategy for enhancing the performance of neural operators: correcting the prediction of a trained neural operator by solving a linear variational problem based on the PDE residual. We show that a trained neural operator with error correction can achieve a quadratic reduction of its approximation error, all while retaining substantial computational speedups of posterior sampling when models are governed by highly nonlinear PDEs. The strategy is applied to two numerical examples of BIPs based on a nonlinear reaction–diffusion problem and deformation of hyperelastic materials. We demonstrate that posterior representations of the two BIPs produced using trained neural operators are greatly and consistently enhanced by error correction.

97 MATHEMATICS AND COMPUTING↗

Plasmaspheric Mass Loss and Refilling as a Result of a Magnetic Storm

Using the sounding measurements from the radio plasma imager on IMAGE, and a plasma density inversion algorithm, we derive the plasma density profiles along the magnetic field in a few L-shells every 14 hours at magnetic local noon before, during, and after the March 31,2001 magnetic storm. An empirical model of the plasmaspheric plasma density distribution is derived as a reference using the measurements before the storm. During the storm, the equatorial plasma was substantially depleted in a range of L-shells. The flux tubes were refilled after the storm. The filling ratio, the equatorial plasma density normalized by its quiet-time value before the storm, is introduced to assess the time evolution of the depletion and refilling processes. The depletion, more than two thirds of the quiet time content, appeared to occur rather quickly after the storm onset, as determined by the limited temporal resolution of the measurements. The refilling proceeded, although more slowly than the depletion process, significantly faster than the theoretical prediction of a 3-day time scale. Dynamic structures are observed in situ and confirmed by the extreme ultraviolet imager (EUV) measurements.

Reinisch, B. W.↗

OSOAA: A Vector Radiative Transfer Model of Coupled Atmosphere-Ocean System for a Rough Sea Surface Application to the Estimates of the Directional Variations of the Water Leaving Reflectance to Better Process Multi-angular Satellite Sensors Data Over the Ocean

In this study, we present a radiative transfer model, so-called OSOAA, that is able to predict the radiance and degree of polarization within the coupled atmosphere-ocean system in the presence of a rough sea surface. The OSOAA model solves the radiative transfer equation using the successive orders of scattering method. Comparisons with another operational radiative transfer model showed a satisfactory agreement within 0.8%. The OSOAA model has been designed with a graphical user interface to make it user friendly for the community. The radiance and degree of polarization are provided at any level, from the top of atmosphere to the ocean bottom. An application of the OSOAA model is carried out to quantify the directional variations of the water leaving reflectance and degree of polarization for phytoplankton and mineral-like dominated waters. The difference between the water leaving reflectance at a given geometry and that obtained for the nadir direction could reach 40%, thus questioning the Lambertian assumption of the sea surface that is used by inverse satellite algorithms dedicated to multi-angular sensors. It is shown as well that the directional features of the water leaving reflectance are weakly dependent on wind speed. The quantification of the directional variations of the water leaving reflectance obtained in this study should help to correctly exploit the satellite data that will be acquired by the current or forthcoming multi-angular satellite sensors.

Ocean surface↗

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗

Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from two decades of spectroscopic surveys at the Apache Point Observatory

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Lyα forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, r d , from eight different samples and six measurements of the growth rate parameter, fσ 8 , from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, ΛCDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization, under the same model, the BAO data provide nearly an order of magnitude improvement on curvature constraints relative to primary CMB constraints alone. Independent of distance measurements, the SDSS RSD data complement weak lensing measurements from the Dark Energy Survey (DES) in demonstrating a preference for a flat ΛCDM cosmological model when combined with Planck measurements. The combined BAO and RSD measurements indicate σ 8 =0.85±0.03, implying a growth rate that is consistent with predictions from Planck temperature and polarization data and with General Relativity. When combining the results of SDSS BAO and RSD, Planck, Pantheon Type Ia supernovae (SNe Ia), and DES weak lensing and clustering measurements, all multiple-parameter extensions remain consistent with a ΛCDM model. Regardless of cosmological model, the precision on each of the three parameters, Ω Λ , H 0 , and σ 8 , remains at roughly 1%, showing changes of less than 0.6% in the central values between models. In a model that allows for free curvature and a time-evolving equation of state for dark energy, the combined samples produce a constraint Ω K =-0.0022±0.0022. The dark energy constraints lead to w 0 =-0.909±0.081 and w a =$-0.49_{-0.30}^{+0.35}$, corresponding to an equation of state of w p =-1.018±0.032 at a pivot redshift z p =0.29 and a Dark Energy Task Force Figure of Merit of 94. The inverse distance ladder measurement under this model yields H 0 =68.18±0.79 km s -1 Mpc -1 , remaining in tension with several direct determination methods; the BAO data allow Hubble constant estimates that are robust against the assumption of the cosmological model. In addition, the BAO data allow estimates of H 0 that are independent of the CMB data, with similar central values and precision under a ΛCDM model. Our most constraining combination of data gives the upper limit on the sum of neutrino masses at Σm ν <0.115 eV (95% confidence). Finally, we consider the improvements in cosmology constraints over the last decade by comparing our results to a sample representative of the period 2000-2010. We compute the relative gain across the five dimensions spanned by w, Ω k , Σm ν , H 0 , and σ 8 and find that the SDSS BAO and RSD data reduce the total posterior volume by a factor of 40 relative to the previous generation. Finally, adding again the Planck, DES, and Pantheon SN Ia samples leads to an overall contraction in the five-dimensional posterior volume of 3 orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the Impedance Resonance of Piezoelectric Multi-Fiber Composite Stacks

Multi-Fiber CompositesTM (MFC's) produced by Smart Materials Corp behave essentially like thin planar stacks where each piezoelectric layer is composed of a multitude of fibers. We investigate the suitability of using previously published inversion techniques for the impedance resonances of monolithic co-fired piezoelectric stacks to the MFCTM to determine the complex material constants from the impedance data. The impedance equations examined in this paper are those based on the derivation. The utility of resonance techniques to invert the impedance data to determine the small signal complex material constants are presented for a series of MFC's. The technique was applied to actuators with different geometries and the real coefficients were determined to be similar within changes of the boundary conditions due to change of geometry. The scatter in the imaginary coefficient was found to be larger. The technique was also applied to the same actuator type but manufactured in different batches with some design changes in the non active portion of the actuator and differences in the dielectric and the electromechanical coupling between the two batches were easily measureable. It is interesting to note that strain predicted by small signal impedance analysis is much lower than high field stains. Since the model is based on material properties rather than circuit constants, it could be used for the direct evaluation of specific aging or degradation mechanisms in the actuator as well as batch sorting and adjustment of manufacturing processes.

Multi-fiber composites↗