Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,225 records · Page 68

Time scales of pattern evolution from cross-spectrum analysis of advanced very high resolution radiometer and coastal zone color scanner imagery

We have selected square subareas (110 km on a side) from coastal zone color scanner (CZCS) and advanced very high resolution radiometer (AVHRR) images for 1981 in the California Current region off northern California for which we could identify sequences of cloud-free data over periods of days to weeks. We applied a two-dimensional fast Fourier transformation to images after median filtering, (x, y) plane removal, and cosine tapering. We formed autospectra and coherence spectra as functions of a scalar wavenumber. Coherence estimates between pairs of images were plotted against time separation between images for several wide wavenumber bands to provide a temporal lagged coherence function. The temporal rate of loss of correlation (decorrelation time scale) in surface patterns provides a measure of the rate of pattern change or evolution as a function of spatial dimension. We found that patterns evolved (or lost correlation) approximately twice as rapidly in upwelling jets as in the 'quieter' regions between jets. The rapid evolution of pigment patterns (lifetime of about 1 week or less for scales of 50-100 km) ought to hinder biomass transfer to zooplankton predators compared with phytoplankton patches that persist for longer times. We found no significant differences between the statistics of CZCS and AVHRR images (spectral shape or rate of decorrelation). In addition, in two of the three areas studied, the peak correlation between AVHRR and CZCS images from the same area occurred at zero lag, indicating that the patterns evolved simutaneously. In the third area, maximum coherence between thermal and pigment patterns occurred when pigment images lagged thermal images by 1-2 days, mirroring the expected lag of high pigment behind low temperatures (and high nutrients) in recently upwelled water. We conclude that in dynamic areas such as coastal upwelling systems, the phytoplankton cells (identified by pigment color patterns) behave largely as passive scalars at the mesoscale and that growth, death, and sinking of phytoplankton collectively play at most a mariginal role in determining the spectral statistics of the pigment patterns.

Denman, Kenneth L.↗

Field-Induced Desorption of Cesium from Rhenium

Desorption of cesium atoms from a rhenium point emitter in an electric field has been studied by method of field-desorption microscopy (FDM). Dependence of the desorbing field on degree of emitter surface coverage with Cs atoms is established. The heat of Cs-atom adsorption, which is dependent on the surface concentration of adsorbate and work function of the surface, is estimated in the framework of the model of mirror-image forces for field-induced desorption. Saturation of the rhenium emitter with carbon leads to changes in the character of desorption and in the distribution of desorption zones on the emitter surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Assessment of the dose reduction potential of a model-based iterative reconstruction algorithm using a task-based performance metrology

Purpose: Different computed tomography (CT) reconstruction techniques offer different image quality attributes of resolution and noise, challenging the ability to compare their dose reduction potential against each other. The purpose of this study was to evaluate and compare the task-based imaging performance of CT systems to enable the assessment of the dose performance of a model-based iterative reconstruction (MBIR) to that of an adaptive statistical iterative reconstruction (ASIR) and a filtered back projection (FBP) technique. Methods: The ACR CT phantom (model 464) was imaged across a wide range of mA setting on a 64-slice CT scanner (GE Discovery CT750 HD, Waukesha, WI). Based on previous work, the resolution was evaluated in terms of a task-based modulation transfer function (MTF) using a circular-edge technique and images from the contrast inserts located in the ACR phantom. Noise performance was assessed in terms of the noise-power spectrum (NPS) measured from the uniform section of the phantom. The task-based MTF and NPS were combined with a task function to yield a task-based estimate of imaging performance, the detectability index (d′). The detectability index was computed as a function of dose for two imaging tasks corresponding to the detection of a relatively small and a relatively large feature (1.5 and 25 mm, respectively). The performance of MBIR in terms of the d′ was compared with that of ASIR and FBP to assess its dose reduction potential. Results: Results indicated that MBIR exhibits a variability spatial resolution with respect to object contrast and noise while significantly reducing image noise. The NPS measurements for MBIR indicated a noise texture with a low-pass quality compared to the typical midpass noise found in FBP-based CT images. At comparable dose, the d′ for MBIR was higher than those of FBP and ASIR by at least 61% and 19% for the small feature and the large feature tasks, respectively. Compared to FBP and ASIR, MBIR indicated a 46%–84% dose reduction potential, depending on task, without compromising the modeled detection performance. Conclusions: The presented methodology based on ACR phantom measurements extends current possibilities for the assessment of CT image quality under the complex resolution and noise characteristics exhibited with statistical and iterative reconstruction algorithms. The findings further suggest that MBIR can potentially make better use of the projections data to reduce CT dose by approximately a factor of 2. Alternatively, if the dose held unchanged, it can improve image quality by different levels for different tasks.

