Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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,783 records · Page 99

Large-scale Gravitational Lens Modeling with Bayesian Neural Networks for Accurate and Precise Inference of the Hubble Constant

We investigate the use of approximate Bayesian neural networks (BNNs) in modeling hundreds of time delay gravitational lenses for Hubble constant (H 0 ) determination. Here, our BNN was trained on synthetic Hubble Space Telescope quality images of strongly lensed active galactic nuclei with lens galaxy light included. The BNN can accurately characterize the posterior probability density functions (PDFs) of model parameters governing the elliptical power-law mass profile in an external shear field. We then propagate the BNN-inferred posterior PDFs into an ensemble H 0 inference, using simulated time delay measurements from a plausible dedicated monitoring campaign. Assuming well-measured time delays and a reasonable set of priors on the environment of the lens, we achieve a median precision of 9.3% per lens in the inferred H 0 . A simple combination of a set of 200 test lenses results in a precision of 0.5 km s –1 Mpc –1 (0.7%), with no detectable bias in this H 0 recovery test. The computation time for the entire pipeline—including the generation of the training set, BNN training and H0 inference—translates to 9 minutes per lens on average for 200 lenses and converges to 6 minutes per lens as the sample size is increased. Being fully automated and efficient, our pipeline is a promising tool for exploring ensemble-level systematics in lens modeling for H 0 inference.

79 ASTRONOMY AND ASTROPHYSICS↗

A rapid-pressure correlation representation consistent with the Taylor-Proudman theorem materially-frame-indifferent in the 2D limit

A nonlinear representation for the rapid-pressure correlation appearing in the Reynolds stress equations, consistent with the Taylor-Proudman theorem, is presented. The representation insures that the modeled second-order equations are frame-invariant with respect to rotation when the flow is two-dimensional in planes perpendicular to the axis of rotation. The representation satisfies realizability in a new way: a special ansatz is used to obtain analytically, the values of coefficients valid away from the realizability limit: the model coefficients are functions of the state of the turbulence that are valid for all states of the mechanical turbulence attaining their constant limiting values only when the limit state is achieved. Utilization of all the mathematical constraints are not enough to specify all the coefficients in the model. The unspecified coefficients appear as free parameters which are used to insure that the representation is asymptotically consistent with the known equilibrium states of a homogeneous sheared turbulence. This is done by insuring that the modeled evolution equations have the same fixed points as those obtained from computer and laboratory experiments for the homogeneous shear. Results of computations of the homogeneous shear, with and without rotation, and with stabilizing and destabilizing curvature, are shown. Results are consistently better, in a wide class of flows which the model not been calibrated, than those obtained with other nonlinear models.

Ristorcelli, J. R.↗

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis↗

Segmented Mirror Image Degradation Due to Surface Dust, Alignment and Figure

In 1996 an algorithm was developed to include the effects of surface roughness in the calculation of the point spread function of a telescope mirror. This algorithm has been extended to include the effects of alignment errors and figure errors for the individual elements, and an overall contamination by surface dust. The final algorithm builds an array for a guard-banded pupil function of a mirror that may or may not have a central hole, a central reflecting segment, or an outer ring of segments. The central hole, central reflecting segment, and outer ring may be circular or polygonal, and the outer segments may have trimmed comers. The modeled point spread functions show that x-tilt and y-tilt, or the corresponding R-tilt and theta-tilt for a segment in an outer ring, is readily apparent for maximum wavefront errors of 0.1 lambda. A similar sized piston error is also apparent, but integral wavelength piston errors are not. Severe piston error introduces a focus error of the opposite sign, so piston could be adjusted to compensate for segments with varying focal lengths. Dust affects the image principally by decreasing the Strehl ratio, or peak intensity of the image. For an eight-meter telescope a 25% coverage by dust produced a scattered light intensity of 10(exp -9) of the peak intensity, a level well below detectability.

Schreur, Julian J.↗

PSF Estimation in Crowded Astronomical Imagery as a Convolutional Dictionary Learning Problem

In this paper, we present a new algorithm for estimating the Point Spread Function (PSF) in wide-field astronomical images with extreme source crowding. Robust and accurate PSF estimation in crowded astronomical images dramatically improves the fidelity of astrometric and photometric measurements extracted from wide-field sky monitoring imagery. Our radically new approach utilizes convolutional sparse representations to model the continuous functions involved in the image formation. This approach avoids the need to detect and precisely localize individual point sources that is shared by existing methods. In experiments involving simulated astronomical imagery, it significantly outperforms the recent alternative method with which it is compared.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxies in X-ray selected clusters and groups in Dark Energy Survey data – II. Hierarchical Bayesian modelling of the red-sequence galaxy luminosity function

