Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convergence analysis”

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
NASA NTRS1982

Application of a finite element algorithm for the solution of steady transonic Euler equations

A finite element algorithm for the solution of two-dimensional, steady Euler equations is presented which, through a Clebsch-type transformation for the velocity vector, solves the conservation of mass equation with one primary variable and two additional equations for the convection of two new variables. The accuracy and efficiency of this scheme is discussed, and the second-order accuracy attained in the analysis of the convection of the vorticity is demonstrated together with the efficient treatment of rotational and irrotational flow subregions. A sample problem is used to show the accuracy of the numerical scheme and its convergence characteristics.

Akay, H. U.↗
NASA NTRS1989

Rotordynamic coefficients for stepped labyrinth gas seals

The basic equations are derived for compressible flow in a stepped labyrinth gas seal. The flow is assumed to be completely turbulent in the circumferential direction where the friction factor is determined by the Blasius relation. Linearized zeroth and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variables solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are presented in the form of a parametric study, since there are no known experimental data for the rotordynamic coefficients of stepped labyrinth gas seals. The parametric study investigates the relative rotordynamic stability of convergent, straight and divergent stepped labyrinth gas seals. The results show that, generally, the divergent seal is more stable, rotordynamically, than the straight or convergent seals. The results also show that the teeth-on-stator seals are not always more stable, rotordynamically, then the teeth-on-rotor seals as was shown by experiment by Childs and Scharrer (1986b) for a 15 tooth seal.

Scharrer, Joseph K.↗
NASA NTRS1997

Digital storage and analysis of color Doppler echocardiograms

Color Doppler flow mapping has played an important role in clinical echocardiography. Most of the clinical work, however, has been primarily qualitative. Although qualitative information is very valuable, there is considerable quantitative information stored within the velocity map that has not been extensively exploited so far. Recently, many researchers have shown interest in using the encoded velocities to address the clinical problems such as quantification of valvular regurgitation, calculation of cardiac output, and characterization of ventricular filling. In this article, we review some basic physics and engineering aspects of color Doppler echocardiography, as well as drawbacks of trying to retrieve velocities from video tape data. Digital storage, which plays a critical role in performing quantitative analysis, is discussed in some detail with special attention to velocity encoding in DICOM 3.0 (medical image storage standard) and the use of digital compression. Lossy compression can considerably reduce file size with minimal loss of information (mostly redundant); this is critical for digital storage because of the enormous amount of data generated (a 10 minute study could require 18 Gigabytes of storage capacity). Lossy JPEG compression and its impact on quantitative analysis has been studied, showing that images compressed at 27:1 using the JPEG algorithm compares favorably with directly digitized video images, the current goldstandard. Some potential applications of these velocities in analyzing the proximal convergence zones, mitral inflow, and some areas of future development are also discussed in the article.

NASA Discipline Cardiopulmonary↗
NASA NTRS1999

Analysis of the ST-T complex of the electrocardiogram using the Karhunen--Loeve transform: adaptive monitoring and alternans detection

The Karhunen-Loeve transform (KLT) is applied to study the ventricular repolarisation period as reflected in the ST-T complex of the surface ECG. The KLT coefficients provide a sensitive means of quantitating ST-T shapes. A training set of ST-T complexes is used to derive a set of KLT basis vectors that permits representation of 90% of the signal energy using four KLT coefficients. As a truncated KLT expansion tends to favor representation of the signal over any additive noise, a time series of KLT coefficients obtained from successive ST-T complexes is better suited for representation of both medium-term variations (such as ischemic changes) and short-term variations (such as ST-T alternans) than discrete parameters such as the ST level or other local indices. For analysis of ischemic changes, an adaptive filter is described that can be used to estimate the KLT coefficient, yielding an increase in the signal-to-noise ratio of 10 dB (u = 0.1), with a convergence time of about three beats. A beat spectrum of the unfiltered KLT coefficient series is used for detection of ST-T alterans. These methods are illustrated with examples from the European ST-T Database. About 20% of records revealed quasi-periodic salvos of ischemic ST-T change episodes and another 20% exhibit repetitive, but not clearly periodic patterns of ST-T change episodes. About 5% of ischemic episodes were associated with ST-T alterans.

NASA Discipline Cardiopulmonary↗
DOE OSTI2023

An Initial Study of the Convergence Rate of Griffin’s Pebble Bed Reactors Algorithm

