SEARCH · Engineering Papers

Results for “numerical approximation & 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.

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Modifying the Asynchronous Jacobi Method for Data Corruption Resilience

Moving scientific computation from high-performance computing (HPC) and cloud computing (CC) environments to devices on the edge, i.e., physically near instruments of interest, has received tremendous interest in recent years. Such edge computing environments can operate on data in situ, offering enticing benefits over data aggregation to HPC and CC facilities that include avoiding costs of transmission, increased data privacy, and real-time data analysis. Because of the inherent unreliability of edge computing environments, new fault-tolerant approaches must be developed before the benefits of edge computing can be realized. Motivated by algorithm-based fault tolerance, a variant of the asynchronous Jacobi (ASJ) method is developed that achieves resilience to data corruption by rejecting solution approximations from neighbor devices according to a bound derived from convergence theory. Numerical results on a two-dimensional Poisson problem show that the new rejection criterion, along with a novel approximation to the shortest path length on which the criterion depends, restores convergence for the ASJ variant in the presence of certain types data corruption. Numerical results are obtained for when the singular values in the analytic bound are approximated. Additional linear systems are also explored, one with a more dense sparsity pattern and one that includes advection. All results indicate that successful resilience to data corruption depends on whether the bound tightens fast enough to reject corrupted data before the iteration evolution deviates significantly from that predicted by the convergence theory defining the bound. This observation generalizes to future work on algorithm-based fault tolerance for other asynchronous algorithms, including upcoming approaches that leverage Krylov subspaces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Code Coverage Status of the ARC Code RCT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. REBUS plays a pivotal role in the ARC system as the primary fuel cycle analysis capability for fast reactor problems. Over its 60 year history, ARC software usage with REBUS has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. The RCT code is a later addition and uses the REBUS restart files to define its input. The RCT code was built to provide pin depletion details on EBR-II models and thus many features of RCT were specifically tailored to the needs of EBR-II models. Additional approximations were invoked which are likely only valid for the EBR-II reactor and the particular fuel management that was done for it. The purpose of the present work is to identify a set of test problems for RCT and assess the code coverage for those test problems. The goal is to document what parts of the existing RCT code are touched by the set of test problems and which are not. Because no detailed verification work has been done on RCT, the existing regression testing suite was chosen for the code coverage assessment. The code coverage analysis of RCT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of RCT. The detailed coverage tables are given for each part of RCT. As will be discussed and shown, some parts of the RCT capability that are known to be used by the EBR-II analysis work are not tested by the regression testing suite. These aspects should be resolved before major source code changes are taken for the RCT software. Because REBUS and DIF3D are not subroutines of RCT, the coverage changes in both of those codes is not altered by RCT. The same is true for all of the modules of DIF3D that are used by RCT such as SYSLIB and SEGLIB.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1974

Post-flight differential correction analysis using Vinti's spheroidal method for the small astronomy satellite orbit

The results of an intensive analysis of a differential orbit improvement method utilizing observational data for a 550-kilometer altitude, near-circular, near-equatorial satellite orbit are presented. Observations of the Small Astronomy Satellite (SAS-1) are in the form of direction cosines as measured at two ground interferometer tracking stations near the Equator during the first 22 orbital revolutions (approximately 37 hours) after launch of the spacecraft. Numerical results, in both tabular and graphical form, are displayed for numerous iterated fittings of various observational arcs by differential correction of the orbital elements. Parameters varied in these comparative cases include the time duration of the observational data block, the number of pairs of direction cosine data and the number of tracking station passes included in the solution, the distribution of such passes between the two available tracking stations, and the acceptance criterion for the observational residuals in the least squares fitting procedure. It was found that three observational pairs of direction cosine data, the minimum number possible for a uniquely determined solution in theory, are sufficient to promote covergence to an accurate solution, if properly selected.

Walden, H.↗

NASA NTRS1987

An evaluation of higher-order modal methods for calculating transient structural response

The present study evaluates a higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗

NASA NTRS1999

Evidence for Methane Segregation at the Surface of Pluto