60 APPLIED LIFE SCIENCES↗

Poster - 31: Predicting IQ and hearing loss following radiotherapy in pediatric brain tumors: proton vs photon

Purpose: The increased sparing of normal tissues in intensity modulated proton therapy (IMPT) in pediatric brain tumor treatments should translate into improved neurocognitive outcomes. Models were used to estimate the intelligence quotient (IQ) and the risk of hearing loss 5 years post radiotherapy and to compare outcomes of proton against photon in pediatric brain tumors. Methods: Patients who had received intensity modulated radiotherapy (IMRT) were randomly selected from our retrospective database. The existing planning CT and contours were used to generate IMPT plans. The RBE-corrected dose was calculated for both IMPT and IMRT. For each patient, the IQ was estimated via a Monte Carlo technique, whereas the reported incidence of hearing loss as a function of cochlear dose was used to estimate the probability of occurrence. Results: The integrated brain dose was reduced in all IMPT plans, translating into a gain of 2 IQ points on average for protons for the whole cohort at 5 years post-treatment. In terms of specific diseases, the gains in IQ ranged from 0.8 points for medulloblastoma, to 2.7 points for craniopharyngioma. Hearing loss probability was evaluated on a per-ear-basis and was found to be systematically less for proton versus photon: overall 2.9% versus 7.2%. Conclusions: A method was developed to predict IQ and hearing outcomes in pediatric brain tumor patients on a case-by-case basis. A modest gain was systematically observed for proton in all patients. Given the uncertainties within the model used and our reinterpretation, these gains may be underestimated.

60 APPLIED LIFE SCIENCES↗

Functional specializations in human cerebral cortex analyzed using the Visible Man surface-based atlas

We used surface-based representations to analyze functional specializations in the human cerebral cortex. A computerized reconstruction of the cortical surface of the Visible Man digital atlas was generated and transformed to the Talairach coordinate system. This surface was also flattened and used to establish a surface-based coordinate system that respects the topology of the cortical sheet. The linkage between two-dimensional and three-dimensional representations allows the locations of published neuroimaging activation foci to be stereotaxically projected onto the Visible Man cortical flat map. An analysis of two activation studies related to the hearing and reading of music and of words illustrates how this approach permits the systematic estimation of the degree of functional segregation and of potential functional overlap for different aspects of sensory processing.

Non-NASA Center↗

Observing System Simulations for Small Satellite Formations Estimating Bidirectional Reflectance

The bidirectional reflectance distribution function (BRDF) gives the reflectance of a target as a function of illumination geometry and viewing geometry, hence carries information about the anisotropy of the surface. BRDF is needed in remote sensing for the correction of view and illumination angle effects (for example in image standardization and mosaicing), for deriving albedo, for land cover classification, for cloud detection, for atmospheric correction, and other applications. However, current spaceborne instruments provide sparse angular sampling of BRDF and airborne instruments are limited in the spatial and temporal coverage. To fill the gaps in angular coverage within spatial, spectral and temporal requirements, we propose a new measurement technique: Use of small satellites in formation flight, each satellite with a VNIR (visible and near infrared) imaging spectrometer, to make multi-spectral, near-simultaneous measurements of every ground spot in the swath at multiple angles. This paper describes an observing system simulation experiment (OSSE) to evaluate the proposed concept and select the optimal formation architecture that minimizes BRDF uncertainties. The variables of the OSSE are identified; number of satellites, measurement spread in the view zenith and relative azimuth with respect to solar plane, solar zenith angle, BRDF models and wavelength of reflection. Analyzing the sensitivity of BRDF estimation errors to the variables allow simplification of the OSSE, to enable its use to rapidly evaluate formation architectures. A 6-satellite formation is shown to produce lower BRDF estimation errors, purely in terms of angular sampling as evaluated by the OSSE, than a single spacecraft with 9 forward-aft sensors. We demonstrate the ability to use OSSEs to design small satellite formations as complements to flagship mission data. The formations can fill angular sampling gaps and enable better BRDF products than currently possible.

small satellite↗

THE MASS DISTRIBUTION FUNCTION OF PLANETS

The distribution of orbital period ratios of adjacent planets in extrasolar planetary systems discovered by the Kepler space telescope exhibits a peak near ∼1.5–2, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below ∼1.5. We find from these data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill’s criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed 10{sup −3.2}m{sub ∗}, or about two-thirds of Jupiter’s mass for solar-mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill’s criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of logm/M{sub ⊕} of ∼(0.6−1.0) and ∼(1.1−1.2), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

