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,009 records · Page 56

Combining Static Analysis and Model Checking for Software Analysis

We present an iterative technique in which model checking and static analysis are combined to verify large software systems. The role of the static analysis is to compute partial order information which the model checker uses to reduce the state space. During exploration, the model checker also computes aliasing information that it gives to the static analyzer which can then refine its analysis. The result of this refined analysis is then fed back to the model checker which updates its partial order reduction. At each step of this iterative process, the static analysis computes optimistic information which results in an unsafe reduction of the state space. However we show that the process converges to a fired point at which time the partial order information is safe and the whole state space is explored.

Brat, Guillaume↗

Comparison of Artificial Compressibility Methods

Various artificial compressibility methods for calculating three-dimensional, steady and unsteady, laminar and turbulent, incompressible Navier-Stokes equations are compared in this work. Each method is described in detail along with appropriate physical and numerical boundary conditions. Analysis of well-posedness and numerical solutions to test problems for each method are provided. A comparison based on convergence behavior, accuracy, stability and robustness is used to establish the relative positive and negative characteristics of each method.

Kiris, Cetin↗

Assessment of Radiative Heating Uncertainty for Hyperbolic Earth Entry

This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.

Johnston, Christopher O.↗

Evaluating ESA CCI Soil Moisture in East Africa

To assess growing season conditions where ground based observations are limited or unavailable, food security and agricultural drought monitoring analysts rely on publicly available remotely sensed rainfall and vegetation greenness. There are also remotely sensed soil moisture observations from missions like the European Space Agency (ESA) Soil Moisture and Ocean Salinity (SMOS) and NASAs Soil Moisture Active Passive (SMAP), however these time series are still too short to conduct studies that demonstrate the utility of these data for operational applications, or to provide historical context for extreme wet or dry events. To promote the use of remotely sensed soil moisture in agricultural drought and food security monitoring, we use East Africa as a case study to evaluate the quality of a 30+ year time series of merged active-passive microwave soil moisture from the ESA Climate Change Initiative (CCI-SM). Compared to the Normalized Difference Vegetation index (NDVI) and modeled soil moisture products, we found substantial spatial and temporal gaps in the early part of the CCI-SM record, with adequate data coverage beginning in 1992. From this point forward, growing season CCI-SM anomalies were well correlated (R greater than 0.5) with modeled, seasonal soil moisture, and in some regions, NDVI. We use correlation analysis and qualitative comparisons at seasonal time scales to show that remotely sensed soil moisture can add information to a convergence of evidence framework that traditionally relies on rainfall and NDVI in moderately vegetated regions.

Amy McNally↗

Design of optimally smoothing multi-stage schemes for the Euler equations

In this paper, a method is developed for designing multi-stage schemes that give optimal damping of high-frequencies for a given spatial-differencing operator. The objective of the method is to design schemes that combine well with multi-grid acceleration. The schemes are tested on a nonlinear scalar equation, and compared to Runge-Kutta schemes with the maximum stable time-step. The optimally smoothing schemes perform better than the Runge-Kutta schemes, even on a single grid. The analysis is extended to the Euler equations in one space-dimension by use of 'characteristic time-stepping', which preconditions the equations, removing stiffness due to variations among characteristic speeds. Convergence rates independent of the number of cells in the finest grid are achieved for transonic flow with and without a shock. Characteristic time-stepping is shown to be preferable to local time-stepping, although use of the optimally damping schemes appears to enhance the performance of local time-stepping. The extension of the analysis to the two-dimensional Euler equations is hampered by the lack of a model for characteristic time-stepping in two dimensions. Some results for local time-stepping are presented.

Van Leer, Bram↗

Square plate with clamped edges under normal pressure producing large deflections

A theoretical analysis is given for the stresses and deflections of a square plate with clamped edges under normal pressure producing large deflections. Values of the bending stress and membrane stress at the center of the plate and at the midpoint of the edge are given for center deflections up to 1.9 times the plate thickness. The shape of the deflected surface is given for low pressures and for the highest pressure considered. Convergence of the solution is considered and it is estimated that the possible error is less than 2 percent. The results are compared with the only previous approximate analysis known to the author and agree within 5 percent. They are also shown to compare favorably with the known exact solutions for the long rectangular plate and the circular plate.

