Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,495 records · Page 83

A synoptic analysis of the first AVE-SESAME '79 period

Key features of a severe convection observed during April 10-11, 1979 as part of the Atmospheric Variability Experiment-Severe Environmental Storms and Mesoscale Experiment (AVE-SESAME) are examined. Three-hourly rawinsonde readings from 23 stations were taken, and vertical motion and divergence parameters are considered. The data were converted into a 127 km grid at the surface, and at 50 mb intervals from 900 mb to 100 mb by an objective analysis scheme, while a kinematic method was used to compute vertical motion. A weak upper tropospheric short wave trough was found to propagate from New Mexico into the Texas panhandle, while a jet maximum propagated eastward. The development of a strong wind maximum over Oklahoma and Kansas was associated with a rapid increase in upper-level divergence and the development of a small-scale pressure perturbation in the Texas panhandle, as well as a low-level jet and convergence, which led to rapid changes over the Red River Valley, where stability was decreased.

Moore, J. T.↗

Evaluation and modeling of autonomous attitude thrust control for the Geostation Operational Environmental Satellite (GOES)-8 orbit determination

The Geostationary Operational Environmental Satellite (GOES)-8 experienced a series of orbital perturbations from autonomous attitude control thrusting before perigee raising maneuvers. These perturbations influenced differential correction orbital state solutions determined by the Goddard Space Flight Center (GSFC) Goddard Trajectory Determination System (GTDS). The maneuvers induced significant variations in the converged state vector for solutions using increasingly longer tracking data spans. These solutions were used for planning perigee maneuvers as well as initial estimates for orbit solutions used to evaluate the effectiveness of the perigee raising maneuvers. This paper discusses models for the incorporation of attitude thrust effects into the orbit determination process. Results from definitive attitude solutions are modeled as impulsive thrusts in orbit determination solutions created for GOES-8 mission support. Due to the attitude orientation of GOES-8, analysis results are presented that attempt to absorb the effects of attitude thrusting by including a solution for the coefficient of reflectivity, C(R). Models to represent the attitude maneuvers are tested against orbit determination solutions generated during real-time support of the GOES-8 mission. The modeling techniques discussed in this investigation offer benefits to the remaining missions in the GOES NEXT series. Similar missions with large autonomous attitude control thrusting, such as the Solar and Heliospheric Observatory (SOHO) spacecraft and the INTELSAT series, may also benefit from these results.

Forcey, W.↗

Characteristics of the field-aligned current system in the nighttime sector during auroral substorms

Fujii et al. (1994) obtained characteristics of the electrodynamic parameters, that is, field-aligned currents, electric fields, and electron precipitation, which are associted with auroral substorm events in the nighttime sector, through a unique analysis that places the ionospheric measurements of these parameters into the context of a generic substorm determined from global auroral images. In this paper we investigate in considerably more detail the characteristics of the field-aligned currents using data from the same set of passes as the previuos study. We show for the first time that the net upward field-aligned currents throughout the surge and surge horn are sufficient to account for most if not all of the converging currents of the auroral electrojets. Current densities are largest in the surge and surge horn. Current region continuity does not appear to exist across the substorm bulge region. Much of the auroral substorm field-aligned current is composed of filamentary currents and finite current segments at large angles to each other. The westward electrojet may contain large gradients in intensity both in local time and latitude due to sets of localized field-aligned currents. The net downward current for several hours to the west of the surge is insufficient to account for the eastward electrojet, consistent with the concept that this electrojet originates primarily on the dayside. Our pattern of field-aligned currents associated with the surge has common features and also differs significantly from the patterns previously derived from data from radars and ground-based magnetometer arrays. Our pattern is considerably more complex, probably due to the much higher resolution in latitude of the satellite data. It is also larger in area, since our average substorm is much larger than those pertaining to the previous patterns, giving a substorm wedge considerably wider than that obtained from the radar and array data.

Hoffman, R. A.↗

Analysis of hourglass instabilities and control in underintegrated finite element methods