This paper presents an initial study of the convergence properties of an iterative algorithm for computing the burnup distribution in a pebble bed reactor (PBR) in its equilibrium core condition. The algorithm is implemented in the Griffin code. Griffin is a reactor multiphysics analysis application jointly developed by Idaho National Laboratory (INL) and Argonne National Laboratory (ANL). Griffin’s PBR algorithm is discussed and simulation data are presented. An alternative matrix formulation of the algorithm is presented that facilitates analysis of the iterative algorithm. The dependence of the spectral radius of the iterative algorithm on operational and discretization parameters is investigated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Uncertainty Underprediction in Monte Carlo Eigenvalue Calculations

It is well-known that statistical estimates obtained from Monte Carlo criticality simulations can be adversely affected by cycle-to-cycle correlations in the fission source, which can lead to estimates of statistical uncertainties that are lower than the true uncertainty by a factor of 5 or more. However, several other more fundamental issues such as adequate source sampling over the fissionable regions and source convergence can have a significant impact on the uncertainties for the calculated eigenvalue and localized tally means, and these issues may be mistaken for effects resulting from cycle-to-cycle correlations. In worst-case scenarios, the uncertainty may be underpredicted by a factor of 40 or more. Since Monte Carlo methods are widely used in criticality safety applications and are increasingly being used for benchmarking reactor analyses, an in-depth understanding of the effects of these issues must be developed in order to support the practical use of Monte Carlo software packages. A rigorous statistical analysis of eigenvalue and localized tally results in Monte Carlo criticality calculations is presented using the SCALE/KENO-VI (continuous-energy version) and MCNP codes. The purpose of this analysis is to investigate the underprediction of uncertainty and its sensitivity to problem characteristics and calculational parameters using two of the most widely used Monte Carlo criticality codes. For the problems considered here, which are fuel rod and fuel assembly problems with reflecting boundary conditions on all four horizontal sides, we show that adequate source convergence along with proper specification of Monte Carlo parameters can reduce the magnitude of uncertainty underprediction to reasonable levels, below a factor of 2 in most cases.

Mervin, Brenden T↗
NASA NTRS2014

Compositional Solution Space Quantification for Probabilistic Software Analysis

Probabilistic software analysis aims at quantifying how likely a target event is to occur during program execution. Current approaches rely on symbolic execution to identify the conditions to reach the target event and try to quantify the fraction of the input domain satisfying these conditions. Precise quantification is usually limited to linear constraints, while only approximate solutions can be provided in general through statistical approaches. However, statistical approaches may fail to converge to an acceptable accuracy within a reasonable time. We present a compositional statistical approach for the efficient quantification of solution spaces for arbitrarily complex constraints over bounded floating-point domains. The approach leverages interval constraint propagation to improve the accuracy of the estimation by focusing the sampling on the regions of the input domain containing the sought solutions. Preliminary experiments show significant improvement on previous approaches both in results accuracy and analysis time.

Monte Carlo Methods↗
DOE OSTI2022

Numerical error analysis of SOLPS-ITER simulations of EAST

Abstract Plasma edge simulations with codes like SOLPS-ITER are widely employed to interpret fusion experiments. However, numerical errors appearing in such simulations are rarely investigated, despite their potential large impact on simulation results. These errors consist of the statistical error and the bias, both resulting from the finite number of employed EIRENE Monte Carlo particles and incomplete convergence, and the discretization error due to the finite resolution of the computational grids. In this contribution, the resulting numerical errors on simulations of pure deuterium and neon seeded H-mode EAST discharges are examined. The statistical error can be kept small compared to other numerical error contributions by averaging the plasma profiles. This allows investigating the bias and discretization errors using Richardson extrapolation. It is shown that grid refinement and the number of employed Monte Carlo particles have the largest influence on the result, in agreement with similar studies of an ITER deuterium case. For the first time, numerical error bars on the entire simulated target profiles are determined showing that the largest numerical error is 17.9%, mainly due to the plasma grid discretization. On top, also numerical errors on simulated neutral pressures are investigated in detail, for which the statistical error is dominant. The analysis demonstrates which setup is needed to keep numerical errors limited: the SOLPS-ITER averaging procedure should be employed including enough EIRENE particles, and the involved grids should be sufficiently refined to reduce discretization errors.

Boeyaert, Dieter (ORCID:0000000309208660)↗
NASA NTRS1995

Precision Cleaning Verification of Fluid Components by Air/Water Impingement and Total Carbon Analysis