In May 1995, a set of spectrophotometric curves of the system Pluto-Charon was recorded with the UKIRT telescope equipped with the spectrometer CGS4. As for the previous observations, the spectra cover a part of the near infrared range, between 1.4 and 2.55 micrometers, but with a higher resolution of approximately 700. In both the 1992 and 1995 data, the existence of solid methane is confirmed by numerous absorption bands, and the carbon monoxide and the nitrogen ices are identified by their respective signatures at 2.35 and 2.15 um. The solid nitrogen seems to be the principal icy component and forms a matrix in which the CH4 and CO molecules are diluted. However a spectroscopic analysis of the 1995 observations indicates that pure methane may coexist with its diluted phase in N2. In order to derive the horizontal and vertical distribution of these different species and to obtain some quantitative information about their characteristics, we have modeled the spectrum of May 15 that corresponds to the maximum of Pluto's visible light curve. This was achieved by means of a radiative transfer algorithm dealing with compact and stratified media. Among the various representations we have tested to describe the surface of Pluto, only a geographical mixture of three distinct units explains all the significant structures of the analyzed spectrum. The first unit is a thin granular layer of pure CH4 covering a compact polycrystalline substratum of N2-CH4-CO, which are in a molecular mixture (concentrations of and CO of the order of 0.45%, 0.1-0.2% respectively). It covers about 70% of the observed area and corresponds to volatile deposits that are sublimating under solar illumination. The second unit is either (a) a single thick layer of pure granular methane or (b) a unit similar to the first unit but with the two components inverted (i.e. with CH4 forming a substratum and the N2-CH4-CO mixture a superficial layer of fine grains). Covering 20% of the surface, it represents some old surfaces that have been sublimated for a long time, and eventually recovered later by very small amounts of fresh deposits of the molecular mixture N2-CH4-CO. Finally, the third unit may result from the condensation of very fine grains of nearly pure N2. It covers the remainder of the surface (about 10%). All these results allow a better understanding of the processes of deposition, metamorphism, sublimation and transport affecting the different ices detected on Pluto during its climatic cycles.

Doute, S.↗

NASA NTRS1996

Theoretical and Numerical Investigation of Radiative Extinction of Diffusion Flames

The influence of soot radiation on diffusion flames was investigated using both analytical and numerical techniques. Soot generated in diffusion flames dominate the flame radiation over gaseous combustion products and can significantly lower the temperature of the flame. In low gravity situations there can be significant accumulation of soot and combustion products in the vicinity of the primary reaction zone owing to the absence of any convective buoyant flow. Such situations may result in substantial suppression of chemical activities in a flame, and the possibility of a radiative extinction may also be anticipated. The purpose of this work was to not only investigate the possibility of radiative extinction of a diffusion flame but also to qualitatively and quantitatively analyze the influence of soot radiation on a diffusion flame. In this study, first a hypothetical radiative loss profile of the form of a sech(sup 2) was assumed to influence a pure diffusion flame. It was observed that the reaction zone can, under certain circumstances, move through the radiative loss zone and locate itself on the fuel side of the loss zone contrary to our initial postulate. On increasing the intensity and/or width of the loss zone it was possible to extinguish the flame, and extinction plots were generated. In the presence of a convective flow, however, the movement of the temperature and reaction rate peaks indicated that the flame behavior is more complicated compared to a pure diffusional flame. A comprehensive model of soot formation, oxidation and radiation was used in a more involved analysis. The soot model of Syed, Stewart and Moss was used for soot nucleation and growth and the model of Nagle and Strickland-Constable was used for soot oxidation. The soot radiation was considered in the optically thin limit. An analysis of the flame structure revealed that the radiative loss term is countered both by the reaction term and the diffusion term. The essential balance for the soot volume fraction was found to be between the processes of soot convection and soot growth. Such a balance yielded to analytical treatment and the soot volume fraction could be expressed in the form of an integral. The integral was evaluated using two approximate methods and the results agreed very well with the numerical solutions for all cases examined.

Ray, Anjan↗

DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

NASA NTRS1995

The analysis of control trajectories using symbolic and database computing

This final report comprises the formal semi-annual status reports for this grant for the periods June 30-December 31, 1993, January 1-June 30, 1994, and June 1-December 31, 1994. The research supported by this grant is broadly concerned with the symbolic computation, mixed numeric-symbolic computation, and database computation of trajectories of dynamical systems, especially control systems. A review of work during the report period covers: trajectories and approximating series, the Cayley algebra of trees, actions of differential operators, geometrically stable integration algorithms, hybrid systems, trajectory stores, PTool, and other activities. A list of publications written during the report period is attached.

Grossman, Robert↗

DOE OSTI2020

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗

NASA NTRS1988

Analysis of long term heart rate variability: methods, 1/f scaling and implications

The use of spectral techniques to quantify short term heart rate fluctuations on the order of seconds to minutes has helped define the autonomic contributions to beat-to-beat control of heart rate. We used similar techniques to quantify the entire spectrum (0.00003-1.0 Hz) of heart rate variability during 24 hour ambulatory ECG monitoring. The ECG from standard Holter monitor recordings from normal subjects was sampled with the use of a phase locked loop, and a heart rate time series was constructed at 3 Hz. Frequency analysis of the heart rate signal was performed after a nonlinear filtering algorithm was used to eliminate artifacts. A power spectrum of the entire 24 hour record revealed power that was inversely proportional to frequency, 1/f, over 4 decades from 0.00003 to 0.1 Hz (period approximately 10 hours to 10 seconds). Displaying consecutive spectra calculated at 5 minute intervals revealed marked variability in the peaks at all frequencies throughout the 24 hours, probably accounting for the lack of distinct peaks in the spectra of the entire records.

Non-NASA Center↗

DOE OSTI2023

NOW-23: The 2023 National Offshore Wind Data Set

17 WIND ENERGY↗

DOE OSTI2023

NOW-23: the 2023 National Offshore Wind Data Set

