Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate deconvolution”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Regularizing the linearly extrapolated BDF2 scheme for incompressible flows with time relaxation

This paper presents a highly-efficient finite element scheme for the time relaxation model (TRM). The efficiency is achieved through the second-order BDF2 time-stepping scheme with linear extrapolation (BDF2LE). The accuracy of the scheme is also greatly enhanced through the use of the divergence-free Scott-Vogeulis finite elements, and van Cittert approximate deconvolution. A complete finite element analysis is provided, which includes rigorous proofs for the stability, well-possessedness, and convergence of both velocity and pressure solutions. Furthermore, we also demonstrate that the inclusion of the linear time relaxation term preserves the long-time stability of the unregularized BDF2LE scheme. Finally, numerical experiments are presented that demonstrate the added stability and accuracy that time relaxation can provide.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
NASA NTRS2002

Properties of the Residual Stress of the Temporally Filtered Navier-Stokes Equations

The development of a unifying framework among direct numerical simulations, large-eddy simulations, and statistically averaged formulations of the Navier-Stokes equations, is of current interest. Toward that goal, the properties of the residual (subgrid-scale) stress of the temporally filtered Navier-Stokes equations are carefully examined. Causal time-domain filters, parameterized by a temporal filter width 0 less than Delta less than infinity, are considered. For several reasons, the differential forms of such filters are preferred to their corresponding integral forms; among these, storage requirements for differential forms are typically much less than for integral forms and, for some filters, are independent of Delta. The behavior of the residual stress in the limits of both vanishing and in infinite filter widths is examined. It is shown analytically that, in the limit Delta to 0, the residual stress vanishes, in which case the Navier-Stokes equations are recovered from the temporally filtered equations. Alternately, in the limit Delta to infinity, the residual stress is equivalent to the long-time averaged stress, and the Reynolds-averaged Navier-Stokes equations are recovered from the temporally filtered equations. The predicted behavior at the asymptotic limits of filter width is further validated by numerical simulations of the temporally filtered forced, viscous Burger's equation. Finally, finite filter widths are also considered, and a priori analyses of temporal similarity and temporal approximate deconvolution models of the residual stress are conducted.

Pruett, C. D.↗
NASA NTRS2004

Temporal Large-Eddy Simulation

In 1999, Stolz and Adams unveiled a subgrid-scale model for LES based upon approximately inverting (defiltering) the spatial grid-filter operator and termed .the approximate deconvolution model (ADM). Subsequently, the utility and accuracy of the ADM were demonstrated in a posteriori analyses of flows as diverse as incompressible plane-channel flow and supersonic compression-ramp flow. In a prelude to the current paper, a parameterized temporal ADM (TADM) was developed and demonstrated in both a priori and a posteriori analyses for forced, viscous Burger's flow. The development of a time-filtered variant of the ADM was motivated-primarily by the desire for a unifying theoretical and computational context to encompass direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaged Navier-Stokes simulation (RANS). The resultant methodology was termed temporal LES (TLES). To permit exploration of the parameter space, however, previous analyses of the TADM were restricted to Burger's flow, and it has remained to demonstrate the TADM and TLES methodology for three-dimensional flow. For several reasons, plane-channel flow presents an ideal test case for the TADM. Among these reasons, channel flow is anisotropic, yet it lends itself to highly efficient and accurate spectral numerical methods. Moreover, channel-flow has been investigated extensively by DNS, and a highly accurate data base of Moser et.al. exists. In the present paper, we develop a fully anisotropic TADM model and demonstrate its utility in simulating incompressible plane-channel flow at nominal values of Re(sub tau) = 180 and Re(sub tau) = 590 by the TLES method. The TADM model is shown to perform nearly as well as the ADM at equivalent resolution, thereby establishing TLES as a viable alternative to LES. Moreover, as the current model is suboptimal is some respects, there is considerable room to improve TLES.

Pruett, C. D.↗
NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗
NASA NTRS1992

An isotopic biogeochemical study of the Green River oil shale

Thirty-five different samples from three different sulfur cycles were examined in this stratigraphically oriented study of the Shell 22x-l well (U.S.G.S. C177 core) in the Piceance Basin, Colorado. Carbon isotopic compositions of constituents of Green River bitumens indicate mixing of three main components: products of primary photoautotrophs and their immediate consumers (delta approximately -30% vs PDB), products of methanotrophic bacteria (delta approximately -85%), and products of unknown bacteria (delta approximately -40%). For individual compounds synthesized by primary producers, delta-values ranged from -28 to -32%. 13C contents of individual primary products (beta-carotane, steranes, acyclic isoprenoids, tricyclic triterpenoids) were not closely correlated, suggesting diverse origins for these materials. 13C contents of numerous hopanoids were inversely related to sulfur abundance, indicating that they derived both from methanotrophs and from other bacteria, with abundances of methanotrophs depressed when sulfur was plentiful in the paleoenvironment. gamma-Cerane coeluted with 3 beta(CH3),17 alpha(H),21 beta(H)-hopane, but delta-values could be determined after deconvolution. gamma-Cerane (delta approximately -25%) probably derives from a eukaryotic heterotroph grazing on primary materials, the latter compound (delta approximately -90%) must derive from methanotrophic organisms. 13C contents of n-alkanes in bitumen differed markedly from those of paraffins generated pyrolytically. Isotopic and quantitative relationships suggest that alkanes released by pyrolysis derived from a resistant biopolymer of eukaryotic origin and that this was a dominant constituent of total organic carbon.