Using ~100 X-ray selected clusters in the Dark Energy Survey Science Verification data, in this work we constrain the luminosity function (LF) of cluster red-sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the LF simultaneously in redshift (0.1 < z < 1.05) and cluster mass ($13.5 \le \rm {log_{10}}(M_{200crit}) \sim \lt 15.0$). We pay particular attention to completeness issues and the detection limit of the galaxy sample. We then apply a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schechter function, including its characteristic break (m*) to a faint end power-law slope (α). Our approach allows us to avoid known issues in similar analyses based on stacking or binning the clusters. We find weak and statistically insignificant (~1.9σ) evolution in the faint end slope α versus redshift. We also find no dependence in α or m* with the X-ray inferred cluster masses. Yet, the amplitude of the LF as a function of cluster mass is constrained to ${\sim } 20{{\ \rm per\ cent}}$ precision. As a by-product of our algorithm, we utilize the correlation between the LF and cluster mass to provide an improved estimate of the individual cluster masses as well as the scatter in true mass given the X-ray inferred masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, greatly improving the sensitivity of the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Hollow Cathode Thermal Modelling and Self-Consistent Plasma Solution: Work Function Evaluation for a LaB6 Cathode.

Self-heating hollow cathodes are central components in modern electric thrusters. The plasma discharge inside these devices heats the internal components, thus maintaining the temperatures required for electron emission. Precise knowledge of the physical phenomena governing hollow cathode operation is key to predict their lifetime, specifically, their thermionic emission characteristics. A simulation platform has been built to couple plasma and thermal models of the self-heating hollow cathode to produce a self-consistent solution. A self-consistent solution has been found for a LaB6 hollow cathode operating at 25A and 13 sccm where the work function is assumed to be spatially uniform along the emitter with a value which is allowed to vary as the coupled model iterates to a self-consistent solution. The emitter temperature from the converged solution does not agree with experimental temperature measurements, however. The results of a sensitivity analysis suggest that none of the tolerances in the measurement are responsible for the discrepancy. We hypothesize that either the work function needs to be a function of position along the emitting surfaces or the heat fluxes have been overestimated in the plasma solver.

Meiron, Daniel I.↗

Measurement-based quantum simulation of Abelian lattice gauge theories

Numerical simulation of lattice gauge theories is an indispensable tool in high energy physics, and their quantum simulation is expected to become a major application of quantum computers in the future. In this work, for an Abelian lattice gauge theory in d d spacetime dimensions, we define an entangled resource state (generalized cluster state) that reflects the spacetime structure of the gauge theory. We show that sequential single-qubit measurements with the bases adapted according to the former measurement outcomes induce a deterministic Hamiltonian quantum simulation of the gauge theory on the boundary. Our construction includes the (2+1) ( 2 + 1 ) -dimensional Abelian lattice gauge theory simulated on three-dimensional cluster state as an example, and generalizes to the simulation of Wegner’s lattice models M_{(d,n)} M ( d , n ) that involve higher-form Abelian gauge fields. We demonstrate that the generalized cluster state has a symmetry-protected topological order with respect to generalized global symmetries that are related to the symmetries of the simulated gauge theories on the boundary.Our procedure can be generalized to the simulation of Kitaev’s Majorana chain on a fermionic resource state. We also study the imaginary-time quantum simulation with two-qubit measurements and post-selections, and a classical-quantum correspondence, where the statistical partition function of the model M_{(d,n)} M ( d , n ) is written as the overlap between the product of two-qubit measurement bases and the wave function of the generalized cluster state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Embracing fine-root system complexity in terrestrial ecosystem modeling

Projecting the dynamics and functioning of the biosphere requires a holistic consideration of whole-ecosystem processes. However, biases toward leaf, canopy, and soil modeling since the 1970s have constantly left fine-root systems being rudimentarily treated. As accelerated empirical advances in the last two decades establish clearly functional differentiation conferred by the hierarchical structure of fine-root orders and associations with mycorrhizal fungi, a need emerges to embrace this complexity to bridge the data-model gap in still extremely uncertain models. In this work, we propose a three-pool structure comprising transport and absorptive fine roots with mycorrhizal fungi (TAM) to model vertically resolved fine-root systems across organizational and spatial–temporal scales. Emerging from a conceptual shift away from arbitrary homogenization, TAM builds upon theoretical and empirical foundations as an effective and efficient approximation that balances realism and simplicity. A proof-of-concept demonstration of TAM in a big-leaf model both conservatively and radically shows robust impacts of differentiation within fine-root systems on simulating carbon cycling in temperate forests. Theoretical and quantitative support warrants exploiting its rich potentials across ecosystems and models to confront uncertainties and challenges for a predictive understanding of the biosphere. Echoing a broad trend of embracing ecological complexity in integrative ecosystem modeling, TAM may offer a consistent framework where modelers and empiricists can work together toward this grand goal.

