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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Characterisation of the spatial variability of material properties of Gilsocarbon and NBG-18 using random fields

Graphite is a candidate material for Generation IV concepts and is used as a moderator in Advanced Gas-cooled Reactors (AGR) in the UK. Spatial material variability is present within billets causing different material property values between different components. Variations in material properties and irradiation effects can produce stress concentrations and diverse mechanical responses in a nuclear reactor graphite core. In order to characterise the material variability, geostatistical techniques called variography and random field theory were adapted for studying the density and Young's modulus of a billet of Gilsocarbon and NBG-18 graphite grades. Variography is a technique for estimating the distance over which material property values have significant spatial correlation, known as the scale of fluctuation or spatial correlation length. The paper uses random field theory to create models that mimic the original spatial and statistical distributions of the original data set. This study found different values of correlation length for density and Young's modulus around the edges of a Gilsocarbon billet, while in the case of NBG-18, similar correlation lengths where found across the billet. Furthermore examples of several random fields are given to reproduce the spatial patterns and values found in the original data.

36 MATERIALS SCIENCE↗

NASA NTRS1972

On the use of the Hill variables in artificial satellite theory - Brouwer's theory.

By reformulating Brouwer's first-order satellite theory in terms of the Hill variables, it is demonstrated that these variables possess several advantages over the Delaunay variables. The final algorithm is simpler and more compact, there is no singularity at zero eccentricity, and the process of deriving the perturbations is considerably simplified. An algorithm for the computation of position and velocity is presented which is significantly simpler and faster in use than Brouwer's algorithm. Moreover, whereas the latter algorithm breaks down for circular and equatorial orbits, the present one holds for all eccentricities less than one and for all inclinations, with the exception of the critical inclination.

Aksnes, K.↗

DOE OSTI2017

Design principles from multiscale simulations to predict nanostructure in self-assembling ionic liquids

Molecular dynamics simulations (up to the nanoscale) were performed on the 3-methyl-1-pentylimidazolium ionic liquid cation paired with three anions; chloride, nitrate, and thiocyanate as aqueous mixtures, using the effective fragment potential (EFP) method, a computationally inexpensive way of modeling intermolecular interactions. The simulations provided insight (preferred geometries, radial distribution functions and theoretical proton NMR resonances) into the interactions within the ionic domain and are validated against 1 H NMR spectroscopy and small- and wide-angle X-ray scattering experiments on 1-decyl-3-methylimidazolium. Ionic liquids containing thiocyanate typically resist gelation and form poorly ordered lamellar structures upon mixing with water. Conversely, chloride, a strongly coordinating anion, normally forms strong physical gels and produces well-ordered nanostructures adopting a variety of structural motifs over a very wide range of water compositions. Nitrate is intermediate in character, whereby upon dispersal in water it displays a range of viscosities and self-assembles into nanostructures with considerable variability in the fidelity of ordering and symmetry, as a function of water content in the binary mixtures. The observed changes in the macro and nanoscale characteristics were directly correlated to ionic domain structures and intermolecular interactions as theoretically predicted by the analysis of MD trajectories and calculated RDFs. Specifically, both chloride and nitrate are positioned in the plane of the cation. Anion to cation proximity is dependent on water content. Thiocyanate is more susceptible to water insertion into the second solvent shell. Experimental 1 H NMR chemical shifts monitor the site-specific competition dependence with water content in the binary mixtures. As a result, thiocyanate preferentially sits above and below the aromatic ring plane, a state disallowing interaction with the protons on the imidazolium ring.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1982

Asymptotic Expansions and Estimates for the Capillary Problem

The asymptotic properties for the small Bond number B of the equilibrium capillary interface interior to a circular cylindrical tube vertically dipped in an infinite reservoir of liquid are discussed. (The Bond number B is a dimensionless parameter which is the ratio of gravitational to capillary forces.) The formal expansion in powers of B of the solution to the differential equation describing the equilibrium surface (as can be obtained by standard perturbation methods) is proved to be truly asymptotic -- to all orders and uniformly in the variable and parameter gamma, the contact angle. Sequences of general estimates, in closed form, from above and from below, are also given for the solution and related functions. The M-th term in these sequences are asymptotically exact to order m. An idiosyncrasy of the problem, crucial in obtaining these estimates, is the absolute monotonicity of the structural function of the system in integral form.

