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,279 records · Page 71

Design project: LONGBOW supersonic interceptor

A recent white paper entitled 'From the Sea' has spotlighted the need for Naval Aviation to provide overland support to joint operations. The base for this support, the Aircraft Carrier (CVN), will frequently be unable to operate within close range of the battleground because of littoral land-based air and subsurface threats. A high speed, long range, carrier capable aircraft would allow the CVN to provide timely support to distant battleground operations. Such an aircraft, operating as a Deck-Launched Interceptor (DLI), would also be an excellent counter to Next Generation Russian Naval Aviation (NGRNA) threats consisting of supersonic bombers, such as the Backfire, equipped with the next generation of high-speed, long-range missiles. Additionally, it would serve as an excellent high speed Reconnaissance airplane, capable of providing Battle Force commanders with timely, accurate pre-mission targeting information and post-mission Bomb Damage Assessment (BDA). Recent advances in computational hypersonic airflow modeling has produced a method of defining aircraft shapes that fit a conical shock flow model to maximize the efficiency of the vehicle. This 'Waverider' concept provides one means of achieving long ranges at high speeds. A Request for Proposal (RFP) was issued by Professor Conrad Newberry that contained design requirements for an aircraft to accomplish the above stated missions, utilizing Waverider technology.

Stoney, Robert↗

Modelling atmospheric scatterers using spacecraft observations

Voyager images of Triton indicate considerable spatial variability in the concentration of at least two different scattering components in the atmosphere. Data from high phase angle limb scans were fit to Mie scattering models to derive mean particle sizes, number densities, and vertical extent for both types of scattering material at ten different locations between 15 deg S and 70 deg S. These fits reveal a thin haze at latitudes equatorward of 25-30 deg S. The imaging data can be fit reasonably well by both conservatively scattering and absorbing hazes with particle sizes near 0.18 micron and optical depths of order 0.001-0.01. Rayleigh scattering haze fits the imaging data somewhat less well, and can be totally ruled out by combining the imaging and UVS measurements. At high southern latitudes, Triton displays clouds below an altitude of approximately 8 km, as well as the haze at higher altitudes. The clouds have particle sizes which may range from 0.7-2.0 microns, or may be near 0.25 micron. The atmospheric optical depth poleward of 30 deg S must be generally greater than 0.1, but need not be more than 0.3. Horizontal inhomogeneities are quite noticeable, especially at longitudes east of (i.e., higher than) 180 deg.

Rages, Kathy A.↗

From Coexpression to Coregulation: An Approach to Inferring Transcriptional Regulation Among Gene Classes from Large-Scale Expression Data

We provide preliminary evidence that existing algorithms for inferring small-scale gene regulation networks from gene expression data can be adapted to large-scale gene expression data coming from hybridization microarrays. The essential steps are (I) clustering many genes by their expression time-course data into a minimal set of clusters of co-expressed genes, (2) theoretically modeling the various conditions under which the time-courses are measured using a continuous-time analog recurrent neural network for the cluster mean time-courses, (3) fitting such a regulatory model to the cluster mean time courses by simulated annealing with weight decay, and (4) analysing several such fits for commonalities in the circuit parameter sets including the connection matrices. This procedure can be used to assess the adequacy of existing and future gene expression time-course data sets for determining transcriptional regulatory relationships such as coregulation.

Mjolsness, Eric↗

Prototype Global Analysis of LISA Data with Multiple Source Types

The novel data analysis challenges posed by the Laser Interferometer Space Antenna (LISA) arise from the overwhelmingly large number of astrophysical sources in the measurement band and the density with which they are found in the data. Robust detection and characterization of the numerous gravitational wave sources in LISA data can not be done sequentially, but rather through a simultaneous global fit of a data model containing the full suite of astrophysical and instrumental features present in the data. While previous analyses have focused on individual source types in isolation, here we present the first demonstration of a LISA global fit analysis containing combined astrophysical populations. The prototype pipeline uses a blocked Metropolis Hastings algorithm to alternatingly fit to a population of ultra compact galactic binaries, known “verification binaries” already identified by electromagnetic observations, a population of massive black hole mergers, and an instrument noise model. The Global LISA Analysis Software Suite (GLASS) is assembled from independently developed samplers for the different model components. The modular design enables flexibility to future development by defining standard interfaces for adding new, or updating additional, components to the global fit without being overly prescriptive for how those modules must be internally designed. The GLASS pipeline is demonstrated on data simulated for the LISA Data Challenge 2b. Results of the analysis and a road-map for continued development are described in detail.

Gravitational Waves↗

Investigation of the Saturn dust environment from the analysis of energetic charged particle measurements