ECO2N V2.0

ECO2N V2.0 is a fluid property module for the TOUGH2 simulator (Version 2.1) that was designed for applications to geologic sequestration of CO2 in saline aquifers and enhanced geothermal reservoirs. ECO2N V2.0 is an enhanced version of the previous ECO2N V1.0 module (Pruess, 2005). It expands the temperature range up to about 300oC whereas V1.0 can only be used for temperatures below about 110oC. V2.0 includes a comprehensive description of the thermodynamic and thermophysical properties of H2O - NaCl - CO2 mixtures, that reproduces fluid properties largely within experimental error for the temperature, pressure and salinity conditions 10 °C < T < 300 °C, P < 600 bar, and salinity up to halite saturation. This includes density, viscosity, and specific enthalpy of fluid phases as functions of temperature, pressure, and composition, as well as partitioning of mass components H2O, NaCl and CO2 among the different phases. In particular, V2.0 accounts for the effects of water on the thermophysical properties of the CO2-rich phase, which was ignored in V1.0, using a model consistent with the solubility models developed by Spycher and Pruess (2005, 2010). In terms of solubility models, V2.0 uses the same model for partitioning of mass components among the different phases (Spycher and Pruess, 2005) as V1.0 for the low temperature range (<99oC) but uses a new model (Spycher and Pruess, 2010) for the high temperature range (>109oC). In the transition range (99-109oC), a smooth interpolation is applied to estimate the partitioning as a function of the temperature. Flow processes can be modeled isothermally or non-isothermally, and phase conditions represented may include a single (aqueous or CO2-rich) phase, as well as two-phase (brine-CO2) mixtures. Fluid phases may appear or disappear in the course of a simulation, and solid salt may precipitate or dissolve. Note that the model cannot be applied to subcritical conditions that involves both liquid and gaseous CO2 unless thermol process is ignored (i.e.,isothermal run). For those cases, a user may use the fluid property module ECO2M (Pruess, 2011) instead.

Pan, Lehua↗

Near Rectilinear Halo Orbit Determination with Simulated DSN Observations

This paper presents the results of a high-fidelity simulation of spacecraft orbit determination in a near rectilinear halo orbit (NRHO). Others in the literature have examined this problem with linear covariance analysis, but the highly-nonlinear dynamics of this orbit challenge the assumptions underlying such analyses. The present work builds on similar analysis performed by other authors to contribute a fuller understanding of the operational requirements for NRHO navigation. The present work serves as a check to the assumptions of previous studies and an independent verification of those results. The results from the literature are extended by quantifying the space of orbital states from which a spacecraft with given control authority can safely return to the nominal path. Spacecraft state uncertainty estimates are evaluated as a function of time. Simulated range and range-rate measurements with the Deep Space Network (DSN) ground stations are used to model orbit determination accuracy. Orbit maintenance maneuvers are performed using both short-horizon and long-horizon stationkeeping targeting. Monte Carlo analysis of orbit determination and stationkeeping is performed. This paper quantifies the achievable state uncertainty with deep space network (DSN)-only range and range-rate observations. This paper also addresses requirements on the frequency of DSN observation periods and correlates ground contact frequency with navigation accuracy. The results of several related studies are presented and discussed: the effect of missing ground station passes, the effect of missing stationkeeping maneuvers, the sensitivity of the spacecraft state estimate to realistic error sources, and stationkeeping propellant budget.

Nathan L Parrish↗

THE QUEST FOR DUSTY STAR-FORMING GALAXIES AT HIGH REDSHIFT z ≳ 4

We exploit the continuity equation approach and “main-sequence” star formation timescales to show that the observed high abundance of galaxies with stellar masses ≳ a few 10{sup 10} M {sub ⊙} at redshift z ≳ 4 implies the existence of a galaxy population featuring large star formation rates (SFRs) ψ ≳ 10{sup 2} M {sub ⊙} yr{sup −1} in heavily dust-obscured conditions. These galaxies constitute the high-redshift counterparts of the dusty star-forming population already surveyed for z ≲ 3 in the far-IR band by the Herschel Space Observatory . We work out specific predictions for the evolution of the corresponding stellar mass and SFR functions out to z ∼ 10, determining that the number density at z ≲ 8 for SFRs ψ ≳ 30 M {sub ⊙} yr{sup −1} cannot be estimated relying on the UV luminosity function alone, even when standard corrections for dust extinction based on the UV slope are applied. We compute the number counts and redshift distributions (including galaxy-scale gravitational lensing) of this galaxy population, and show that current data from the AzTEC - LABOCA , SCUBA-2 , and ALMA - SPT surveys are already addressing it. We demonstrate how an observational strategy based on color preselection in the far-IR or (sub-)millimeter band with Herschel and SCUBA-2 , supplemented by photometric data from on-source observations with ALMA , can allow us to reconstruct the bright end of the SFR functions out to z ≲ 8. In parallel, such a challenging task can be managed by exploiting current UV surveys in combination with (sub-)millimeter observations by ALMA and NIKA2 and/or radio observations by SKA and its precursors.

