Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Multi-Step Adaptive Interpolation Approach to Mitigating the Impact of Nonlinearity on Dynamic State Estimation

Accurate estimation of dynamic states (such as synchronous machine rotor angle and speed) is critical for monitoring and controlling transient stability. Extended Kalman filter (EKF)--based approaches have been developed for estimating dynamic states. In order to improve the EKF's performance in estimating dynamic states of a synchronous machine, this paper proposes a multi-step adaptive interpolation (MSAI) approach to achieve balance between estimation accuracy and computational time. This approach consists of three major steps. First, two indexes are calculated to quantify the nonlinearity of the state transition function and measurement function, respectively. Second, based on the nonlinearity indexes, the interpolation factor is determined using a finite state machine model. And finally, to mitigate the negative impact of nonlinearity on the estimation accuracy, pseudo-measurements are added between consecutive measurements through linear interpolation. A simple example is used to validate the proposed nonlinear indexes. The two-area four-machine system and 16-machine 68-bus system are used to evaluate the effectiveness of the proposed MSAI approach. It is shown through the Monte-Carlo method that the estimation accuracy can be improved through interpolation. In addition, a good trade-off between estimation accuracy and computational time can be achieved effectively through the proposed MSAI approach.

Dynamic State Estimation, extended Kalman filter (↗
DOE OSTI2023

Calibrating Constant Elasticity of Substitution Technologies to Bottom-up Cost Estimates

We propose a method for calibrating an industry-level technology to engineering (bottom-up) estimates with a particular focus on abatement opportunities. As a demonstration, substitution elasticities across inputs are adjusted in the nested cost function for the electricity sector to best fit a target marginal abatement cost (MAC) curve derived from engineering assessments of available technologies. Elasticities are optimized over an entire relevant range of the MAC, whereas current techniques use local point estimates under little or no abatement. In the context of fitting to a given MAC we evaluate alternative nesting structures and find that, while complexity in nesting improves the fit, even relatively simple nesting structures can reasonably approximate the target MAC. In our example, focused on the electricity sector, we find standard elasticities adopted in top-down models moderately overstate abatement costs relative to the engineering targets. In our preferred specification the most important adjustment is to escalate the substitution elasticity between energy and value-added inputs. This is consistent with an argument that the current set of point estimates fail to properly account for new capital-based technologies. These conclusions, however, are sensitive to our assumption about output-intensity abatement and consumer price responsiveness, both of which are not delineated in engineering estimates.

abatement cost↗
DOE OSTI2015

Vertical deformation at western part of Sumatra

This research tries to make advancement in GPS signal processing to estimate the interseismic vertical deformation field at western part of Sumatra Island. The data derived by Continuous Global Positioning System (CGPS) from Badan Informasi Geospasial (BIG) between 2010 and 2012. GPS Analyze at Massachusetts Institute of Technology (GAMIT) software and Global Kalman Filter (GLOBK) software are used to process the GPS signal to estimate the vertical velocities of the CGPS station. In order to minimize noise due to atmospheric delay, Vienna Mapping Function 1 (VMF1) is used as atmospheric parameter model and include daily IONEX file provided by the Center for Orbit Determination in Europe (CODE) as well. It improves GAMIT daily position accuracy up to 0.8 mm. In a second step of processing, the GLOBK is used in order to estimate site positions and velocities in the ITRF08 reference frame. The result shows that the uncertainties of estimated displacement velocity at all CGPS stations are smaller than 1.5 mm/yr. The subsided deformation patterns are seen at the northern and southern part of west Sumatra. The vertical deformation at northern part of west Sumatra indicates postseismic phase associated with the 2010 and 2012 Northern Sumatra earthquakes and also the long-term postseismic associated with the 2004 and 2005 Northern Sumatra earthquakes. The uplifted deformation patterns are seen from Bukit Tinggi to Seblat which indicate a long-term interseismic phase after the 2007 Bengkulu earthquake and 2010 Mentawai earthquake. GANO station shows a subsidence at rate 12.25 mm/yr, indicating the overriding Indo-Australia Plate which is dragged down by the subducting Southeast Asian Plate.

58 GEOSCIENCES↗
DOE OSTI2020

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗
DOE OSTI2016

Molecular transformations of phenolic SOA during photochemical aging in the aqueous phase: competition among oligomerization, functionalization, and fragmentation

Abstract. Organic aerosol is formed and transformed in atmospheric aqueous phases (e.g., cloud and fog droplets and deliquesced airborne particles containing small amounts of water) through a multitude of chemical reactions. Understanding these reactions is important for a predictive understanding of atmospheric aging of aerosols and their impacts on climate, air quality, and human health. In this study, we investigate the chemical evolution of aqueous secondary organic aerosol (aqSOA) formed during reactions of phenolic compounds with two oxidants – the triplet excited state of an aromatic carbonyl (3C∗) and hydroxyl radical (•OH). Changes in the molecular composition of aqSOA as a function of aging time are characterized using an offline nanospray desorption electrospray ionization mass spectrometer (nano-DESI MS) whereas the real-time evolution of SOA mass, elemental ratios, and average carbon oxidation state (OSC) are monitored using an online aerosol mass spectrometer (AMS). Our results indicate that oligomerization is an important aqueous reaction pathway for phenols, especially during the initial stage of photooxidation equivalent to ∼ 2 h irradiation under midday winter solstice sunlight in Northern California. At later reaction times functionalization (i.e., adding polar oxygenated functional groups to the molecule) and fragmentation (i.e., breaking of covalent bonds) become more important processes, forming a large variety of functionalized aromatic and open-ring products with higher OSC values. Fragmentation reactions eventually dominate the photochemical evolution of phenolic aqSOA, forming a large number of highly oxygenated ring-opening molecules with carbon numbers (nC) below 6. The average nC of phenolic aqSOA decreases while average OSC increases over the course of photochemical aging. In addition, the saturation vapor pressures (C∗) of dozens of the most abundant phenolic aqSOA molecules are estimated. A wide range of C∗ values is observed, varying from < 10−20 µg m−3 for functionalized phenolic oligomers to > 10 µg m−3 for small open-ring species. The detection of abundant extremely low-volatile organic compounds (ELVOC) indicates that aqueous reactions of phenolic compounds are likely an important source of ELVOC in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Novel angular velocity estimation technique for plasma filaments

Magnetic field aligned filaments such as blobs and edge localized mode filaments carry significant amounts of heat and particles to the plasma facing components and they decrease their lifetime. The dynamics of these filaments determine at least a part of the heat and particle loads. These dynamics can be characterized by their translation and rotation. In this paper, we present an analysis method novel for fusion plasmas, which can estimate the angular velocity of the filaments on frame-by-frame time resolution. After pre-processing, the frames are two-dimensional (2D) Fourier-transformed, then the resulting 2D Fourier magnitude spectra are transformed to log-polar coordinates, and finally the 2D cross-correlation coefficient function (CCCF) is calculated between the consecutive frames. The displacement of the CCCF’s peak along the angular coordinate estimates the angle of rotation of the most intense structure in the frame. Further, the proposed angular velocity estimation method is tested and validated for its accuracy and robustness by applying it to rotating Gaussian-structures. The method is also applied to gas-puff imaging measurements of filaments in National Spherical Torus Experiment plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

CHERENCUBE: Concept definition and implementation challenges of a Cherenkov-based detector block for PET

Purpose: A new concept for a depth-of-interaction (DOI) capable time-of-flight (TOF) PET detector is defined, based only on the detection of Cherenkov photons. The proposed “CHERENCUBE” consists of a cubic Cherenkov radiator with position-sensitive photodetectors covering each crystal face. By means of the spatial distribution of the detected photons and their time of arrival, the point of interaction of the gamma-ray in the crystal can be determined. This study analyzes through theoretical calculations and Monte Carlo simulations the potential advantages of the concept toward reaching a Cherenkov-only detector for TOF-PET with DOI capability. Furthermore, an algorithm for the DOI estimation is presented and the requirements for a practical implementation of the proposed concept are defined. Methods: The Monte Carlo simulations consisted of a cubic crystal with one photodetector coupled to each one of the faces of the cube. The sensitive area of the detector matched exactly the crystal size, which was varied in 1 mm steps between 1 × 1 × 1 mm{sup 3} and 10 × 10 × 10 mm{sup 3}. For each size, five independent simulations of ten thousand 511 keV gamma-rays were triggered at a fixed distance of 10 mm. The crystal chosen was PbWO{sub 4}. Its scintillation properties were simulated, but only Cherenkov photons were analyzed. Photodetectors were simulated having perfect photodetection efficiency and infinite time resolution. For every generated particle, the analysis considered its creation process, parent and daughter particles, energy, origin coordinates, trajectory, and time and position of detection. The DOI determination is based on the distribution of the emission time of all photons per event. These values are calculated as a function of the coordinates of detection and origin for every photon. The common origin is estimated by finding the distribution with the most similar emission time-points. Results: Detection efficiency increases with crystal size from 8.2% (1 × 1 × 1 mm{sup 3}) to 58.6% (10 × 10 × 10 mm{sup 3}) and decreases applying a photon detection threshold of 5/10/20 photons to 6.3%/4.3%/0.7% and 49.3%/30.4%/2.8%, respectively. The detection rate in the six photodetectors is uniform due to the nearly isotropic cone emission. Most cones originated after a photoelectric effect interaction, with two dominating peaks for the kinetic energy of the electron at 422.99 and 441.47 keV. The detection distance between same-event photons defines the spatial resolution of the detector required for individual photon recognition, with 20% of the detected photons having their closest neighbor within a distance of 5% of the length of the cube. Same-event photons are detected within a time window whose width is determined by the crystal size, with values of 30 and 150 ps for a 1 × 1 × 1 mm{sup 3} and a 10 × 10 × 10 mm{sup 3} cube, respectively. The DOI reconstruction has an accuracy of approximately 23% of the length of the cube, with an average value of 2.2 mm for a 10 × 10 × 10 mm{sup 3} CHERENCUBE. Conclusions: The proposed concept requires a detector with high photodetection efficiency. The structure of the sensitive surface of the detector should be a two dimensional array of microcells, able to provide individual detection coordinates and time stamps. The microcell size determines the ability to recognize individual photons, influencing detection efficiency. The 3D DOI recognition relies on the accuracy of the time stamps and detection coordinates, without the need for a recognition of the projected patterns of photons. The refractive index of the material defines a detector intrinsic energy-based rejection of scattered PET events at the cost of reduced sensitivity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

The image quality of ion computed tomography at clinical imaging dose levels

Purpose: Accurately predicting the range of radiotherapy ions in vivo is important for the precise delivery of dose in particle therapy. Range uncertainty is currently the single largest contribution to the dose margins used in planning and leads to a higher dose to normal tissue. The use of ion CT has been proposed as a method to improve the range uncertainty and thereby reduce dose to normal tissue of the patient. A wide variety of ions have been proposed and studied for this purpose, but no studies evaluate the image quality obtained with different ions in a consistent manner. However, imaging doses ion CT is a concern which may limit the obtainable image quality. In addition, the imaging doses reported have not been directly comparable with x-ray CT doses due to the different biological impacts of ion radiation. The purpose of this work is to develop a robust methodology for comparing the image quality of ion CT with respect to particle therapy, taking into account different reconstruction methods and ion species. Methods: A comparison of different ions and energies was made. Ion CT projections were simulated for five different scenarios: Protons at 230 and 330 MeV, helium ions at 230 MeV/u, and carbon ions at 430 MeV/u. Maps of the water equivalent stopping power were reconstructed using a weighted least squares method. The dose was evaluated via a quality factor weighted CT dose index called the CT dose equivalent index (CTDEI). Spatial resolution was measured by the modulation transfer function. This was done by a noise-robust fit to the edge spread function. Second, the image quality as a function of the number of scanning angles was evaluated for protons at 230 MeV. In the resolution study, the CTDEI was fixed to 10 mSv, similar to a typical x-ray CT scan. Finally, scans at a range of CTDEI’s were done, to evaluate dose influence on reconstruction error. Results: All ions yielded accurate stopping power estimates, none of which were statistically different from the ground truth image. Resolution (as defined by the modulation transfer function = 10% point) was the best for the helium ions (18.21 line pairs/cm) and worst for the lower energy protons (9.37 line pairs/cm). The weighted quality factor for the different ions ranged from 1.23 for helium to 2.35 for carbon ions. For the angle study, a sharp increase in absolute error was observed below 45 distinct angles, giving the impression of a threshold, rather than smooth, limit to the number of angles. Conclusions: The method presented for comparing various ion CT modalities is feasible for practical use. While all studied ions would improve upon x-ray CT for particle range estimation, helium appears to give the best results and deserves further study for imaging.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Derivation of the Two-Exponential Probability Density Function for Rossi-Alpha Measurements of Reflected Assemblies and Validation for the Special Case of Shielded Measurements

Rossi-alpha measurements of fissionable assemblies are used to estimate the prompt neutron decay constant, α, of a system, which can be used to infer reactivity. Previous works have demonstrated that two-exponential fits for Rossi-alpha measurements of reflected assemblies are better than one exponential fits; however, the two-exponential probability density function that is needed to obtain α from the fit parameters has not been derived. This paper derives the two-region point kinetics model (two-exponential fit) for Rossi-alpha measurements, a generalization of the current, one region model (one-exponential fit). The new model is validated with measurements of reflected, fissionable assemblies (highly enriched uranium with k eff ≈ 0.95 and weapons-grade plutonium with k eff > 0.77) with α held constant. Here, the results show that the twoexponential model can (i) predict the constant α, and (ii) deconvolve α and the time a neutron spends in the reflector, neither of which is possible with the one-exponential model.

42 ENGINEERING↗
DOE OSTI2016

Lensing-induced morphology changes in CMB temperature maps in modified gravity theories

Lensing of the Cosmic Microwave Background (CMB) changes the morphology of pattern of temperature fluctuations, so topological descriptors such as Minkowski Functionals can probe the gravity model responsible for the lensing. We show how the recently introduced two-to-two and three-to-one kurt-spectra (and their associated correlation functions), which depend on the power spectrum of the lensing potential, can be used to probe modified gravity theories such as f ( R ) theories of gravity and quintessence models. We also investigate models based on effective field theory, which include the constant-Ω model, and low-energy Hořava theories. Estimates of the cumulative signal-to-noise for detection of lensing-induced morphology changes, reaches O(10{sup 3}) for the future planned CMB polarization mission COrE{sup +}. Assuming foreground removal is possible to ℓ{sub max}=3000, we show that many modified gravity theories can be rejected with a high level of significance, making this technique comparable in power to galaxy weak lensing or redshift surveys. These topological estimators are also useful in distinguishing lensing from other scattering secondaries at the level of the four-point function or trispectrum. Examples include the kinetic Sunyaev-Zel'dovich (kSZ) effect which shares, with lensing, a lack of spectral distortion. We also discuss the complication of foreground contamination from unsubtracted point sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Van der Waals coefficients beyond the classical shell model

Van der Waals (vdW) coefficients can be accurately generated and understood by modelling the dynamic multipole polarizability of each interacting object. Accurate static polarizabilities are the key to accurate dynamic polarizabilities and vdW coefficients. In this work, we present and study in detail a hollow-sphere model for the dynamic multipole polarizability proposed recently by two of the present authors (JT and JPP) to simulate the vdW coefficients for inhomogeneous systems that allow for a cavity. The inputs to this model are the accurate static multipole polarizabilities and the electron density. A simplification of the full hollow-sphere model, the single-frequency approximation (SFA), circumvents the need for a detailed electron density and for a double numerical integration over space. We find that the hollow-sphere model in SFA is not only accurate for nanoclusters and cage molecules (e.g., fullerenes) but also yields vdW coefficients among atoms, fullerenes, and small clusters in good agreement with expensive time-dependent density functional calculations. However, the classical shell model (CSM), which inputs the static dipole polarizabilities and estimates the static higher-order multipole polarizabilities therefrom, is accurate for the higher-order vdW coefficients only when the interacting objects are large. For the lowest-order vdW coefficient C{sub 6}, SFA and CSM are exactly the same. The higher-order (C{sub 8} and C{sub 10}) terms of the vdW expansion can be almost as important as the C{sub 6} term in molecular crystals. Application to a variety of clusters shows that there is strong non-additivity of the long-range vdW interactions between nanoclusters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

COPSS II: the molecular gas content of ten million cubic megaparsecs at redshift z ∼ 3

We present a measurement of the abundance of carbon monoxide in the early universe, utilizing the final results from the CO Power Spectrum Survey (COPSS). Between 2013 and 2015, we performed observations with the Sunyaev–Zel’dovich Array to measure aggregate CO emission from z∼3 galaxies with the intensity mapping technique. Data were collected on 19 fields, covering an area of 0.7 square degrees, over the frequency range 27--35 GHz. With these data, along with data analyzed in COPSS I, we are able to observe the CO(1–0) transition within the redshift range z=2.3--3.3 for spatial frequencies between k=0.5--10 h Mpc{sup −1}, spanning a comoving volume of 4.9×10{sup 6} h{sup −3} Mpc{sup 3}. We present estimates of contributions from continuum sources and ground illumination within our measurement. We constrain the amplitude of the CO power spectrum to P{sub CO}=3.0{sub −1.3}{sup +1.3}×10{sup 3} μK{sup 2}(h{sup −1}Mpc){sup 3}, or Δ{sub CO}{sup 2}(k=1 h Mpc{sup −1})=1.5{sub −0.7}{sup +0.7}×10{sup 3} μK{sup 2}, at 68% confidence, and P{sub CO}>0 at 98.9% confidence. These results are a factor of 10 improvement in sensitivity compared to those of COPSS I. With this measurement, we constrain on the CO(1–0) galaxy luminosity function at z∼3. Assuming that CO emission is proportional to halo mass and using theoretical estimates of the scatter in this relationship, we constrain the ratio of CO(1--0) luminosity to halo mass to A{sub CO}=6.3{sub −2.1}{sup +1.4}×10{sup −7} L{sub ⊙} M{sub ⊙}{sup −1}. Assuming a Milky Way-like linear relationship between CO luminosity and molecular gas mass, we estimate a mass fraction of molecular gas of f{sub H{sub 2}}=5.5{sub −2.2}{sup +3.4}×10{sup −2} for halos with masses of ∼10{sup 12}M{sub ⊙}. Using theoretical estimates for the scaling of molecular gas mass fraction and halo mass, we estimate the cosmic molecular gas density to be ρ{sub z∼3}(H{sub 2})=1.1{sub −0.4}{sup +0.7}×10{sup 8} M{sub ⊙} Mpc{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Moment Magnitudes of Local/Regional Events from 1D Coda Calibrations in the Broader Middle East Region [Moment Magnitudes of Small to Moderate Size Regional Events from Coda Calibration in the Broader Middle East]

Reliable moment magnitude estimates for seismic events in the Middle East region can be difficult to obtain due to the uneven distribution of stations, the complex tectonic structure, and regions of high attenuation. In this paper, we take advantage of the many new broadband seismic stations that have become available through improved national networks and numerous temporary deployments. We make coda envelope-amplitude measurements for 2247 events recorded by 68 stations over 13 narrow frequency bands ranging between 0.03 and 8 Hz. The absolute scaling of these spectra was calculated based on independent waveform modeling solutions of the moment magnitudes for a subset of these events to avoid circularity. Using our 1D path calibrations, we determined coda-based magnitudes for a majority of the events. We obtain fairly good agreement with waveform-modeled seismic moments for the larger events (Mw>4.5) at low frequencies (<0.7 Hz). As expected, the coda-derived source spectra become increasingly scattered at higher frequencies (>0.7 Hz) because of unaccounted 2D path effects, as well as mixing of both Sn coda and Lg coda, which have different attenuation behavior. This scatter leads to increased variance in the magnitudes estimated for smaller events in which low-frequency amplitudes are below the noise levels and the higher frequencies are the only signals available. We quantify the expected variance in coda envelope amplitudes as a function of frequency using interstation scatter as our metric. Finally, the net results of this study provide thousands of new 1D coda magnitude estimates for events in the broad region, as well as the necessary initial starting model for use in a new related 2D coda study (Pasyanos et al., 2016).

58 GEOSCIENCES↗
DOE OSTI2015

Uncertainty quantification and propagation in nuclear density functional theory

Nuclear density functional theory (DFT) is one of the main theoretical tools used to study the properties of heavy and superheavy elements, or to describe the structure of nuclei far from stability. While on-going eff orts seek to better root nuclear DFT in the theory of nuclear forces, energy functionals remain semi-phenomenological constructions that depend on a set of parameters adjusted to experimental data in fi nite nuclei. In this study, we review recent eff orts to quantify the related uncertainties, and propagate them to model predictions. In particular, we cover the topics of parameter estimation for inverse problems, statistical analysis of model uncertainties and Bayesian inference methods. Illustrative examples are taken from the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Early arrival waveform inversion using data uncertainties and matching filters with application to near-surface seismic refraction data

We develop an early arrival waveform inversion (EAWI) technique for high-resolution near-surface velocity estimation by iteratively updating the P-wave velocity model to minimize the difference between the observed and calculated seismic refraction data. Traditional EAWI uses a least-squares penalty function and an acoustic forward-modeling engine. Conventional least-squares error is sensitive to data with low signal-to-noise ratio (S/N) and iterations of EAWI stop at a local-minimum data misfit or at the preassigned maximum number of iterations. These stopping criteria can result in overfitting the data. In addition, fitting the elastic field data with an acoustic modeling engine can introduce artifacts in velocity estimation, especially in land data with significant elastic effects. To overcome these challenges, we develop a robust EAWI (REAWI) method by (1) incorporating the data uncertainties into the penalty function and (2) mitigating the elastic effects using a matching filter workflow. The data uncertainties are estimated from waveform reciprocal errors. When full-waveform reciprocity is not available, trace interpolation is applied. The proposed method prevents closely fitting data with low S/N, avoids overall overfitting by stopping the iterations when a normalized chi-square ([Formula: see text]) waveform misfit of one is achieved, and is less affected by elastic effects. Numerical examples and application to near-surface refraction data at a groundwater contamination site suggest that the final REAWI models are more accurate than the corresponding EAWI models, at the same level of misfit. This is the first known application of a matching filter workflow to real land data. The final REAWI models satisfy an appropriate misfit between the real data and predicted elastic P-wave data, making this approach in this respect equivalent to elastic waveform inversion. We also develop a method to analyze model constraint by examining the energy of the wavefield Fréchet derivative thereby avoiding the influence of the data residuals in traditional Fréchet kernels.

Geochemistry & Geophysics↗
DOE OSTI2017

Entrainment in Superfluid Neutron-Star Crusts: Hydrodynamic Description and Microscopic Origin

In spite of the absence of viscous drag, the neutron superfluid permeating the inner crust of a neutron star cannot flow freely and is entrained by the nuclear lattice similarly to laboratory superfluid atomic gases in optical lattices. The role of entrainment on the neutron superfluid dynamics is reviewed. For this purpose, a minimal hydrodynamical model of superfluidity in neutron-star crusts is presented. This model relies on a fully 4-dimensionally covariant action principle. The equivalence of this formulation with the more traditional approach is demonstrated. In addition, the different treatments of entrainment in terms of dynamical effective masses or superfluid density are clarified. The nuclear energy density functional theory employed for the calculations of all the necessary microscopic inputs is also reviewed, focusing on superfluid properties. In particular, the microscopic origin of entrainment and the different methods to estimate its importance are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