PLATES-PRESSURE↗

Analysis of physics-based preconditioning for single-phase subchannel equations

The (single-phase) subchannel approximations are used throughout nuclear engineering to provide an efficient flow simulation because the computational burden is much smaller than for computational fluid dynamics (CFD) simulations, and empirical relations have been developed and validated to provide accurate solutions in appropriate flow regimes. Here, the subchannel equations have been recast in a residual form suitable for a multi-physics framework. The Eigen spectrum of the Jacobian matrix, along with several potential physics-based preconditioning approaches, are evaluated, and the the potential for improved convergence from preconditioning is assessed. The physics-based preconditioner options include several forms of reduced equations that decouple the subchannels by neglecting crossflow, conduction, and/or both turbulent momentum and energy exchange between subchannels. Eigen-scopy analysis shows that preconditioning moves clusters of eigenvalues away from zero and toward one. A test problem is run with and without preconditioning. Without preconditioning, the solution failed to converge using GMRES, but application of any of the preconditioners allowed the solution to converge. (authors)

97 MATHEMATICS AND COMPUTING↗

Lanczos algorithm for lattice QCD matrix elements

Recent work [M. L. Wagman, Lanczos, the transfer matrix, and the signal-to-noise problem, .] found that an analysis formalism based on the Lanczos algorithm allows energy levels to be extracted from Euclidean correlation functions with faster ground-state convergence than effective masses, convergent estimators for multiple states from a single correlator, and two-sided error bounds. After filtering out spurious eigenvalues and using outlier-robust estimators within a nested bootstrap framework, Lanczos estimators behave more like multistate fit results than effective masses—but without involving statistical fitting. We extend this formalism to the determination of matrix elements from three-point correlation functions and provide a physical picture of “spurious-state filtering” involving restriction to a Hermitian subspace. We demonstrate similar advantages for matrix elements as for spectroscopy through example applications to noiseless mock-data and (bare) forward matrix elements of the strange scalar current between both ground and excited states with the quantum numbers of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flux stabilization in neutron problems with fixed sources

Although critical core calculations are the most common in design and safety analysis, fixed source calculations are needed for specific applications, e.g. to compute ex-core detector response functions, to develop new methodologies for dilution and reload error accidents and more in general for all situations involving sub-critical shut-down states. It is well known that the source problem becomes difficult to be solved with core configuration close to criticality, i.e. with the multiplication factor approaching unity, for the occurrence of numerical ill-conditioning and very high number of iterations, possibly leading to failure in the flux convergence. In this work, the Wielandt Eigen-shift technique used in iterative methods of critical problems is developed for source problems too, in order to stabilize the solution. The mathematical basis and the proof of the convergence are discussed. Compared to the existing methods, this technique allows also more control to avoid singular behavior at inner iterations. Numerical tests with a 1D analytical benchmark are reported to prove the robustness of the technique. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On the rate of convergence as t→+∞ of the distributions of solutions to the stationary measure for the stochastic system of the Lorenz model describing a baroclinic atmosphere

The paper is concerned with a nonlinear system of partial differential equations with parameters which describes the two-layer quasi- solenoidal Lorenz model for a baroclinic atmosphere on a rotating two- dimensional sphere. The right-hand side of the system is perturbed by white noise. A unique stationary measure for the Markov semigroup defined by the solutions of the Cauchy problem for this problem is considered. An estimate for the rate of convergence of the distributions of all solutions in a certain class of this system to the unique stationary measure as t→+∞ is proposed. A similar result is obtained for the equation of a barotropic atmosphere and the two-dimensional Navier-Stokes equation. A comparative analysis with some of the available related results is given for the latter. Bibliography: 39 titles. (paper)

97 MATHEMATICS AND COMPUTING↗

Consistency and convergence for numerical radiation conditions

The problem of imposing radiation conditions at artificial boundaries for the numerical simulation of wave propagation is considered. Emphasis is on the behavior and analysis of the error which results from the restriction of the domain. The theory of error estimation is briefly outlined for boundary conditions. Use is made of the asymptotic analysis of propagating wave groups to derive and analyze boundary operators. For dissipative problems this leads to local, accurate conditions, but falls short in the hyperbolic case. A numerical experiment on the solution of the wave equation with cylindrical symmetry is described. A unified presentation of a number of conditions which have been proposed in the literature is given and the time dependence of the error which results from their use is displayed. The results are in qualitative agreement with theoretical considerations. It was found, however, that for this model problem it is particularly difficult to force the error to decay rapidly in time.