Non-NASA Center↗
NASA NTRS1994

Hubble Space Telescope imaging of the Seyfert 2 galaxy NGC 2110

Hubble Space Telescope (HST) images of the Seyfert 2 galaxy NGC 2110 in the emission lines of (O III) lambda lambda-4959, 5007 and H-alpha + (N II) lambda lambda-6548, 6584 and the nearby continuum are presented. After deconvolution, the resolution of these images is approximately 0.1 sec, corresponding to a linear scale of approximately 23 pc at the assumed distance (46.8 Mpc) to this galaxy. The dominant line emission is associated with a strongly curved, narrow, 1 sec (230 pc) long feature emanating from the nucleus to the north. The morphology of this feature suggests that it may represent gas flowing outward from the nucleus (perhaps as a 'jet'). While this optical emission-line feature and the jetlike radio emission extend along a similar P.A., their detailed brightness distributions are loosely anticorrelated. In particular, there is no line emission associated with the off-nuclear 6 cm peaks and only weak radio emission along the optical 'jet'. On a larger scale, there is an S-shaped distribution of ionized gas (seen best in H-alpha) extending approximately 4 sec north and south of the nucleus. Ground-based kinematical data suggest this larger scale component represents normally rotating, ambient, interstellar gas photoionized by the central source. A continuum color map derived from HST observations shows evidence for an unresolved region (scale less than 23 pc) of very strong reddening at the location of the optical and radio nucleus. The location and size of the reddened region are consistent with obscuration by a dusty torus, but the observed reddening provides only a lower limit (A(sub V) greater than or equal to 4.6(sub -1.7 sup +1.2) mag) on the amount of obscuration to the nucleus. Nevertheless, this obscuration is sufficient to account for the 'energy deficit' previously reported in this object.

Mulchaey, John S.↗
NASA NTRS1992

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction.

Abgrall, Remi↗
NASA NTRS1993

Focus history of the Hubble Space Telescope: Launch to May 1993

Since the launch of the Hubble Space Telescope (HST) the secondary mirror of the telescope has been moved several times in order to collimate the telescope and also to define a position of best focus. In addition to these moves the focus position changes over time because of water desorption by the graphite epoxy in the metering truss. The authors report here the focus history of the telescope based on a knowledge of the mirror moves made and an analysis of desorption monitoring data obtained by the Faint Object Camera (FOC) in the F/96 mode and of the routine calibration data obtained by the Wide Field and Planetary Cameras. Focus values are extracted using two different methods. In the first method the distance between the center of the point spread function (PSF) and the shadows of the pads supporting the HST primary mirror are related to the focus error. In the second method an analytical formula for the PSF with variable aberration coefficients is fitted to the data. Focus positions derived from the two methods show good agreement. The data show that a desorption of about 83 microns has taken place since Aug. 16, 1990. The desorption has clearly not leveled off as expected from the trend of the earlier data. Long term variations of the secondary mirror position of approximately 3-15 microns from the 'best' focus position have been observed. Variations of the order of 2-5 microns over an orbital period have also been noted. Focus changes resulting from secondary mirror movements greater than approximately 5 microns changes the point spread function significantly and makes deconvolution and quantitative measurements difficult.

Hasan, H.↗
NASA NTRS1976

Energy spectra of cosmic ray nuclei: 4z26 and .3E2 GeV/amu

Energy spectra of cosmic ray nuclei in the charge range 5 is less than or equal to z less than or equal to 26 have been derived from the response of an acrylic plastic Cerenkov detector. Data were obtained using a balloon borne detector and cover the energy range 320 is approximately less than e approximately less than 2200 MeV. amu. Spectra are derived from a formal deconvolution using the method of Lezniak (1975). Relative spectra of different elements are compared by observing charge ratios. Secondary primary ratios are observed to decrease with increasing energy, consistent with the effect previously observed at higher energy. Primary to primary ratios are constant for 6 is less than or equal to z less than or equal to 26 and 14 is less than or equal to z less than or equal to 26 but vary for 10 is less than or equal to z less than or equal to 14. This data is found to be consistent with existing data where comparable and lends strong support ot the idea of two separate source populations contributing to the cosmic ray composition.