In order to assist the Cassini project in evaluating risks of collisions with particulate matter in the rotational equatorial plane, available Pioneer 11 and Voyager energetic charged particle data were reinvestigated to constrain the column mass density of absorbing material within several radial ranges. Within the orbit of Mimas, CRAND proton phase space densities maximize near 2.67 R sub s and exhibit secondary maximum at 2.43 R sub s. From the condition that sources must exceed losses near these maxima and using available theoretical models for CRAND proton production rates, upper bounds on the column mass density at these two radial locations are calculated. Detailed fits of radial diffusion models to the observed flux maxima yield somewhat more restrictive upper limits. The upper limits compare to a lower limit on the column mass density, estimated from previous model calculations by Van Allen. Aside from continuous rings, longitudinally limited, low optical depth clouds of particulates may exist in orbit with several of the inner satellites including Mimas and Enceladus. A brief review of Voyager energetic particle microsignatures that suggest the presence of material co-orbiting with these two satellites is presented. Finally, Pioneer 11 and Voyager measurements of low energy electron fluxes exhibit minima near the location of the tenuous E Ring centered on approx. 4 R sub s. Pioneer 11 pitch angle distributions appear to support the possibility that direct absorption by Ring E particulates produced the observed flux decreases.

Hood, L. L.↗

Low Reynolds number two-equation modeling of turbulent flows

A k-epsilon model that accounts for viscous and wall effects is presented. The proposed formulation does not contain the local wall distance thereby making very simple the application to complex geometries. The formulation is based on an existing k-epsilon model that proved to fit very well with the results of direct numerical simulation. The new form is compared with nine different two-equation models and with direct numerical simulation for a fully developed channel flow at Re = 3300. The simple flow configuration allows a comparison free from numerical inaccuracies. The computed results prove that few of the considered forms exhibit a satisfactory agreement with the channel flow data. The model shows an improvement with respect to the existing formulations.

Michelassi, V.↗

Accelerated Aging of Lead-Free Propellant

Following higher than expected 2-NDPA depletion rates in a lead-free doublebase formulation (RPD-422), an accelerated aging study was conducted to verify the depletion rates. A test plan was prepared to compare the aging characteristics of lead-free propellant and NOSIH-AA2. The study was also designed to determine which lead-free ballistic modifiers accelerated 2-NDPA depletion. The increased depletion rate occurred in propellants containing monobasic copper salicylate. Four lead-free propellants were then formulated to improved aging characteristics over previous lead-free propellant formulations. The new formulations reduced or replaced the monobasic copper salicylate. The new formulations had improved aging characteristics. Their burn rates, however, were unacceptable for use in a 2.75 inch rocket. To compare aging characteristics, stabilizer depletion rates of RPD-422, AA2, M28, and RLC 470/6A were measured or taken from the literature. The data were fit to a kinetic model. The model contained first and zero order terms which allowed the stabilizer concentration to go to zero. In the model, only the concentration of the primary stabilizer was considered. Derivatives beyond the first nitrated or nitroso derivative of 2-NPDA were not considered. The rate constants were fit to the Arrhenius equation and extrapolated to lower temperatures. The time to complete stabilizer depletion was estimated using the kinetic model. The four propellants were compared and the RPD-422 depleted faster at 45 C than both A22 and M28. These types of predictions depend on the validity of the model and on confidence in the Arrhenius relationship holding at lower temperatures. At 45 C, the zero order portion of the model dominates the depletion rate.

Furrow, Keith W.↗

Constraints on the Neutron Star and Inner Accretion Flow in Serpens X-1 Using Nustar

We report on an observation of the neutron star low-mass X-ray binary Serpens X-1, made with NuSTAR. The extraordinary sensitivity afforded by NuSTAR facilitated the detection of a clear, robust, relativistic Fe K emission line from the inner disk. A relativistic profile is required over a single Gaussian line from any charge state of Fe at the 5 sigma level of confidence, and any two Gaussians of equal width at the same confidence. The Compton back-scattering "hump" peaking in the 10-20 keV band is detected for the first time in a neutron star X-ray binary. Fits with relativistically blurred disk reflection models suggest that the disk likely extends close to the innermost stable circular orbit (ISCO) or stellar surface. The best-fit blurred reflection models constrain the gravitational redshift from the stellar surface to be ZnS (is) greater than 0.16. The data are broadly compatible with the disk extending to the ISCO; in that case,ZnS(is) greater than 0.22 and RNS (is) less than12.6 km (assuming MnS = 1.4 solar mass and a = 0, where a = cJ/GM2). If the star is as large or larger than its ISCO, or if the effective reflecting disk leaks across the ISCO to the surface, the redshift constraints become measurements. We discuss our results in the context of efforts to measure fundamental properties of neutron stars, and models for accretion onto compact objects.