Hagstrom, Thomas↗

Study on the well-posedness, convergence and the stability of the semi-implicit upwind numerical solver for the multi-fluid model

The well-posedness, convergence and the stability of the two-fluid code has been studied for a long time. Most of the investigations concern the semi-implicit upwind solution scheme for the six equation two-fluid model such as used in RELAP5 3 or TRACE 21. Since the system code, SPACE 2, adopts one more field, a droplet field, it consists of nine equations (3 mass, 3 momentum and 3 energy balance equations) and thus more involved investigations are necessary to confirm the stability and convergence. For this objective, the old issue of the well-posedness, convergence and the stability is revisited and some general guidelines to develop a well-posed numerical multi-fluid model are derived as follows; (1) Hyperbolicity of the corresponding system of partial differential equations is not a necessary condition for the development of a numerical model for multi-phase flow, but whether or not it is hyperbolic can provide guidance relative to initial conditions, boundary conditions, and expected high frequency behavior of the model. (2) A necessary condition for a well-posed numerical model is stability in the von Neumann sense, i.e. growth factor less than 1.0 for the shortest wave-length, 2{Delta}x. (3) The smallest node size used for convergence studies should be of the order of the characteristic dimension of the average description, i.e. smaller nodes can be used so long as they do not result in unphysical growth of wave-lengths less than the characteristic dimension. The usual mathematical definition of convergence i.e. the behavior of the solution as the node size approaches zero, is not appropriate for the discrete averaged numerical model, since there is diminished physical meaning to behavior at wavelengths less than the characteristic dimension of the average description. Under these guidelines, dispersion analysis and von Neumann stability analysis are performed for the three field multi-fluid, semi-implicit, upwind numerical model to show that the necessary conditions for well-posedness are met. To study the non-linear stability and the convergence, various runs with the torus problem using the SPACE code are utilized to confirm the frequency cascading effects that augment non-linear stability. A robust mechanism for the flow regime change is also a very important factor for developing non-linearly stable and convergent code. A open pipe flow problem is also simulated to investigate the non-linear stability effects in a flow geometry more typical of real safety code applications. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Imploding and exploding shocks in negative ion degenerate plasmas

Imploding and exploding shocks are studied in nonplanar geometries for negative ion degenerate plasma. Deformed Korteweg de Vries Burgers (DKdVB) equation is derived by using reductive perturbation method. Two level finite difference scheme is used for numerical analysis of DKdVB. It is observed that compressive and rarefactive shocks are observed depending on the value of quantum parameter. The effects of temperature, kinematic viscosity, mass ratio of negative to positive ions and quantum parameter on diverging and converging shocks are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Applications of Sea-Ice Motion and Deformation Derived from Satellite Data

QuikSCAT backscatter, AMSR and DMSP SSM/I radiance data have been used to derive sea ice motion for both the Arctic and Antarctic region using the wavelet analysis tracking method. All results from QuikSCAT, AMSR and SSM5 are compatible with buoys and can then be merged by data fusion method to generate composite sea ice motion maps for more complete coverage. Furthermore, based on this merged data set daily sea-ice deformation (shear, and convergence) maps have been produced and show consistent spatial and temporal patterns. Temporal correlation maps between ice deformation and SSM/II ice concentration show interesting results in the Arctic, especially the coast area. In this study, principal component analysis for both the merged ice tracking result from satellite data and pressure field from buoy have also been examined for the relationship between the principal components and eigenvectors from these two data sets. While the result shows that principal components of modes 1 and 2 from two data sets are highly correlated which confirms that wind forcing is a major factor driving the ice drift, it also reveals that other high energy modes are not highly correlated which maybe caused by coastal effects. Principal component analysis of Arctic sea-ice motion during fall/winter period in different years shows the reverse of dominant modes or patterns is related to the Arctic Oscillation.

Liu, Antony K.↗

OPE convergence in non-relativistic conformal field theories