Belytschko et al. (1981, 1984) has developed stabilization methods for the treatment of underintegrated FEM problems; these methods involve the computation of an underintegrated stiffness matrix, which is rank-deficient, and the addition of a stabilization matrix which effectively eliminates the spurious modes. An attempt is presently made to give this a priori stabilization method a mathematical means of support. Attention is also given to an a posteriority stabilization method for hourglass control, in which an approximate solution of the underintegrated system is obtained and then subjected to a special projection in order to eliminate the hourglass modes. A proof is obtained for the convergence of this stabilized underintegrated approximation to the exact solution of a model problem at almost the same rate (as the mesh is refined) as the fully integrated solutions.

Jacquotte, O.-P.↗

Development of Renormalization Group Analysis of Turbulence

The renormalization group (RG) procedure for nonlinear, dissipative systems is now quite standard, and its applications to the problem of hydrodynamic turbulence are becoming well known. In summary, the RG method isolates self similar behavior and provides a systematic procedure to describe scale invariant dynamics in terms of large scale variables only. The parameterization of the small scales in a self consistent manner has important implications for sub-grid modeling. This paper develops the homogeneous, isotropic turbulence and addresses the meaning and consequence of epsilon-expansion. The theory is then extended to include a weak mean flow and application of the RG method to a sequence of models is shown to converge to the Navier-Stokes equations.

L M Smith↗

Applications of singular value analysis and partial-step algorithm for nonlinear orbit determination

An adaptive method in which cruise and nonlinear orbit determination problems can be solved using a single program is presented. It involves singular value decomposition augmented with an extended partial step algorithm. The extended partial step algorithm constrains the size of the correction to the spacecraft state and other solve-for parameters. The correction is controlled by an a priori covariance and a user-supplied bounds parameter. The extended partial step method is an extension of the update portion of the singular value decomposition algorithm. It thus preserves the numerical stability of the singular value decomposition method, while extending the region over which it converges. In linear cases, this method reduces to the singular value decomposition algorithm with the full rank solution. Two examples are presented to illustrate the method's utility.

Ryne, Mark S.↗

A mathematical analysis of the theory of interplanetary scintillation in the weak scattering approximation

A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.

Mitchell, D. G.↗

BeyondPlanck: XII. Cosmological parameter constraints with end-to-end error propagation

We present cosmological parameter constraints estimated using the Bayesian BEYONDPLANCK analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data onto final cosmological parameters. As a first demonstration of the method, we analyzed time-ordered Planck LFI observations, combined with selected external data (WMAP 33–61 GHz, Planck HFI DR4 353 and 857 GHz, and Haslam 408 MHz) in the form of pixelized maps that are used to break critical astrophysical degeneracies. Overall, all the results are generally in good agreement with previously reported values from Planck 2018 and WMAP, with the largest relative difference for any parameter amounting about 1σ when considering only temperature multipoles between 30 ≤ ℓ ≤ 600. In cases where there are differences, we note that the BEYONDPLANCK results are generally slightly closer to the high-ℓ HFI-dominated Planck 2018 results than previous analyses, suggesting slightly less tension between low and high multipoles. Using low-ℓ polarization information from LFI and WMAP, we find a best-fit value of τ = 0.066 ± 0.013, which is higher than the low value of τ = 0.052 ± 0.008 derived from Planck 2018 and slightly lower than the value of 0.069 ± 0.011 derived from the joint analysis of official LFI and WMAP products. Most importantly, however, we find that the uncertainty derived in the BEYONDPLANCK processing is about 30 % greater than when analyzing the official products, after taking into account the different sky coverage. We argue that this uncertainty is due to a marginalization over a more complete model of instrumental and astrophysical parameters, which results in more reliable and more rigorously defined uncertainties. We find that about 2000 Monte Carlo samples are required to achieve a robust convergence for a low-resolution cosmic microwave background (CMB) covariance matrix with 225 independent modes, and producing these samples takes about eight weeks on a modest computing cluster with 256 cores.

79 ASTRONOMY AND ASTROPHYSICS↗

SASS wind ambiguity removal by direct minimization

