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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Reanalysis Activities at the NASA Global Modeling and Assimilation Office

This talk presents an overview of recent reanalysis activities at the NASA Global Modeling and Assimilation Office (GMAO) as part of a multi-faceted strategy towards an Integrated Earth System retrospective analysis, coupling components of the atmosphere, ocean, chemistry, land, and ice. While elements of the atmosphere-ocean coupled Goddard Earth Observing System (GEOS) model and data assimilation are being actively developed, a suite of reanalysis products is designed to provide further understanding of key aspects of Earth system coupling in a reanalysis context: The baseline atmospheric reanalysis, the GEOS Retrospective analysis for the early 21st Century (GEOS-R21C), features recent advances in the GEOS model and data assimilation, and targets the NASA’s Earth Observing System EOS and post-EOS satellite observations; GEOS-IT, a user-tailored low-resolution atmospheric reanalysis, serves as a second baseline to the NASA Instrument Teams for validation and calibration and drives a one-way coupled ocean reanalysis, GEOSIT-Ocean; PolarMERRA, a high-resolution downscaled product for the polar regions, focuses on improving the representation of polar atmospheric processes with an assessment of current cryospheric biases, and targeted improvements to surface sea ice and glacier conditions; Finally, R21C-Chem, an off-line atmospheric chemistry and composition reanalysis, includes both tropospheric and stratospheric trace gases. The diversity of these reanalysis activities presents unique opportunities for collaborations cross-teams/institutions, with new commercial data partners, and with end-user groups. This talk will discuss these opportunities and explore leveraging the lessons learned along the way on key drivers in Earth system interactions as we converge towards the next generation of the Modern-Era Retrospective analysis for Research and Applications (MERRA) suite.

Amal El Akkraoui↗
NASA NTRS1991

Acceleration of convergence by shifting the spectrum of implicit finite difference operators associated with the equations of gas dynamics

Eigensystem analysis techniques are applied to finite difference formulations of the Navier-Stokes equations in one dimension. Spectra of the resulting implicit difference operators are computed. The largest eigenvalues are calculated by using a combination of the Frechet derivative of the operators and Arnoldi's method. The accuracy of Arnoldi's method is tested by comparing the rate of convergence of the iterative method with the dominant eigenvalue of the original iteration matrix. On the basis of the pattern of eigenvalue distributions for various flow configurations, a shifting of the implicit operators in question is devised. This procedure has improved the rates of convergence of CFD codes by 20 - 50 percent.

Cheer, A.↗
NASA NTRS1978

Application of a trigonometric finite difference procedure to numerical analysis of compressive and shear buckling of orthotropic panels

A numerical analysis developed for the buckling of rectangular orthotropic layered panels under combined shear and compression is described. This analysis uses a central finite difference procedure based on trigonometric functions instead of using the conventional finite differences which are based on polynomial functions. Inasmuch as the buckle mode shape is usually trigonometric in nature, the analysis using trigonometric finite differences can be made to exhibit a much faster convergence rate than that using conventional differences. Also, the trigonometric finite difference procedure leads to difference equations having the same form as conventional finite differences; thereby allowing available conventional finite difference formulations to be converted readily to trigonometric form. For two-dimensional problems, the procedure introduces two numerical parameters into the analysis. Engineering approaches for the selection of these parameters are presented and the analysis procedure is demonstrated by application to several isotropic and orthotropic panel buckling problems. Among these problems is the shear buckling of stiffened isotropic and filamentary composite panels in which the stiffener is broken. Results indicate that a break may degrade the effect of the stiffener to the extent that the panel will not carry much more load than if the stiffener were absent.

Stein, M.↗
DOE OSTI2017

Numerical Uncertainties vs. Model Uncertainties in Two-phase Flow Simulations