59 BASIC BIOLOGICAL SCIENCES↗

Lineage-based functional types: characterising functional diversity to enhance the representation of ecological behaviour in Land Surface Models

Process-based vegetation models attempt to represent the wide range of trait variation in biomes by grouping ecologically similar species into plant functional types (PFTs). This approach has been successful in representing many aspects of plant physiology and biophysics but struggles to capture biogeographic history and ecological dynamics that determine biome boundaries and plant distributions. Grass-dominated ecosystems are broadly distributed across all vegetated continents and harbour large functional diversity, yet most Land Surface Models (LSMs) summarise grasses into two generic PFTs based primarily on differences between temperate C 3 grasses and (sub)tropical C 4 grasses. Incorporation of species-level trait variation is an active area of research to enhance the ecological realism of PFTs, which form the basis for vegetation processes and dynamics in LSMs. Using reported measurements, we developed grass functional trait values (physiological, structural, biochemical, anatomical, phenological, and disturbance-related) of dominant lineages to improve LSM representations. In this work, our method is fundamentally different from previous efforts, as it uses phylogenetic relatedness to create lineage-based functional types (LFTs), situated between species-level trait data and PFT-level abstractions, thus providing a realistic representation of functional diversity and opening the door to the development of new vegetation models.

54 ENVIRONMENTAL SCIENCES↗

Unusual Exchange Couplings and Intermediate Temperature Weyl State in Co 3 Sn 2 S 2

Understanding magnetism and its possible correlations to topological properties has emerged to the forefront as a difficult topic in studying magnetic Weyl semimetals. Co 3 Sn 2 S 2 is a newly discovered magnetic Weyl semimetal with a kagome lattice of cobalt ions and has triggered intense interest for rich fantastic phenomena. In this study, we report the magnetic exchange couplings of Co 3 Sn 2 S 2 using inelastic neutron scattering and two density functional theory (DFT) based methods: constrained magnetism and multiple-scattering Green’s function methods. Co 3 Sn 2 S 2 exhibits highly anisotropic magnon dispersions and linewidths below T C , and paramagnetic excitations above T C . The spin-wave spectra in the ferromagnetic ground state is well described by the dominant third-neighbor “across-hexagon” J d model. Our density functional theory calculations reveal that both the symmetry-allowed 120° antiferromagnetic orders support Weyl points in the intermediate temperature region, with distinct numbers and the locations of Weyl points. Our study highlights the important role Co 3 Sn 2 S 2 can play in advancing our understanding of kagome physics and exploring the interplay between magnetism and band topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Application of Reflected GPS Signals to Ocean Remote Sensing

The L-band broadcast signal from the Global Positioning System (GPS) which has reflected off of the sea surface is under study for use as a ocean, coastal and wetlands remote sensing tool. The reflected signal from a given GPS satellite is cross-correlated with the pseudorandom noise code uniquely identifying that satellite. The shape of this cross-correlation, ordinarily a very sharp triangle when tracking a direct line of sight signal, becomes wider and smoother as the mean square slope of the reflecting surface increases. It is proposed that the surface wind speed can be determined by matching the recorded shape of this cross-correlation to that predicted by theoretical models as a function of wind speed and direction. The significance of these effects increases with altitude of the receiver. Experimental data have been collected using a specially modified GPS receiver on aircraft and on a balloon at altitudes of up to 25 km. These data compare favorably with predictions of analytical models and demonstrate the dependence of the waveform shape on surface wind speed and receiver altitude. The advantages that this measurement technique has over conventional scatterometers is the small size, low cost and simplicity of the receiver hardware, no requirement for a transmitter, and the ability to simultaneously collect data from usually 10 or more points (from a low earth orbiting satellite). This number could if the Russian GLONASS (Global Navigation Satellite System) satellites are also considered as additional sources of radiation. Furthermore, the bistatic scattering geometry is complementary to the backscatter used by conventional scatterometers.

Garrison, James L.↗

Oxy-Combustion Modeling for Direct-Fired Supercritical CO 2 Power Cycles