An objective analysis procedure is presented which combines Seasat-A satellite scatterometer (SASS) data with other available data on wind speeds by minimizing an objective function of gridded wind speed values. The functions are defined as the loss functions for the SASS velocity data, the forecast, the SASS velocity magnitude data, and conventional wind speed data. Only aliases closest to the analysis were included, and a method for improving the first guess while using a minimization technique and slowly changing the parameters of the problem is introduced. The model is employed to predict the wind field for the North Atlantic on Sept. 10, 1978. Dealiased SASS data is compared with available ship readings, showing good agreement between the SASS dealiased winds and the winds measured at the surface. Expansion of the model to take in low-level cloud measurements, pressure data, and convergence and cloud level data correlations is discussed.

Hoffman, R. N.↗

The tectonic setting of the Seychelles, Mascarene and Amirante Plateaus in the Western Equatorial Indian Ocean

A system of marine plateaus occurs in the western equatorial Indian Ocean, forming an arcuate series of wide and shallow banks with small islands in places. The oceanic basins that surround the Seychelles - Amirante region are of various ages and reflect a complex seafloor spreading pattern. The structural analysis of the Seychelle - Amirante - Mascarene region reflects the tectonic evolution of the western equatorial Indian Ocean. It is suggested that due to the seafloor spreading during a tectonic stage, the Seychelles continental block drifted southwestwards to collide with the oceanic crust of the Mascarene Basin, forming an elongated folded structure at first, and then a subduction zone. The morphological similarity, the lithological variability and the different origin of the Seychelles Bank, the Mascarene Plateau and the Amirante Arc emphasizes the significant convergent effects of various plate tectonic processes on the development of marine plateaus.

Mart, Y.↗

Aeroelastic stability analyses of two counter rotating propfan designs for a cruise missile model

Aeroelastic stability analyses were performed to insure structural integrity of two counterrotating propfan blade designs for a NAVY/Air Force/NASA cruise missile model wind tunnel test. This analysis predicted if the propfan designs would be flutter free at the operating conditions of the wind tunnel test. Calculated stability results are presented for the two blade designs with rotational speed and Mach number as the parameters. A aeroelastic analysis code ASTROP2 (Aeroelastic Stability and Response of Propulsion Systems - 2 Dimensional Analysis), developed at LeRC, was used in this project. The aeroelastic analysis is a modal method and uses the combination of a finite element structural model and two dimensional steady and unsteady cascade aerodynamic models. This code was developed to analyze single rotation propfans but was modified and applied to counterrotating propfans for the present work. Modifications were made to transform the geometry and rotation of the aft rotor to the same reference frame as the forward rotor, to input a non-uniform inflow into the rotor being analyzed, and to automatically converge to the least stable aeroelastic mode.

Mahajan, Aparajit J.↗

Single-engine tail interference model

The data presented in this contribution were obtained in the NASA Langley 16-Foot Transonic Tunnel. Multiple test entries were completed and the results have been completely reported in five NASA reports. The objective of the initial investigation was to determine the effect of empennage (tail) interference on the drag characteristics of an axisymmetric model with a single engine fighter aft-end with convergent divergent nozzles. Two nozzle power settings, dry and maximum afterburning, were investigated. Several empennage arrangements and afterbody modifications were investigated during the initial investigation. Subsequent investigations were used to determine the effects of other model variables including tail incidence, tail span, and nozzle shape. For the final investigation, extensive surface pressure instrumentation was added to the model in order to develop and understanding of the flow interactions associated with afterbody/empennage integration and also to provide data for code validation. Extensive computational analysis has been conducted on the staggered empennage configuration at a Mach number of 0.6 utilizing a three-dimensional Navier Stokes code. Most of the investigations were conducted at Mach numbers from 0.60 to 1.20 and at ratios of jet total pressure to free stream static pressure (nozzle pressure ratio) from 0.1 (jet off) to 8.0. Some angle of attack variation was obtained at jet off conditions.

Bobby L. Berrier↗

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗

Aeroelastic Calculations Based on Three-Dimensional Euler Analysis

