Search NASA⌕ Search

SEARCH · Search NASA

Results for “Exponential integrator”

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 163 records · Page 9

Attitude/attitude-rate estimation from GPS differential phase measurements using integrated-rate parameters

A sequential filtering algorithm is presented for attitude and attitude-rate estimation from Global Positioning System (GPS) differential carrier phase measurements. A third-order, minimal-parameter method for solving the attitude matrix kinematic equation is used to parameterize the filter's state, which renders the resulting estimator computationally efficient. Borrowing from tracking theory concepts, the angular acceleration is modeled as an exponentially autocorrelated stochastic process, thus avoiding the use of the uncertain spacecraft dynamic model. The new formulation facilitates the use of aiding vector observations in a unified filtering algorithm, which can enhance the method's robustness and accuracy. Numerical examples are used to demonstrate the performance of the method.

Oshman, Yaakov↗

CREKID: A computer code for transient, gas-phase combustion of kinetics

A new algorithm was developed for fast, automatic integration of chemical kinetic rate equations describing homogeneous, gas-phase combustion at constant pressure. Particular attention is paid to the distinguishing physical and computational characteristics of the induction, heat-release and equilibration regimes. The two-part predictor-corrector algorithm, based on an exponentially-fitted trapezoidal rule, includes filtering of ill-posed initial conditions, automatic selection of Newton-Jacobi or Newton iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm was found to compare favorably with LSODE on two representative test problems drawn from combustion kinetics.

Pratt, D. T.↗

The cosmic-ray source composition

The cosmic-ray source composition is shown to be related to the observed cosmic-ray abundances by a set of linear equations which are valid for either ad hoc path length distributions or quantitative solutions of galactic propagation models. In the limit of no ionization energy loss, this method is exact for all path length distributions and can replace more detailed numerical integrations; it provides a simple way to calculate directly the source composition implied by a particular set of cross sections, observed abundances, and propagation model. This method is used to determine the source abundances of the galactic cosmic rays in a model-independent manner. Compared to a simple exponential, the path length distribution required by observations is deficient in the contributions from short (not greater than 0.5 g/sq cm) paths at energies of a few GeV per nucleon.

Margolis, S. H.↗

Simultaneous Retrieval of Surface Roughness Parameters for Bare Soils From Combined Active–Passive Microwave SMAP Observations

An active–passive microwave retrieval algorithm for simultaneous determination of soil surface roughness parameters [vertical root-mean-square (RMS) height ( s ) and horizontal correlation length ( l )] is presented for bare soils. The algorithm is based on active–passive microwave covariation, including the improved Integral Equation Method (I(exp 2) EM), and is tested with global soil moisture active passive (SMAP) observations. The estimated retrieval results for s and l are overall consistent with values in the literature, indicating the validity of the proposed algorithm. Sensitivity analyses showed that the developed roughness retrieval algorithm is independent of permittivity for E(sub s) > 10 [-]. Furthermore, the physical model basis of this approach (I(exp 2) EM) allows the application of different autocorrelation functions (ACF), such as Gaussian and exponential ACFs. Global roughness retrieval results confirm bare areas in deserts such as Sahara or Gobi. However, the type of ACF used within roughness parameter estimation is important. Retrieval results for the Gaussian ACF describe a rougher surface than retrieval results for the exponential ACF. No correlations were found between roughness results and the amount of precipitation or the soil texture, which could be due to the coarse spatial resolution of the SMAP data. The extension of this approach to vegetated soils is planned as an add-on study.

Correlation length↗

Functional methods for waves in random media

Some basic ideas in functional methods for waves in random media are illustrated through a simple random differential equation. These methods are then generalized to solve certain random parabolic equations via an exponential representation given by the Feynman-Kac formula. It is shown that these functional methods are applicable to a number of problems in random wave propagation. They include the forward-scattering approximation in Gaussian white-noise media; the solution of the optical beam propagation problem by a phase-integral method; the high-frequency scattering by bounded random media, and a derivation of approximate moment equations from the functional integral representation.

Chow, P. L.↗

Functional methods for waves in random media

Some basic ideas in functional methods for waves in random media are illustrated through a simple random differential equation. These methods are then generalized to solve certain random parabolic equations via an exponential representation given by the Feynman-Kac formula. It is shown that these functional methods are applicable to a number of problems in random wave propagation. They include the forward-scattering approximation in Gaussian white-noise media; the solution of the optical beam propagation problem by a phase-integral method; the high-frequency scattering by bounded random media; and a derivation of approximate moment equations from the functional integral representation.

Chow, P. L.↗

SMM observation of a cosmic gamma-ray burst from 20 keV to 100 MeV