In recent years, verification, validation and uncertainty quantification (VVUQ) have become a common practice in thermal-hydraulics analysis. In general, these activities deal with propagations of uncertainties in computer code simulations, e.g., through system analysis codes. However, most existing such activities in thermal-hydraulics analysis have been primarily focused on model uncertainties, while numerical uncertainties were largely overlooked. Numerical uncertainties (or numerical errors) can appear in many different forms, e.g., round-off error, statistical sampling error, and discretization error. As pointed out by Oberkampf and Roy, discretization error is usually the largest however the most difficult type of error to estimate reliably. In thermal-hydraulics analysis, especially two-phase flow simulations commonly encountered in reactor safety analysis, the lack of consideration on discretization error is mainly due to the difficulty in estimating them. Accurate estimations on discretization error require continuous mesh refinement and/or implementation of high-order numerical schemes in system analysis code, both of which are unfortunately difficult to achieve for existing codes. At first, it should be emphasized that mesh refinement (or mesh sensitivity study) is a well-understood concept in the thermal-hydraulics field. However, it is not commonly practiced for many reasons, e.g., lack of understanding on numerical errors by code users, consideration on computational cost, deteriorated numerical stability performance with refined meshes, and etc. Secondly, high-order (e.g., second-order) numerical schemes are simply unavailable in almost any existing system analysis code. Implementation of high-order numerical schemes in these codes is difficult, if not impossible, given the complexity of these codes. On the other hand, people tend to argue that numerical uncertainties (or discretization error as we have emphasized) are not as important as the commonly large uncertainties in two-phase flow models, e.g., closure correlations. As a result, nodalization is sometimes played unintentionally, or maybe intentionally, to match experimental data, which is referred as 'user effects' or 'compensating errors' in thermal-hydraulics community. Such a practice is highly criticized by Levy. In this work, we will at first build a computer code that incorporates both first-order and second-order numerical methods to solve the two-phase flow problems. The first-order method resembles the one used in many existing system analysis codes; and the second-order method works as the reference to estimate numerical errors. Numerical verifications on spatial discretization schemes will be presented in the form of mesh convergence study. It will also be demonstrated that, in practical scenarios, discretization errors can be as large as, or even larger than, model uncertainties via case studies.

42 ENGINEERING↗
NASA NTRS1993

Analysis of microwave-heated rocket engines for space propulsion

Microwave-heated plasmas in convergent nozzles are analyzed using a coupled Maxwell and Navier-Stokes solver to examine relevant issues associated with microwave thermal propulsion. Parametric studies are conducted to understand the effect of power, pressure, and plasma location with respect to the nozzle throat. For nozzles in the 0.5 to 3 N range with helium flow, results show that specific impulses up to 550-650 sec are possible, with further increases being limited by severe wall-heating. Coupling efficiencies of over 90 percent are consistently obtained, with overall efficiencies ranging from 40 percent to 80 percent. Size scale-up studies (done by scaling the frequency from 2.45 GHz to 0.91 GHz) indicate that plasma migration toward the walls occurs more frequently for the lower frequency. Increasing the cavity aspect ratio and detuning the cavity are found to be effective ways of keeping the plasma on-axis.

Schwer, D. A.↗
NASA NTRS2017

Benchmark Shock Tube Experiments for Radiative Heating Relevant to Earth Re-Entry

Detailed spectrally and spatially resolved radiance has been measured in the Electric Arc Shock Tube (EAST) facility for conditions relevant to high speed entry into a variety of atmospheres, including Earth, Venus, Titan, Mars and the Outer Planets. The tests that measured radiation relevant for Earth re-entry are the focus of this work and are taken from campaigns 47, 50, 52 and 57. These tests covered conditions from 8 km/s to 15.5 km/s at initial pressures ranging from 0.05 Torr to 1 Torr, of which shots at 0.1 and 0.2 Torr are analyzed in this paper. These conditions cover a range of points of interest for potential fight missions, including return from Low Earth Orbit, the Moon and Mars. The large volume of testing available from EAST is useful for statistical analysis of radiation data, but is problematic for identifying representative experiments for performing detailed analysis. Therefore, the intent of this paper is to select a subset of benchmark test data that can be considered for further detailed study. These benchmark shots are intended to provide more accessible data sets for future code validation studies and facility-to-facility comparisons. The shots that have been selected as benchmark data are the ones in closest agreement to a line of best fit through all of the EAST results, whilst also showing the best experimental characteristics, such as test time and convergence to equilibrium. The EAST data are presented in different formats for analysis. These data include the spectral radiance at equilibrium, the spatial dependence of radiance over defined wavelength ranges and the mean non-equilibrium spectral radiance (so-called 'spectral non-equilibrium metric'). All the information needed to simulate each experimental trace, including free-stream conditions, shock time of arrival (i.e. x-t) relation, and the spectral and spatial resolution functions, are provided.