This paper presents representative results from an aeroelastic code (TURBO-AE) based on an Euler/Navier-Stokes unsteady aerodynamic code (TURBO). Unsteady pressure, lift, and moment distributions are presented for a helical fan test configuration which is used to verify the code by comparison to two-dimensional linear potential (flat plate) theory. The results are for pitching and plunging motions over a range of phase angles, Good agreement with linear theory is seen for all phase angles except those near acoustic resonances. The agreement is better for pitching motions than for plunging motions. The reason for this difference is not understood at present. Numerical checks have been performed to ensure that solutions are independent of time step, converged to periodicity, and linearly dependent on amplitude of blade motion. The paper concludes with an evaluation of the current state of development of the TURBO-AE code and presents some plans for further development and validation of the TURBO-AE code.

Bakhle, Milind A.↗

A Collaborative Analysis Tool for Integrated Hypersonic Aerodynamics, Thermal Protection Systems, and RBCC Engine Performance for Single Stage to Orbit Vehicles

Presented is a computer-based tool that connects several disciplines that are needed in the complex and integrated design of high performance reusable single stage to orbit (SSTO) vehicles. Every system is linked to every other system, as is the case of SSTO vehicles with air breathing propulsion, which is currently being studied by NASA. An RBCC propulsion system integrates airbreathing and rocket propulsion into a single engine assembly enclosed within a cowl or duct. A typical RBCC propulsion system operates as a ducted rocket up to approximately Mach 3. Then there is a transition to a ramjet mode for supersonic-to-hypersonic acceleration. Around Mach 8 the engine transitions to a scramjet mode. During the ramjet and scramjet modes, the integral rockets operate as fuel injectors. Around Mach 10-12 (the actual value depends on vehicle and mission requirements), the inlet is physically closed and the engine transitions to an integral rocket mode for orbit insertion. A common feature of RBCC propelled vehicles is the high degree of integration between the propulsion system and airframe. At high speeds the vehicle forebody is fundamentally part of the engine inlet, providing a compression surface for air flowing into the engine. The compressed air is mixed with fuel and burned. The combusted mixture must be expanded to an area larger than the incoming stream to provide thrust. Since a conventional nozzle would be too large, the entire lower after body of the vehicle is used as an expansion surface. Because of the high external temperatures seen during atmospheric flight, the design of an airbreathing SSTO vehicle requires delicate tradeoffs between engine design, vehicle shape, and thermal protection system (TPS) sizing in order to produce an optimum system in terms of weight (and cost) and maximum performance. To adequately determine the performance of the engine/vehicle, the Hypersonic Flight Inlet Model (HYFIM) module was designed to interface with the RBCC engine model. HYFIM performs the aerodynamic analysis of forebodies and inlet characteristics of RBCC powered SSTO launch vehicles. HYFIM is applicable to the analysis of the ramjet/scramjet engine operations modes (Mach 3-12), and provides estimates of parameters such as air capture area, shock-on-lip Mach number, design Mach number, compression ratio, etc., based on a basic geometry routine for modeling axisymmetric cones, 2-D wedge geometries. HYFIM also estimates the variation of shock layer properties normal to the forebody surface. The thermal protection system (TPS) is directly linked to determination of the vehicle moldline and the shaping of the trajectory. Thermal protection systems to maintain the structural integrity of the vehicle must be able to mitigate the heat transfer to the structure and be lightweight. Herein lies the interdependency, in that as the vehicle's speed increases, the TPS requirements are increased. And as TPS masses increase the effect on the propulsion system and all other systems is compounded. The need to analyze vehicle forebody and engine inlet is critical to be able to design the RBCC vehicle. To adequately determine insulation masses for an RBCC vehicle, the hypersonic aerodynamic environment and aeroheating loads must be calculated and the TPS thicknesses must be calculated for the entire vehicle. To accomplish this an ascent or reentry trajectory is obtained using the computer code Program to Optimize Simulated Trajectories (POST). The trajectory is then used to calculate the convective heat rates on several locations on the vehicles using the Miniature Version of the JA70 Aerodynamic Heating Computer Program (MINIVER). Once the heat rates are defined for each body point on the vehicle, then insulation thicknesses that are required to maintain the vehicle within structural limits are calculated using Systems Improved Numerical Differencing Analyzer (SINDA) models. If the TPS masses are too heavy for the performance of the vehicle the process may be repeated altering the trajectory or some other input to reduce the TPS mass. E-PSURBCC is an "engine performance" model and requires the specification of inlet air static temperature and pressure as well as Mach number (which it pulls from the HYFIM and POST trajectory files), and calculates the corresponding stagnation properties. The engine air flow path geometry includes inlet, a constant area section where the rocket is positioned, a subsonic diffuser, a constant area afterburner, and either a converging nozzle or a converging-diverging nozzle. The current capabilities of E-PSURBCC ejector and ramjet mode treatment indicated that various complex flow phenomena including multiple choking and internal shocks can occur for combinations of geometry/flow conditions. For a given input deck defining geometry/flow conditions, the program first goes through a series of checks to establish whether the input parameters are sound in terms of a solution path. If the vehicle/engine performance fails mission goals, the engineer is able to collaboratively alter the vehicle moldline to change aerodynamics, or trajectory, or some other input to achieve orbit. The problem described is an example of the need for collaborative design and analysis. RECIPE is a cross-platform application capable of hosting a number of engineers and designers across the Internet for distributed and collaborative engineering environments. Such integrated system design environments allow for collaborative team design analysis for performing individual or reduced team studies. To facilitate the larger number of potential runs that may need to be made, RECIPE connects the computer codes that calculate the trajectory data, aerodynamic data based on vehicle geometry, heat rate data, TPS masses, and vehicle and engine performance, so that the output from each tool is easily transferred to the model input files that need it.