Brulois, F. P.↗

DOE OSTI2019

Arbitrary-order time-accurate semi-Lagrangian spectral approximations of the Vlasov–Poisson system

The Vlasov–Poisson system, modeling the evolution of non-collisional plasmas in the electrostatic limit, is approximated by a semi-Lagrangian technique. Spectral methods of periodic type are implemented through a collocation approach. Groups of particles are represented by the Fourier Lagrangian basis and evolve, for a single timestep, along a high-order accurate representation of the local characteristic lines. The time-advancing technique is based on truncated Taylor series that can be, in principle, of any order of accuracy. A variant is obtained by coupling the phase space discretization with high-order accurate Backward Differentiation Formulas (BDF). At each timestep, particle displacements are reinterpolated and expressed in the original basis to guarantee the order of accuracy in all the variables at relatively low costs. Thus, these techniques combine the excellent features of spectral approximations with high-order time integration. The resulting method has excellent conservation properties. Indeed, it can be proven that the total number of particles, proportional to the total mass and charge, is conserved up to the machine precision. Series of numerical experiments are performed in order to assess the real performance. In particular, comparisons with standard benchmarks are examined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1993

Eddy mixing coefficient upper limit derived from the photochemical balance of O2

This work is based on the study of the photochemical balance of molecular oxygen in the martian atmosphere by using a one-dimensional model of photochemical reactions involving species derived from CO2 and H2O. The model is basically similar to one used previously for the study of the regulation of CO on a global scale, but the chemical rates are taken from another source. In the present scheme, the regulation of molecular oxygen is studied over timescales of the order of its photochemical lifetime (approximately equals 30 yr), which is much shorter than typical escape timescales. Thus, the escape fluxes are fixed to the values given by 3 and 4. We examine the calculated equilibrium abundances of O2 for given thermal, eddy diffusion coefficients and H2O profiles. The thermal profile is taken from in the lower atmosphere. At higher levels, in order to include the diurnal and seasonal thermal profile variability, we have also used the IRTM data. In order to study the influence of both temperature and pressure profiles on the O2 mixing ratio, we have made several tests corresponding to different martian seasons. The results show that the influence of pressure and temperature is quantitatively weak compared to the one of K and of the water vapor density (H2O). Thus, in the following we have fixed the pressure at the surface to a value of 7 mbar and we have used unique standard thermal profile corresponding to a profile roughly averaged over the year, the season, and the day: T equal 205 K at 0 km altitude, 175 K at 25 km, and 145 K at 50 km.

Rosenqvist, J.↗

NASA NTRS2001

Approach for Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

NASA NTRS1995

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗

NASA NTRS1984

A modified least squares formulation for a system of first order equations

Second order equations in terms of auxiliary variables similar to potential and stream functions are obtained by applying a weighted least squares formulation to a first order system. The additional boundary conditions which are necessary to solve the higher order equations are determined and numerical results are presented for the Cauchy-Riemann equations.

Hafez, M. M.↗

NASA NTRS1985

A modified least squares formulation for a system of first-order equations

Hafez, M. M.↗

NASA NTRS2001

Biogeochemical Response to Mesoscale Physical Forcing in the California Current System

In the first part of the project, we investigated the local response of the coastal ocean ecosystems (changes in chlorophyll, concentration and chlorophyll, fluorescence quantum yield) to physical forcing by developing and deploying Autonomous Drifting Ocean Stations (ADOS) within several mesoscale features along the U.S. west coast. Also, we compared the temporal and spatial variability registered by sensors mounted in the drifters to that registered by the sensors mounted in the satellites in order to assess the scales of variability that are not resolved by the ocean color satellite. The second part of the project used the existing WOCE SVP Surface Lagrangian drifters to track individual water parcels through time. The individual drifter tracks were used to generate multivariate time series by interpolating/extracting the biological and physical data fields retrieved by remote sensors (ocean color, SST, wind speed and direction, wind stress curl, and sea level topography). The individual time series of the physical data (AVHRR, TOPEX, NCEP) were analyzed against the ocean color (SeaWiFS) time-series to determine the time scale of biological response to the physical forcing. The results from this part of the research is being used to compare the decorrelation scales of chlorophyll from a Lagrangian and Eulerian framework. The results from both parts of this research augmented the necessary time series data needed to investigate the interactions between the ocean mesoscale features, wind, and the biogeochemical processes. Using the historical Lagrangian data sets, we have completed a comparison of the decorrelation scales in both the Eulerian and Lagrangian reference frame for the SeaWiFS data set. We are continuing to investigate how these results might be used in objective mapping efforts.

