Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,297 records · Page 72

Models of the quasar population. I - A new luminosity function

A new functional form for the quasar luminosity function is tested using recent observational results for both bright and faint quasar count and redshift distributions. The form is of a fairly general type based on three free parameters and allows for quasars to undergo a combination of luminosity evolution and luminosity-dependent density evolution; an advantage to this approach is that it does not constrain quasars to follow a single type of evolution. Models of pure luminosity evolution or luminosity-dependent density evolution can be constructed, but the apparent magnitude distribution of observed quasars is best fitted by a combination model. The combination model also gives the correct redshift distribution for quasars with redshifts less than three and predicts that quasars brighter than B = 22 provide a 2-10 keV X-ray flux that is equal to 32 percent of the observed X-ray background. However, the model is flawed in that it predicts more high-redshift quasars than are observed.

Heisler, Julia↗

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗

Comparison of distributed acceleration and standard models of cosmic-ray transport

Recent cosmic-ray abundance measurements for elements in the range 3 less than or equal to Z less than or equal to 28 and energies 10 MeV/n less than or equal to E less than or equal to 1 TeV/n have been analyzed with computer transport modeling. About 500 elemental and isotopic measurements have been explored in this analysis. The transport code includes the effects of ionization losses, nuclear spallation reactions (including those of secondaries), all nuclear decay modes, stripping and attachment of electrons, escape from the Galaxy, weak reacceleration and solar modulation. Four models of reacceleration (with several submodels of various reacceleration strengths) were explored. A chi (exp 2) analysis show that the reacceleration models yield at least equally good fits to the data as the standard propagation model. However, with reacceleration, the ad hoc assumptions of the standard model regarding discontinuities in the energy dependence of the mean path length traversed by cosmic rays, and in the momentum spectrum of the cosmic-ray source spectrum are eliminated. Futhermore, the difficulty between rigidity dependent leakage and energy independent anisotropy below energies of 10(exp 14) eV is alleviated.

Letaw, J. R.↗

Fiber Breakage Model for Carbon Composite Stress Rupture Phenomenon: Theoretical Development and Applications

Stress rupture failure of Carbon Composite Overwrapped Pressure Vessels (COPVs) is of serious concern to Science Mission and Constellation programs since there are a number of COPVs on board space vehicles with stored gases under high pressure for long durations of time. It has become customary to establish the reliability of these vessels using the so called classic models. The classical models are based on Weibull statistics fitted to observed stress rupture data. These stochastic models cannot account for any additional damage due to the complex pressure-time histories characteristic of COPVs being supplied for NASA missions. In particular, it is suspected that the effects of proof test could significantly reduce the stress rupture lifetime of COPVs. The focus of this paper is to present an analytical appraisal of a model that incorporates damage due to proof test. The model examined in the current paper is based on physical mechanisms such as micromechanics based load sharing concepts coupled with creep rupture and Weibull statistics. For example, the classic model cannot accommodate for damage due to proof testing which every flight vessel undergoes. The paper compares current model to the classic model with a number of examples. In addition, several applications of the model to current ISS and Constellation program issues are also examined.

Murthy, Pappu L. N.↗

Dealing with Ion LET Uncertainties: An Application of Generalized Linear Models