Maehl, R. C.↗
NASA NTRS1983

Comparison of measurements by different instruments of the far-infrared reflectance of rough, optically black coatings

Two comparisons of available specular infrared measurements using identical or very similar coatings reveal that specular BRDF (bidirectional reflectance distribution function) measurements of rough coatings are comparable but that simple measured reflectances are not. This is seen as indicating that the diffuse component of the measured specular reflectance is significantly larger than that of the specular component. It is noted that the diffuse component of the measurement is a convolution of the intrinsic directional reflectivity of the surface with the instrument function and that because the value of that function increases at the wavelengths under consideration, the diffuse component is larger than would otherwise be expected. The BRDF normalization is found to cause a primitive deconvolution of the instrument function and to produce a quantity that is approximately equal to the mean intrinsic directional reflectivity very near the specular direction. The analysis and comparison of data reveal that nonspecular BRDF measurements at angles well outside the instrument function can be calibrated as well by an image of the source upon the detector as by reference to a known surface.

Smith, S. M.↗
NASA NTRS1992

First MEM task IRMEO in IRAF

The first task for image restoration using the Maximum Entropy Method (MEM) in IRAF, called IRMEO, is described. The underlining algorithm is the approximate Newton method for optimization. The basic input images and parameters for deconvolution are described in some detail. Results of preliminary tests, including the number of iterations, required CPU time on a variety of computers, and deconvolved images are reported and compared with those from other deconvolution methods. The merits and limitations of this task are pointed out. The possible development of better MEM tasks on the basis of IRMEO is also discussed.

Wu, Nailong↗
NASA NTRS2010

Mineralogy and Thermal Properties of V-Type Asteroid 956 Elisa: Evidence for Diogenitic Material from the Spitzer IRS (5-35 Micrometers) Spectrum

