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 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗

NASA NTRSDate not supplied

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

Computational Icing↗

NASA NTRS1978

Implementation of natural frequency analysis and optimality criterion design

A description is presented of an effective optimality criterion computer design approach for member size selection to improve frequency characteristics for moderately large structure models. It is shown that the implementation of the simultaneous iteration method within a natural frequency structural design optimization provides a method which is more efficient in isolating the lowest natural frequency modes than the frequently applied Stodola method. Additional computational advantages are derived by using previously converged eigenvectors at the start of the iterations during the second and the following design cycles. Vectors with random components can be used at the first design cycle, which, in relation to the entire computer time for the design program, results in only a moderate computational penalty.

Levy, R.↗

NASA NTRSDate not supplied

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems↗

NASA NTRS2001

Interfacial Stability of Converging Plasma Jets for Magnetized Target Fusion

The merging of a spherical distribution of plasma jets to dynamically form a gaseous liner has been proposed for use in magnetized target fusion propulsion. In this paper, a study is made of the interfacial stability of the interaction of these jets. Specifically, the Orr-Sommerfeld equation is integrated to obtain the growth rate of a perturbation to the primary flow at the interface between the colliding jets. The analysis lead to an estimate on the tolerances on the relative flow velocities of the merging plasma jets to form a stable, imploding liner. The results show that during the merging of the jets to form a liner and before contact with the target plasma the growth of the perturbed flow at the jet interface is not likely to destabilize the liner. These data suggest that, as far as the stability of the interface between the merging jets is concerned, the formation of liner can withstand velocity variation up to 50% between the neighboring jets over the density and temperature ranges investigated.

Cassibry, J. T.↗

NASA NTRS1981

Transonic flow computations in cascades using finite element method

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and 1 deg.

Akay, H. U.↗

NASA NTRS1987

Transient response for interaction of two dynamic bodies

During the launch sequence of any space vehicle complicated boundary interactions occur between the vehicle and the launch stand. At the start of the sequence large forces exist between the two; contact is then broken in a short but finite time which depends on the release mechanism. The resulting vehicle response produces loads which are very high and often form the design case. It is known that the treatment of the launch pad as a second dynamic body is significant for an accurate prediction of launch response. A technique was developed for obtaining loads generated by the launch transient with the effect of pad dynamics included. The method solves uncoupled vehicle and pad equations of motion. The use of uncoupled models allows the simulation of vehicle launch in a single computer run. Modal formulation allows a closed-form solution to be written, eliminating any need for a numerical integration algorithm. When the vehicle is on the pad the uncoupled pad and vehicle equations have to be modified to account for the constraints they impose on each other. This necessitates the use of an iterative procedure to converge to a solution, using Lagrange multipliers to apply the required constraints. As the vehicle lifts off the pad the coupling between the vehicle and the pad is eliminated point by point until the vehicle flies free. Results obtained by this method were shown to be in good agreement with observed loads and other analysis methods. The resulting computer program is general, and was used without modification to solve a variety of contact problems.

Prabhakar, A.↗

NASA NTRS1993

Moist wind relationships

Equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are introduced. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term, while the other arises from temporal and spatial changes in the pressure gradient term. The influence of diabatic processes on the large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. From a scale analysis of the moisture flux equations it is argued that reductions in the large-scale horizontal wind speed, observed within major cumulus cloud systems, help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated. This wind modification induces secondary vertical circulations that contribute to the convergence, creating or supporting long-lived mesoscale flows. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.

Raymond, William H.↗

NASA NTRS1992

Flow analysis and design optimization methods for nozzle-afterbody of a hypersonic vehicle

