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.

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Control of the errors of discretization and idealization in finite element analysis

Understanding of the basic principles which control errors of discretization in finite element analysis has increased very substantially since 1980. The main milestones were: (1) development of the theoretical basis of p-extensions (1981); (2) understanding of the proper interplay between mesh design and assignment of polynomial degree to elements. Practical realization of exponential convergence rates, independently of the smoothness of the exact solution (1984); and (3) industrial experience with the new finite element technology known as the p- or hp-version of the finite element method: General Dynamics reported thirty- to forty-fold savings in terms of human time and large savings in computer time (1986). Lockheed reported favorably on their evaluation of error estimation and quality control capabilities of the p-version in industrial settings (1987). The gains in our understanding of how to control the errors of discretization represent only half of the control necessary to ensure that a numerical model is in fact an accurate representation of the corresponding physical system. Control of the errors of idealization is equally important. A brief overview of the main ideas of how to ensure the quality and reliability of mathematical models of structural systems is presented.

Szabo, Barna A.↗

NASA NTRS1990

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems, task 1: Ducted propfan analysis

The time-dependent three-dimensional Euler equations of gas dynamics were solved numerically to study the steady compressible transonic flow about ducted propfan propulsion systems. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. An implicit residual smoothing operator was used to aid convergence. Two calculation grids were employed in this study. The first grid utilized an H-type mesh network with a branch cut opening to represent the axisymmetric cowl. The second grid utilized a multiple-block mesh system with a C-type grid about the cowl. The individual blocks were numerically coupled in the Euler solver. Grid systems were generated by a combined algebraic/elliptic algortihm developed specifically for ducted propfans. Numerical calculations were initially performed for unducted propfans to verify the accuracy of the three-dimensional Euler formulation. The Euler analyses were then applied for the calculation of ducted propfan flows, and predicted results were compared with experimental data for two cases. The three-dimensional Euler analyses displayed exceptional accuracy, although certain parameters were observed to be very sensitive to geometric deflections. Both solution schemes were found to be very robust and demonstrated nearly equal efficiency and accuracy, although it was observed that the multi-block C-grid formulation provided somewhat better resolution of the cowl leading edge region.

Hall, Edward J.↗

NASA NTRS2000

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗

NASA NTRS1996

Overset Grid Methods for Multidisciplinary Applications in Rotorcraft Problems

A methodology for the coupling of an advanced computational fluid dynamics method based on an overset grid flow-solver and an advanced computational structural dynamics method based on a finite element analysis is presented. Various procedures for the fluid-structure interactions modeling along with their limitations are also discussed. The flight test data for the four-bladed UH-60A Blackhawk helicopter rotor is chosen for the validation of the results. Convergence and accuracy are tested by numerical experiments with a single-bladed rotor. A comparison of airload predictions with flight test data as well as with a rigid blade case is presented. Grid and interpolation related issues for this aeroelastic application are described.

Ahmad, J. U.↗

DOE OSTI2017

Evaluation of Rapid for a UNF cask benchmark problem

To ensure criticality safety and maintain a high level of security of Used Nuclear Fuel (UNF), it is necessary to perform detailed Particle Transport simulations. The commonly used technique is the 3-D, continuous energy Monte Carlo method that is time-consuming and, for eigenvalue calculations, may have difficulties with the source convergence. Wenner and Haghighat and Dufek demonstrated that in criticality problems with high dominance ratio and/or under-sampling possibilities, the fission matrix (FM) methodology should be used in order to be able to obtain unbiased solution. The RAPID (Real-time Analysis for Particle transport and In-situ Detection) code system is based on the Multi-stage Response-function Transport (MRT) methodology, and uses the FM formulation to determine eigenvalue, subcritical multiplication and fission density distribution of nuclear systems. This paper uses a benchmark problem similar to the GBC-32 Cask to evaluate the accuracy and performance of RAPID against the results of a reference MCNP solution. We performed a detailed sensitivity analysis on a single-assembly model in order to determine the MCNP eigenvalue parameters necessary to obtain a reliable reference solution. These results provided a guide for the selection of an appropriate set of parameters for the eigenvalue calculation for the whole cask. The paper is organized as follows: initially, the RAPID code system is briefly discussed; then, the benchmark problem is introduced. The methodology used to analyze results is proposed and some of the results of different cases are analyzed to identify the best set of eigenvalue parameters for a reference MCNP calculation. The selected case is then compared with the RAPID calculation for both a single assembly and a full cask model. In the end some conclusions and future work perspectives are drawn.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2015

Spectral analysis of shielded gamma ray sources using precalculated library data

