Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter continuation method”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Improved Guarantees for Optimal Nash Equilibrium Seeking and Bilevel Variational Inequalities

We consider a class of hierarchical variational inequality (VI) problems that subsumes VI-constrained optimization and several other problem classes, including the optimal solution selection problem and the optimal Nash equilibrium (NE) seeking problem. Our main contribution is threefold. (i) We consider bilevel VIs with monotone and Lipschitz continuous mappings and devise a single-timescale iteratively regularized extragradient method, named IR-EG 𝚖,𝚖 . We improve the existing iteration complexity results for addressing both bilevel VI and VI-constrained convex optimization problems. (ii) Under the strong monotonicity of the outer-level mapping, we develop a method named IR-EG 𝚜,𝚖 and derive faster guarantees than those in (i). We also study the iteration complexity of this method under a constant regularization parameter. These results appear to be new for both bilevel VIs and VI-constrained optimization. (iii) To our knowledge, complexity guarantees for computing the optimal NE in nonconvex settings do not exist. Motivated by this lacuna, we consider VI-constrained nonconvex optimization problems and devise an inexactly projected gradient method, named IPR-EG, where the projection onto the unknown set of equilibria is performed using IR-EG 𝚜,𝚖 with a prescribed termination criterion and an adaptive regularization parameter. We obtain new complexity guarantees in terms of a residual map and an infeasibility metric for computing a stationary point. Here, we validate the theoretical findings using preliminary numerical experiments for computing the best and the worst NEs.

bilevel optimization↗
NASA NTRS1989

Preliminary development of an intelligent computer assistant for engine monitoring

As part of the F-18 high-angle-of-attack vehicle program, an AI method has been developed for the real-time monitoring of the propulsion system and for the identification of recovery procedures for the F404 engine. The aim of the development program is to provide enhanced flight safety and to reduce the duties of the propulsion engineers. As telemetry data is received, the results are continually displayed in a number of different color graphical formats. The system makes possible the monitoring of the engine state and the individual parameters. Anomaly information is immediately displayed to the engineer.

Disbrow, James D.↗
NASA NTRS2018

Assessment of Axial Wave Number and Mean Flow Uncertainty on Acoustic Liner Impedance Education

A key parameter in designing and assessing advanced broadband acoustic liners to achieve the current and future noise reduction goals is the acoustic impedance presented by the liner. This parameter, intrinsic to a specific liner configuration, is dependent on sound pressure level and grazing flow velocity. Current impedance eduction approaches have, in general, provided excellent results and continue to be employed throughout the acoustic liner community. However, some recent applications have indicated a possible dependence of the educed impedance on the direction of incident waves relative to the mean flow. The purpose of the current study is to investigate this unexpected behavior for various impedance eduction methods based on the Pridmore-Brown and convected Helmholtz equations. Specifically, the effects of flow profile and axial wavenumber uncertainties on educed impedances for upstream and downstream sources are investigated. The uniform flow results demonstrate the importance of setting a correct Mach number value in obtaining consistent educed impedances for upstream and downstream sources. In fact, the consistency of results over the two source locations was greatly improved by a slight modification of the uniform flow Mach number. In addition, uncertainty in educed axial wavenumber was also illustrated to correlate well with differences in the educed impedances, even with modified uniform flow Mach number. Finally, while less straightforward than in the uniform flow case, it appears that modification of the mean flow profile may also improve consistency of results for upstream and downstream results when shear flow is included.

Nark, Douglas M.↗
NASA NTRS2024

Double Power-law Formation by Sequential Particle Acceleration

Spectral double power laws are common in solar high-energy phenomena such as flares and interplanetary energetic-electron events. However, the physical mechanism that produces the changes in power-law index within a single spectrum is unclear. We developed a fully analytical method of forming single power-law spectra from sequential acceleration of particles orbiting inside and hopping between simulated large-scale magnetic islands formed by flare reconnection. Here, we extend the analytical method to the formation of double power-law spectra by assuming sequential acceleration in two successive regions with different acceleration and particle-transport rates. The resulting spectral distribution is continuous and smooth, with a flattening at low energies, two power-law regions at mid-energies, and a steep rollover at high energies. The model provides analytical expressions for the spectral indices, all energy breaks, and normalization constants as functions of just three physical parameters of each acceleration region: (1) the energy gain in each accelerator, (2) the percentage of particles transferred between accelerators, and (3) the number of accelerators visited. One of the most salient predictions of our work is that the spectral index at high (low) energies is determined by the parameters of the first "seed" (second) acceleration region. By constructing the spectral distribution through an iterative analytical process, the evolution toward a double power law is easily characterized and explained. Our analytical model provides tools to interpret space- and ground-based observations from RHESSI, FOXSI, NuSTAR, Solar Orbiter/STIX, EOVSA, and future high-energy missions.

