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,063 records · Page 59

All sources NASA NTRS DOE OSTI

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.↗

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 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 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↗

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.↗

NASA NTRS1992

Unstructured grid methods for compressible flows

The implementation of the finite element method on unstructured triangular grids is described and the development of centered finite element schemes for the solution of the compressible Euler equation on general triangular and tetrahedral grids is discussed. Explicit and implicit Lax-Wendroff type methods and a method based upon the use of explicit multistep timestepping are considered. In the latter case, the convergence behavior of the method is accelerated by the incorporation of a fully unstructured multigrid procedure. The advancing front method for generating unstructured grids of triangles and tetrahedra is described and the application of adaptive mesh techniques to both steady and transient flow analysis is illustrated.

Morgan, K.↗

NASA NTRS1992

Approximation methods for stochastic petri nets

Stochastic Marked Graphs are a concurrent decision free formalism provided with a powerful synchronization mechanism generalizing conventional Fork Join Queueing Networks. In some particular cases the analysis of the throughput can be done analytically. Otherwise the analysis suffers from the classical state explosion problem. Embedded in the divide and conquer paradigm, approximation techniques are introduced for the analysis of stochastic marked graphs and Macroplace/Macrotransition-nets (MPMT-nets), a new subclass introduced herein. MPMT-nets are a subclass of Petri nets that allow limited choice, concurrency and sharing of resources. The modeling power of MPMT is much larger than that of marked graphs, e.g., MPMT-nets can model manufacturing flow lines with unreliable machines and dataflow graphs where choice and synchronization occur. The basic idea leads to the notion of a cut to split the original net system into two subnets. The cuts lead to two aggregated net systems where one of the subnets is reduced to a single transition. A further reduction leads to a basic skeleton. The generalization of the idea leads to multiple cuts, where single cuts can be applied recursively leading to a hierarchical decomposition. Based on the decomposition, a response time approximation technique for the performance analysis is introduced. Also, delay equivalence, which has previously been introduced in the context of marked graphs by Woodside et al., Marie's method and flow equivalent aggregation are applied to the aggregated net systems. The experimental results show that response time approximation converges quickly and shows reasonable accuracy in most cases. The convergence of Marie's method and flow equivalent aggregation are applied to the aggregated net systems. The experimental results show that response time approximation converges quickly and shows reasonable accuracy in most cases. The convergence of Marie's is slower, but the accuracy is generally better. Delay equivalence often fails to converge, while flow equivalent aggregation can lead to potentially bad results if a strong dependence of the mean completion time on the interarrival process exists.

Jungnitz, Hauke Joerg↗

NASA NTRS1996

Design of Robust Adaptive Unbalance Response Controllers for Rotors with Magnetic Bearings

Experimental results have recently demonstrated that an adaptive open loop control strategy can be highly effective in the suppression of unbalance induced vibration on rotors supported in active magnetic bearings. This algorithm, however, relies upon a predetermined gain matrix. Typically, this matrix is determined by an optimal control formulation resulting in the choice of the pseudo-inverse of the nominal influence coefficient matrix as the gain matrix. This solution may result in problems with stability and performance robustness since the estimated influence coefficient matrix is not equal to the actual influence coefficient matrix. Recently, analysis tools have been developed to examine the robustness of this control algorithm with respect to structured uncertainty. Herein, these tools are extended to produce a design procedure for determining the adaptive law's gain matrix. The resulting control algorithm has a guaranteed convergence rate and steady state performance in spite of the uncertainty in the rotor system. Several examples are presented which demonstrate the effectiveness of this approach and its advantages over the standard optimal control formulation.

Knospe, Carl R.↗

NASA NTRS2002

Evidence of an Intrinsic Intraseasonal Oscillation over Tropical South America During Austral Summer