In this paper, an approach has been developed for determining the intensity of a shielded source by first determining the thicknesses of three different shielding materials from a passively collected gamma-ray spectrum by making comparisons with predetermined shielded spectra. These evaluations are dependent on the accuracy and validity of the predetermined library spectra which were created by changing the thicknesses of the three chosen materials lead, aluminum and wood that are used to simulate any actual shielding. Each of the spectra produced was generated using MCNP5 with a sufficiently large number of histories to ensure a low relative error at each channel. The materials were held in the same respective order from source to detector, where each material consisted of three individual thicknesses and a null condition. This then produced two separate data sets of 27 total shielding material situations and subsequent predetermined libraries that were created for each radionuclide source used. The technique used to calculate the thicknesses of the materials implements a Levenberg–Marquardt nonlinear search that employs a tri-linear interpolation with the respective predetermined libraries within each channel for the supplied input unknown spectrum. Given that the nonlinear parameters require an initial guess for the calculations, the approach demonstrates first that when the correct values are input, the correct thicknesses are found. It then demonstrates that when multiple trials of random values are input for each of the nonlinear parameters, the average of the calculated solutions that successfully converges also produced the correct thicknesses. Under situations with sufficient information known about the detection situation at hand, the method was shown to behave in a manner that produces reasonable results and can serve as a good preliminary solution. Finally, this technique has the capability to be used in a variety of full spectrum inverse analysis problems including homeland security issues.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1995

A hybrid neurocomputing/numerical strategy for nonlinear structural analysis

A hybrid neurocomputing/numerical strategy is presented for geometrically nonlinear analysis of structures. The strategy combines model-free data processing capabilities of computational neural networks with a Pade approximants-based perturbation technique to predict partial information about the nonlinear response of structures. In the hybrid strategy, multilayer feedforward neural networks are used to extend the validity of solutions by using training samples produced by Pade approximations to the Taylor series expansion of the response function. The range of validity of the training samples is taken to be the radius of convergence of Pade approximants and is estimated by setting a tolerance on the diverging approximants. The norm of residual vector of unbalanced forces in a given element is used as a measure to assess the quality of network predictions. To further increase the accuracy and the range of network predictions, additional training data are generated by either applying linear regression to weight matrices or expanding the training data by using predicted coefficients in a Taylor series. The effectiveness of the hybrid strategy is assessed by performing large-deflection analysis of a doubly-curved composite panel with a circular cutout, and postbuckling analyses of stiffened composite panels subjected to an in-plane edge shear load. In all the problems considered, the hybrid strategy is used to predict selective information about the structural response, namely the total strain energy and the maximum displacement components only.

Szewczyk, Z. Peter↗

NASA NTRS1994

Program Aids Analysis And Optimization Of Design

NETS/ PROSSS (NETS Coupled With Programming System for Structural Synthesis) computer program developed to provide system for combining NETS (MSC-21588), neural-network application program and CONMIN (Constrained Function Minimization, ARC-10836), optimization program. Enables user to reach nearly optimal design. Design then used as starting point in normal optimization process, possibly enabling user to converge to optimal solution in significantly fewer iterations. NEWT/PROSSS written in C language and FORTRAN 77.

Rogers, James L., Jr.↗

NASA NTRS2014

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

pulse detonation engines↗

NASA NTRS2014

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

propulsion↗

DOE OSTI2016

Preparatory Verification of a New Monte Carlo Code cosRMC based on VERA Core Physics Benchmark

To promote nuclear power plant system design and simulation software self-reliance development, DOE of China approved a R and D key project to develop an integrated engine for pressurized water reactor core and system design and analysis, referred as COSINE hereafter. The scope of this code package includes reactor physics, thermal hydraulics, safety analysis, severe accidents, probability safety assessment and etc. Many universities, research institutes and nuclear related corporations in China are involved in the code development, among which the State Nuclear Power Software Development Center (SNPSDC) plays the leading role. With the requirement of accurate three-dimensional modeling of the advanced, special or new complex core in reactor neutronics and shielding analysis, Monte Carlo method is becoming a more powerful tool for reactor analysis. Thus, a new Monte Carlo transport code named cosRMC (abbreviation for COSINE Reactor Monte Carlo) has been developed jointly by SNPSDC and Department of Engineering Physics, Tsinghua University, the main developer of RMC code. The cosRMC now has many functions as criticality calculation, fixed-source calculation, burnup calculation, and etc. Some techniques for geometry treatment, new burnup algorithm, source convergence acceleration, massive tally and parallel calculation have been implemented in cosRMC to improve the efficiency. In order to ensure its reliability, several types of V and V work have been and are being carried out for cosRMC. In this paper, some preparatory verification results based on VERA Core Physics Benchmark are presented. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗

DOE OSTI2016

Released air during vapor and air cavitation