We present the thermal infrared (5-35 micrometer) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph ("IRS"; Houck,1.R. et .11. [20041. Astrophys, 1. SuppL 154, 18-24) together with new ground-based lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 +/- 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (Hv) at that rotational phase to be 12.58 +/- 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 +/- 0.4 km with a visible albedo pv = 0.142+/- 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 +/- 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 +/- 2.8 K and beaming parameter eta = 1.16 +/- 0.05. Thermophysical modeling places a lower limit of 20 J m(exp -2)K(exp -1)s(exp -1/2) on the thermal inertia of the asteroid's surface layer (if the surface is very smooth) but more likely values fall between 30 and 150 J m(exp -2)K(exp -1)s(exp -1/2) depending on the sense of rotation. The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 micrometer reststrahlen band, the 15-16.5 micrometer Si-O-Si stretching region, and the 16-25 micrometer reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range W0(sub 2+/-1)En(sub 74+/-2)Fs(sub 24+/-1). Spectral deconvolution of the 9-12 micrometer reststrahlen features indicates that up to approximately 20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component. Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J,M., Pieters, C.M., Pratt, S.F.[1990]. J. Geophys. Res, 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite- like composition, which has very rarely been observed among asteroids.

Lim, Lucy F.↗
NASA NTRS1994

Deep HST imaging of distant weak radio and field galaxies

We present deep Hubble Space Telescope (HST) Wide-Field Camera (WFC) V- and I-band images of three distant weak radio galaxies with z = 0.311-2.390 and seven field galaxies with z = 0.131-0.58. The images were deconvolved with both the Lucy and multiresolution CLEAN methods, which yield a restoring Full Width at Half Maximum (FWHM) of less than or equal to 0.2 sec, (nearly) preserve photons and signal-to-noise ratio at low spatial frequencies, and produce consistent light profiles down to our 2 sigma surface brightness sensitivity limit of V approximately 27.2 and I approximately 25.9 mag/sq arcsec. Multi-component image modeling was used to provide deconvolution-independent estimates of structural parameters for symmetric galaxies. We present 12-band (m(sub 2750) UBVRIgriJHK) photometry for a subset of the galaxies and bootstrap the unknown FOC/48 zero point at 2750 A in three independent ways (yielding m(sub 2750) = 21.34 +/- 0.09 mag for 1.0 e(-)/s). Two radio galaxies with z = 0.311 and 0.528, as well as one field galaxy with z = 0.58, have the colors and spectra of early-type galaxies, and a(exp 1/4)-like light profiles in the HST images. The two at z greater than 0.5 have little or no color gradients in V - I and are likely giant ellipticals, while the z = 0.311 radio galaxy has a dim exponential disk and is likely an S0. Six of the seven field galaxies have light profiles that indicate (small) inner bulges following a(exp 1/4) laws and outer exponential disks, both with little or no color gradients. These are (early-type) spiral galaxies with z = 0.131-0.528. About half have faint companions or bars. One shows lumpy structure, possibly a merger. The compact narrow-line galaxy 53W002 at z = 2.390 has less than or = 30% +/- 10% of its HST V and I flux in the central kiloparsec (due to its weak Active Galactic Nucleus (AGN)). Most of its light (V approximately equal to 23.3) occurs in a symmetric envelope with a regular a(exp 1/4)-like profile of effective radius a approximately equal to 1.1 sec (approximately equal to 12 kpc for H(sub 0) = 50, q(sub 0) = 0.1. Its (HST) V - I color varies at most from approximately 0.3 mag at a approximately equal to 0.2 sec to approximately 1.2 mag at a approximately greater than 0.4 sec, and possibly to approximately greater than 2.2 mag at a approximately greater than 1.2 sec. Together with its I - K color (approximately equal to 2.5 mag for a approximately greater than 1.0 sec-2.0 sec), this is consistent with an aging stellar population approximately 0.3-0.5 Gyr old in the galaxy center (a approx. less than 2 kpc radius), and possibly approximately 0.5-1.0 Gyr old at a approximately greater than 10 kpc radius. While its outer part may thus have started to collapse at z = 2.5-4, its inner part still is aligned with its redshifted Ly(alpha) cloud and its radio axis, possibly caused by star formation associated with the radio jet, or by reflection from its AGN cone.

Windhorst, R. A.↗
NASA NTRS1973

Doppler broadening effect on collision cross section functions - Deconvolution of the thermal averaging

The surprising feature of the Doppler problem in threshold determination is the 'amplification effect' of the target's thermal energy spread. The small thermal energy spread of the target molecules results in a large dispersion in relative kinetic energy. The Doppler broadening effect in connection with thermal energy beam experiments is discussed, and a procedure is recommended for the deconvolution of molecular scattering cross-section functions whose dominant dependence upon relative velocity is approximately that of the standard low-energy form.

Bernstein, R. B.↗
NASA NTRS1989

High spatial resolution far-infrared scans of W3(OH)

The strong far-infrared continuum source associated with the OH maser-ultracompact H II region W3(OH) was observed from the Kuiper Airborne Observatory in order to determine its size at 50 and 100 microns. Maximum entropy deconvolutions of the slit scans made in orthogonal directions across the source show it to be approximately circular with FWHM sizes of 8 arcsec and 13 arcsec at 50 and 100 microns, respectively. For these sizes, W3(OH) is optically thick in the far-infrared. A spherically symmetric, dusty molecular cloud model is developed to fit the data.

Campbell, M. F.↗
NASA NTRS1978

The application of a time-domain deconvolution technique for identification of experimental acoustic-emission signals

A method is presented for the signature analysis of pulses by reconstructing in the time domain the shape of the pulse prior to its passing through the measurement system. This deconvolution technique is first evaluated using an idealized system and analytical pulse models and is shown to provide improved results. An experimental situation is then treated; system-component models are developed for the digitizer, tape recorder, filter, transducer and mechanical structure. To accommodate both calibration results and manufacturer's data, and to provide stable mathematical models entails considerable effort: some 30 parameters must be identified to model this system - which is still a substantial approximation - albeit of very high order. Experimental pulses generated by a ball drop, spark discharge and a tearing crack are then deconvoluted 'back through' the system as modeled, using this technique. These results are compared and indicate (a) that consistent shapes may be expected from a given type of source and (b) that some sources can be identified with greater clarity using the deconvolution approach.

Houghton, J. R.↗
DOE OSTI2022

Differentiable Predictions for Large Scale Structure with SHAMNet

In simulation-based models of the galaxy-halo connection, theoretical predictions for galaxy clustering and lensing are typically made based on Monte Carlo realizations of a mock universe. In this paper, we use Subhalo Abundance Matching (SHAM) as a toy model to introduce an alternative to stochastic predictions based on mock population, demonstrating how to make simulation-based predictions for clustering and lensing that are both exact and differentiable with respect to the parameters of the model. Conventional implementations of SHAM are based on iterative algorithms such as Richardson-Lucy deconvolution; here we use the JAX library for automatic differentiation to train SHAMNet, a neural network that accurately approximates the stellar-to-halo mass relation (SMHM) defined by abundance matching. In our approach to making differentiable predictions of large-scale structure, we map parameterized PDFs onto each simulated halo, and calculate gradients of summary statistics of the galaxy distribution by using autodiff to propagate the gradients of the SMHM through the statistical estimators used to measure one- and two-point functions. Our techniques are quite general, and we conclude with an overview of how they can be applied in tandem with more complex, higher-dimensional models, creating the capability to make differentiable predictions for the multi-wavelength universe of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