In this report, we present the latest wind resource data specifically tailored for offshore regions in the United States. The data set, known as the 2023 National Offshore Wind data set (NOW-23), has been developed by the National Renewable Energy Laboratory (NREL) and its partners, and surpasses the previous resource data set, the Wind Integration National Dataset (WIND) Toolkit, which was released approximately ten years ago, in its offshore component. The WIND Toolkit has been widely utilized by stakeholders involved in wind resource assessments across the continental United States. However, with significant advancements in numerical weather prediction modeling over the past decade, the NOW-23 data set incorporates the latest research and development progress to provide stakeholders with an updated and cutting-edge resource for offshore wind analysis. The NOW-23 data set is created using the Weather Research and Forecasting (WRF) model and its output is available, as for its predecessor, at 5-minute time resolution and 2-kilometer horizontal spatial resolution. However, the NOW-23 data set improves upon the WIND Toolkit through: 1. A modeling period of at least 20 years (and as long as 23 years in selected regions), starting in 2000 (compared to the 7-year 2007–2013 modeling period in the WIND Toolkit). 2. For several offshore regions, a region-specific sensitivity analysis, driven by an ensemble of WRF simulations, to assess the most adequate region-specific WRF setup. 3. An updated WRF model, from Version 3.4 used in the WIND Toolkit to Version 4.2.1 used for the NOW-23 data set, which incorporates significant research advancements. 4. The use of the state-of-the-art reanalysis product ERA5 (which supersedes the older ERA-Interim used in the WIND Toolkit) to provide atmospheric forcing at the WRF domain boundaries. Figure 1 shows the NOW-23 mean wind speed at 160 m above sea level (asl), which we use as proxy for hub-height of a commercial offshore wind turbine in this report across all regions.

17 WIND ENERGY↗

NASA NTRS1980

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗

NASA NTRS1982

Heat generation in aircraft tires under free rolling conditions

A method was developed for calculating the internal temperature distribution in an aircraft tire while free rolling under load. The method uses an approximate stress analysis of each point in the tire as it rolls through the contact patch, and from this stress change the mechanical work done on each volume element may be obtained and converted into a heat release rate through a knowledge of material characteristics. The tire cross-section is then considered as a body with internal heat generation, and the diffusion equation is solved numerically with appropriate boundary conditions of the wheel and runway surface. Comparison with data obtained with buried thermocouples in tires shows good agreement.

Clark, S. K.↗

NASA NTRS1983

Heat generation in aircraft tires

Clark, S. K.↗

NASA NTRS1984

Optimum design sensitivity based on approximation concepts and dual methods

An optimum design sensitivity analysis capability is reported which exploits the approximation concepts-dual method formulation of the minimum weight structural sizing problem. An efficient iterative solution technique is used to facilitate determination of sensitivity derivatives for both primal and dual variables. Estimates on the useful range of parameter perturbations, over which the optimum design sensitivity projections can be expected to yield satisfactory revised optimum designs, are also obtained. Numerical results for several example problems will be presented to illustrate the effectiveness of the capability reported.

Schmit, L. A.↗

NASA NTRS1985

Heat generation in aircraft tires

Clark, S. K.↗

NASA NTRS1994

Preferred latitudes of the intertropical convergence zone

The latitude preference of the intertropical convergence zone (ITCZ) is examined on the basis of observations, theory, and a modeling analysis. Observations show that convection is enhanced at latitudes of about 4 deg to 10 deg relative to the equator, even in regions where the sea surface temperature (SST) is maximum on the equator. Both linear shallow-water theory and a moist primitive equation model suggest a new explanation for the off-equatorial latitude preference of the ITCZ that requires neither the existence of zonally propagating disturbances nor an off-equatorial maximum in SST. The shallow-water theory indicates that a finite-width, zonally oriented, midtropospheric heat source (i.e., an ITCZ) produces the greatest local low-level convergence when placed a finite distance away from the equator. This result suggests that an ITCZ is most likely to be supported via low-level convergence of moist energy when located at these "preferred" latitudes away from the equator. For a plausible range of heating widths and damping parameters, the theoretically predicted latitude is approximately equal to the observed position (s) of the ITCZ (s). Analysis with an axially symmetric, moist, primitive equation model indicates that when the latent heating field is allowed to be determined internally, a positive feedback develops between the midtropospheric latent heating and the low-level convergence, with the effect of enhancing the organization of convection at latitudes of about 4 deg to 12 deg. Numerical experiments show that (1) two peaks in convective precipitation develop straddling the equator when the SST maximum is located on the equator; (2) steady ITCZ-like structures form only when the SST maximum is located away from the equator; and (3) peaks in convection can develop away from the maximum in SST, with a particular preference for latitudes of about 4 deg to 12 deg, even in the ('cold') hemisphere without the SST maximum. The relationship between this mechanism and earlier theories is discussed, as are implications for the coupled ocean-atmosphere system and the roles played by midlevel latent heating and SST gradients in forcing the low-level atmospheric circulation in the tropics.

Waliser, Duane E.↗