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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The arbitrary-order virtual element method for linear elastodynamics models. Convergence, stability and dispersion-dissipation analysis.

We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under $\textit{h}$-refinement in both the energy and the $L^2$ norms, and optimal error estimates under $\textit{p}$-refinement in the energy norm. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally observed under p-refinement. Finally, we present a dispersion-dissipation analysis for both the semi-discrete and fully-discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion-dissipation properties.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2020

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method↗
DOE OSTI2023

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
NASA NTRS1976

Convergence of an iterative procedure for large-scale static analysis of structural components

The paper proves convergence of an iterative procedure for calculating the deflections of built-up component structures which can be represented as consisting of a dominant, relatively stiff primary structure and a less stiff secondary structure, which may be composed of one or more substructures that are not connected to one another but are all connected to the primary structure. The iteration consists in estimating the deformation of the primary structure in the absence of the secondary structure on the assumption that all mechanical loads are applied directly to the primary structure. The j-th iterate primary structure deflections at the interface are imposed on the secondary structure, and the boundary loads required to produce these deflections are computed. The cycle is completed by applying the interface reaction to the primary structure and computing its updated deflections. It is shown that the mathematical condition for convergence of this procedure is that the maximum eigenvalue of the equation relating primary-structure deflection to imposed secondary-structure deflection be less than unity, which is shown to correspond with the physical requirement that the secondary structure be more flexible at the interface boundary.

Austin, F.↗
NASA NTRS1977

Moisture convergence from a combined mesoscale moisture analysis and wind field for 24 April 1975

Precipitable water values inferred from the Vertical Temperature Profile Radiometer data of the polar orbiting NOAA-4 satellite are used in conjunction with wind-field analyses obtained from Synchronous Meteorological Satellite visible-channel data to study the moisture convergence in the boundary layer immediately preceding a storm. This combination of data simulates the information that will be available from the Visible and Infrared Spin-Scan Radiometer on board the GOES-D satellite, which is scheduled to begin operation in the 1980s. Serviceable representations of boundary layer flow are developed through analysis of the satellite infrared cumulus velocities, although the flow representations are not exactly located in the vertical.

Negri, A. J.↗
NASA NTRS2008

Deflection Analysis of the Space Shuttle External Tank Door Drive Mechanism

Upon observing an abnormal closure of the Space Shuttle s External Tank Doors (ETD), a dynamic model was created in MSC/ADAMS to conduct deflection analyses of the Door Drive Mechanism (DDM). For a similar analysis, the traditional approach would be to construct a full finite element model of the mechanism. The purpose of this paper is to describe an alternative approach that models the flexibility of the DDM using a lumped parameter approximation to capture the compliance of individual parts within the drive linkage. This approach allows for rapid construction of a dynamic model in a time-critical setting, while still retaining the appropriate equivalent stiffness of each linkage component. As a validation of these equivalent stiffnesses, finite element analysis (FEA) was used to iteratively update the model towards convergence. Following this analysis, deflections recovered from the dynamic model can be used to calculate stress and classify each component s deformation as either elastic or plastic. Based on the modeling assumptions used in this analysis and the maximum input forcing condition, two components in the DDM show a factor of safety less than or equal to 0.5. However, to accurately evaluate the induced stresses, additional mechanism rigging information would be necessary to characterize the input forcing conditions. This information would also allow for the classification of stresses as either elastic or plastic.

Tosto, Michael A.↗
DOE OSTI2025

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1984

A two-dimensional analysis of laser heat addition in converging nozzles

The two-dimensional equations of motion describing the interaction between a laser beam and a flowing gas are considered. An implicit numerical scheme is used to solve these equations for unchoked flow through a converging-diverging nozzle. Separate grids are used for the fluid dynamics and the radiation equations. The effects of beam focusing and cross-beam intensity profiles are included. The calculations are based upon real gas properties for all quantities except the gas absorptivity, which is taken as a constant. The solutions contain the expected hot central core region with cool gas near the walls. This results in steep temperature gradients in both the streamwise and cross-stream directions. The absorption zone acts as a blockage in the nozzle causing a nonuniform velocity profile at the inlet and an overall decrease in mass flow. The absorption region also forces the streamlines to move away from the axis of symmetry, although this effect is not strong.

Molvik, G. A.↗
DOE OSTI2012

Acceleration of Monte Carlo Criticality Calculations Using Deterministic-Based Starting Sources

