Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model fitting”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

On-the-fly Doppler Broadening with Probability Table Interpolation for Unresolved Resonance Region in RMC

With the rapid development of computational power and parallel algorithms, Monte Carlo method has been widely investigated and used in neutron transport simulations due to its advantages in geometry modelling and usage of continuous energy point-wise cross sections. Nowadays the road has been paved to couple Monte Carlo codes with thermal-hydraulics codes to make multi-physics calculations. One of the most important aspects of thermal-hydraulics feedback is the detailed temperature distribution, which demands temperature dependent cross sections in neutron transport simulations. The conventional pre-generated cross section libraries at each temperature point are no longer feasible for detailed temperature distributions because of tremendous memory footprint. Hence, different on-the-fly Doppler broadening methods have been proposed for resolved resonance energy region and thermal energy region, such as the interpolation, regress model fitting, multipole representation (MPR), target motion sampling (TMS) and different integration algorithms based on Doppler broadening equation. Compared to resolved resonance energy region and thermal energy region, the representation of cross sections in unresolved resonance region is slightly different. As a result of the limitation of experiment measurements, the overlapped resonances cannot be distinguished from each other and the statistic average resonance parameters are given in the evaluated nuclear data file. In early applications, infinite-dilute pointwise cross sections were obtained from average parameters and were directly used in neutron simulations. However, it was proven that these infinite-dilute pointwise cross sections cannot capture the energy self-shielding effect, thus leading to an overestimate of reaction rates and usually an underestimate of k{sub eff} owing to the increased absorption rates and increased scattering probability to energy below unresolved resonance region where the capture-to-fission ratio is typically higher. This phenomenon is of great importance especially for problems where there is a considerable flux in unresolved resonance region. Then the probability table method was put forward by Levitt to take energy self-shielding effect into account, which provides probability distribution function of cross sections in unresolved resonance region instead of specific cross section values. The probability distribution function is presented in the form of tables and the mean values are equal to the infinite dilution cross sections. To obtain the temperature dependent cross sections in unresolved energy region, there are two different methods. One is to calculate cross sections directly from average resonance parameters and their statistical distributions , which is effective and precise but time-consuming. The other is to interpolate between pre-generated probability table files, which is much more efficient. In this paper, on-the-fly Doppler broadening of cross sections in unresolved resonance regions based on probability tables interpolation was implemented in RMC codes as a supplement to on-the-fly Doppler broadening for resolved resonance energy region and thermal region, which is based on the integration of Doppler broadening equation. The results of criticality calculation and depletion calculation show that on-the-fly Doppler broadening for the whole energy region is accurate and efficient. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Constraining GREA, an alternative theory accounting for the present cosmic acceleration

The origin of the Universe's late-time accelerated expansion remains unknown. The General Relativistic Entropic Acceleration (GREA) theory offers a compelling alternative to $Λ$CDM, attributing cosmic acceleration to entropy growth associated with cosmic and black hole horizons, without invoking a cosmological constant. We test GREA against the latest DESI DR2 Baryon Acoustic Oscillations (BAO), multiple Type Ia supernova compilations (Union3, Pantheon$\texttt{+}$, DES-SN5YR), and cosmic microwave background (CMB) distance measurements. While GREA is not nested within $Λ$CDM, it achieves a comparable goodness-of-fit, highlighting its potential as a theoretically motivated framework that circumvents some of the fine-tuning issues of the standard $Λ$CDM cosmology. We find that the best-fit model features a transient phantom crossing at $z \lesssim 2$, with $w_a\equiv \mathrm{d} w(a=1)/\mathrm{d}a \simeq-0.3$, in good agreement with observations. However, its present-day value $w_0\equiv w(z=0)$ is tightly constrained to be $w_0\simeq-1$. Upcoming low-redshift (i.e. $z < 1$) cosmological probes, from both background and perturbations, will offer promising avenues for further exploring the viability of the GREA theory.