accretion↗

Effects of curvilinear motion in large-amplitude bending of C3

The geometry of the bending of a linear triatomic molecule is analyzed, and an expression for the average rotational constant is derived. A harmonic oscillator model of C3 is fitted to the observed rotational constant within 0.6%. The bond distance between atoms at zero bending is 1.287 A according to this model; this is noticeably larger than the average internuclear distance of 1.277 A for the vibrational ground state. The first order perturbation solutions for the vibrational energy levels, taking into account the effect of a quartic perturbation potential, closely match observed levels. For a square well potential model of C3, the effect of curvilinear motion in bending is similar to that found for the harmonic oscillator model, though the decreases in energy are about twice as large. In both models, the average energy decrease is relatively constant at approximately 10% over a wide range of vibrational quantum number.

Hansen, C. F.↗

Inhomogeneous chemical evolution of the Galactic disk

We present analytical models for inhomogeneous chemical evolution (ICE) of systems in which the star formation history resembles a series of bursts, localized in space and/or time, with intermittent periods of remixing. The additional parameter of this model is the metallicity increment of bursting subsystems, but this parameter is constrained by the spread in the age-metallicity relation. We apply this model to the solar annulus in the Galactic disk and show that ICE models yield an improved fit to the observed shape of the stellar abundance distribution function (ADF). The G-dwarf problem can be alleviated with ICE models, but infall of metal poor gas and/or some preenrichment of the disk during the epoch of protogalactic evolution is still required to explain the paucity of low-metallicity dwarfs. ICE models also suggest an explanation of the reduced frequency of metal-rich G-dwarfs relative to the predictions of the simple model. It does not seem likely that chemical evolution of the solar annulus proceeded in a medium that was well-mixed at all times.

Malinie, Guy↗

The Mathematics of Navigating the Solar System

In navigating spacecraft throughout the solar system, the space navigator relies on three academic disciplines - optimization, estimation, and control - that work on mathematical models of the real world. Thus, the navigator determines the flight path that will consume propellant and other resources in an efficient manner, determines where the craft is and predicts where it will go, and transfers it onto the optimal trajectory that meets operational and mission constraints. Mission requirements, for example, demand that observational measurements be made with sufficient precision that relativity must be modeled in collecting and fitting (the estimation process) the data, and propagating the trajectory. Thousands of parameters are now determined in near real-time to model the gravitational forces acting on a spacecraft in the vicinity of an irregularly shaped body. Completing these tasks requires mathematical models, analyses, and processing techniques. Newton, Gauss, Lambert, Legendre, and others are justly famous for their contributions to the mathematics of these tasks. More recently, graduate students participated in research to update the gravity model of the Saturnian system, including higher order gravity harmonics, tidal effects, and the influence of the rings. This investigation was conducted for the Cassini project to incorporate new trajectory modeling features in the navigation software. The resulting trajectory model will be used in navigating the 4-year tour of the Saturnian satellites. Also, undergraduate students are determining the ephemerides (locations versus time) of asteroids that will be used as reference objects in navigating the New Millennium's Deep Space 1 spacecraft autonomously.

Hintz, Gerald↗

AXAF VETA-I mirror encircled energy measurements and data reduction

The AXAF VETA-I mirror encircled energy was measured with a series of apertures and two flow gas proportional counters at five X-ray energies ranging from 0.28 to 2.3 keV. The proportional counter has a thin plastic window with an opaque wire mesh supporting grid. Depending on the counter position, this mesh can cause the X-ray transmission to vary as much as +/-9 percent, which directly translates into an error in the encircled energy. In order to correct this wire mesh effect, window scan measurements were made, in which the counter was scanned in both horizontal (Y) and vertical (Z) directions with the aperture fixed. Post VETA measurement of the VXDS setup were made to determine the exact geometry and position of the mesh grid. Computer models of the window mesh were developed to simulate the X-ray transmission based on this measurement. The window scan data were fitted to such mesh models and corrections were made. After this study, the mesh effect was well understood and the final results of the encircled energy were obtained with an uncertainty of less than 0.8 percent.

Zhao, Ping↗

AXAF VETA-I mirror encircled energy measurements and data reduction