The Solar Maximum Mission gamma-ray spectrometer has detected an intense gamma-ray burst that occurred on August 5, 1984. The burst originated from a source in the constellation Hydra and lasted about 45 s. Its integral fluence at 20 keV was 0.003 erg/sq cm. Spectral evolution similar to other bursts detected by SMM was observed. The overall shape of the spectrum from 20 keV to 100 MeV, on timescales as short as 2 s, is relatively constant. This shape can be fitted by the sum of an exponential-type function and a power law. There is no evidence for narrow or broadened emission lines.

Share, G. H.↗

Evidence for fast dynamo action in a chaotic web

The evolution of a magnetic field in a chaotic web is studied. The model flow possessing the web is closely related to the nearly integrable ABC flow with A = B and C much less than 1. The magnetic diffusivity is taken to be zero and the field is followed using the Cauchy solution. It is found that the flow folds the magnetic field constructively, in the sense that the average magnetic field in a chaotic region grows exponentially in time. This is suggestive of fast dynamo action, although the effect of diffusion of the strong streamwise magnetic field remains to be assessed.

Gilbert, A. D.↗

Bayesian estimation of life parameters in the Weibull distribution.

Development of a Bayesian analysis of the scale and shape parameters in the Weibull distribution and the corresponding reliability function with respect to the usual life-testing procedures. For the scale parameter theta, Bayesian estimates of theta and reliability are obtained for the uniform, exponential, and inverted gamma prior probability densities. Bhattacharya's results (1967) for the one-parameter exponential life-testing distribution are reduced to a special case of these results. A fully Bayesian analysis of both the scale and shape parameters is developed by assuming independent prior distributions; since in the latter case, analytical tractability is not possible, Bayesian estimates are obtained through a conjunction of Monte Carlo simulation and numerical-integration techniques. In both cases, a computer simulation is carried out, and a comparison is made between the Bayesian and the corresponding minimum-variance unbiased, or maximum likelihood, estimates. As expected, the Bayesian estimates are superior.

Canavos, G. C.↗

Control logic to track the outputs of a command generator or randomly forced target

A procedure is presented for synthesizing time-invariant control logic to cause the outputs of a linear plant to track the outputs of an unforced (or randomly forced) linear dynamic system. The control logic uses feed-forward of the reference system state variables and feedback of the plant state variables. The feed-forward gains are obtained from the solution of a linear algebraic matrix equation of the Liapunov type. The feedback gains are the usual regulator gains, determined to stabilize (or augment the stability of) the plant, possibly including integral control. The method is applied here to the design of control logic for a second-order servomechanism to follow a linearly increasing (ramp) signal, an unstable third-order system with two controls to track two separate ramp signals, and a sixth-order system with two controls to track a constant signal and an exponentially decreasing signal (aircraft landing-flare or glide-slope-capture with constant velocity).

Trankle, T. L.↗

Variation in light intensity with height and time from subsequent lightning return strokes

Photographic measurements of relative light intensity as a function of height and time have been conducted for seven return strokes in two lightning flashes at 7.8 and 8.7 km ranges, using film which possesses an approximately constant spectral response in the 300-670 nm range. The amplitude of the initial light peak is noted to decrease exponentially with height, with a decay constant of 0.6-0.8 km. The logarithm of the peak light intensity near the ground is found to be approximately proportional to the initial peak electric field intensity, implying that the current decrease with height may be much slower than the light decrease. Absolute light intensity is presently estimated through the integration of the photographic signals from individual channel segments, in order to simulate the calibrated, all-sky photoelectric data of Guo and Krider (1982).

Jordan, D. M.↗

A Diagnostic Analysis of the Kennedy Space Center LDAR Network: Data Characteristics - 1

An analytical framework is developed in which to analyze climatological VHF (66 MHz) radiation measurements taken by the Kennedy Space Center Lightning Detection and Ranging (LDAR) network. A 19-month noise-filtered sample of LDAR observations is examined using this framework. It is found that the climatological VHF source density as observed by LDAR falls off approximately 10 dB every 71 km of ground range away from the network centroid (a 31 km e-folding scale). From this framework it is inferred that the underlying source distribution is likely inverse exponential in amplitude, with a log-slope of 126-147 mW(sup -1/2). The underlying vertical distribution of VHF sources is approximately normally distributed with a mean altitude of 9 km and a standard deviation of 2.7 km; this implies that the loss of below-horizon sources has a negligible effect on column-integrated source densities within 200 km ground range. At medium to far ranges, location errors are primarily radial and have a slightly asymmetric distribution whose first moment increases as range squared. Error moments estimated from observed lightning are significantly higher than those from aircraft-based signal generator or analytic solution estimates, suggesting that timing errors arising from poor signal identification and discrimination may dominate over timing errors arising from nominal sensor resolution. The VHF source properties of individual LDAR-observed flashes are computed and an analytic expression for flash detection efficiency vs. range is derived. This reveals nearly constant flash detection efficiency to 80-94 km range from the network centroid.