Calderon, R. [Prague, Inst. Phys.] (ORCID:00000002↗
NASA NTRS2007

Fitting a Two-Component Scattering Model to Polarimetric SAR Data from Forests

Two simple scattering mechanisms are fitted to polarimetric synthetic aperture radar (SAR) observations of forests. The mechanisms are canopy scatter from a reciprocal medium with azimuthal symmetry and a ground scatter term that can represent double-bounce scatter from a pair of orthogonal surfaces with different dielectric constants or Bragg scatter from a moderately rough surface, which is seen through a layer of vertically oriented scatterers. The model is shown to represent the behavior of polarimetric backscatter from a tropical forest and two temperate forest sites by applying it to data from the National Aeronautic and Space Agency/Jet Propulsion Laboratory's Airborne SAR (AIRSAR) system. Scattering contributions from the two basic scattering mechanisms are estimated for clusters of pixels in polarimetric SAR images. The solution involves the estimation of four parameters from four separate equations. This model fit approach is justified as a simplification of more complicated scattering models, which require many inputs to solve the forward scattering problem. The model is used to develop an understanding of the ground-trunk double-bounce scattering that is present in the data, which is seen to vary considerably as a function of incidence angle. Two parameters in the model fit appear to exhibit sensitivity to vegetation canopy structure, which is worth further exploration. Results from the model fit for the ground scattering term are compared with estimates from a forward model and shown to be in good agreement. The behavior of the scattering from the ground-trunk interaction is consistent with the presence of a pseudo-Brewster angle effect for the air-trunk scattering interface. If the Brewster angle is known, it is possible to directly estimate the real part of the dielectric constant of the trunks, a key variable in forward modeling of backscatter from forests. It is also shown how, with a priori knowledge of the forest height, an estimate for the attenuation coefficient of the canopy can be obtained directly from the multi-incidence-angle polarimetric observations. This attenuation coefficient is another key variable in forward models and is generally related to the canopy density.

temperate forests↗
NASA NTRS1992

Fitting a three-component scattering model to polarimetric SAR data

A new technique for fitting a three-component scattering mechanism model to the polarimetric synthetic aperture radar (SAR) data itself, without utilizing any ground truth measurements, is presented. The three scattering mechanism components included in the model are volume scatter from randomly oriented dipoles, first-order Bragg surface scatter and a dihedral scattering mechanism for two surfaces with different dielectric constants. The model fit yields an estimate of the contribution to the total backscatter of each of the three components. The backscatter contributions can also be compared to give the relative percentage weight of each. The model fit has an equal number of input parameters (the polarimetric radar backscatter measurements) and output parameters (the backscatter parameters describing them). The model can be applied to entire images or to small areas within an image to give a first-order estimate of the relevant scattering mechanisms. The model was applied to many C-, L- and P-band Airborne SAR (AIRSAR) images of different types of terrain. Results were presented at the workshop.

Freeman, A.↗
NASA NTRS1994

Analysis of helium-rich subdwarf O stars. 1: NLTE models, methods, and fits for 21 Palomar Green survey sdOs

Atmospheric parameters for 21 helium-rich hot subdwarf O stars from the Palomar Green survey are found from fits of non-local thermodynamic equilibrium (NLTE) models to optical spectra. About 250 new NLTE models in the parameter range T(sub eff) from 35,000 to 65,000 K, log (g) from 4.0 to 6.5, and epsilon(He) from 50% He to 99% He have been calculated. A fit for each object is presented. Estimated distances and luminosities are calculated, assuming a mass of 0.5 solar masses. Large distances above the Galactic plane are found, and this implies that the majority of the sdO stars belong to a different stellar population than the sdB, planetary nebulae and white dwarf stars. Possibilities may include the halo and a 'thick disk.' Kinematical data from the literature are also discussed in view of the found distances. Five of the stars have been observed, modeled, and analyzed by Dreizler et al., and significant differences exist between their results and ours for these stars. The reason for these differences is not yet known, but the problem does not alter the conclusions about population membership that we present.

Thejll, P.↗
DOE OSTI2016

A method for simultaneous linear optics and coupling correction for storage rings with turn-by-turn beam position monitor data

We propose a method to simultaneously correct linear optics errors and linear coupling for storage rings using turn-by-turn (TbT) beam position monitor (BPM) data. The independent component analysis (ICA) method is used to isolate the betatron normal modes from the measured TbT BPM data. The betatron amplitudes and phase advances of the projections of the normal modes on the horizontal and vertical planes are then extracted, which, combined with dispersion measurement, are used to fit the lattice model. The fitting results are used for lattice correction. The method has been successfully demonstrated on the NSLS-II storage ring.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A Redox-Electrodialysis Model with Zero Fitting Parameters: Insights into Process Limitations, Design, and Material Interventions

Redox-Electrodialysis (r-ED) is an electrochemical desalination cell architecture that has recently received considerable interest, due to its low energy demand relative to electrochemical desalination technologies that rely on electrode-based ion removal. To further improve the energy efficiency of r-ED, we developed a lumped mathematical model with no adjustable parameters to investigate the various sources of overpotential within the cell. Existing models of electrodialysis and r-ED cells either do not accurately incorporate all phenomena contributing to the overpotential or utilize empirical fitting parameters. The model developed here indicates that ohmic overpotentials, especially in the diluate chamber, are the most significant contributors to energy losses. Based on this insight, we hypothesized that adding an ion exchange resin wafer in the diluate compartment would increase the ionic conductivity and decrease the energy demand. Experimental results showed an 18% reduction in specific energy use while achieving the same degree of salt removal (20 mM to 12 mM). Furthermore, the resin wafer enabled complete desalination to potable drinking levels at a current density previously unachievable within practical operating voltage limits (4.93 mA cm -2 ). We also expanded the model to explore differences in r-ED energy use between configurations using multiple cells and a single cell with increased area.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

RUDR: A New Paradigm for Model Extraction and Fitting

This paper describes a new algorithmic paradigm for solving problems where a model is extracted from or it to data. This paradigm has numerous applications, particularly in computer vision.

algorithmic computer vision RUDR↗
DOE OSTI2014

MODELING THE NONLINEAR CLUSTERING IN MODIFIED GRAVITY MODELS. I. A FITTING FORMULA FOR THE MATTER POWER SPECTRUM OF f(R) GRAVITY

Based on a suite of N-body simulations of the Hu-Sawicki model of f(R) gravity with different sets of model and cosmological parameters, we develop a new fitting formula with a numeric code, MGHalofit, to calculate the nonlinear matter power spectrum P(k) for the Hu-Sawicki model. We compare the MGHalofit predictions at various redshifts (z ≤ 1) to the f(R) simulations and find that the relative error of the MGHalofit fitting formula of P(k) is no larger than 6% at k ≤ 1 h Mpc{sup –1} and 12% at k in (1, 10] h Mpc{sup –1}, respectively. Based on a sensitivity study of an ongoing and a future spectroscopic survey, we estimate the detectability of a signal of modified gravity described by the Hu-Sawicki model using the power spectrum up to quasi-nonlinear scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Rejection Sampling with Autodifferentiation -- Case study: Fitting a Hadronization Model

We present an autodifferentiable rejection sampling algorithm termed Rejection Sampling with Autodifferentiation (RSA). In conjunction with reweighting, we show that RSA can be used for efficient parameter estimation and model exploration. Additionally, this approach facilitates the use of unbinned machine-learning-based observables, allowing for more precise, data-driven fits. To showcase these capabilities, we apply an RSA-based parameter fit to a simplified hadronization model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Mapping Incidence and Prevalence Peak Data for SIR Modeling Applications

Infectious disease modeling and forecasting have played a key role in helping assess and respond to epidemics and pandemics. Recent work has leveraged data on disease peak infection and peak hospital incidence to fit compartmental models for the purpose of forecasting and describing the dynamics of a disease outbreak. Incorporating these data can greatly stabilize a compartmental model fit on early observations, where slight perturbations in the data may lead to model fits that forecast wildly unrealistic peak infection. We introduce a new method for incorporating historic data on the value and time of peak incidence of hospitalization into the fit for a Susceptible-Infectious-Recovered (SIR) model by formulating the relationship between an SIR model’s starting parameters and peak incidence as a system of two equations that can be solved computationally. We demonstrate how to calculate SIR parameter estimates – which describe disease dynamics such as transmission and recovery rates – using this method, and determine that there is a noticeable loss in accuracy whenever prevalence data is misspecified as incidence data. To exhibit the modeling potential, we update the Dirichlet-Beta State Space modeling framework to use hospital incidence data, as this framework was previously formulated to incorporate only data on total infections. This approach is assessed for practicality in terms of accuracy and speed of computation via simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Spectral observations of Ellerman bombs and fitting with a two-cloud model

We study the Hα and Ca II 8542 Å line spectra of four typical Ellerman bombs (EBs) in the active region NOAA 11765 on 2013 June 6, observed with the Fast Imaging Solar Spectrograph installed at the 1.6 m New Solar Telescope at Big Bear Solar Observatory. Considering that EBs may occur in a restricted region in the lower atmosphere, and that their spectral lines show particular features, we propose a two-cloud model to fit the observed line profiles. The lower cloud can account for the wing emission, and the upper cloud is mainly responsible for the absorption at line center. After choosing carefully the free parameters, we get satisfactory fitting results. As expected, the lower cloud shows an increase of the source function, corresponding to a temperature increase of 400-1000 K in EBs relative to the quiet Sun. This is consistent with previous results deduced from semi-empirical models and confirms that local heating occurs in the lower atmosphere during the appearance of EBs. We also find that the optical depths can increase to some extent in both the lower and upper clouds, which may result from either direct heating in the lower cloud, or illumination by an enhanced radiation on the upper cloud. The velocities derived from this method, however, are different from those obtained using the traditional bisector method, implying that one should be cautious when interpreting this parameter. The two-cloud model can thus be used as an efficient method to deduce the basic physical parameters of EBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE LEECH EXOPLANET IMAGING SURVEY: CHARACTERIZATION OF THE COLDEST DIRECTLY IMAGED EXOPLANET, GJ 504 b, AND EVIDENCE FOR SUPERSTELLAR METALLICITY

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly imaged exoplanets were all L type. Recently, Kuzuhara et al. announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ∼500 K temperature that bridges the gap between the first directly imaged planets (∼1000 K) and our own solar system's Jupiter (∼130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 μm), spanning the red end of the broad methane fundamental absorption feature (3.3 μm) as part of the LBTI Exozodi Exoplanet Common Hunt (LEECH) exoplanet imaging survey. By comparing our new photometry and literature photometry with a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well fit by models with the following parameters: T{sub eff} = 544 ± 10 K, g < 600 m s{sup −2}, [M/H] = 0.60 ± 0.12, cloud opacity parameter of f{sub sed} = 2–5, R = 0.96 ± 0.07 R{sub Jup}, and log(L) = −6.13 ± 0.03 L{sub ⊙}, implying a hot start mass of 3–30 M{sub jup} for a conservative age range of 0.1–6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a superstellar metallicity. Since planet formation can create objects with nonstellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

TOWARD PRECISION BLACK HOLE MASSES WITH ALMA: NGC 1332 AS A CASE STUDY IN MOLECULAR DISK DYNAMICS

We present first results from a program of Atacama Large Millimeter/submillimeter Array (ALMA) CO(2–1) observations of circumnuclear gas disks in early-type galaxies. The program was designed with the goal of detecting gas within the gravitational sphere of influence of the central black holes (BHs). In NGC 1332, the 0.″3-resolution ALMA data reveal CO emission from the highly inclined (i≈83{sup ∘}) circumnuclear disk, spatially coincident with the dust disk seen in Hubble Space Telescope images. The disk exhibits a central upturn in maximum line-of-sight velocity, reaching ±500 km s{sup −1} relative to the systemic velocity, consistent with the expected signature of rapid rotation around a supermassive BH. Rotational broadening and beam smearing produce complex and asymmetric line profiles near the disk center. We constructed dynamical models for the rotating disk and fitted the modeled CO line profiles directly to the ALMA data cube. Degeneracy between rotation and turbulent velocity dispersion in the inner disk precludes the derivation of strong constraints on the BH mass, but model fits allowing for a plausible range in the magnitude of the turbulent dispersion imply a central mass in the range of ∼(4–8) × 10{sup 8} M{sub ⊙}. We argue that gas-kinematic observations resolving the BH’s projected radius of influence along the disk’s minor axis will have the capability to yield BH mass measurements that are largely insensitive to systematic uncertainties in turbulence or in the stellar mass profile. For highly inclined disks, this is a much more stringent requirement than the usual sphere-of-influence criterion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Petascale Hierarchical Modeling VIA Parallel Execution

The research allows more effective model building. By allowing researchers to fit complex models to large datasets in a scalable manner, our algorithms and software enable more effective scientific research. In the new area of “big data,” it is often necessary to fit “big models” to adjust for systematic differences between sample and population. For this task, scalable and efficient model-fitting tools are needed, and these have been achieved with our new Hamiltonian Monte Carlo algorithm, the no-U-turn sampler, and our new C++ program, Stan. In layman’s terms, our research enables researchers to create improved mathematical modes for large and complex systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2015

Atmospheric Turbulence Modeling for Aero Vehicles: Fractional Order Fits

Atmospheric turbulence models are necessary for the design of both inlet/engine and flight controls, as well as for studying coupling between the propulsion and the vehicle structural dynamics for supersonic vehicles. Models based on the Kolmogorov spectrum have been previously utilized to model atmospheric turbulence. In this paper, a more accurate model is developed in its representative fractional order form, typical of atmospheric disturbances. This is accomplished by first scaling the Kolmogorov spectral to convert them into finite energy von Karman forms and then by deriving an explicit fractional circuit-filter type analog for this model. This circuit model is utilized to develop a generalized formulation in frequency domain to approximate the fractional order with the products of first order transfer functions, which enables accurate time domain simulations. The objective of this work is as follows. Given the parameters describing the conditions of atmospheric disturbances, and utilizing the derived formulations, directly compute the transfer function poles and zeros describing these disturbances for acoustic velocity, temperature, pressure, and density. Time domain simulations of representative atmospheric turbulence can then be developed by utilizing these computed transfer functions together with the disturbance frequencies of interest.

Atmospheric Turbulence Modeling↗
NASA NTRS2002

The conical fit approach to modeling ionospheric total electron content

The Global Positioning System (GPS) can be used to measure the integrated electron density along raypaths between satellites and receivers. Such measurements may, in turn, be used to construct regional and global maps of the ionospheric total electron content (TEC). Maps are generated by fitting measurements to an assumed ionospheric model.

ionosphere WAAS GPS TEC↗
NASA NTRS1989

Aeroelastic modeling for the FIT (Functional Integration Technology) team F/A-18 simulation

As part of Langley Research Center's commitment to developing multidisciplinary integration methods to improve aerospace systems, the Functional Integration Technology (FIT) team was established to perform dynamics integration research using an existing aircraft configuration, the F/A-18. An essential part of this effort has been the development of a comprehensive simulation modeling capability that includes structural, control, and propulsion dynamics as well as steady and unsteady aerodynamics. The structural and unsteady aerodynamics contributions come from an aeroelastic mode. Some details of the aeroelastic modeling done for the Functional Integration Technology (FIT) team research are presented. Particular attention is given to work done in the area of correction factors to unsteady aerodynamics data.

Zeiler, Thomas A.↗