Stanley, Thomas Troy↗

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

Determination of trace gases in learjet and U-2 whole air samples collected during the intertropical convergence zone experiment

The contribution of Washington State University (WSU) to the Intertropical Convergence Zone (ITCZ) study, conducted in the Panama Canal Zone in July 1977 is reported. The ITCZ program had the objective to study the meteorology and air chemistry to determine, if possible, the magnitude of injection of tropospheric air into the stratosphere in the ITCZ. WSU was responsible for two phases of the air chemistry program for the ITCZ study: (1) the collection of whole-air samples from ground level to 13.7 km using a Learjet as a sampling platform, and the analysis of the samples for selected halocarbons, C2 hydrocarbons, N2O, and SF6; and (2) the analysis for selected halocarbon species and N2O of low pressure whole-air samples, collected between 13.7 km and 21.3 km by a U-2 aircraft. The report is divided into separate sections describing the Learjet and the U-2 sampling methods, analytical methods, and results. Complete data tables for all samples analyzed are given.

Cronn, D. R.↗

Human Tracking Performance in Uncoupled and Coupled Two-Axis Systems

This report presents tile results of an experimental and analytical study of human performance in uncoupled and coupled control systems. Human pilot performance in single and two-axis systems was mathematically modeled by linear second-order describing functions. Model parameters were determined using model matching techniques. Analysis of the models showed that the amplitude ratio and phase lead of the describing function increased with training indicating an increase in open loop bandwidth. The phase margin also decreased with training. Increasing the plant lag time constant resulted in an increase in the model lead time constant and a decrease in the zero frequency gain. No significant difference was found to exist in the normalized tracking error per axis between the two-axis tasks and the single-axis tasks. However tile model lead time constant was significantly greater in two-axis tracking. Manual tracking of two-axis systems with cross-coupling was studied experimentally and analytically. Approximate methods for modeling two-axis performance were developed and checked using a precise spectral analysis approach. Coupled and uncoupled, symmetrical and asymmetrical two-axis performance was compared. The results show that modeling of cross-coupled systems is feasible and that trained subjects are capable of decoupling the axes of some systems. A methodology study compared the identification performance of continuous, iterative, and extrapolation model matching techniques. An iterative technique employing sensitivity equations for the generation of influence coefficients was found to be the best technique due to its excellent identification accuracy and ease of implementation. Convergence in iterative techniques can be improved substantially by equalizing the parameter adjustment rates and limiting the maximum correction per iteration.

Todosiev, E. P.↗