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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Cluster decomposition of full configuration interaction wave functions: A tool for chemical interpretation of systems with strong correlation

Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size, thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e., those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI algorithm. Furthermore, the cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the π space of polyacenes, and (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with an increasing rank for the first three systems, even connected octuple excitations still appear important in Cr 2 , suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry.

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

Dose Rate Evaluation for the ES-3100 Package with HEU Content Using MCNP, ADVANTG, Monaco, and MAVRIC

Here, this paper presents a comparative study of dose rate calculations for the ES-3100 package with highly enriched uranium (HEU) content for different source configurations using computer codes: (1) MCNP; (2) Automated Variance Reduction Generator (ADVANTG)/MCNP; (3) Monaco; and (4) Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC). The Model ES-3100 package was developed at the Y-12 National Security Complex for domestic and international transportation of Type B fissile radioactive material. In this study, six different source configurations (i.e., solid cylinder, cylindrical hemi-shell, cylindrical shell, rectangular plate, cylindrical rod, or cylindrical segment form) having 36 kg of HEU metal inside the package containment vessel (based on configurations in the ES-3100/HEU safety analysis report for packaging) are evaluated. Dose rates at the package surfaces and 1 m from the package surfaces are calculated for these different source configurations. MCNP and Monaco cases are run without any biasing options to accelerate the convergence. Consistent Adjoint Driven Importance Sampling (CADIS) and Forward-Weighted CADIS (FW-CADIS) methods developed at the Oak Ridge National Laboratory are implemented in ADVANTG/MCNP and MAVRIC codes to accelerate the convergence. ADVANTG generates variance reduction parameters using Denovo code, and MCNP is used with the variance reduction parameters to accelerate the convergence. MAVRIC uses Denovo code to construct an importance map and a biased source distribution that are supplied to Monaco to accelerate the Monte Carlo simulation. The FW-CADIS option in ADVANTG and MAVRIC is used to accelerate the convergence in this study. The accelerated convergence cases (ADVANTG/MCNP and MAVRIC) are about 100 times faster with 100 times less particle simulation than those cases run without biasing options (analog MCNP and analog Monaco). Lastly, the MCNP, ADVANTG/MCNP, Monaco, and MAVRIC calculated dose rates at the package surfaces and at 1 m from the package surfaces for the different source configurations are compared and are found to be in general agreement.

42 ENGINEERING↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2012

Computational Challenges for Power System Operation

As the power grid technology evolution and information technology revolution converge, power grids are witnessing a revolutionary transition, represented by emerging grid technologies and large scale deployment of new sensors and meters in networks. This transition brings opportunities, as well as computational challenges in the field of power grid analysis and operation. This paper presents some research outcomes in the areas of parallel state estimation using the preconditioned conjugated gradient method, parallel contingency analysis with a dynamic load balancing scheme and distributed system architecture. Based on this research, three types of computational challenges are identified: highly coupled applications, loosely coupled applications, and centralized and distributed applications. Recommendations for future work for power grid applications are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Application of Detectability Analysis for Power System Dynamic State Estimation

To guide the measurement placement and selection for dynamic state estimation in power systems, it is important to determine the existence of an observer that can estimate the dynamic states. Past studies have been focused on observability analysis, which determines whether the initial states can be uniquely determined. In this paper, detectability analysis is proposed to determine whether the current states can be estimated as time involving. Through the proposed detectability analysis, the states of an unobservable system can be categorized into the detectable states and undetectable states. It is shown that if the eigenvalues corresponding to the unobservable states are stable, the states are detectable and their estimates will converge to the true ones asymptotically. The study can be used to guide measurement placement and selection.

Detectability, dynamic state estimation↗
DOE OSTI2014

TH-E-BRE-03: A Novel Method to Account for Ion Chamber Volume Averaging Effect in a Commercial Treatment Planning System Through Convolution