The intraseasonal variation (ISV) in the 30-60 day band, also known as Madden-Julian oscillation (MJO), has been studied for decades. Madden and Julian showed that the oscillation originated from the western Indian Ocean, propagated eastward, got enhanced over the maritime continent and weakened after passing over the dateline. Composite studies showed evidences of a signal in upper and lower level zonal wind propagating around the globe during an oscillation. Theoretical studies pointed out that the interaction with the warm ocean surface and the coupling with the convective and radiative processes in the atmosphere could manifest the oscillation, which propagates eastward via mutual feedbacks between the wave motions and the cumulus heating. Over tropical South America, no independent 30-60 day oscillation has been reported so far, despite that Amazon is the most distinct tropical convection center over the western hemisphere and the fluxes from its surface of tropical rainforests are close to that from the warm tropical ocean. Liebmann et al. showed a distinct spectral peak of 40-50 day oscillation in outgoing longwave radiation (OLR) over tropical South America and considered that was manifested by the MJO propagation. Nogues-Paegle et al. (2000) focused on a dipole pattern of the OLR anomaly with centers of action over the South Atlantic Convergence Zone (SACZ) and the subtropical plain. They used the regional 10-90 day filtered data and demonstrated this pattern could be represented by the fifth mode of the rotated empirical orthogonal function. Its principal component was further analyzed using the singular spectrum analysis. Their result showed two oscillatory modes with periods of 36-40 days and 22-28 days, of which the former was related to the MJO influence and the latter linked to the remote forcing over southwest of Australia, which produced a wave train propagating southeastward, rounding the southern tip of South America and returning back toward the northeast. The 22-28 day mode has distinct impact on SACZ, responsible for the regional seesaw pattern of alternating dry and wet conditions. In this study we will focus on the 30-60-day spectral band and investigate whether the independent oscillation source over tropical South America is existed. First, we will show the seasonal dependence of the tropical South American ISV in Section 3. Then, the leading principal modes of 30-60 day bandpass filtered 850-hPa velocity potential (VP850) will be computed to distinguish the stationary ISV over tropical South America (SISA) from the propagating MJO in the austral summertime in Section 4. The importance of SISA in representing the regional ISV over South America will be discussed. In Section 5, we will demonstrate the mass oscillation regime of SISA, which is well separated from that of MJO by the Andes, and the convective coupling with rainfall. The dynamical response of SISA and the impact on the South American summer monsoon (SASM) will be presented. Finally, we will give the concluding remarks.

Zhou, Jiayu↗

NASA NTRS2023

The Hazard Components of Representative Key Risks. The Physical Climate Perspective

The framework of Representative Key Risks (RKRs) has been adopted by the Intergovernmental Panel on Climate Change Working Group II (WGII) to categorize, assess and communicate a wide range of regional and sectoral key risks from climate change. These are risks expected to become severe due to the potentially detrimental convergence of changing climate conditions with the exposure and vulnerability of human and natural systems. Other papers in this special issue treat each of eight RKRs holistically by assessing their current status and future evolution as a result of this convergence. However, in these papers, such assessment cannot always be organized according to a systematic gradation of climatic changes. Often the big-picture evolution of risk has to be extrapolated from either qualitative effects of “low”, “medium” and “high” warming, or limited/focused analysis of the consequences of particular mitigation choices (e.g., benefits of limiting warming to 1.5 or 2C), together with consideration of the socio-economic context and possible adaptation choices. In this study we offer a representation – as systematic as possible given current literature and assessments – of the future evolution of the hazard components of RKRs. We identify the relevant hazards for each RKR, based upon the WGII authors’ assessment, and we report on their current state and expected future changes in magnitude, intensity and/or frequency, linking these changes to Global Warming Levels (GWLs) to the extent possible. We draw on the assessment of changes in climatic impact-drivers relevant to RKRs described in the 6th Assessment Report by Working Group I supplemented when needed by more recent literature. For some of these quantities - like regional trends in oceanic and atmospheric temperature and precipitation, some heat and precipitation extremes, permafrost thaw and Northern Hemisphere snow cover - a strong and quantitative relationship with increasing GWLs has been identified. For others - like frequency and intensity of tropical cyclones and extra-tropical storms, and fire weather - that link can only be described qualitatively. For some processes - like the behavior of ice sheets, or changes in circulation dynamics - large uncertainties about the effects of different GWLs remain, and for a few others - like ocean pH and air pollution - the composition of the scenario of anthropogenic emissions is most relevant, rather than the warming reached. In almost all cases, however, the basic message remains that every small increment in CO2 concentration in the atmosphere and associated warming will bring changes in climate phenomena that will contribute to increasing risk of impacts on human and natural systems, in the absence of compensating changes in these systems’ exposure and vulnerability, and in the absence of effective adaptation. Our picture of the evolution of RKR-relevant climatic impact-drivers complements and enriches the treatment of RKRs in the other papers in at least two ways: by filling in their often only cursory or limited representation of the physical climate aspects driving impacts, and by providing a fuller representation of their future potential evolution, an important component – if never the only one – of the future evolution of risk severity.

Representative Key Risks↗

NASA NTRS1986

Comment on 'Convergence characteristics of a vortex-lattice method for nonlinear configuration aerodynamics'

It is proposed that the study of Rusak et al. (1985), which reports numerical modeling sensitivities on longitudinal force/moment properties for a vortex-lattice method incorporating free vortex filaments to represent the leading-edge vortex separation, employs a formula that is strongly affected by the particular points of analysis chosen. This results in a narrowly applicable curve fit, where numerical sensitivities of the theory are inappropriately traded off against physical effects that are not modeled in that theory. Attention is also given to questionable drag estimate computations.

Luckring, J. M.↗