A new automatic approach that uses approximate deterministic solutions for providing the starting fission source for Monte Carlo eigenvalue calculations was evaluated in this analysis. By accelerating the Monte Carlo source convergence and decreasing the number of cycles that has to be skipped before the tally's estimation, this approach was found to increase the efficiency of the overall simulation, even with the inclusion of the extra computational time required by the deterministic calculation. This approach was also found to increase the reliability of the Monte Carlo criticality calculations of loosely coupled systems because the use of the better starting source reduces the likelihood of producing an undersampled k eff due to the inadequate source convergence. The efficiency improvement was demonstrated using two of the standard test problems devised by the OECD/NEA Expert Group on Source Convergence in Criticality-Safety Analysis to measure the source convergence in Monte Carlo criticality calculations. For a fixed uncertainty objective, this approach increased the efficiency of the overall simulation by factors between 1.2 and 3 depending on the difficulty of the source convergence in these problems. The reliability improvement was demonstrated in a modified version of the 'k eff of the world' problem that was specifically designed to demonstrate the limitations of the current Monte Carlo power iteration techniques. For this problem, the probability of obtaining a clearly undersampled k eff decreased from 5% with a uniform starting source to zero with a deterministic starting source when batch sizes with more than 15,000 neutron/cycle were used.

42 ENGINEERING↗
DOE OSTI2022

A deterministic gradient-based approach to avoid saddle points

Abstract Loss functions with a large number of saddle points are one of the major obstacles for training modern machine learning (ML) models efficiently. First-order methods such as gradient descent (GD) are usually the methods of choice for training ML models. However, these methods converge to saddle points for certain choices of initial guesses. In this paper, we propose a modification of the recently proposed Laplacian smoothing gradient descent (LSGD) [Osher et al., arXiv:1806.06317 ], called modified LSGD (mLSGD), and demonstrate its potential to avoid saddle points without sacrificing the convergence rate. Our analysis is based on the attraction region, formed by all starting points for which the considered numerical scheme converges to a saddle point. We investigate the attraction region’s dimension both analytically and numerically. For a canonical class of quadratic functions, we show that the dimension of the attraction region for mLSGD is $\lfloor (n-1)/2\rfloor$ , and hence it is significantly smaller than that of GD whose dimension is $n-1$ .

Mathematics↗
DOE OSTI2025

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Cosmology from clustering, cosmic shear, CMB lensing, and cross correlations: combining Rubin observatory and Simons Observatory

ABSTRACT In the near future, the overlap of the Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) will present an ideal opportunity for joint cosmological data set analyses. In this paper, we simulate the joint likelihood analysis of these two experiments using six two-point functions derived from galaxy position, galaxy shear, and CMB lensing convergence fields. Our analysis focuses on realistic noise and systematics models and we find that the dark energy Figure of merit (FoM) increases by 53 per cent (92 per cent) from LSST-only to LSST+SO in Year 1 (Year 6). We also investigate the benefits of using the same galaxy sample for both clustering and lensing analyses, and find the choice improves the overall signal to noise by ${\sim}30{-}40{{\ \rm per\ cent}}$, which significantly improves the photo-z calibration and mildly improves the cosmological constraints. Finally, we explore the effects of catastrophic photo-z outliers finding that they cause significant parameter biases when ignored. We develop a new mitigation approach termed ‘island model’, which corrects a large fraction of the biases with only a few parameters while preserving the constraining power.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Statistical Symbolic Execution with Informed Sampling

Symbolic execution techniques have been proposed recently for the probabilistic analysis of programs. These techniques seek to quantify the likelihood of reaching program events of interest, e.g., assert violations. They have many promising applications but have scalability issues due to high computational demand. To address this challenge, we propose a statistical symbolic execution technique that performs Monte Carlo sampling of the symbolic program paths and uses the obtained information for Bayesian estimation and hypothesis testing with respect to the probability of reaching the target events. To speed up the convergence of the statistical analysis, we propose Informed Sampling, an iterative symbolic execution that first explores the paths that have high statistical significance, prunes them from the state space and guides the execution towards less likely paths. The technique combines Bayesian estimation with a partial exact analysis for the pruned paths leading to provably improved convergence of the statistical analysis. We have implemented statistical symbolic execution with in- formed sampling in the Symbolic PathFinder tool. We show experimentally that the informed sampling obtains more precise results and converges faster than a purely statistical analysis and may also be more efficient than an exact symbolic analysis. When the latter does not terminate symbolic execution with informed sampling can give meaningful results under the same time and memory limits.