Although most SEE rate estimation methods presume a fit to SEE cross section vs. LET, fitting SEE data is challenging because the data are not compatible with the assumptions of many common fitting techniques (e.g. linear regression. The difficulty of fitting such data is compounded when the LET of the ion responsible for an SEE is uncertain. We modify a Generalized Linear Model SEE data fitting methodology to accommodate uncertain LET and apply the method to the problem of backside heavy-ion SEE testing to demonstrate the utility of the method, explore the dependence of systematic errors that arise from improper treatment of LET uncertainty and develop guidelines for minimizing such systematic errors when proper treatment is not possible. Additional applications are suggested and assessed for suitability of treatment by the model.

Single-event effects↗

Error enhancement in geomagnetic models derived from scalar data

Models of the main geomagnetic field are generally represented by a scalar potential gamma expanded in a finite number of spherical harmonics. Very accurate observations of F were used, but indications exist that the accuracy of models derived from them is considerably lower. One problem is that F does not always characterize gamma uniquely. It is not clear whether such ambiguity can be encountered in deriving gamma from F in geomagnetic surveys, but there exists a connection, due to the fact that the counterexamples of Backus are related to the dipole field, while the geomagnetic field is dominated by its dipole component. If the models are recovered with a finite error (i.e. they cannot completely fit the data and consequently have a small spurious component), this connection allows the error in certain sequences of harmonic terms in gamma to be enhanced without unduly large effects on the fit of F to the model.

Stern, D. P.↗

What causes the seasonal cycle of stationary waves in the southern stratosphere?

Stationary planetary waves in the southern stratosphere display a characteristic seasonal cycle. Previous research based on a one-dimensional model suggests that this behavior is mainly determined by seasonally varying transmission properties of the atmosphere with respect to wave propagation. The issue is investigated with the help of a hemispheric, linear, quasi-geostrophic model. It reproduces well some of the observed qualitative features and is internally consistent in the sense that its seasonal wave cycle can be explained in terms of varying wave transmission properties of the mean circulation. On the other hand, the model does not yield the observed seasonal cycle. Despite considerable sensitivity to modifications in the basic-state wind and dissipation parameterization, the model could not be reasonably fit to reproduce the observed seasonal cycle. Possible reasons for the model deficiency are put forward. In summary, even though suggestive, the present study is not entirely conclusive about the degree to which the observed cycle is determined by wave transmission properties alone.

Wirth, Volkmar↗

Determination of the Charon/Pluto Mass Ratio from Center-of-Light Astrometry

The Charon/Pluto mass ratio is a fundamental but poorly known parameter of the two-body system. Previous values for the mass ratio have ranged from 0.0837 plus or minus 0.0147 (Null et al., 1993, Astron. J. 105, 2319-2335) to 0.1566 plus or minus 0.0035 (Young et al., 1994, Icarus 108,186-199). We report here a new determination of the Charon/Pluto mass ratio, using five sets of groundbased images taken at four sites in support of Pluto occultation predictions. Unlike the Null et al. and Young et A determinations, where the centers of light for Pluto and Charon could be determined separately, this technique examines the motion of the center of light of the blended Pluto-Charon image. We compute the offsets of the observed center-of-light position of Pluto-Charon from the ephemeris position of the system and fit these offsets to a model of the Pluto-Charon system. The least-squares fits to the five data sets agree within their errors, and the weighted mean mass ratio is 0.117 plus or minus 0.006. The effects of errors in the Charon light fraction, semimajor axis, and ephemeris have been examined and are equal to only a small fraction of the formal error from the fit. This result is intermediate between those of Null et al., and Young et al. and matches a new value of 0.124 plus or minus 0.008 by Null and Owen (1996, Astron. J. 111, 1368-1381). The mass ratio and resulting individual masses and densities of Pluto and Charon are consistent with a collisional origin for the Pluto-Charon system.

Foust, Jeffrey A.↗

Radial Variations in the Io Plasma Torus during the Cassini Era

A radial scan through the midnight sector of the Io plasma torus was made by the Cassini Ultraviolet Imaging Spectrograph on 14 January 2001, shortly after closest approach to Jupiter. From these data, Steffl et al. (2004a) derived electron temperature, plasma composition (ion mixing ratios), and electron column density as a function of radius from L = 6 to 0 as well as the total luminosity. We have advanced our homogeneous model of torus physical chemistry (Delamere and Bagenal, 2003) to include latitudinal and radial variations in a manner similar to the two-dimensional model by Schreier et al. (1998). The model variables include: (1) neutral source rate, (2) radial transport coefficient, (3) the hot electron fraction, (4) hot electron temperature, and (5) the neutral O/S ratio. The radial variation of parameters 1-4 are described by simple power laws, making a total of nine parameters. We have explored the sensitivity of the model results to variations in these parameters and compared the best fit with previous Voyager era models (schreier et al., 1998), galileo data (Crary et al., 1998), and Cassini observations (steffl et al., 2004a). We find that radial variations during the Cassini era are consistent with a neutral source rate of 700-1200 kg/s, an integrated transport time from L = 6 to 9 of 100-200 days, and that the core electron temperature is largely determined by a spatially and temporally varying superthermal electron population.

Delamere, P. A.↗

Challenges Beyond the Standard Model: An Examination and Evaluation of Alternatives with Respect to Unsolved Problems in Energetics and Signaling

The standard model of particle physics is regarded as one of the greatest achievements of modern physics. It has successfully predicted the existence of numerous hadrons; e.g. the Higgs Boson, the W and Z Bosons, B mesons, etc. Moreover it accurately models many phenomena such as Beta and Radioactive decay, pair production, hadron mass splitting, and others. In spite of this, the standard model has many shortfalls. For example it fails to account for or explain the presence of dark matter and dark energy. Many alternative models of the substructure of matter etc. have been proposed. Some of these are examined. In addition the strongly related standard model of cosmology, that is the big bang theory, accounts for all observed celestial phenomena. What it cannot do is satisfactorily justify all astrophysical events. As noted, dark matter and dark energy are not well understood yet they are core elements of the big bang model. Many elements of this standard model of cosmology are challenged regularly, yet somehow the community always fits the problems with the model into another mysterious element. This model and others are examined for consistency and reconciliation. Evidence of particles that move faster than the speed of light has long been sought, but remains elusive. Some of the known space-like observations of physics are examined. Potential consequences of such particles are likewise investigated. Some conclusions are drawn with respect to the framework of tachyons in a rational but non-rigorous sense. Finally, some of my ideas are presented, concerning all of the above material. The framework of theories of everything is examined, considering what such a theory requires for validity. Sociological notes on the structure of research in physics are given.

Coxe, Alexander M.↗

Axisymmetric model of the ionized gas in the Orion Nebula

New ionization and thermal equilibrium models for the ionized gas in the Orion Nebula with an axisymmetric two-dimensional 'blister' geometry/density distribution are presented. The HII region is represented more realistically than in previous models, while the physical detail of the microphysics and radiative transfer of the earlier spherical modeling is maintained. The predicted surface brightnesses are compared with observations for a large set of lines at different positions to determine the best-fitting physical parameters. The model explains the strong singly ionized line emission along the lines of sight near the Trapezium.

Rubin, R. H.↗

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗

Strong-field Gravitational Effects in the Timing and Spectroscopy of X-Ray Binaries

Using the ray-tracing code described, we developed additional features to the hot spot model to explain the broadening of the quasi-periodic oscillation (QPO) peaks as well as the damping of higher frequency harmonics in the power spectrum. We derived a number of analytic results that agree closely with more detailed numerical calculations. Three primary pieces are developed: the superposition of multiple hot spots with finite lifetimes and random phases, a width in the radial distribution of geodesic orbits, and the scattering of photons from the hot spot through a corona around the black hole. The complete model is used to fit the observed power spectra of both type A and type B QPOs seen in XTE 51550-564, giving confidence limits on each of the model parameters. These results were presented in a conference paper and in a refereed publication. We also explored the possibility of using higher-order statistics, particularly the bispectrum and bicoherence, to distinguish between competing QPO models. These methods were applied to the hot spot model and a method was proposed to map out the spacetime metric around Kerr black holes in Ref. 3 below.

Bertschinger, Edmund↗

A Detailed Study of Two Narrow-Line Seyfert 1s

Spectral modeling and variability studies are reported in this paper, including fitting with relativistically blurred ionized-disk models designed to reproduce the radiation spectrum from the inner accretion disk around a black hole. PG 1404+226 is highly variable in flux and has a steep X-ray spectrum that varies over the course of the observation. Its strong observed soft X-ray excess can be well fitted with an ionized-disk model. The parameters derived from this model indicate that the central black hole is rotating, since the inner radius of the accretion disk lies within the last stable orbit for a non-rotating black hole. The harder spectrum of PG 1244426 also shows a soft X-ray excess, but it is better represented by a more conventional two-component phenomenonological model of a power-law and thermal emission.

Brandt, Niel↗

High Tide or Riptide on the Cosmic Shoreline? A Water-Rich Atmosphere or Stellar Contamination for the Warm Super-Earth GJ 486b from JWST Observations

Planets orbiting M-dwarf stars are prime targets in the search for rocky exoplanet atmospheres. The small size of M dwarfs renders their planets exceptional targets for transmission spectroscopy, facilitating atmospheric characterization. However, it remains unknown whether their host stars’ highly variable extreme-UV radiation environments allow atmospheres to persist. With JWST, we have begun to determine whether or not the most favorable rocky worlds orbiting M dwarfs have detectable atmospheres. Here, we present a 2.8—5.2 µm JWST NIRSpec/G395H transmission spectrum of the warm (700 K, 40.3× Earth’s insolation) super-Earth GJ 486b (1.3 R⊕ and 3.0 M⊕). The measured spectrum from our two transits of GJ 486b deviates from a flat line at 2.2 − 3.3σ, based on three independent reductions. Through a combination of forward and retrieval models, we determine that GJ 486b either has a water-rich atmosphere (with the most stringent constraint on the retrieved water abundance of H2O > 10% to 2σ) or the transmission spectrum is contaminated by water present in cool unocculted starspots. We also find that the measured stellar spectrum is best fit by a stellar model with cool starspots and hot faculae. While both retrieval scenarios provide equal quality fits (χ2ν = 1.0) to our NIRSpec/G395H observations, shorter wavelength observations can break this degeneracy and reveal if GJ 486b sustains a water-rich atmosphere.

JWST↗

Exploring the Limits of the Data-Model-Theory Synergy: “Hot” MW Transitions for Rovibrational IR Studies

In order to further improve the accuracy of rovibrational IR line lists generated from the “Best Theory +Reliable High-resolution Experiment” (BTRHE) strategy from 0.01-0.05 cm-1, or 300-1500 MHz, to ~10 MHz, we explore the current limits of the Data-Model-Theory synergy by examining the accuracy and consistency of existing data, then propose that “hot” bands in microwave (MW) spectra is the solution we need for future enhancements. The Ames SO2 J=0-20 rovibrational energy levels computed on the semi-empirically refined Ames-2 potential energy surface (PES) are fit to the Effective Hamiltonian (EH) model regularly used in the experimental infrared (IR) analysis for SO2 isotopologues. In the fitted EH(Ames) model, the rotational constants A/B/C and all 5 quartic centrifugal distortion constants display clear, systematic, and consistent patterns along the vibrational state energy or quanta. Such consistent patterns may facilitate the vibrational assignments for MW hot bands and extract more information from high temperature MW spectra. Some EH(Expt) analyses were carried out with the lowest order Coriolis Coupling term, C1. Their constants should not be directly compared with other EH(Expt) and EH(Ames) results. After excluding them, our  = EH(Ames)- EH(Expt) analyses for 5 isotopologues (626, 636, 646, 628 and 828) indicates some loss of accuracy and consistency starting from vibrational states as low as 22 or 1000 cm-1. Some EH parameters, e.g. K, may have relative deviations as large as 50-100% and totally lose any recognizable patterns. This simply means that current EH(Expt) models do not have the system-wide consistency we need to further refine the EH(Ames) and Ames rovibrational IR line lists. A large part of such defects are probably inherited from the limited precision of experimental line positions, i.e. 1E-3 ~ 1E-4 cm-1, or 3-30 MHz. This is confirmed in a series of truncation tests using the Ames data. Although the EH(Ames) consistency may help identify unreliable rovibrational EH(Expt) parameters, and make reliable predictions for minor isotopologues and unobserved vibrational bands, we believe only the highly accurate “hot” MW transitions can provide real enhancements for EH(Expt) accuracy and consistency. “Hot” MW spectra should play a more significant role in the future synergy of Data, Model, and Theory in the field of rovibrational IR studies.

Xinchuan Huang↗

RXTE Observation of the Tycho Supernova Remnant

SN1006 [4] and Cas A [1, 9] supernova remnants have been shown convincingly to have a hard X-ray power-law continuum. This continuum is thought to be the synchrotron radiation from accelerated electrons of approx. 100 TeV at the shock fronts. Our goal of AO2 RXTE observation is to detect the hard X-ray continuum and to determine the nature of the continuum from Tycho SNR. A detection of a power-law continuum from Tycho SNR can strongly argue for SNRs are the source of cosmic rays with the first order Fermi acceleration as the energizing process. We report the results of our AO2 RXTE 1 x 10(exp 5) sec observation of Tycho SNR. We detect two components of the X-ray spectrum from Tycho SNR both at better than 3 omega confidence. The best two component models are: bremsstrahlung (kT=2.67 +/- 0.13 keV) + bremsstrahlung (kT=7.07 +/- 2.21/1.72 keV) or bremsstrahlung (kT=2.36 +/- 0.21/0.57 keV) + power-law (gamma=2.58 +/- 0.12/0.09 ). This result is an improvement compaxed with the previous most sensitive X-ray measurements by Ginga which shows Tycho's observed X-ray continuum requires a two-component model to yield acceptable fits with the hard component parameters being highly uncertain. Our RXTE measurements constrain all parameter within 3o, ranges. However, we cannot yet distinguish between thermal and nonthermal models for the hard component. In the followings, we describe what we accomplished in the period covered by the grant proposal.

The, Lih-Sin↗

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: from Theory to Fermi Observations

We compute the patterns of gamma-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed gamma-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model gamma-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and gamma-ray emission on the gamma-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the gamma-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

stars↗