Purpose: Fourier-based deconvolution approaches used to eliminate ion chamber volume averaging effect (VAE) suffer from measurement noise. This work aims to investigate a novel method to account for ion chamber VAE through convolution in a commercial treatment planning system (TPS). Methods: Beam profiles of various field sizes and depths of an Elekta Synergy were collected with a finite size ion chamber (CC13) to derive a clinically acceptable beam model for a commercial TPS (Pinnacle{sup 3}), following the vendor-recommended modeling process. The TPS-calculated profiles were then externally convolved with a Gaussian function representing the chamber (σ = chamber radius). The agreement between the convolved profiles and measured profiles was evaluated with a one dimensional Gamma analysis (1%/1mm) as an objective function for optimization. TPS beam model parameters for focal and extra-focal sources were optimized and loaded back into the TPS for new calculation. This process was repeated until the objective function converged using a Simplex optimization method. Planar dose of 30 IMRT beams were calculated with both the clinical and the re-optimized beam models and compared with MapCHEC™ measurements to evaluate the new beam model. Results: After re-optimization, the two orthogonal source sizes for the focal source reduced from 0.20/0.16 cm to 0.01/0.01 cm, which were the minimal allowed values in Pinnacle. No significant change in the parameters for the extra-focal source was observed. With the re-optimized beam model, average Gamma passing rate for the 30 IMRT beams increased from 92.1% to 99.5% with a 3%/3mm criterion and from 82.6% to 97.2% with a 2%/2mm criterion. Conclusion: We proposed a novel method to account for ion chamber VAE in a commercial TPS through convolution. The reoptimized beam model, with VAE accounted for through a reliable and easy-to-implement convolution and optimization approach, outperforms the original beam model in standard IMRT QA verification.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

Dose Rate Evaluation for the ES-3100 Package with HEU Content Using MCNP, ADVANTG, Monaco, and MAVRIC

This paper presents a comparative study of dose rate calculations for the ES-3100 package with highly enriched uranium (HEU) content for different source configurations using computer codes: (1) MCNP; (2) Automated Variance Reduction Generator (ADVANTG)/MCNP; (3) Monaco; and (4) Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC). The Model ES-3100 package was developed at the Y-12 National Security Complex for domestic and international transportation of Type B fissile radioactive material. In this study, six different source configurations (i.e., solid cylinder, cylindrical hemi-shell, cylindrical shell, rectangular plate, cylindrical rod, or cylindrical segment form) having 36 kg of HEU metal inside the package containment vessel (based on configurations in the ES-3100/HEU safety analysis report for packaging) are evaluated. Dose rates at the package surfaces and 1 m from the package surfaces are calculated for these different source configurations. MCNP and Monaco cases are run without any biasing options to accelerate the convergence. Consistent Adjoint Driven Importance Sampling (CADIS) and Forward-Weighted CADIS (FW-CADIS) methods developed at the Oak Ridge National Laboratory are implemented in ADVANTG/MCNP and MAVRIC codes to accelerate the convergence. ADVANTG generates variance reduction parameters using Denovo code, and MCNP is used with the variance reduction parameters to accelerate the convergence. MAVRIC uses Denovo code to construct an importance map and a biased source distribution that are supplied to Monaco to accelerate the Monte Carlo simulation. The FW-CADIS option in ADVANTG and MAVRIC is used to accelerate the convergence in this study. The accelerated convergence cases (ADVANTG/MCNP and MAVRIC) are about 100 times faster with 100 times less particle simulation than those cases run without biasing options (analog MCNP and analog Monaco). The MCNP, ADVANTG/MCNP, Monaco, and MAVRIC calculated dose rates at the package surfaces and at 1 m from the package surfaces for the different source configurations are compared and are found to be in general agreement.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Power density calculation of spontaneous radiation emission from a finite emittance electron beam in planar undulators

We extend the angular power density formulae for spontaneous radiation emission from planar undulators to include finite emittance and angular misalignment of the electron beam. Then we calculate and compare power densities estimated from integral approach with Gaussian beam distribution and summation approach with macro-particle allocation covering particle beam phase space. After showing that both approaches converge to each other, we apply the macro-particle approach to study power absorbed and transmitted by various apertures at front end of 9-ID beamline at National Synchrotron Light Source-II for misaligned electron beam. Our analysis indicate that the electron beam misalignment could lead to unwarranted power deposition at tighter apertures that would have been otherwise difficult to account for during the aperture design/choice process.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Optimal Power Dispatch of DGs in Radial and Mesh AC Grids: A Hybrid Solution Methodology between the Salps Swarm Algorithm and Successive Approximation Power Flow Method