79 ASTRONOMY AND ASTROPHYSICS↗

Approximate effect of parameter pseudonoise intensity on rate of convergence for EKF parameter estimators

When using parameter estimation methods based on extended Kalman filter (EKF) theory, it is common practice to assume that the unknown parameter values behave like a random process, such as a random walk, in order to guarantee their identifiability by the filter. The present work is the result of an ongoing effort to quantitatively describe the effect that the assumption of a fictitious noise (called pseudonoise) driving the unknown parameter values has on the parameter estimate convergence rate in filter-based parameter estimators. The initial approach is to examine a first-order system described by one state variable with one parameter to be estimated. The intent is to derive analytical results for this simple system that might offer insight into the effect of the pseudonoise assumption for more complex systems. Such results would make it possible to predict the estimator error convergence behavior as a function of the assumed pseudonoise intensity, and this leads to the natural application of the results to the design of filter-based parameter estimators. The results obtained show that the analytical description of the convergence behavior is very difficult.

Hill, Bryon K.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM (atmospheric radiation measurement) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth observing system) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of 'radiative smoothing' of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. 'Offbeam' Lidar returns are now being used to estimate the cloud 'radiative Green's function' (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the 'Independent Pixel Approximation' or IPA. This and other measurements of 3D (three dimensional) transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international 'Intercomparison of 3D Radiation Codes' (I3RC) program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

ECO2N V. 2.0: A New TOUGH2 Fluid Property Module for Mixtures of Water, NaCl, and CO 2

ECO2N V2.0 is a fluid property module for the TOUGH2 simulator (Version 2.1) that was designed for applications to geologic sequestration of CO 2 in saline aquifers and enhanced geothermal reservoirs. ECO2N V2.0 is an enhanced version of the previous ECO2N V1.0 module (Pruess, 2005). It expands the temperature range up to about 300°C whereas V1.0 can only be used for temperatures below about 110°C. V2.0 includes a comprehensive description of the thermodynamics and thermophysical properties of H 2 O - NaCl -CO 2 mixtures, that reproduces fluid properties largely within experimental error for the temperature, pressure and salinity conditions of interest (10 °C < T < 300 °C; P < 600 bar; salinity up to halite saturation). This includes density, viscosity, and specific enthalpy of fluid phases as functions of temperature, pressure, and composition, as well as partitioning of mass components H 2 O, NaCl and CO 2 among the different phases. In particular, V2.0 accounts for the effects of water on the thermophysical properties of the CO 2 -rich phase, which was ignored in V1.0, using a model consistent with the solubility models developed by Spycher and Pruess (2005, 2010). In terms of solubility models, V2.0 uses the same model for partitioning of mass components among the different phases (Spycher and Pruess, 2005) as V1.0 for the low temperature range (<99°C) but uses a new model (Spycher and Pruess, 2010) for the high temperature range (>109°C). In the transition range (99-109°C), a smooth interpolation is applied to estimate the partitioning as a function of the temperature. Flow processes can be modeled isothermally or non-isothermally, and phase conditions represented may include a single (aqueous or CO 2 -rich) phase, as well as two-phase mixtures. Fluid phases may appear or disappear in the course of a simulation, and solid salt may precipitate or dissolve. This report gives technical specifications of ECO2N V2.0 and includes instructions for preparing input data

15 GEOTHERMAL ENERGY↗

Transferability and accuracy by combining dispersionless density functional and incremental post-Hartree-Fock theories: Noble gases adsorption on coronene/graphene/graphite surfaces