NASA personnel at Kennedy Space Center's Material Science Laboratory have developed new environmentally sound precision cleaning and verification techniques for systems and components found at the center. This technology is required to replace existing methods traditionally employing CFC-113. The new patent-pending technique of precision cleaning verification is for large components of cryogenic fluid systems. These are stainless steel, sand cast valve bodies with internal surface areas ranging from 0.2 to 0.9 m(exp 2). Extrapolation of this technique to components of even larger sizes (by orders of magnitude) is planned. Currently, the verification process is completely manual. In the new technique, a high velocity, low volume water stream impacts the part to be verified. This process is referred to as Breathing Air/Water Impingement and forms the basis for the Impingement Verification System (IVS). The system is unique in that a gas stream is used to accelerate the water droplets to high speeds. Water is injected into the gas stream in a small, continuous amount. The air/water mixture is then passed through a converging-diverging nozzle where the gas is accelerated to supersonic velocities. These droplets impart sufficient energy to the precision cleaned surface to place non-volatile residue (NVR) contaminants into suspension in the water. The sample water is collected and its NVR level is determined by total organic carbon (TOC) analysis at 880 C. The TOC, in ppm carbon, is used to establish the NVR level. A correlation between the present gravimetric CFC-113 NVR and the IVS NVR is found from experimental sensitivity factors measured for various contaminants. The sensitivity has the units of ppm of carbon per mg-ft(exp 2) of contaminant. In this paper, the equipment is described and data are presented showing the development of the sensitivity factors from a test set including four NVR's impinged from witness plates of 0.05 to 0.75 m(exp 2).

Barile, Ronald G.↗
NASA NTRS1977

Analysis of satellite derived winds for April 24, 1975

Cumulus cloud tracking based on short-interval GOES satellite imagery provided a representation of the low-level wind fields existing prior to the outbreak of severe thunderstorms. Through use of the satellite-derived winds and surface mixing ratios, moisture convergence fields were estimated. The effect of the spatial distribution of trackable clouds on the computed dynamic fields was also investigated.

Maddox, R. A.↗
NASA NTRS1991

Time domain convergence properties of Lyapunov stable penalty methods

Linear hyperbolic partial differential equations are analyzed using standard techniques to show that a sequence of solutions generated by the Liapunov stable penalty equations approaches the solution of the differential-algebraic equations governing the dynamics of multibody problems arising in linear vibrations. The analysis does not require that the system be conservative and does not impose any specific integration scheme. Variational statements are derived which bound the error in approximation by the norm of the constraint violation obtained in the approximate solutions.

Kurdila, A. J.↗
NASA NTRS1993

Analysis of the Command and Control Segment (CCS) attitude estimation algorithm

This paper categorizes the qualitative behavior of the Command and Control Segment (CCS) differential correction algorithm as applied to attitude estimation using simultaneous spin axis sun angle and Earth cord length measurements. The categories of interest are the domains of convergence, divergence, and their boundaries. Three series of plots are discussed that show the dependence of the estimation algorithm on the vehicle radius, the sun/Earth angle, and the spacecraft attitude. Common qualitative dynamics to all three series are tabulated and discussed. Out-of-limits conditions for the estimation algorithm are identified and discussed.

Stockwell, Catherine↗
NASA NTRSDate not supplied

Exploring the Geochemistry of Tholeiitic Basalts and Hyaloclastites Formed in Submarine Volcanoes for Comparison With Venus Lavas

Underwater basaltic lavas exhibit distinct rock textures, mineralogy and morphological differences compared to those formed on land, primarily due to rapid cooling caused by contact with cold seawater instead of air and to high pressures – similar to Venus. The geomorphological similarities between terrestrial submarine volcanoes and some morphological classes of Venusian volcanoes have been investigated, but a comprehensive comparison in geochemical terms has yet to be conducted. Only three missions to Venus (Venera-13, Venera-14, and Vega-2) were able to report bulk composition of Venusian rocks, while other four missions (Venera-8, Venera-9, Venera-10, and Vega-1) provided limited geochemical information. Using X-Ray fluorescent (XRF) analysis, we saw that rocks analyzed on the surface of Venus are tholeiitic and alkalic basalts. This study, as part of the “Analogs for VENus's GEologically Recent Surfaces” (AVENGERS) initiative, aims to examine the possibility of seamounts/submarine volcanoes lavas serving as analogues to Venus lava geochemistry and morphology/formational processes. By sampling and analyzing submarine lavas from three distinct tectonic settings on Earth - hot-spot, convergent margins, and divergent margins - we will assess in which tectonic context the analogy with Venus is most comparable. Fieldwork will be conducted in further detail for submarine volcanoes located in Sicily, Hawaiʻi, and the Southern Pacific (Fiji, Vanuatu, Tonga, Samoa); we will use previously collected samples where these submarine areas are not accessible. Sampling efforts will target areas likely to yield submarine tholeiitic lavas and/or hyaloclastites. Samples collected during fieldwork will undergo detailed whole-rock geochemical analysis, including XRF, X-Ray diffraction (XRD) and Inductively Coupled Plasma Mass Spectrometry (ICPMS). Thin sections will be used to constrain differences in morphology/textures. The obtained geochemical data will then be compared with existing data on Venusian lava geochemistry, allowing for a comprehensive evaluation of potential analogies between submarine lavas from Earth and volcanic rocks on Venus. By elucidating these similarities and pinpointing specific tectonic contexts, our study has the potential to enhance our understanding of Venusian geology and the possibility of ephemeral plate tectonics on the planet.