In this paper, we address the problem of the optimal power dispatch of Distributed Generators (DGs) in Alternating Current (AC) networks, better known as the Optimal Power Flow (OPF) problem. We used, as the objective function, the minimization of power losses (P loss ) associated with energy transport, which are subject to the set of constraints that compose AC networks in an environment of distributed generation. To validate the effectiveness of the proposed methodology in solving the OPF problem in any network topology, we employed one 10-node mesh test system and three radial text systems: 10, 33, and 69 nodes. In each test system, DGs were allowed to inject 20%, 40%, and 60% of the power supplied by the slack generator in the base case. To solve the OPF problem, we used a master–slave methodology that integrates the optimization method Salps Swarm Algorithm (SSA) and the load flow technique based on the Successive Approximation (SA) method. Moreover, for comparison purposes, we employed some of the algorithms reported in the specialized literature to solve the OPF problem (the continuous genetic algorithm, the particle swarm optimization algorithm, the black hole algorithm, the antlion optimization algorithm, and the Multi-Verse Optimizer algorithm), which were selected because of their excellent results in solving such problems. The results obtained by the proposed solution methodology demonstrate its superiority and convergence capacity in terms of minimization of P loss in both radial and mesh systems. It provided the best reduction in minimum P loss in short processing times and showed excellent repeatability in each test system and scenario under analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1989

A three-node C deg element for analysis of laminated composite sandwich shells

A three-node flat shell element with C deg rotation fields has been developed for analysis of arbitrary composite shells. The element may consist of any number of orthotropic layers, each layer having different material properties and angular orientation. The formulation includes coupling between bending and extension, which is essential for analysis of unsymmetric laminates. Shearing deflections are included, since laminated and sandwich construction frequently results in shear stiffness much smaller than bending stiffness. Formulation of the element is straightforward, and calculation of its stiffness matrix is simple and fast. Convergence of solutions with mesh refinement is uniform for both thin and thick shells and is insensitive to element shape, although not as rapid as some other elements that lack one or more capabilities of the newly developed element. An experimental verification of the shall element is reported in the appendix.

Martin, C. Wayne↗
NASA NTRS1994

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗
NASA NTRS2004

On the Hilbert-Huang Transform Theoretical Foundation

The Hilbert-Huang Transform [HHT] is a novel empirical method for spectrum analysis of non-linear and non-stationary signals. The HHT is a recent development and much remains to be done to establish the theoretical foundation of the HHT algorithms. This paper develops the theoretical foundation for the convergence of the HHT sifting algorithm and it proves that the finest spectrum scale will always be the first generated by the HHT Empirical Mode Decomposition (EMD) algorithm. The theoretical foundation for cutting an extrema data points set into two parts is also developed. This then allows parallel signal processing for the HHT computationally complex sifting algorithm and its optimization in hardware.

Kizhner, Semion↗
NASA NTRS2004

A Subsonic Aircraft Design Optimization With Neural Network and Regression Approximators

The Flight-Optimization-System (FLOPS) code encountered difficulty in analyzing a subsonic aircraft. The limitation made the design optimization problematic. The deficiencies have been alleviated through use of neural network and regression approximations. The insight gained from using the approximators is discussed in this paper. The FLOPS code is reviewed. Analysis models are developed and validated for each approximator. The regression method appears to hug the data points, while the neural network approximation follows a mean path. For an analysis cycle, the approximate model required milliseconds of central processing unit (CPU) time versus seconds by the FLOPS code. Performance of the approximators was satisfactory for aircraft analysis. A design optimization capability has been created by coupling the derived analyzers to the optimization test bed CometBoards. The approximators were efficient reanalysis tools in the aircraft design optimization. Instability encountered in the FLOPS analyzer was eliminated. The convergence characteristics were improved for the design optimization. The CPU time required to calculate the optimum solution, measured in hours with the FLOPS code was reduced to minutes with the neural network approximation and to seconds with the regression method. Generation of the approximators required the manipulation of a very large quantity of data. Design sensitivity with respect to the bounds of aircraft constraints is easily generated.

Patnaik, Surya N.↗
DOE OSTI2012

Equilibrium cycle pin by pin transport depletion calculations with DeCART