We report that supercritical CO 2 power cycles for fossil energy power generation will likely employ oxy-combustion at very high pressures, possibly exceeding 300 bar. At these high pressures, a direct fired oxy-combustor is more likely to behave like a rocket engine than any type of conventional gas turbine combustor. Issues such as injector design, wall heat transfer, and combustion dynamics may play a challenging role in combustor design. Computational fluid dynamics modeling will not only be useful, but may be a necessity in the combustor design process. To accurately model turbulent reacting flows, combustion submodels appropriate for the conditions of interest as defined by the turbulent time and length scales as well as chemical kinetic time scales are necessary. Here, this paper presents a comparison of various turbulence–chemistry interaction (TCI) modeling approaches on a canonical, single injector, direct-fired sCO 2 combustor. Large eddy simulation is used to model the turbulent combustion process with varying levels of injector oxygen concentration while comparing the effect of the combustion submodel on CO emissions and flame shape. While experimental data are not yet available to validate the simulations, the sensitivity of CO production and flame shape can be studied as a function of combustion modeling approach and oxygen concentration in an effort to better understand how to approach combustion modeling at these unique conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Curie–Weiss model with Complex Temperature: Phase Transitions

We study the partition function of the Curie–Weiss model with complex temperature, and partially describe its phase transitions. As a consequence, we obtain information on the locations of zeros of the partition function (the Fisher zeros).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multipurpose model of Hermes-Columbus docking mechanism

One of the foreseen missions of the HERMES spacevehicle is the servicing to the Columbus Free Flying Laboratory (MTFF). Docking between the two spacecraft is a critical operation in which the Docking Mechanism (DM) has a major role. In order to analyze and assess robustness of initially selected concepts and to identify suitable implementation solutions, through the investigation of main parameters involved in the docking functions, a multipurpose model of DM was developed and tested. This paper describes the main design features as well as the process of calibrating and testing.

Gonzalez-Vallejo, J. J.↗

Constraints on the Optical Depth to Reionization from Balloon-borne Cosmic Microwave Background Measurements

We assess the uncertainty with which a balloon-borne experiment, nominally called Tau Surveyor (τS), can measure the optical depth to reionization σ(τ) with given realistic constraints of instrument noise and foreground emissions. Using a τS fiducial design with six frequency bands between 150 and 380 GHz, with white and uniform map noise of 7 μK arcmin, achievable with a single midlatitude flight, and including Planck's 30 and 44 GHz data, we assess the error σ(τ) obtained with three foreground models and as a function of sky fraction f sky between 40% and 54%. We carry out the analysis using both parametric and blind foreground separation techniques. We compare the σ(τ) values to those obtained with low-frequency and high-frequency versions of the experiment called τS-lf and τS-hf, which have only four and up to eight frequency bands with narrower and wider frequency coverage, respectively. We find that with τS, the lowest constraint is σ(τ) = 0.0034, obtained for one of the foreground models with f sky = 54%. σ(τ) is larger, in some cases by more than a factor of 2, for smaller sky fractions, with τS-lf, or as a function of foreground model. The τS-hf configuration does not lead to significantly tighter constraints. The exclusion of the 30 and 44 GHz data, which give information about synchrotron emission, leads to significant τ misestimates. Decreasing noise by an ambitious factor of 10, while keeping f sky = 40%, gives σ(τ) = 0.0031. The combination of σ(τ) = 0.0034, baryon acoustic oscillation data from DESI, and future cosmic microwave background B-mode lensing data from the CMB-S3/CMB-S4 experiments could give σ(∑m v ) = 17 meV.

cosmic microwave background radiation↗

Adaptive Reweighted Variance Estimation for Monte Carlo Eigenvalue Simulations

Monte Carlo (MC) simulation is used to solve the eigenvalue form of the Boltzmann transport equation to estimate various parameters such as fuel pin flux distributions that are crucial for the safe and efficient operation of nuclear systems (e.g., a nuclear reactor). Monte Carlo eigenvalue simulation uses a sample mean over many stationary cycles (iterations) to estimate various parameters important to nuclear systems. A variance estimate of the sample mean is often used for calculating the confidence intervals. However, MC eigenvalue simulation variance estimators that ignore the intercycle correlation underestimate the true variance of the estimated quantity. This paper presents novel data-adaptive approaches based on a simple autoregressive (AR) model and sigmoid functions to improve MC variance estimation. The standard MC sample-based variance estimator (or naïve estimator) and the spectral density–based MC variance estimator are enhanced by adding data-adaptive components that reduce their bias and improve performance. By investigating the frequency pattern of the AR(1) (order 1) model, two adaptive spectral estimators and one adaptive naïve estimator are proposed. The proposed estimators manifest superior performance when applied to three test problems compared to the standard spectral density–based estimator previously introduced by the authors. In conclusion, these new estimators are straightforward, as they use online algorithms and do not require storage of tallies from all active cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