The AXAF VETA-I mirror encircled energy was measured with a series of apertures and two flow gas proportional counters at five X-ray energies ranging from 0.28 to 2.3 keV. The proportional counter has a thin plastic window with an opaque wire mesh supporting grid. Depending on the counter position, this mesh can cause the X-ray transmission to vary as much as +/- 9 percent, which directly translates into an error in the encircled energy. In order to correct this wire mesh effect, window scan measurements were made, in which the counter was scanned in both horizontal (Y) and vertical (Z) directions with the aperture fixed. Post VETA measurement of the VXDS setup were made to determine the exact geometry and position of the mesh grid. Computer models of the window mesh were developed to simulate the X-ray transmission based on this measurement. The window scan data were fitted to such mesh models and corrections were made. After this study, the mesh effect was well understood and the final results of the encircled energy were obtained with an uncertainty of less than 0.8 percent.

Zhao, Ping↗

Interannual Variability of Martian Global Dust Storms: Simulations with a Low-Order Model of the General Circulation

We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good. Interannual variability of GDSs is possible in the IBM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up its a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.

Pankine, A. A.↗

AG Channel Measurement and Modeling Results for Over-Sea Conditions

This report describes results from flight tests conducted in an over-sea environment, for the purpose of characterizing the air-to-ground (AG) channel, for future unmanned aircraft system (UAS) communication system analysis and design. These results are for the first of a set of several flight tests conducted in different ground site (GS) environments. An ultimate aim of all these tests is the development of models for the AG channel that can be used in communication system evaluation. In this report we provide measured results for propagation path loss, root-mean square delay spread (RMS-DS), and the correlation coefficient of the primary received signal components on the four antennas (two antennas for C-band, two for L-band). For path loss, the curved-earth two-ray model provides a reasonable fit to the measured data, altered by several dB at the shortest link distances by aircraft antenna pattern effects. This two-ray model also accounts for the majority of measured RMS-DS results of a few tens of nanoseconds, except for the occasional intermittent reflections from surface objects. These intermittent reflections yield RMS-DS values up to several hundred nanoseconds. For portions of the flight path that were over a harbor area highly populated with boats, the channel was found to be more "continuously dispersive," with RMS-DS reaching approximately 250 ns. A separate model will be developed for this over-harbor setting. The correlation coefficient results are still undergoing analysis; preliminary observations are that correlation between signals on the same-band antennas is generally large (>0.6) for the C-band straight flight paths, whereas for the L-band signals and for the oval-shaped flight paths the correlation is generally small (below 0.4). Inter-band correlations are typically very small, and are well modeled as zero-mean Gaussian in distribution, with a standard deviation less than 0.2. Hence the over-sea channel effects in the two bands can be considered uncorrelated, which will allow for good diversity gains in dual-band systems. We describe initial modeling approaches for the over-sea channel; complete models for this and the over-harbor setting will appear in a subsequent report.

Aircraft Communications↗

A Gaussian Process Enhancement to Linear Parameter Varying Models

Simulation and analysis for modern engineering systems now routinely requires the merging of multiple disciplines, physical-domains, time-scales, and data sets — all at ever increasing levels. These capabilities are especially needed in the domain of Advanced Air Mobility, where rapidly emerging vehicle designs are significantly more complex, while having to be both cost-effective and safe. To meet these engineering challenges, machine learning methods are an attractive option for merging models and data across multiple areas while providing uncertainty quantification and maintaining computational efficiency. This paper examines the use of Gaussian process machine learning to generalize and enhance the commonly used class of quasi-Linear Parameter Varying models for fast full-envelope simulation while also supporting control system design and analysis with model uncertainty. Gaussian process machine learning is selected because it: can fuse multiple data sets, enables an easy trade-off between data fitting and smoothing, provides model uncertainty quantification, scales well with increasing complexity, and does not generally require starting from a large training data set. To demonstrate the benefits of the approach, a robust stability analysis with Gaussian process uncertainty is shown for a NASA reference design of an electric quad-rotor air-taxi concept vehicle with motor parameter uncertainty.

Gaussian Process↗

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗

A physical model of Titan's clouds

A physical model of the formation and growth of aerosols in the atmosphere of Titan has been constructed in light of the observed correlation between variations in Titan's albedo and the sunspot cycle. The model was developed to fit spectral observations of deep methane bands, pressures, temperature distributions, and cloud structure, and is based on a one-dimensional physical-chemical model developed to simulate the earth's stratospheric aerosol layer. Sensitivity tests reveal the model parameters to be relatively insensitive to particle shape but sensitive to particle density, with high particle densities requiring larger aerosol mass production rates to produce compatible clouds. Solution of the aerosol continuity equations for particles of sizes 13 A to about 3 microns indicates the importance of a warm upper atmosphere and a high-altitude mass injection layer, and the production of aerosols at very low aerosol optical depths. Limits are obtained for the chemical production of aerosol mass and the eddy diffusion coefficient, and it is found that an increase in mass input causes a decrease in mean particle size.

Toon, O. B.↗