Niiler, Pearn P.↗

DOE OSTI2022

Variable Ion Compositions of Solar Energetic Particle Events in the Inner Heliosphere: A Field Line Braiding Model with Compound Injections

We propose a model for interpreting highly variable ion composition ratios in solar energetic particle (SEP) events recently observed by the Parker Solar Probe (PSP) at 0.3–0.45 au. We use numerical simulations to calculate SEP propagation in a turbulent interplanetary magnetic field with a Kolmogorov power spectrum from large scales down to the gyration scale of energetic particles. We show that when the source regions of different species are offset by a distance comparable to the size of the source regions, the observed energetic particle composition He/H can be strongly variable over more than two orders of magnitude, even if the source ratio is at the nominal value. Assuming a 3 He/ 4 He source ratio of 10% in impulsive 3 He-rich events and the same spatial offset of the source regions, the 3 He/ 4 He ratio at observation sites also vary considerably. The variability of the ion composition ratios depends on the radial distance, which can be tested by observations made at different radial locations. We discuss the implications of these results on the variability of ion composition of impulsive events and on further PSP and Solar Orbiter observations close to the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2008

The a(3) Scheme--A Fourth-Order Space-Time Flux-Conserving and Neutrally Stable CESE Solver

The CESE development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To initiate a systematic CESE development of high order schemes, in this paper we provide a thorough discussion on the structure, consistency, stability, phase error, and accuracy of a new 4th-order space-time flux-conserving and neutrally stable CESE solver of an 1D scalar advection equation. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables (the numerical analogues of the dependent variable and its 1st-order and 2ndorder spatial derivatives, respectively) and three equations per mesh point, the new scheme is referred to as the a(3) scheme. Through the von Neumann analysis, it is shown that the a(3) scheme is stable if and only if the Courant number is less than 0.5. Moreover, it is established numerically that the a(3) scheme is 4th-order accurate.

Chang, Sin-Chung↗

NASA NTRS1991

Variable sensitivity moire interferometry

This paper demonstrates that optical microlithography can be used to produce a crossed grating which diffracts light into multiple orders sufficient to record moire interferograms with sensitivities ranging from 2.0 to 0.285 micron/fringe. The grating profile produced by the method is analyzed to establish the diffraction efficiency in each diffraction order, and generalized expressions are given for variable sensitivity moire interferometry. Experimental tests are conducted to verify analytical arguments. In one of these tests, two different diffraction order pairs are used simultaneously to verify that surface displacement can be measured at different sensitivities.

Johnson-Cole, Helen↗

NASA NTRS1982

Validation of zero-order feedback strategies for medium range air-to-air interception in a horizontal plane

A zero order feedback solution of a variable speed interception game between two aircraft in the horizontal plane, obtained by using the method of forced singular perturbation (FSP), is compared with the exact open loop solution. The comparison indicates that for initial distances of separation larger than eight turning radii of the evader, the accuracy of the feedback approximation is better than one percent. The result validates the zero order FSP approximation for medium range air combat analysis.

Shinar, J.↗

DOE OSTI2019

Long-Term Measurements of Methane Ebullition From Thaw Ponds

Arctic regions are experiencing rapid warming, leading to permafrost thaw and formation of numerous water bodies. Although small ponds in particular are considered hot spots for methane (CH 4 ) release, long-term studies of CH 4 efflux from these surfaces are rare. We have collected an extensive data set of CH 4 ebullition (bubbling) measurements from eight small thaw ponds (<0.001 km 2 ) with different physical and hydrological characteristics over four summer seasons, the longest set of observations from thaw ponds to date. The measured fluxes were highly variable with an average of 20.0 mg CH 4 · m -2 · day -1 (median: 4.1 mg CH 4 · m -2 · day -1 , n = 2,063) which is higher than that of most nearby lakes. The ponds were categorized into four types based on clear and significant differences in bubble flux. We found that the amount of CH 4 released as bubbles from ponds was very weakly correlated with environmental variables, like air temperature and atmospheric pressure, and was potentially more related to differences in physical characteristics of the ponds. Using our measured average daily bubble flux plus the available literature, we estimate circumpolar thaw ponds <0.001 km 2 in size to emit between 0.2 and 1.0 Tg of CH 4 through ebullition. Our findings exemplify the importance of high-frequency measurements over long study periods in order to adequately capture the variability of these water bodies. Through the expansion of current spatial and temporal monitoring efforts, we can increase our ability to estimate CH 4 emissions from permafrost pond ecosystems now and in the future.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2003