Motivated by applications to the study of ultracold atomic gases near the unitarity limit, we investigate the structure of the operator product expansion (OPE) in non-relativistic conformal field theories (NRCFTs). The main tool used in our analysis is the representation theory of charged (i.e. non-zero particle number) operators in the NR-CFT, in particular the mapping between operators and states in a non-relativistic “radial quantization” Hilbert space. Our results include: a determination of the OPE coefficients of descendant operators in terms of those of the underlying primary state, a demonstration of convergence of the (imaginary time) OPE in certain kinematic limits, and an estimate of the decay rate of the OPE tail inside matrix elements which, as in relativistic CFTs, depends exponentially on operator dimensions. To illustrate our results we consider several examples, including a strongly interacting field theory of bosons tuned to the unitarity limit, as well as a class of holographic models. Given the similarity with known statements about the OPE in SO(2, d) invariant field theories, our results suggest the existence of a bootstrap approach to constraining NRCFTs, with applications to bound state spectra and interactions. We briefly comment on a possible implementation of this non-relativistic conformal bootstrap program.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Differentiable but exact formulation of density-functional theory

The universal density functional F of density-functional theory is a complicated and ill-behaved function of the density—in particular, F is not differentiable, making many formal manipulations more complicated. While F has been well characterized in terms of convex analysis as forming a conjugate pair (E, F) with the ground-state energy E via the Hohenberg–Kohn and Lieb variation principles, F is nondifferentiable and subdifferentiable only on a small (but dense) subset of its domain. In this article, we apply a tool from convex analysis, Moreau–Yosida regularization, to construct, for any ε > 0, pairs of conjugate functionals ({sup ε}E, {sup ε}F) that converge to (E, F) pointwise everywhere as ε → 0{sup +}, and such that {sup ε}F is (Fréchet) differentiable. For technical reasons, we limit our attention to molecular electronic systems in a finite but large box. It is noteworthy that no information is lost in the Moreau–Yosida regularization: the physical ground-state energy E(v) is exactly recoverable from the regularized ground-state energy {sup ε}E(v) in a simple way. All concepts and results pertaining to the original (E, F) pair have direct counterparts in results for ({sup ε}E, {sup ε}F). The Moreau–Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of {sup ε}F, a rigorous formulation of Kohn–Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn–Sham theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparing temperature convergence of shocked thin films of tin and iron to a bulk temperature source

An outstanding challenge in developing a complete equation of state for materials at elevated pressure and temperature is a robust method of determining the bulk temperature state under dynamic conditions. In metals, the determination of bulk temperature states by optical pyrometry is complicated by the small optical depth and thermal conduction effects. These effects lead to observed temperatures differing by 20% or more from the bulk temperature state. In this work, we show the presence of thermal conduction effects in temperature measurements of tin and iron coatings during dynamic compression experiments. We demonstrate that tin, in contrast to iron, coatings can fail to converge to a bulk temperature source over the time scale of the experiment, requiring the experimenter to modify assumptions, design, or analysis. This work bounds thermal transport at shocked conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Average-atom treatment of relaxation time in x-ray Thomson scattering from warm dense matter

Here, the influence of finite relaxation times on Thomson scattering from warm dense plasmas is examined within the framework of the average-atom approximation. Presently most calculations use the collision-free Lindhard dielectric function to evaluate the free-electron contribution to the Thomson cross section. In this work, we use the Mermin dielectric function, which includes relaxation time explicitly. The relaxation time is evaluated by treating the average atom as an impurity in a uniform electron gas and depends critically on the transport cross section. The calculated relaxation rates agree well with values inferred from the Ziman formula for the static conductivity and also with rates inferred from a fit to the frequency-dependent conductivity. Transport cross sections determined by the phase-shift analysis in the average-atom potential are compared with those evaluated in the commonly used Born approximation. The Born approximation converges to the exact cross sections at high energies; however, differences that occur at low energies lead to corresponding differences in relaxation rates. The relative importance of including relaxation time when modeling x-ray Thomson scattering spectra is examined by comparing calculations of the free-electron dynamic structure function for Thomson scattering using Lindhard and Mermin dielectric functions. Applications are given to warm dense Be plasmas, with temperatures ranging from 2 to 32 eV and densities ranging from 2 to 64 g/cc.

74 ATOMIC AND MOLECULAR PHYSICS↗