This report summarizes the methods developed for the aerodynamic analysis and the shape optimization of the nozzle-afterbody section of a hypersonic vehicle. Initially, exhaust gases were assumed to be air. Internal-external flows around a single scramjet module were analyzed by solving the 3D Navier-Stokes equations. Then, exhaust gases were simulated by a cold mixture of Freon and Ar. Two different models were used to compute these multispecies flows as they mixed with the hypersonic airflow. Surface and off-surface properties were successfully compared with the experimental data. The Aerodynamic Design Optimization with Sensitivity analysis was then developed. Pre- and postoptimization sensitivity coefficients were derived and used in this quasi-analytical method. These coefficients were also used to predict inexpensively the flow field around a changed shape when the flow field of an unchanged shape was given. Starting with totally arbitrary initial afterbody shapes, independent computations were converged to the same optimum shape, which rendered the maximum axial thrust.

Baysal, O.↗

NASA NTRS2016

FUN3D Analyses in Support of the Second Aeroelastic Prediction Workshop

This paper presents the computational aeroelastic results generated in support of the second Aeroelastic Prediction Workshop for the Benchmark Supercritical Wing (BSCW) configurations and compares them to the experimental data. The computational results are obtained using FUN3D, an unstructured grid Reynolds- Averaged Navier-Stokes solver developed at NASA Langley Research Center. The analysis results include aerodynamic coefficients and surface pressures obtained for steady-state, static aeroelastic equilibrium, and unsteady flow due to a pitching wing or flutter prediction. Frequency response functions of the pressure coefficients with respect to the angular displacement are computed and compared with the experimental data. The effects of spatial and temporal convergence on the computational results are examined.

Chwalowski, Pawel↗

NASA NTRS1982

Analytical and experimental investigation of stator endwall contouring in a small axial-flow turbine. 1: Stator performance

Three stator configurations were studied to determine the effect of stator outer endwall contouring on stator performance. One configuration was a cylindrical stator design. One contoured stator configuration had an S-shaped outer endwall, the other had a conical-convergent outer endwall. The experimental investigation consisted of annular surveys of stator exit total pressure and flow angle for each stator configuration over a range of stator pressure ratio. Radial variations in stator loss and aftermixed flow conditions were obtained when these data were compared with the analytical results to assess the validity of the analysis, good agreement was found.

Haas, J. E.↗

NASA NTRS1985

Vectorizable multigrid algorithms for transonic flow calculations

The analysis and incorporation into a multigrid scheme of several vectorizable algorithms are discussed. Von Neumann analyses of vertical line, horizontal line, and alternating direction ZEBRA algorithms were performed; and the results were used to predict their multigrid damping rates. The algorithms were then successfully implemented in a transonic conservative full-potential computer program. The convergence acceleration effect of multiple grids is shown and the convergence rates of the vectorizable algorithms are compared to the convergence rates of standard successive line overrelaxation (SLOR) algorithms.

Melson, N. D.↗

NASA NTRS1986

Solution of elliptic PDEs by fast Poisson solvers using a local relaxation factor