Single-Vector Calibration of Wind-Tunnel Force Balances

An improved method of calibrating a wind-tunnel force balance involves the use of a unique load application system integrated with formal experimental design methodology. The Single-Vector Force Balance Calibration System (SVS) overcomes the productivity and accuracy limitations of prior calibration methods. A force balance is a complex structural spring element instrumented with strain gauges for measuring three orthogonal components of aerodynamic force (normal, axial, and side force) and three orthogonal components of aerodynamic torque (rolling, pitching, and yawing moments). Force balances remain as the state-of-the-art instrument that provide these measurements on a scale model of an aircraft during wind tunnel testing. Ideally, each electrical channel of the balance would respond only to its respective component of load, and it would have no response to other components of load. This is not entirely possible even though balance designs are optimized to minimize these undesirable interaction effects. Ultimately, a calibration experiment is performed to obtain the necessary data to generate a mathematical model and determine the force measurement accuracy. In order to set the independent variables of applied load for the calibration 24 NASA Tech Briefs, October 2003 experiment, a high-precision mechanical system is required. Manual deadweight systems have been in use at Langley Research Center (LaRC) since the 1940s. These simple methodologies produce high confidence results, but the process is mechanically complex and labor-intensive, requiring three to four weeks to complete. Over the past decade, automated balance calibration systems have been developed. In general, these systems were designed to automate the tedious manual calibration process resulting in an even more complex system which deteriorates load application quality. The current calibration approach relies on a one-factor-at-a-time (OFAT) methodology, where each independent variable is incremented individually throughout its full-scale range, while all other variables are held at a constant magnitude. This OFAT approach has been widely accepted because of its inherent simplicity and intuitive appeal to the balance engineer. LaRC has been conducting research in a "modern design of experiments" (MDOE) approach to force balance calibration. Formal experimental design techniques provide an integrated view to the entire calibration process covering all three major aspects of an experiment; the design of the experiment, the execution of the experiment, and the statistical analyses of the data. In order to overcome the weaknesses in the available mechanical systems and to apply formal experimental techniques, a new mechanical system was required. The SVS enables the complete calibration of a six-component force balance with a series of single force vectors.

Parker, P. A.↗

DOE OSTI2018

Consequences of neglecting the interannual variability of the solar resource: A case study of photovoltaic power among the Hawaiian Islands

The interannual variability of the solar irradiance and meteorological conditions are often ignored in favor of single-year data sets for modeling power generation and evaluating the economic value of photovoltaic (PV) power systems. Yet interannual variability significantly impacts the generation from one year to another of renewable power systems such as wind and PV. Consequently, the interannual variability of power generation corresponds to the interannual variability of capital returns on investment. The penetration of PV systems within the Hawaiian Electric Companies' portfolio has rapidly accelerated in recent years and is expected to continue to increase given the state's energy objectives laid out by the Hawaii Clean Energy Initiative. We use the National Solar Radiation Database (1998-2015) to characterize the interannual variability of the solar irradiance and meteorological conditions across the State of Hawaii. These data sets are passed to the National Renewable Energy Laboratory's System Advisory Model (SAM) to calculate an 18-year PV power generation data set to characterize the variability of PV power generation. We calculate the interannual coefficient of variability (COV) for annual average global horizontal irradiance (GHI) on the order of 2% and COV for annual capacity factor on the order of 3% across the Hawaiian archipelago. Regarding the interannual variability of seasonal trends, we calculate the COV for monthly average GHI values on the order of 5% and COV for monthly capacity factor on the order of 10%. We model residential-scale and utility-scale PV systems and calculate the economic returns of each system via the payback period and the net present value. We demonstrate that studies based on single-year data sets for economic evaluations reach conclusions that deviate from the true values realized by accounting for interannual variability.

14 SOLAR ENERGY↗