Cavitation today is a very important problem that is solved by means of experimental and mathematical methods. The article deals with the generation of cavitation in convergent divergent nozzle of rectangular cross section. Measurement of pressure, flow rate, temperature, amount of dissolved air in the liquid and visualization of cavitation area using high-speed camera was performed for different flow rates. The measurement results were generalized by dimensionless analysis, which allows easy detection of cavitation in the nozzle. For numerical simulation the multiphase mathematical model of cavitation consisting of water and vapor was created. During verification the disagreement with the measurements for higher flow rates was proved, therefore the model was extended to multiphase mathematical model (water, vapor and air), due to release of dissolved air. For the mathematical modeling the multiphase turbulence RNG k-ε model for low Reynolds number flow with vapor and air cavitation was used. Subsequently the sizes of the cavitation area were verified. In article the inlet pressure and loss coefficient depending on the amount of air added to the mathematical model are evaluated. On the basis of the approach it may be create a methodology to estimate the amount of released air added at the inlet to the modeled area.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Flow-based likelihoods for non-Gaussian inference

We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We call these flow-based likelihoods (FBL). We analyze the accuracy and precision of the reconstructed likelihoods on mock Gaussian data, and show that simply gauging the quality of samples drawn from the trained model is not a sufficient indicator that the true likelihood has been learned. We nevertheless demonstrate that the likelihood can be reconstructed to a precision equal to that of sampling error due to a finite sample size. We then apply FBLs to mock weak lensing convergence power spectra, a cosmological observable that is significantly non-Gaussian (NG). We find that the FBL captures the NG signatures in the data extremely well, while other commonly used data-driven likelihoods, such as Gaussian mixture models and independent component analysis, fail to do so. This suggests that works that have found small posterior shifts in NG data with data-driven likelihoods such as these could be underestimating the impact of non-Gaussianity in parameter constraints. By introducing a suite of tests that can capture different levels of NG in the data, we show that the success or failure of traditional data-driven likelihoods can be tied back to the structure of the NG in the data. Here, unlike other methods, the flexibility of the FBL makes it successful at tackling different types of NG simultaneously. Because of this, and consequently their likely applicability across datasets and domains, we encourage their use for inference when sufficient mock data are available for training.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Crystallographic data processing for free-electron laser sources

A processing pipeline for diffraction data acquired using the ‘serial crystallography’ methodology with a free-electron laser source is described with reference to the crystallographic analysis suite CrystFEL and the pre-processing program Cheetah. A processing pipeline for diffraction data acquired using the ‘serial crystallography’ methodology with a free-electron laser source is described with reference to the crystallographic analysis suite CrystFEL and the pre-processing program Cheetah. A detailed analysis of the nature and impact of indexing ambiguities is presented. Simulations of the Monte Carlo integration scheme, which accounts for the partially recorded nature of the diffraction intensities, are presented and show that the integration of partial reflections could be made to converge more quickly if the bandwidth of the X-rays were to be increased by a small amount or if a slight convergence angle were introduced into the incident beam.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1984

Performance of the h, p and h-p versions

There are three basic versions of the finite element method, called the h, p and h-p versions. They are characterized by the way in which the finite element meshes and polynomial degree of elements are chosen. They differ in computer implementation (program architecture) and mathematical analysis. The manner in which the meshes and polynomial degree of the elements affect the accuracy of finite element solutions is examined. The approach is to fix certain parameters and increase the number of degrees of freedom so that the finie element solutions converge to the exact solution. Such a systematic increase of the number of degrees of freedom is called extension because it can be interpreted as a systematic extension of finite element spaces.

Babuska, I.↗

NASA NTRS1986

Non-reflecting boundary conditions for the compressible Navier-Stokes equations

A small perturbation analysis, in the long wavelength regime, is used to obtain the downstream boundary condition for the pressure for the flow over a flat plate. The methodology is extendable to other geometries. Numerical results for high Reynolds number laminar flows show great improvement in convergence rate to steady state as well as the quality of the results.

Abarbanel, S.↗

NASA NTRS1991

An introduction to the global positioning system and some geological applications

The fundamental principles of the global positioning system (GPS) are reviewed, with consideration given to geological and geophysical applications and related accuracy requirements. Recent improvements are emphasized which relate to areas such as equipment cost, limitations in the GPS satellite constellation, data analysis, uncertainties in satellite orbits and propagation delays, and problems in resolving carrier phase cycle ambiguities. Earthquake processes and near-fault crustal deformation monitoring have been facilitated by advances in GPS data acquisition and analysis. Horizontal positioning capability has been improved by new satellite constellation, better models, and global tracking networks. New classes of tectonic problems may now be studied through GPS, such as kinematic descriptions of crustal deformation and the measurement of relative plate motion at convergent boundaries. Continued improvements in the GPS are foreseen.

Dixon, T. H.↗