Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Estimated capital costs of fish exclusion technologies for hydropower facilities

Hydropower is a reliable source of renewable energy, and its future expansion is likely to be in the form of either smaller new stream development (NSD) projects or powering existing non-powered dams. Thresholds for entrainment risk to fish and the requirements for fish exclusion at hydropower facilities often differ depending on the species involved, the characteristics of the facility, and the goals of stakeholders, but little quantitative information is present within the literature regarding the specific costs of fish exclusion measures. Cost data associated with protection, mitigation, and enhancement (PM&E) measures related to positive barrier screening were identified using keyword searches of an existing environmental mitigation cost data set and manual extraction from regulatory licensing documents available in the Federal Energy Regulatory Commission (FERC) eLibrary. This approach yielded a total of 50 p.m.&E mitigation measures with estimated capital construction costs pertaining to positive barrier screens and represented <10% of the 171 total FERC project dockets available in the data set. These data were highly skewed toward conventional relicensing projects, as <7% were associated with NSD projects. Results indicate highly variable costs are associated with fish screening, with flow-normalized costs one to two orders of magnitude higher for screening with the highest exclusion capability (≤0.09 in. spacing) compared with coarser screening (1–2 in.). These data provide an initial baseline for estimating exclusion costs for hydropower development and may help developers consider options for more fish-friendly generation technologies, though gaps remain relating to a lack of data, particularly for NSD projects.

13 HYDRO ENERGY↗
DOE OSTI2020

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. Here, these observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Influence of uniaxial single-ion anisotropy on the magnetic and thermal properties of Heisenberg antiferromagnets within unified molecular field theory

Here, the influence of uniaxial single-ion anisotropy –DS 2 z on the magnetic and thermal properties of Heisenberg antiferromagnets (AFMs) is investigated. The uniaxial anisotropy is treated exactly and the Heisenberg interactions are treated within unified molecular field theory (MFT), where thermodynamic variables are expressed in terms of directly measurable parameters. The properties of collinear AFMs with ordering along the z axis (D>0) in applied field H z = 0 are calculated versus D and temperature T, including the ordered moment μ, the Néel temperature T N , the magnetic entropy, internal energy, heat capacity, and the anisotropic magnetic susceptibilities χ ∥ and χ ⊥ in the paramagnetic (PM) and AFM states. The high-field average magnetization per spin μ z (H z ,D,T) is found, and the critical field H c (D,T) is derived at which the second-order AFM to PM phase transition occurs. The magnetic properties of the spin-flop (SF) phase are calculated, including the zero-field properties T N (D) and μ(D,T). The high-field μ z (H z ,D,T) is determined, together with the associated spin-flop field H SF (D,T) at which a second-order SF to PM phase transition occurs. The free energies of the AFM, SF, and PM phases are derived from which H z –T phase diagrams are constructed. For f J =–1 and –0.75, where f J = θ pJ /T NJ and θ pJ and T NJ are the Weiss temperature in the Curie-Weiss law and the Néel temperature due to exchange interactions alone, respectively, phase diagrams in the H z –T plane similar to previous results are obtained. However, for f J = 0 we find a topologically different phase diagram where a spin-flop bubble with PM and AFM boundaries occurs at finite H z and T. Also calculated are properties arising from a perpendicular magnetic field, including the perpendicular susceptibility χ ⊥ (D,T), the associated effective torque at low fields arising from the –DS 2 z term in the Hamiltonian, the high-field perpendicular magnetization μ ⊥ , and the perpendicular critical field H c⊥ at which the second-order AFM to PM phase transition occurs. In addition to the above results for D > 0, the T N (D) and ordered moment μ(T,D) for collinear AFM ordering along the x axis with D < 0 are determined. In order to compare the properties of the above spin systems with those of noninteracting systems with –DS 2 z uniaxial anisotropy with either sign of D, Supplemental Material is provided in which results for the thermal and magnetic properties of such noninteracting spin systems are given.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Secondary-Side Active Rectifier Synchronization and Control in LCC-LC Compensated Inductive Power Transfer Systems

This paper proposes a communication-less synchronization and output power control strategy for LCC-LC compensated inductive power transfer systems with a secondaryside active rectifier. By replacing the passive rectifier with an active bridge, the proposed architecture eliminates the need for a separate secondary-side DC-DC converter and enables direct regulation of battery power. The key challenge in such a system is achieving robust synchronization of the secondary active bridge with the primary inverter without access to primary-side signals. To address this issue, the secondary-side variables are analyzed and the current through the secondary series inductor is identified as a load-independent synchronization variable with a fixed phase relationship to the primary excitation. A fixed-frequency second order generalized (SOGI)-based phase locked loop (PLL) is used to estimate the switching phase and frequency, and a phase-shift-based output power controller is developed using a phasor-transformer model of the active bridges. A modified control formulation is further introduced to remove phase-dependent loop-gain variation and simplify controller design. Simulation results validate the proposed method under both fixed and varying switching-frequency conditions, demonstrating successful frequency tracking from 82 kHz to 88 kHz and accurate output power regulation over a wide operating range. The proposed method offers a compact and effective solution for high-power IPT systems by enabling secondary-side synchronization and control without wireless communication or an additional DC-DC conversion stage.