As the Advanced Fuel Cycle Initiative (AFCI) program has matured it has become more important to utilize more advanced simulation methods. The work reported here was performed as part of the AFCI fellowship program to develop and demonstrate the capability of performing high fidelity equilibrium cycle calculations. As part of the work here, a new multi-cycle analysis capability was implemented in the DeCART code which included modifying the depletion modules to perform nuclide decay calculations, implementing an assembly shuffling pattern description, and modifying iteration schemes. During the work, stability issues were uncovered with respect to converging simultaneously the neutron flux, isotopics, and fluid density and temperature distributions in 3-D. Relaxation factors were implemented which considerably improved the stability of the convergence. To demonstrate the capability two core designs were utilized, a reference UOX core and a CORAIL core. Full core equilibrium cycle calculations were performed on both cores and the discharge isotopics were compared. From this comparison it was noted that the improved modeling capability was not drastically different in its prediction of the discharge isotopics when compared to 2-D single assembly or 2-D core models. For fissile isotopes such as U-235, Pu-239, and Pu-241 the relative differences were 1.91%, 1.88%, and 0.59%), respectively. While this difference may not seem large it translates to mass differences on the order of tens of grams per assembly, which may be significant for the purposes of accounting of special nuclear material. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Numerically approximating centrality for graph ranking guarantees

Many real-world datasets can be represented as graphs. Using iterative solvers to approximate graph centrality measures allows us to obtain a ranking vector on the nodes of the graph, consisting of a number for each vertex in the graph identifying its relative importance. In this study the centrality measures we use are Katz Centrality and PageRank. Given an approximate solution, we use the residual to accurately estimate how much of the ranking matches the ranking given by the exact solution. Using probabilistic matrix norms, we obtain bounds on the accuracy of the approximation compared to the exact solution with respect to the highly ranked nodes and apply numerical analysis to the computation of centrality with iterative methods. This relates the numerical accuracy of the linear solver to the data analysis accuracy of finding the correct ranking. In particular, we answer the question of which pairwise rankings are reliable given an approximate solution to the linear system. Experiments on many real-world undirected and directed networks up to several million vertices and several hundred million edges validate our theory and show that we are able to accurately estimate large portions of the approximation. We also analyze the difference between global centrality scores and personalized scores (w.r.t. specific seed vertices). By analyzing convergence error, we develop confidence in the ranking schemes of data mining. We show we are able to accurately guarantee ranking of vertices with an approximation to centrality metrics faster than current methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2007

Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.

Swanson, Roy C., Jr.↗
DOE OSTI2026

Asymptotic inconsistency of the cumulative algorithm for laser-induced damage probability analysis

The “cumulative algorithm” is a data analysis method that has been proposed to provide an objective, nonparametric determination of laser-induced damage probability as a function of fluence from experimental data that contain both damaged sites and undamaged sites (i.e., 1-on-1 or S-on-1 testing protocols). In this work, the limitations of this approach are explored by considering the asymptotic limit of a large number of test sites. It is shown that the cumulative algorithm does not converge to the true probability distribution and significantly underestimates the damage probability near the damage onset. Here, based on the results of this work, the cumulative algorithm is not recommended for accurate estimation of damage probability.

Computational methods↗
DOE OSTI2020

The Role of Multiscale Interaction in Tropical Cyclogenesis and Its Predictability in Near-Global Aquaplanet Cloud-Resolving Simulations

Tropical cyclogenesis (TCG) is a multiscale process that involves interactions between large-scale circulation and small-scale convection. A near-global aquaplanet cloud-resolving model (NGAqua) with 4-km horizontal grid spacing that produces tropical cyclones (TCs) is used to investigate TCG and its predictability. This study analyzes an ensemble of three 20-day NGAqua simulations, with initial white-noise perturbations of low-level humidity. TCs develop spontaneously from the northern edge of the intertropical convergence zone (ITCZ), where large-scale flows and tropical convection provide necessary conditions for barotropic instability. Zonal bands of positive low-level absolute vorticity organize into cyclonic vortices, some of which develop into TCs. A new algorithm is developed to track the cyclonic vortices. A vortex-following framework analysis of the low-level vorticity budget shows that vertical stretching of absolute vorticity due to convective heating contributes positively to the vorticity spinup of the TCs. A case study and composite analyses suggest that sufficient humidity is key for convective development. TCG in these three NGAqua simulations undergoes the same series of interactions. The locations of cyclonic vortices are broadly predetermined by planetary-scale circulation and humidity patterns associated with ITCZ breakdown, which are predictable up to 10 days. Whether and when the cyclonic vortices become TCs depend on the somewhat more random feedback between convection and vorticity.

54 ENVIRONMENTAL SCIENCES↗