Shock Tube↗
DOE OSTI2015

Dominating Controls for Wetter South Asian Summer Monsoon in the Twenty-First Century

This study analyzes a suite of global climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) archives to understand the mechanisms behind a net increase in the South Asian summer monsoon precipitation in response to enhanced radiative forcing during the twenty-first century. An increase in radiative forcing fuels an increase in the atmospheric moisture content through warmer temperatures, which overwhelms the weakening of monsoon circulation and results in an increase of moisture convergence and therefore summer monsoon precipitation over South Asia. Moisture source analysis suggests that both regional (local recycling, the Arabian Sea, the Bay of Bengal) and remote (including the south Indian Ocean) sources contribute to the moisture supply for precipitation over South Asia during the summer season that is facilitated by the monsoon dynamics. For regional moisture sources, the effect of excessive atmospheric moisture is offset by weaker monsoon circulation and uncertainty in the response of the evapotranspiration over land, so anomalies in their contribution to the total moisture supply are either mixed or muted. In contrast, weakening of the monsoon dynamics has less influence on the moisture supply from remote sources that not only is a dominant moisture contributor in the historical period but is also the net driver of the positive summer monsoon precipitation response in the twenty-first century. Finally, the results also indicate that historic measures of the monsoon dynamics may not be well suited to predict the nonstationary moisture-driven South Asian summer monsoon precipitation response in the twenty-first century.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Genesis of Pre-Hurricane Felix (2007): The Role of the Easterly Wave Critical Layer - Part 1

The formation of pre Hurricane Felix (2007) in a tropical easterly wave is examined in a two-part study using the Weather Research and Forecasting (WRF) model with a high-resolution nested grid configuration that permits the representation of cloud system processes. The simulation commences during the wave stage of the precursor African easterly-wave disturbance. Here the simulated and observed developments are compared, while in Part II of the study various large-scale analyses, physical parameterizations, and initialization times are explored to document model sensitivities. In this first part the authors focus on the wave/vortex morphology, its interaction with the adjacent intertropical convergence zone complex, and the vorticity balance in the neighborhood of the developing storm. Analysis of the model simulation points to a bottom-up development process within the wave critical layer and supports the three new hypotheses of tropical cyclone formation proposed recently by Dunkerton, Montgomery, and Wang. It is shown also that low-level convergence associated with the ITCZ helps to enhance the wave signal and extend the "wave pouch" from the jet level to the top of the atmospheric boundary layer. The region of a quasi-closed Lagrangian circulation within the wave pouch provides a focal point for diabatic merger of convective vortices and their vortical remnants. The wave pouch serves also to protect the moist air inside from dry air intrusion, providing a favorable environment for sustained deep convection. Consistent with the authors' earlier findings, the tropical storm forms near the center of the wave pouch via system-scale convergence in the lower troposphere and vorticity aggregation. Components of the vorticity balance are shown to be scale dependent, with the immediate effects of cloud processes confined more closely to the storm center than the overturning Eliassen circulation induced by diabatic heating, the influence of which extends to larger radii.

Wang, Zhuo↗
DOE OSTI2020

Constraining physical models at gigabar pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative infor mation. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. As a result, the goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Preferred latitudes of the intertropical convergence zone

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

Waliser, Duane E.↗
DOE OSTI2016

Controlling excited-state contamination in nucleon matrix elements

We present a detailed analysis of methods to reduce statistical errors and excited-state contamination in the calculation of matrix elements of quark bilinear operators in nucleon states. All the calculations were done on a 2+1-flavor ensemble with lattices of size 32 3 × 64 generated using the rational hybrid Monte Carlo algorithm at 𝑎 = 0.081 fm and with 𝑀 𝜋 = 312 MeV. The statistical precision of the data is improved using the all-mode-averaging method. We compare two methods for reducing excited-state contamination: a variational analysis and a 2-state fit to data at multiple values of the source-sink separation 𝑡 sep . We show that both methods can be tuned to significantly reduce excited-state contamination and discuss their relative advantages and cost effectiveness. Here a detailed analysis of the size of source smearing used in the calculation of quark propagators and the range of values of 𝑡 sep needed to demonstrate convergence of the isovector charges of the nucleon to the 𝑡 sep → ∞ estimates is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2000