Solar flares↗
NASA NTRS2000

Continuous Processing With Mars Gases

Current Martian missions call for the production of oxygen for breathing, and fuel and oxygen for propulsion to be produced from atmospheric carbon dioxide (CO2). Adsorption and freezing are the two methods considered for capturing CO2 from the atmosphere. However, the nitrogen (N2) and argon (Ar), which make up less than 5 percent of the atmosphere, cause difficulties with both of these processes by blocking the CO2. This results in the capture process rapidly changing from a pressure driven process to a diffusion controlled process. To increase the CO2 capture rates, some type of mechanical pump is usually proposed to remove the N2 and Ar. The N2 and Ar are useful and have been proposed for blanketing and pressurizing fuel tanks and as buffer gas for breathing air for manned missions. Separation of the Martian gases with the required purity can be accomplished with a combination of membranes. These membrane systems do not require a high feed pressure and provide suitable separation. Therefore, by use of the appropriate membrane combination with the Martian atmosphere supplied by a compressor a continuous Supply Of CO2 for fuel and oxygen production can be supplied. This phase of our program has focused on the selection of the membrane system. Since permeation data for membranes did not exist for Martian atmospheric pressures and temperatures, this information had to be compiled. The general trend as the temperature was lowered was for the membranes to become more selective. In addition, the relative permeation rates between the three gases changed with temperature. The end result was to provide design parameters that could be used to separate CO2 from N2 and Ar. This paper will present the membrane data, provide the design requirements for a compressor, and compare the results with adsorption and freezer methods.

Parrish, Clyde↗
NASA NTRS2001

Continuous Processing with Mars Gases

Current Martian missions call for the production of oxygen for breathing, and fuel and oxygen for propulsion to be produced from atmospheric carbon dioxide (CO2). Adsorption and freezing are the two methods considered for capturing CO, from the atmosphere. However, the nitrogen (N2) and argon (Ar), which make up less than 5 percent of the atmosphere, cause difficulties with both of these processes by blocking the CO2, This results in the capture process rapidly changing from a pressure driven process to a diffusion controlled process. To increase the CO, capture rates, some type of mechanical pump is usually proposed to remove the N2 and Ar. The N2 and Ar are useful and have been proposed for blanketing and pressurizing fuel tanks and as buffer gas for breathing air for manned missions. Separation of the Martian gases with the required purity can be accomplished with a combination of membranes. These membrane systems do not require a high feed pressure and provide suitable separation. Therefore, by use of the appropriate membrane combination with the Martian atmosphere supplied by a compressor a continuous supply of CO2 for fuel and oxygen production can be supplied. This phase of our program has focused on the selection of the membrane system. Since permeation data for membranes did not exist for Martian atmospheric pressures and temperatures, this information had to be compiled. The general trend as the temperature was lowered was for the membranes to become more selective. In addition, the relative permeation rates between the three gases changed with temperature. The end result was to provide design parameters that could be used to separate CO2 from N2 and Ar. This paper will present the membrane data, provide the design requirements for a compressor, and compare the results with adsorption and freezer methods.

Parrish, Clyde↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1983

Technical and investigative support for high density digital satellite recording systems

Recent results of dropout measurements and defect analysis conducted on one reel of Ampex 721 which was submitted for evaluation by the manufacturer are described. The results or status of other tape evaluation activities are also reviewed. Several changes in test interpretations and applications are recommended. In some cases, deficiencies in test methods or equipment became apparent during continued work on this project and other IITRI tape evaluation projects. Techniques and equipment for future tasks such as tape qualification are also recommended and discussed. Project effort and expenditures were kept at a relatively low level. This rate provided added development time and experience with the IITRI Dropout Measurement System, which is approaching its potential as a computer based dropout analysis tool. Another benefit is the expanded data base on critical parameters that can be achieved from tests on different tape types and lots as they become available. More consideration and effort was directed toward identification of critical parameters, development of meaningful repeatable test procedures, and tape procurement strategy.