Boccippio, D. J.↗

Exploring Alternate Parameterizations for Snowfall with Validation from Satellite and Terrestrial Radars

Increases in computational resources have allowed operational forecast centers to pursue experimental, high resolution simulations that resolve the microphysical characteristics of clouds and precipitation. These experiments are motivated by a desire to improve the representation of weather and climate, but will also benefit current and future satellite campaigns, which often use forecast model output to guide the retrieval process. Aircraft, surface and radar data from the Canadian CloudSat/CALIPSO Validation Project are used to check the validity of size distribution and density characteristics for snowfall simulated by the NASA Goddard six-class, single-moment bulk water microphysics scheme, currently available within the Weather Research and Forecast (WRF) Model. Widespread snowfall developed across the region on January 22, 2007, forced by the passing of a midlatitude cyclone, and was observed by the dual-polarimetric, C-band radar King City, Ontario, as well as the NASA 94 GHz CloudSat Cloud Profiling Radar. Combined, these data sets provide key metrics for validating model output: estimates of size distribution parameters fit to the inverse-exponential equations prescribed within the model, bulk density and crystal habit characteristics sampled by the aircraft, and representation of size characteristics as inferred by the radar reflectivity at C- and W-band. Specified constants for distribution intercept and density differ significantly from observations throughout much of the cloud depth. Alternate parameterizations are explored, using column-integrated values of vapor excess to avoid problems encountered with temperature-based parameterizations in an environment where inversions and isothermal layers are present. Simulation of CloudSat reflectivity is performed by adopting the discrete-dipole parameterizations and databases provided in literature, and demonstrate an improved capability in simulating radar reflectivity at W-band versus Mie scattering assumptions.

Molthan, Andrew L.↗

Exploring Alternative Parameterizations for Snowfall with Validation from Satellite and Terrestrial Radars

Increases in computational resources have allowed operational forecast centers to pursue experimental, high resolution simulations that resolve the microphysical characteristics of clouds and precipitation. These experiments are motivated by a desire to improve the representation of weather and climate, but will also benefit current and future satellite campaigns, which often use forecast model output to guide the retrieval process. The combination of reliable cloud microphysics and radar reflectivity may constrain radiative transfer models used in satellite simulators during future missions, including EarthCARE and the NASA Global Precipitation Measurement. Aircraft, surface and radar data from the Canadian CloudSat/CALIPSO Validation Project are used to check the validity of size distribution and density characteristics for snowfall simulated by the NASA Goddard six-class, single moment bulk water microphysics scheme, currently available within the Weather Research and Forecast (WRF) Model. Widespread snowfall developed across the region on January 22, 2007, forced by the passing of a mid latitude cyclone, and was observed by the dual-polarimetric, C-band radar King City, Ontario, as well as the NASA 94 GHz CloudSat Cloud Profiling Radar. Combined, these data sets provide key metrics for validating model output: estimates of size distribution parameters fit to the inverse-exponential equations prescribed within the model, bulk density and crystal habit characteristics sampled by the aircraft, and representation of size characteristics as inferred by the radar reflectivity at C- and W-band. Specified constants for distribution intercept and density differ significantly from observations throughout much of the cloud depth. Alternate parameterizations are explored, using column-integrated values of vapor excess to avoid problems encountered with temperature-based parameterizations in an environment where inversions and isothermal layers are present. Simulation of CloudSat reflectivity is performed by adopting the discrete-dipole parameterizations and databases provided in literature, and demonstrate an improved capability in simulating radar reflectivity at W-band versus Mie scattering assumptions.

Molthan, Andrew L.↗

Simulation Propulsion System and Trajectory Optimization