Session on High Speed Civil Transport Design Capability Using MDO and High Performance Computing

Since the inception of CAS in 1992, NASA Langley has been conducting research into applying multidisciplinary optimization (MDO) and high performance computing toward reducing aircraft design cycle time. The focus of this research has been the development of a series of computational frameworks and associated applications that increased in capability, complexity, and performance over time. The culmination of this effort is an automated high-fidelity analysis capability for a high speed civil transport (HSCT) vehicle installed on a network of heterogeneous computers with a computational framework built using Common Object Request Broker Architecture (CORBA) and Java. The main focus of the research in the early years was the development of the Framework for Interdisciplinary Design Optimization (FIDO) and associated HSCT applications. While the FIDO effort was eventually halted, work continued on HSCT applications of ever increasing complexity. The current application, HSCT4.0, employs high fidelity CFD and FEM analysis codes. For each analysis cycle, the vehicle geometry and computational grids are updated using new values for design variables. Processes for aeroelastic trim, loads convergence, displacement transfer, stress and buckling, and performance have been developed. In all, a total of 70 processes are integrated in the analysis framework. Many of the key processes include automatic differentiation capabilities to provide sensitivity information that can be used in optimization. A software engineering process was developed to manage this large project. Defining the interactions among 70 processes turned out to be an enormous, but essential, task. A formal requirements document was prepared that defined data flow among processes and subprocesses. A design document was then developed that translated the requirements into actual software design. A validation program was defined and implemented to ensure that codes integrated into the framework produced the same results as their standalone counterparts. Finally, a Commercial Off the Shelf (COTS) configuration management system was used to organize the software development. A computational environment, CJOPT, based on the Common Object Request Broker Architecture, CORBA, and the Java programming language has been developed as a framework for multidisciplinary analysis and Optimization. The environment exploits the parallelisms inherent in the application and distributes the constituent disciplines on machines best suited to their needs. In CJOpt, a discipline code is "wrapped" as an object. An interface to the object identifies the functionality (services) provided by the discipline, defined in Interface Definition Language (IDL) and implemented using Java. The results of using the HSCT4.0 capability are described. A summary of lessons learned is also presented. The use of some of the processes, codes, and techniques by industry are highlighted. The application of the methodology developed in this research to other aircraft are described. Finally, we show how the experience gained is being applied to entirely new vehicles, such as the Reusable Space Transportation System. Additional information is contained in the original.

Rehder, Joe↗
DOE OSTI2022

An in-line burn-up algorithm: proposal and preliminary analysis

A new burn-up algorithm is proposed for coupling neutronics and isotopic depletion. The algorithm is analogous to those proposed for enforcing xenon equilibrium and converging neutronics and thermal-hydraulics. During the corrector step of standard predictor-corrector methods, the corrector nuclide density is updated as neutronics converges. This is expected to improve stability while reducing computational expense compared to other implicit algorithms. The stability of the algorithm is investigated by Fourier analysis for a simple case whose results have been successfully compared with numerical experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. The governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models are described in detail.

Towne, Charles E.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗
DOE OSTI2024

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Separation analysis, a tool for analyzing multigrid algorithms

The separation of vectors by multigrid (MG) algorithms is applied to the study of convergence and to the prediction of the performance of MG algorithms. The separation operator for a two level cycle algorithm is derived. It is used to analyze the efficiency of the cycle when mixing of eigenvectors occurs. In particular cases the separation analysis reduces to Fourier type analysis. The separation operator of a two level cycle for a Schridubger eigenvalue problem, is derived and analyzed in a Fourier basis. Separation analysis gives information on how to choose performance relaxations and inter-level transfers. Separation analysis is a tool for analyzing and designing algorithms, and for optimizing their performance.

Costiner, Sorin↗