Reliability↗
DOE OSTI2012

IMPLODING IGNITION WAVES. I. ONE-DIMENSIONAL ANALYSIS

We show that converging spherical and cylindrical shock waves may ignite a detonation wave in a combustible medium, provided the radius at which the shocks become strong exceeds a critical radius, R{sub crit}. An approximate analytic expression for R{sub crit} is derived for an ideal gas equation of state and a simple (power-law-Arrhenius) reaction law, and shown to reproduce the results of numerical solutions. For typical acetylene-air experiments we find R{sub crit} {approx} 100 {mu}m (spherical) and R{sub crit} {approx} 1 mm (cylindrical). We suggest that the deflagration to detonation transition (DDT) observed in these systems may be due to converging shocks produced by the turbulent deflagration flow, which reaches sub- (but near) sonic velocities on scales >>R{sub crit}. Our suggested mechanism differs from that proposed by Zel'dovich et al., in which a fine-tuned spatial gradient in the chemical induction time is required to be maintained within the turbulent deflagration flow. Our analysis may be readily extended to more complicated equations of state and reaction laws. An order of magnitude estimate of R{sub crit} within a white dwarf at the pre-detonation conditions believed to lead to Type Ia supernova explosions is 0.1 km, suggesting that our proposed mechanism may be relevant for DDT initiation in these systems. The relevance of our proposed ignition mechanism to DDT initiation may be tested by both experiments and numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Necessary and sufficient conditions for parameter convergence in adaptive control

Using Generalized Harmonic Analysis, a complete description of parameter convergence in Model Reference Adaptive Control (MRAC) is given in terms of the spectrum of the exogenous reference input signal. Roughly speaking, if the reference signal 'contains enough frequencies' then the parameter vector converges to its correct value. If not, it converges to an easily characterizable subspace in parameter space.

Boyd, Stephen↗
NASA NTRS1989

Multiple boundary condition testing error analysis

Techniques for interpreting data from multiple-boundary-condition (MBC) ground tests of large space structures are developed analytically and demonstrated. The use of MBC testing to validate structures too large to stand alone on the ground is explained; the generalized least-squares mass and stiffness curve-fitting methods typically applied to MBC test data are reviewed; and a detailed error analysis is performed. Consideration is given to sensitivity coefficients, covariance-matrix theory, the correspondence between test and analysis modes, constraints and step sizes, convergence criteria, and factor-analysis theory. Numerical results for a simple beam problem are presented in tables and briefly characterized. The improved error-updating capabilities of MBC testing are confirmed, and it is concluded that reasonably accurate results can be obtained using a diagonal covariance matrix.

Glaser, R. J.↗
NASA NTRS1975

Analysis and testing of two-dimensional slot nozzle ejectors with variable area mixing sections

Finite difference computer techniques have been used to calculate the detailed performance of air-to-air two-dimensional ejectors with symmetric variable area mixing sections and coaxial converging primary nozzles. The analysis of the primary nozzle assumed correct expansion of the flow and is suitable for subsonic and slightly supersonic velocity levels. The variation of the mixing section channel walls is assumed to be gradual so that the static pressure can be assumed uniform on planes perpendicular to the axis. A test program was run to provide two-dimensional ejector test data for verification of the computer analysis. A primary converging nozzle with a discharge geometry of 3.18 mm x 203 mm was supplied with 0.340 kg/sec of air at about 2.43 bar and 356 K. This nozzle was combined with two mixing section geometries with throat sizes of 31.8 mm x 203 mm and 47.6 mm x 203 mm and was tested at a total of 11 operating points.

Gilbert, G. B.↗
NASA NTRS2001

Analysis of Preconditioning and Relaxation Operators for the Discontinuous Galerkin Method Applied to Diffusion

The explicit stability constraint of the discontinuous Galerkin method applied to the diffusion operator decreases dramatically as the order of the method is increased. Block Jacobi and block Gauss-Seidel preconditioner operators are examined for their effectiveness at accelerating convergence. A Fourier analysis for methods of order 2 through 6 reveals that both preconditioner operators bound the eigenvalues of the discrete spatial operator. Additionally, in one dimension, the eigenvalues are grouped into two or three regions that are invariant with order of the method. Local relaxation methods are constructed that rapidly damp high frequencies for arbitrarily large time step.

Atkins, H. L.↗