N Mari↗
DOE OSTI2024

Distributed optimization for multi-commodity urban traffic control

A distributed method for concurrent traffic signal and routing control of traffic networks is proposed. The method is based on the multi-commodity store-and-forward model, in which the destinations are the commodities. The system benefits from the communication between vehicles and infrastructure, providing optimal signal timings to intersections and routes to vehicles on a link-by-link basis. Using the augmented Lagrangian to model the constraints into the objective, the baseline centralized problem is decomposed into a set of objective-coupled subproblems, one for each intersection, enabling the solution to be computed by a distributed- gradient projection algorithm. Further, the intersection agents only need to communicate and coordinate with neighboring intersections to ensure convergence to the optimal solution while tolerating suboptimal iterations that offer more flexibility, unlike other distributed approaches. Through microsimulation, we demonstrate the effectiveness of the proposed algorithm in traffic networks with time-varying demand. Computational analysis shows that the distributed problem is suitable for real-time applications. A robustness analysis show that the distributed formulation enables a graceful degradation of the system in case of failure.

Augmented Lagrangian↗
DOE OSTI2016

Ordered Self-Similar Patterns in Anisotropic Stochastic Growth

Here, we propose an anisotropic stochastic growth model to rationalize the anisotropic self-assembly of supramolecules to form elongated two-dimensional ribbon structures in a recent experiment. The model exhibits distinct growth scenarios that are critically controlled by the ratio of the transverse and the longitudinal growth rate. In the regime of suppressed transverse growth, the model generates the experimentally observed elongated structures through layer-by-layer growing. We further observe the convergence of rough clusters toward smooth regular elliptic patterns by averaging over a number of independent growth processes. Remarkably, these resulting elliptic clusters are self-similar in each instantaneous moment in the growth process. Statistical analysis suggests that the realization of such ordered patterns does not rely on the delicate coordination of different parts in the cluster growth. The self-similarity phenomenon derived from this idealized model may have wider implications, notably in the designed clustering of various elementary building blocks with anisotropic interactions.

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

Tensor network simulations of quasi-GPDs in the massive Schwinger model

Generalized parton distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED 2 ), within the Hamiltonian formulation of lattice field theory. Quasidistributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasiparton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Dynamic Power Distribution System Management With a Locally Connected Communication Network

Coordinated optimization and control of distribution-level assets can enable a reliable and optimal integration of massive amount of distributed energy resources (DERs) and facilitate distribution system management (DSM). Accordingly, the objective is to coordinate the power injection at the DERs to maintain certain quantities across the network, e.g., voltage magnitude, line flows, or line losses, to be close to a desired profile. By and large, the performance of the DSM algorithms has been challenged by two factors: i) the possibly non-strongly connected communication network over DERs that hinders the coordination; ii) the dynamics of the real system caused by the DERs with heterogeneous capabilities, time-varying operating conditions, and real-time measurement mismatches. In this paper, we investigate the modeling and algorithm design and analysis with the consideration of these two factors. In particular, a game theoretic characterization is first proposed to account for a locally connected communication network over DERs, along with the analysis of the existence and uniqueness of the Nash equilibrium (NE) therein. To achieve the equilibrium in a distributed fashion, a projected-gradient-based asynchronous DSM algorithm is then advocated. The algorithm performance, including the convergence speed and the tracking error, is analytically guaranteed under the dynamic setting. In conclusion, extensive numerical tests on both synthetic and realistic cases corroborate the analytical results derived.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

LACEwING: A New Moving Group Analysis Code

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars.

79 ASTRONOMY AND ASTROPHYSICS↗