A large class of two- and three-dimensional, nonseparable elliptic partial differential equations (PDEs) is presently solved by means of novel one-step (D'Yakanov-Gunn) and two-step (accelerated one-step) iterative procedures, using a local, discrete Fourier analysis. In addition to being easily implemented and applicable to a variety of boundary conditions, these procedures are found to be computationally efficient on the basis of the results of numerical comparison with other established methods, which lack the present one's: (1) insensitivity to grid cell size and aspect ratio, and (2) ease of convergence rate estimation by means of the coefficient of the PDE being solved. The two-step procedure is numerically demonstrated to outperform the one-step procedure in the case of PDEs with variable coefficients.

Chang, Sin-Chung↗

NASA NTRS1990

A coordination theory for intelligent machines

A formal model for the coordination level of intelligent machines is established. The framework of the coordination level investigated consists of one dispatcher and a number of coordinators. The model called coordination structure has been used to describe analytically the information structure and information flow for the coordination activities in the coordination level. Specifically, the coordination structure offers a formalism to (1) describe the task translation of the dispatcher and coordinators; (2) represent the individual process within the dispatcher and coordinators; (3) specify the cooperation and connection among the dispatcher and coordinators; (4) perform the process analysis and evaluation; and (5) provide a control and communication mechanism for the real-time monitor or simulation of the coordination process. A simple procedure for the task scheduling in the coordination structure is presented. The task translation is achieved by a stochastic learning algorithm. The learning process is measured with entropy and its convergence is guaranteed. Finally, a case study of the coordination structure with three coordinators and one dispatcher for a simple intelligent manipulator system illustrates the proposed model and the simulation of the task processes performed on the model verifies the soundness of the theory.

Wang, Fei-Yue↗

NASA NTRS1991

Flow analysis and design optimization methods for nozzle afterbody of a hypersonic vehicle

This report summarizes the methods developed for the aerodynamic analysis and the shape optimization of the nozzle-afterbody section of a hypersonic vehicle. Initially, exhaust gases were assumed to be air. Internal-external flows around a single scramjet module were analyzed by solving the three dimensional Navier-Stokes equations. Then, exhaust gases were simulated by a cold mixture of Freon and Argon. Two different models were used to compute these multispecies flows as they mixed with the hypersonic airflow. Surface and off-surface properties were successfully compared with the experimental data. In the second phase of this project, the Aerodynamic Design Optimization with Sensitivity analysis (ADOS) was developed. Pre and post optimization sensitivity coefficients were derived and used in this quasi-analytical method. These coefficients were also used to predict inexpensively the flow field around a changed shape when the flow field of an unchanged shape was given. Starting with totally arbitrary initial afterbody shapes, independent computations were converged to the same optimum shape, which rendered the maximum axial thrust.

Baysal, Oktay↗

NASA NTRS1993

Heating analysis for a Lunar Transfer Vehicle at near-equilibrium flow conditions

A heating analysis for a 15.2 m diameter Lunar Transfer Vehicle (LTV) at 0 and 10.6 deg angle of attack for a nominal trajectory through the earth's atmosphere is described. The analysis utilizes the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) with thin-layer, Navier-Stokes, thermochemical nonequilibrium options. Radiative heating levels are calculated using the Langley Optimized RAdiative Nonequilibrium (LORAN) and the Non-EQuilibrium AIr Radiation (NEQAIR) codes. At peak heating, the shock layer is substantially in equilibrium. Comprehensive spatial and spectral grid convergence studies have been implemented to quantify grid effects on the convective and radiative heating levels. Axisymmetric tests including the coupled effects of radiative energy transfer show negligible change to the convective heating and a 20 percent reduction in the radiative heating.

Gnoffo, Peter A.↗

NASA NTRS1994

An HP Adaptive Discontinuous Galerkin Method for Hyperbolic Conservation Laws

This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.

Bey, Kim S.↗

NASA NTRS2013

Application of Spectroscopic Doppler Velocimetry for Measurement of Streamwise Vorticity

A spectroscopic Doppler velocimetry technique has been developed for measuring two transverse components of velocity and hence streamwise vorticity in free jet flows. The nonintrusive optical measurement system uses Mie scattering from a 200 mW green continuous-wave laser interacting with dust and other tracer particulates naturally present in the air flow to measure the velocities. Scattered light is collected in two opposing directions to provide measurements of two orthogonal velocity components. An air-spaced Fabry-Perot interferometer is used for spectral analysis to determine the optical frequency shift between the incident laser light and the Mie scattered light. This frequency shift is directly proportional to the velocity component in the direction of the bisector of the incident and scattered light wave propagation vectors. Data were acquired for jet Mach numbers of 1.73 and 0.99 using a convergent 1.27-cm diameter round nozzle fitted with a single triangular "delta-tab". The velocity components and the streamwise vorticity calculated from the measurements are presented. The results demonstrate the ability of this novel optical system to obtain velocity and vorticity data without any artificial seeding and using a low power laser system.

Fagan, Amy↗