The accuracy and transferability of the electronic structure approach combining dispersionless density functional theory (DFT) [K. Pernal et al., Phys. Rev. Lett. 103, 263201 (2009)] with the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)], are validated for the interaction between the noble-gas Ne, Ar, Kr, and Xe atoms and coronene/graphene/graphite surfaces. This approach uses the method of increments for surface cluster models to extract intermonomer dispersion-like (2- and 3-body) correlation terms at coupled cluster singles and doubles and perturbative triples level, while periodic dispersionless density functionals calculations are performed to estimate the sum of Hartree-Fock and intramonomer correlation contributions. Dispersion energy contributions are also obtained using DFT-based symmetry-adapted perturbation theory [SAPT(DFT)]. An analysis of the structure of the X/surface (X = Ne, Ar, Kr, and Xe) interaction energies shows the excellent transferability properties of the leading intermonomer correlation contributions across the sequence of noble-gas atoms, which are also discussed using the Drude oscillator model. We further compare these results with van der Waals-(vdW)-corrected DFT-based approaches. As a test of accuracy, the energies of the low-lying nuclear bound states supported by the laterally averaged X/graphite potentials (X = {sup 3}He, {sup 4}He, Ne, Ar, Kr, and Xe) are calculated and compared with the best estimations from experimental measurements and an atom-bond potential model using the ab initio-assisted fine-tuning of semiempirical parameters. The bound-state energies determined differ by less than 6–7 meV (6%) from the atom-bond potential model. The crucial importance of including incremental 3-body dispersion-type terms is clearly demonstrated, showing that the SAPT(DFT) approach effectively account for these terms. With the deviations from the best experimental-based estimations smaller than 2.3 meV (1.9%), the accuracy of the combined DFT and post-HF incremental scheme is established for all the noble-gas atoms. With relative deviations smaller than 4% and 11%, good agreement is also achieved by applying the vdW-corrected DFT treatments PBE-D3 and vdW-DF2 for noble-gas atoms heavier than neon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model Predictive Fault-Tolerant Tracking Control for PDF Control Systems With Packet Losses

In this article, a fault-tolerant tracking control strategy is investigated for nonlinear probability density function (PDF) control systems with the actuator fault, uncertainties, unknown disturbance, and random packet losses. The control input signal dropout and measurement signal dropouts are described as the independent Bernoulli distribution. An adaptive fault diagnosis (FD) observer based on the Lyapunov function is given to simultaneously estimate the fault, disturbance, and state with packet losses. Furthermore, different from the traditional robust fault-tolerant control (FTC), a new active fault-tolerant tracking controller is designed based on the model predictive control framework, which has better adaptive fault-tolerant performance. Finally, the validity of the proposed FTC method has been proved by a simulation study of a papermaking process.

42 ENGINEERING↗

Linear functional minimization for inverse modeling

In this paper, we present a novel inverse modeling strategy to estimate spatially distributed parameters of nonlinear models. The maximum a posteriori (MAP) estimators of these parameters are based on a likelihood functional, which contains spatially discrete measurements of the system parameters and spatiotemporally discrete measurements of the transient system states. The piecewise continuity prior for the parameters is expressed via Total Variation (TV) regularization. The MAP estimator is computed by minimizing a nonquadratic objective equipped with the TV operator. We apply this inversion algorithm to estimate hydraulic conductivity of a synthetic confined aquifer from measurements of conductivity and hydraulic head. The synthetic conductivity field is composed of a low-conductivity heterogeneous intrusion into a high-conductivity heterogeneous medium. Our algorithm accurately reconstructs the location, orientation, and extent of the intrusion from the steady-state data only. Finally, addition of transient measurements of hydraulic head improves the parameter estimation, accurately reconstructing the conductivity field in the vicinity of observation locations.

47 OTHER INSTRUMENTATION↗

Intensity Mapping of Hα, Hβ, [O II], and [O III] Lines at z < 5

Intensity mapping is becoming a useful tool to study the large-scale structure of the universe through spatial variations in the integrated emission from galaxies and the intergalactic medium. We study intensity mapping of the Hα 6563 A-ring , [O iii] 5007 Å, [O ii] 3727 Å, and Hβ 4861 A-ring lines at 0.8⩽z⩽5.2. The mean intensities of these four emission lines are estimated using the observed luminosity functions (LFs), cosmological simulations, and the star formation rate density (SFRD) derived from observations at z≲5. We calculate the intensity power spectra and consider the foreground contamination of other lines at lower redshifts. We use the proposed NASA small explorer SPHEREx (the Spectro-Photometer for the History of the universe, Epoch of Reionization, and Ices Explorer) as a case study for the detectability of the intensity power spectra of the four emission lines. We also investigate the cross-correlation with the 21 cm line probed by the Canadian Hydrogen Intensity Mapping Experiment (CHIME), Tianlai experiment and the Square Kilometer Array (SKA) at 0.8⩽z⩽2.4. We find both the auto and cross power spectra can be well measured for the Hα, [O iii] and [O ii] lines at z≲3, while it is more challenging for the Hβ line. Finally, we estimate the constraint on the SFRD from intensity mapping, and find we can reach an accuracy higher than 7% at z≲4, which is better than with the usual method of measurements using the LFs of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