Schultz, R. A.↗
NASA NTRS1986

A method for screening of plant species for space use

A cost-effective methodology which monitors numerous dynamic aspects of carbon assimilation and allocation kinetics in live, intact plants is discussed. Analogous methods can apply to nitrogen uptake and allocation. This methodology capitalizes on the special properties of the short-lived, positron-gamma emitting isotope C-11 especially when applied as CO2-11 in a special extended square wave (ESW) pattern. The 20.4 minute half-life allows for repeated or continuous experiments on the same plant over periods of minutes, hours, days, or weeks. The steady-state isotope equilibrium approached during the ESW experiments, and the parameters which can be analyzed by this technique are also direct results of that short half-life. Additionally, the paired .511 MeV gamma rays penetrate any amount of tissue and their 180 deg opposite orientation provides good collimation and allows coincidence counting which nearly eliminates background.

Goeschl, J. D.↗
NASA NTRS1986

Boundary perturbation method for free boundary problem in convectively cooled continuous casting

Novel mathematical techniques are here used to treat general boundary conditions in convectively cooled continuous ingot casting problems, applying a constant convective heat transfer coefficient at the ingot sides (so that the temperature distribution along the sides is unknown). The problem is first inverted to treat the physical coordinates as dependent variables of temperature and heat flow, mapping the ingot into a region that is nearly rectangular. A boundary perturbation method is then used to obtain an analytical solution in this near-rectangular region. Solidification interface shapes depend on two dimensionless parameters: one is a function of the casting velocity, while the other is dependent on the convective heat transfer coefficient at the cooled wall.

Siegel, R.↗
NASA NTRS1993

The origin of Venusian channels: Modelling of thermal erosion by lava

Magellan imagery has revealed that channels, apparently volcanic in origin, are abundant on the surface of Venus. There has been much debate about the origin of these channels. Are they the result of erosional (either thermal or mechanical) or constructional processes? A common characteristic of the simple sinuous channels is that they show evidence of erosion near their source and then become purely constructional, forming levees and in some cases roofing over completely. One method of showing that thermal erosion is capable of producing the type of channels seen is to use computer modeling incorporating the physical conditions on Venus and the physical characteristics of the different types of lava that may have been erupted. It is possible to calculate, relatively easily, two channel parameters. The first is the erosion rate, which combined with eruption duration, gives depth. The second is for how long after leaving the source the erupted lava will continue to be capable of thermal erosion before constructional processes dominate. Making assumptions about the rheology of the lava (e.g., assume it behaves as a Bingham plastic) along with the slope angle yields a flow velocity and therefore a distance over which thermal erosion will take place. Due to the resolution (both vertical and horizontal) of the Magellan altimetric data, the distance from the source that the channel is erosional can be much more accurately measured than the depth of the channel. This will remain the case until stereo imagery becomes available for large areas of the planet.

Bussey, D. B. J.↗
NASA NTRS1975

Optimization of structures undergoing harmonic or stochastic excitation

The optimal design was investigated of simple structures subjected to dynamic loads, with constraints on the structures' responses. Optimal designs were examined for one dimensional structures excited by harmonically oscillating loads, similar structures excited by white noise, and a wing in the presence of continuous atmospheric turbulence. The first has constraints on the maximum allowable stress while the last two place bounds on the probability of failure of the structure. Approximations were made to replace the time parameter with a frequency parameter. For the first problem, this involved the steady state response, and in the remaining cases, power spectral techniques were employed to find the root mean square values of the responses. Optimal solutions were found by using computer algorithms which combined finite elements methods with optimization techniques based on mathematical programming. It was found that the inertial loads for these dynamic problems result in optimal structures that are radically different from those obtained for structures loaded statically by forces of comparable magnitude.

Johnson, E. H.↗
NASA NTRS1991

Reliability analysis of continuous fiber composite laminates

This paper describes two methods, the maximum distortion energy (MDE) and the principle of independent action (PIA), developed for the analysis of the reliability of a single continuous composite lamina. It is shown that, for the typical laminated plate structure, the individual lamina reliabilities can be combined in order to produce the upper and the lower bounds of reliability for the laminate, similar in nature to the bounds on properties produced from variational elastic methods. These limits were derived for both the interactive and the model failure considerations. Analytical expressions were also derived for the sensitivity of the reliability limits with respect to changes in the Weibull parameters and in loading conditions.