Fernandes, Arnold Anthony [ORNL] (ORCID:0009000631↗
DOE OSTI2022

Sequential Linearization Method for Bound-Constrained Mathematical Programs with Complementarity Constraints

Here, we propose an algorithm for solving bound-constrained mathematical programs with complementarity constraints on the variables. Each iteration of the algorithm involves solving a linear program with complementarity constraints in order to obtain an estimate of the active set. The algorithm enforces descent on the objective function to promote global convergence to B-stationary points. We provide a convergence analysis and preliminary numerical results on a range of test problems. We also study the effect of fixing the active constraints in a bound-constrained quadratic program that can be solved on each iteration in order to obtain fast convergence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Los Alamos National Laboratory Fall Avian Migration Monitoring Report 2010-2018

During the fall of 2018, biologists at Los Alamos National Laboratory completed the ninth year of monitoring fall migration passerines (songbirds). Songbirds were captured at a mist-netting station located in a wetland/riparian complex in Technical Area 36 on the north side of Pajarito Road in Los Alamos County. Captured birds were identified, measured, and banded with a United States Fish and Wildlife Service (USFWS) migratory bird band. Banding operations took place between August and October, with the completion of a total of 10 mist-netting sessions annually. This project was conducted as part of the implementation of the Biological Resources Management Plan and is in compliance with the 2013 Memorandum of Understanding between the USFWS and the United States Department of Energy/National Nuclear Security Administration and Executive Order 13186. Between 2010 and 2018 the overall number of birds and species captured was variable. We saw significant differences in the bird community composition when we compared earlier years to later years. Although environmental variables including a drought severity index, were not shown to influence the community composition, insectivorous birds and the percentage of hatch year birds were influenced by drought severity. Additionally, Audubon’s Warbler and Virginia’s Warbler numbers decreased from 2010 and showed declining trends in their abundances. The variability in bird populations is likely driven by regional climatic factors, but more data are needed for robust assessments.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Aircraft model prototypes which have specified handling-quality time histories

Several techniques for obtaining linear constant-coefficient airplane models from specified handling-quality time histories are discussed. One technique, the pseudodata method, solves the basic problem, yields specified eigenvalues, and accommodates state-variable transfer-function zero suppression. The method is fully illustrated for a fourth-order stability-axis small-motion model with three lateral handling-quality time histories specified. The FORTRAN program which obtains and verifies the model is included and fully documented.

Johnson, S. H.↗
NASA NTRS1978

A safety margin and flight reference system and display for powered-lift aircraft

A study was conducted to explore the feasibility of a safety margin and flight reference system for those powered-lift aircraft which require a backside piloting technique. The main objective was to display multiple safety margin criteria as a single variable which could be tracked both manually and automatically and which could be monitored in order to derive safety margin status. The study involved a pilot-in-the-loop analysis of several system concepts and a simulator experiment to evaluate those concepts showing promise. A system was ultimately configured which yielded reasonable compromises in controllability, status information content, and the ability to regulate safety margins at some expense of the allowable low speed flight path envelope.

Heffley, R. K.↗
NASA NTRS1979

Stability analysis of a liquid fuel annular combustion chamber

The problems of combustion instability in an annular combustion chamber are investigated. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations. By evaluating these effects, parameters which cause instabilities to occur in the combustion chamber can be ascertained. It is assumed that in the annular combustion chamber, the liquid propellants are injected uniformly across the injector face, the combustion processes are distributed throughout the combustion chamber, and that no time delay occurs in the combustion processes.

Mcdonald, G. H.↗
NASA NTRS1980

A computer program for determining truncation error coefficients for Runge-Kutta methods

The basic structure of a program to generate the truncation error coefficients for Runge-Kutta (RK) methods is reformulated to reduce storage requirements significantly and to accommodate variable dimensioning. This FORTRAN program, SUBROUTINE RKEQ, determines truncation error coefficients for RK algorithms for orders 1 through 10 and extends the order of coefficients through 12 with the 11th- and 12th-order terms determined following the patterns used to establish the lower order coefficients. Both subroutines (the original and RKEQ) are also written to treat RK m-fold methods which utilize m known derivatives of f to increase the order of the algorithm. Setting m = 0 gives the classical RK algorithm.

Horn, M. K.↗
NASA NTRS1981

Computation of the modes and polar moment of inertial of the blades of an HAWT

The coupled differential equations of motion of the blades of a horizontal axis wind turbine are solved numerically, permitting the optimization of the design at relatively low cost. The equation of motion is transformed into a set of first order equations and solved with fourth order Runge-Kutta integrators. This technique is applied to a twisted, tapered blade of variable cross section and stiffness including discontinuities. The first six natural frequencies and mode shapes are obtained. The polar moment of inertia of the blades is obtained as a function of frequency and rotational speed.

Beaulieu, G.↗
NASA NTRS1983

Evidence for temporal variations in polar plumes

The variability of the O VI and Mg X emissions from four plumes is studied in order to place empirical constraints on possible plasma heating processes and acceleration mechanisms for mass flows in polar plumes. The observations comprise a 40 min seqeunce of spectroheliograms acquired with a temporal resolution of approximately 2 min. It is found that the Mg X emission from each plume has a strong 'dc' component and that any variations occurred at approximately the 10 percent (of the mean) level or less. This means that temporal variations in the density and/or temperature were 5 percent or less. Some statistically significant variations in the measured Mg X signals are seen. It is noted that two of the plumes decreased in brightness by approximately 10 percent during the 40 min period of the observations. Besides this, two plumes exhibited short-term (few minute) variations that are seen as statistically significant at a confidence level of more than 90 percent. The latter variations (if real) are thought to be due to small (approximately 20,000 K) fluctuations in temperature caused by fluctuations in the local plasma heating rate.

Withbroe, G. L.↗
NASA NTRS1987

Modeling digital control systems with MA-prefiltered measurements

Three discrete state variable representations are derived for a continuous-time plant driven by a zero-order-hold with a combination of instantaneous measurements and measurements prefiltered by moving-average (MA) digital filters. These representations allow the control system engineer to accurately model plants of this type in a form which permits him to use standard techniques to design digital feedback controllers for them. An example is presented which illustrates how to obtain the coefficient matrices in each representation. Guidelines are presented for choosing the best representation to use for any given point.

Polites, Michael E.↗
NASA NTRS1987

Reexamination of the SAS 2 Cygnus X-3 data

Recent observations of Cygnus X-3 have shown marked variability of the radiation on short time scales. In particular, the bursts lasting on the order of 10 minutes, seen in both the infrared and very high energy (greater than 10 to the 11th eV) gamma-ray regions, and the time-variations on many scales at high energies, have stimulated a reanalysis of the March 6 to 13, 1973 SAS 2 high-energy gamma-ray data. Although a clear periodicity in the E greater 35 MeV gamma radiation is observed at the 4.79 hr period seen in X-rays, there is no evidence for major variations of the radiation from one day to the next, and no statistically significant evidence for bursts on the 10-minute time scale seen in the infrared or very high energy ranges. If the excess observed had been predominantly in the form of ten minute bursts even at a rate as high as two/day, a clearly significant set of bursts would have been seen.

Fichtel, C. E.↗
NASA NTRS1989

Further developments in modeling digital control systems with MA-prefiltered measurements

State variable representations are presented for a continuous-time plant driven by a zero-order-hold with multirate-sampled measurements prefiltered by multi-input/multi-output moving average (MA) processes. These representations have broad application, but are known to be useful in the aerospace field for modeling systems with star trackers and some state-of-the-art rate-gyroscopes and accelerometers.

Polites, Michael E.↗
NASA NTRS1991

State-Variable Representations For Moving-Average Sampling

Two state-variable representations derived for continuous-time plant driven by control algorithm including zero-order hold and measurements sampled at mutliple rates by multiple-input/multiple-output moving-average processes. New representations enhance observability and controllability of plant. Applications include mathematical modeling of navigation systems including star trackers, gyroscopes, and accelerometers.

Polites, Michael E.↗
NASA NTRS1991

A study of the response of nonlinear springs

The various phases to developing a methodology for studying the response of a spring-reinforced arch subjected to a point load are discussed. The arch is simply supported at its ends with both the spring and the point load assumed to be at midspan. The spring is present to off-set the typical snap through behavior normally associated with arches, and to provide a structure that responds with constant resistance over a finite displacement. The various phases discussed consist of the following: (1) development of the closed-form solution for the shallow arch case; (2) development of a finite difference analysis to study (shallow) arches; and (3) development of a finite element analysis for studying more general shallow and nonshallow arches. The two numerical analyses rely on a continuation scheme to move the solution past limit points, and to move onto bifurcated paths, both characteristics being common to the arch problem. An eigenvalue method is used for a continuation scheme. The finite difference analysis is based on a mixed formulation (force and displacement variables) of the governing equations. The governing equations for the mixed formulation are in first order form, making the finite difference implementation convenient. However, the mixed formulation is not well-suited for the eigenvalue continuation scheme. This provided the motivation for the displacement based finite element analysis. Both the finite difference and the finite element analyses are compared with the closed form shallow arch solution. Agreement is excellent, except for the potential problems with the finite difference analysis and the continuation scheme. Agreement between the finite element analysis and another investigator's numerical analysis for deep arches is also good.

Hyer, M. W.↗