A number of new aircraft concepts have recently been proposed which tightly couple the propulsion system design and operation with the overall vehicle design and performance characteristics. These concepts include propulsion technology such as boundary layer ingestion, hybrid electric propulsion systems, distributed propulsion systems and variable cycle engines. Initial studies examining these concepts have typically used a traditional decoupled approach to aircraft design where the aerodynamics and propulsion designs are done a-priori and tabular data is used to provide inexpensive look ups to the trajectory ana-ysis. However the cost of generating the tabular data begins to grow exponentially when newer aircraft concepts require consideration of additional operational parameters such as multiple throttle settings, angle-of-attack effects on the propulsion system, or propulsion throttle setting effects on aerodynamics. This paper proposes a new modeling approach that eliminated the need to generate tabular data, instead allowing an expensive propulsion or aerodynamic analysis to be directly integrated into the trajectory analysis model and the entire design problem optimized in a fully coupled manner. The new method is demonstrated by implementing a canonical optimal control problem, the F-4 minimum time-to-climb trajectory optimization using three relatively new analysis tools: Open M-DAO, PyCycle and Pointer. Pycycle and Pointer both provide analytic derivatives and Open MDAO enables the two tools to be combined into a coupled model that can be run in an efficient parallel manner that helps to cost the increased cost of the more expensive propulsion analysis. Results generated with this model serve as a validation of the tightly coupled design method and guide future studies to examine aircraft concepts with more complex operational dependencies for the aerodynamic and propulsion models.

Optimization↗

Numerical solutions of Navier-Stokes equations for the structure of a trailing vortex

The structure and decay of a trailing vortex were analyzed during the numerical solutions of the full Navier-Stokes equations. Unsteady forms of the governing equations were recast in terms of circulation, vorticity, and stream function as dependent variables, and a second upwind finite difference scheme was used to integrate them with prescribed initial and boundary conditions. The boundary conditions at the outer edge and at the outflow section of the trailing vortex were considered. Different models of the flow were postulated, and solutions were obtained describing the development of the flow as integration proceeds in time. A parametric study was undertaken with a view to understanding the various phenomena that may possibly occur in the trailing vortex. Using the Hoffman and Joubert law of circulation at the inflow section, the results of this investigation were compared with experimental data for a Convair 990 wind model and a rectangular wing. With an exponentially decaying law of circulation at the inflow section and an adverse pressure gradient at the outer edge of the trailing vortex, solutions depict vortex bursting through the sudden expansion of the core and/or through the stagnation and consequent reversal of the flow on the axis. It was found that this bursting takes place at lower values of the swirl ratio as the Reynolds number increases.

Jain, A. C.↗

Suzaku observation of IGR J16318-4848

We report on the first Suzaku observation of IGR J16318-4848, the most extreme example of a new group of highly absorbed X-ray binaries that have recently been discovered by the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The Suzaku observation was carried out between 2006 August 14 and 17, with a net exposure time of 97 ks. The average X-ray spectrum of the source can be well described (chi-square (sub red)= 0.99) with a continuum model typical for neutron stars i.e., a strongly absorbed power law continuum with a photon index of 0.676(42) and an exponential cutoff at 20.5(6) keY. The absorbing column is N(sub H) = 1.95(3) X 10(exp 24)/square cm. Consistent with earlier work, strong fluorescent emission lines of Fe K-alpha, Fe K-beta, and Ni K-alpha are observed. Despite the large N(sub H), no Compton shoulder is seen in the lines, arguing for a non-spherical and inhomogeneous absorber. Seen at an average 5-60 keV absorbed flux of 3.4 x 10(exp -10) erg/square cm/second, the source exhibits significant variability on timescales of hours.

Barragan, Laura↗

Bioreactor Development for CO2-Based In Situ Resource Utilization Manufacturing

Sustainable long-duration manned missions on both the Moon and Mars will require in situ resource utilization (ISRU). Carbon dioxide (CO2) has great potential as a harvestable resource, making up 95% of the atmosphere on Mars and being produced as respiratory waste in spacecraft and future planetary habitats. Through ISRU, biomanufacturing has the capability to produce a near limitless array of products from local space resources. Here, a CO2-based ISRU recombinant protein producing bioreactor and associated biomanufacturing organisms were designed to produce a highly stable carbonic anhydrase (CA). Candidate organisms were selected by growth characterization on acetate and formic acid, carbon substrates that are synthesized via electrochemical conversion of CO2. To improve growth on the CO2 producing substrate formic acid and for direct integration of ISRU CO2, a synthetic Calvin-Benson-Bassam cycle was designed for use in Escherichia coli. Multiplex genetic modification in E. coli was facilitated by a tailored CRISPR/Cas9 and λ red recombineering two-vector system. For expression of CA, a blue light regulated T7 promoter was employed for dynamic and small molecule free induction. Efficient bioproduction through a fed-batch exponential feeding strategy was determined via mass balance analysis from ISRU substrates to biomass and CA yield. Flux balance analysis was used to model ISRU substrate metabolism and metabolic pathway engineering in candidate organisms under cultivation strategy conditions for both metabolism reconstruction and pathway design optimization. Finally, a small-scale, disposable bag bioreactor for use in the NASA Bioculture System infrastructure was designed to enable CO2-based CA biomanufacturing in reduced-gravity environments.

biomanufacturing↗