Thomas, David J.↗
NASA NTRS1991

mu analysis with real parametric uncertainty

The authors give a broad overview, from a LFT (linear fractional transformation)/mu perspective, of some of the theoretical and practical issues associated with robustness in the presence of real parametric uncertainty, with a focus on computation. Recent results on the properties of mu in the mixed case are reviewed, including issues of NP completeness, continuity, computation of bounds, the equivalence of mu and its bounds, and some direct comparisons with Kharitonov-type analysis methods. In addition, some advances in the computational aspects of the problem, including a novel branch and bound algorithm, are briefly presented together with numerical results. The results suggest that while the mixed mu problem may have inherently combinatoric worst-case behavior, practical algorithms with modest computational requirements can be developed for problems of medium size (less than 100 parameters) that are of engineering interest.

Young, Peter M.↗
NASA NTRS2019

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

ice sheet model↗
NASA NTRS1993

A proof of the Woodward-Lawson sampling method for a finite linear array

An extension of the continuous aperture Woodward-Lawson sampling theorem has been developed for a finite linear array of equidistant identical elements with arbitrary excitations. It is shown that by sampling the array factor at a finite number of specified points in the far field, the exact array factor over all space can be efficiently reconstructed in closed form. The specified sample points lie in real space and hence are measurable provided that the interelement spacing is greater than approximately one half of a wavelength. This paper provides insight as to why the length parameter used in the sampling formulas for discrete arrays is larger than the physical span of the lattice points in contrast with the continuous aperture case where the length parameter is precisely the physical aperture length.

Somers, Gary A.↗
NASA NTRS2003

A Computationally Efficient Meshless Local Petrov-Galerkin Method for Axisymmetric Problems

The Meshless Local Petrov-Galerkin (MLPG) method is one of the recently developed element-free methods. The method is convenient and can produce accurate results with continuous secondary variables, but is more computationally expensive than the finite element method. To overcome this disadvantage, a simple Heaviside test function is chosen. The computational effort is significantly reduced by eliminating the domain integral for the axisymmetric potential problems and by simplifying the domain integral for the axisymmetric elasticity problems. The method is evaluated through several patch tests for axisymmetric problems and example problems for which the exact solutions are available. The present method yielded very accurate solutions. The sensitivity of several parameters of the method is also studied.

Raju, I. S.↗
NASA NTRS2010

Drop Axis Ratio Distributions in Stratiform and Convective Rain

A fully calibrated low profile 2D video disdrometer (2DVD) has been recording many different rainfall events in Northern Alabama (USA) since June 2007. An earlier publication reported drop shapes and axis ratio distributions determined for some of the events. For one of the cases examined, a noticeable shift in the 3.5 - 3.75 mm drop axis ratio distribution was noted. In this paper, we extend the earlier work by separating the 2DVD measurements into stratiform and convective rain. The separation is made possible by using the minute-by-minute drop size distribution (DSD) measured by the 2DVD. The 1-minute DSDs are fitted to a gamma distribution, and using a simple indexing technique which involves two of the fitted parameters, periods of convective and stratiform rain are separated for a given event. The output of the DSD indexing technique is qualitatively confirmed by comparing with simultaneous time series observations from a co-located UHF profiler which continuously records height profiles of reflectivity, Doppler mean and spectral width, all of which enable the identification of bright-band periods and, furthermore, periods of moderate and deep convection. Excellent consistency is found between the output of the DSD-based separation method and the profiler observations. Next, we utilize the output of DSD index-based separation method to flag the periods of severe convection for a given event. Drop axis ratios during the flagged periods are derived and compared with those during stratiform rain periods. Five cases have been considered. Axis ratio distributions do not show appreciable differences between stratiform and convective periods for four of the cases. The fifth case (the same case as reported earlier) shows a shift in the 3.5 - 3.75 mm drop axis ratios during a prolonged period of convection. The contoured shapes for these drops determined from the 2DVD camera data indicate the possibility of non-axisymmetric oscillations, compared with the contoured images for other events which fit well to our reference drop shapes. For all of above cases, observations from a C-band polarimetric radar - situated 15 km away are examined. The variations between the co-polar radar reflectivity and the differential reflectivity as well as the specific differential phase are compared with the 2DVD data based scattering calculations for the 5 events. The implications will be discussed.

Thurai, M.↗