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 721 records · Page 40

Simulation of Low-density Nozzle Plumes in Non-zero Ambient Pressures

The direct simulation Monte-Carlo (DSMC) method was applied to the analysis of low-density nitrogen plumes exhausting from a small converging-diverging nozzle into finite ambient pressures. Two cases were considered that simulated actual test conditions in a vacuum facility. The numerical simulations readily captured the complicated flow structure of the overexpanded plumes adjusting to the finite ambient pressures, including Mach disks and barrel shaped shocks. The numerical simulations compared well to experimental data of Rothe.

Chung, Chan-Hong↗

Multidisciplinary design optimization using genetic algorithms

Multidisciplinary design optimization (MDO) is an important step in the conceptual design and evaluation of launch vehicles since it can have a significant impact on performance and life cycle cost. The objective is to search the system design space to determine values of design variables that optimize the performance characteristic subject to system constraints. Gradient-based optimization routines have been used extensively for aerospace design optimization. However, one limitation of gradient based optimizers is their need for gradient information. Therefore, design problems which include discrete variables can not be studied. Such problems are common in launch vehicle design. For example, the number of engines and material choices must be integer values or assume only a few discrete values. In this study, genetic algorithms are investigated as an approach to MDO problems involving discrete variables and discontinuous domains. Optimization by genetic algorithms (GA) uses a search procedure which is fundamentally different from those gradient based methods. Genetic algorithms seek to find good solutions in an efficient and timely manner rather than finding the best solution. GA are designed to mimic evolutionary selection. A population of candidate designs is evaluated at each iteration, and each individual's probability of reproduction (existence in the next generation) depends on its fitness value (related to the value of the objective function). Progress toward the optimum is achieved by the crossover and mutation operations. GA is attractive since it uses only objective function values in the search process, so gradient calculations are avoided. Hence, GA are able to deal with discrete variables. Studies report success in the use of GA for aircraft design optimization studies, trajectory analysis, space structure design and control systems design. In these studies reliable convergence was achieved, but the number of function evaluations was large compared with efficient gradient methods. Applicaiton of GA is underway for a cost optimization study for a launch-vehicle fuel-tank and structural design of a wing. The strengths and limitations of GA for launch vehicle design optimization is studied.

Unal, Resit↗

Status of Axisymmetric CFD of an Eleven Inch Diameter Hybrid Rocket Motor

Current status of a steady state, axisymmetric analysis of an experimental 11 inch diameter hybrid rocket motor internal flow field is given. The objective of this effort is to develop a steady state axisymmetric model of the 11 inch hybrid rocket motor which can be used as a design and/or analytical tool. A test hardware description, modeling approach, and future plans are given. The analysis was performed with FDNS implementing several finite rate chemistry sets. A converged solution for a two equation and five species set on a 'fine' grid is shown.

Ruf, Joseph↗

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗

An Investigation of the Influence of Urban Areas on Rainfall Using a Cloud-Mesoscale Model and the TRMM Satellite

A recent paper by Shepherd and Pierce (conditionally accepted to Journal of Applied Meteorology) used rainfall data from the Precipitation Radar on NASA's Tropical Rainfall Measuring Mission's (TRMM) satellite to identify warm season rainfall anomalies downwind of major urban areas. A convective-mesoscale model with extensive land-surface processes is employed to (a) determine if an urban heat island (UHI) thermal perturbation can induce a dynamic response to affect rainfall processes and (b) quantify the impact of the following three factors on the evolution of rainfall: (1) urban surface roughness, (2) magnitude of the UHI temperature anomaly, and (3) physical size of the UHI temperature anomaly. The sensitivity experiments are achieved by inserting a slab of land with urban properties (e.g. roughness length, albedo, thermal character) within a rural surface environment and varying the appropriate lower boundary condition parameters. Early analysis suggests that urban surface roughness (through turbulence and low-level convergence) may control timing and initial location of UHI-induced convection. The magnitude of the heat island appears to be closely linked to the total rainfall amount with minor impact on timing and location. The physical size of the city may predominantly impact on the location of UHI-induced rainfall anomaly. The UHI factor parameter space will be thoroughly investigated with respect to their effects on rainfall amount, location, and timing. This study extends prior numerical investigations of the impact of urban surfaces on meteorological processes, particularly rainfall development. The work also contains several novel aspects, including the application of a high-resolution (less than I km) cloud-mesoscale model to investigate urban-induce rainfall process; investigation of thermal magnitude of the UHI on rainfall process; and investigation of UHI physical size on rainfall processes.

Shepherd, J. Marshall↗

Overview of the Development for a Suite of Low-Thrust Trajectory Analysis Tools

A NASA intercenter team has developed a suite of low-thrust trajectory analysis tools to make a significant improvement in three major facets of low-thrust trajectory and mission analysis. These are: 1) ease of use, 2) ability to more robustly converge to solutions, and 3) higher fidelity modeling and accuracy of results. Due mostly to the short duration of the development, the team concluded that a suite of tools was preferred over having one integrated tool. This tool-suite, their characteristics, and their applicability will be described. Trajectory analysts can read this paper and determine which tool is most appropriate for their problem.

Kos, Larry D.↗

Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.

Meelan Choudhari↗

Anthropogenic Sulfate Aerosol Pollution in South and East Asia Induces Increased Summer Precipitation Over Arid Central Asia

Precipitation has increased across the arid Central Asia region over recent decades. However, the underlying mechanisms of this trend are poorly understood. Here, we analyze multi-model simulations from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP) to investigate potential drivers of the observed precipitation trend. We find that anthropogenic sulfate aerosols over remote polluted regions in South and East Asia lead to increased summer precipitation, especially convective and extreme precipitation, in arid Central Asia. Elevated concentrations of sulfate aerosols over remote polluted Asia cause an equatorward shift of the Asian Westerly Jet Stream through a fast response to cooling of the local atmosphere at mid-latitudes. This shift favours moisture supply from low-latitudes and moisture flux convergence over arid Central Asia, which is confirmed by a moisture budget analysis. High levels of absorbing black carbon lead to opposing changes in the Asian Westerly Jet Stream and reduced local precipitation, which can mask the impact of sulfate aerosols. This teleconnection between arid Central Asia precipitation and anthropogenic aerosols in remote Asian polluted regions highlights long-range impacts of anthropogenic aerosols on atmospheric circulations and the hydrological cycle.

Atmospheric dynamics↗

The effects of temperatures on the pebble flow in a pebble bed high temperature reactor

The core of a pebble bed high temperature reactor (PBHTR) moves during operation, a feature which leads to better fuel economy (online refueling with no burnable poisons) and lower fuel stress. The pebbles are loaded at the top and trickle to the bottom of the core after which the burnup of each is measured. The pebbles that are not fully burned are recirculated through the core until the target burnup is achieved. The flow pattern of the pebbles through the core is of importance for core simulations because it couples the burnup distribution to the core temperature and power profiles, especially in cores with two or more radial burnup 'zones '. The pebble velocity profile is a strong function of the core geometry and the friction between the pebbles and the surrounding structures (other pebbles or graphite reflector blocks). The friction coefficient for graphite in a helium environment is inversely related to the temperature. The Thorium High Temperature Reactor (THTR) operated in Germany between 1983 and 1989. It featured a two-zone core, an inner core (IC) and outer core (OC), with different fuel mixtures loaded in each zone. The rate at which the IC was refueled relative to the OC in THTR was designed to be 0.56. During its operation, however, this ratio was measured to be 0.76, suggesting the pebbles in the inner core traveled faster than expected. It has been postulated that the positive feedback effect between inner core temperature, burnup, and pebble flow was underestimated in THTR. Because of the power shape, the center of the core in a typical cylindrical PBHTR operates at a higher temperature than the region next to the side reflector. The friction between pebbles in the IC is lower than that in the OC, perhaps causing a higher relative flow rate and lower average burnup, which in turn yield a higher local power density. Furthermore, the pebbles in the center region have higher velocities than the pebbles next to the side reflector due to the interaction between the pebbles and the immobile graphite reflector as well as the geometry of the discharge conus near the bottom of the core. In this paper, the coupling between the temperature profile and the pebble flow dynamics was analyzed by using PEBBED/THERMIX and PEBBLES codes by modeling the HTR-10 reactor in China. Two extreme and opposing velocity profiles are used as a starting point for the iterations. The PEBBED/THERMIX code is used to calculate the burnup, power and temperature profiles with one of the velocity profiles as input. The resulting temperature profile is then passed to PEBBLES code to calculate the updated pebble velocity profile taking the new temperature profile into account. If the aforementioned hypothesis is correct, the strong temperature effect upon the friction coefficients would cause the two cases to converge to different final velocity and temperature profiles. The results of this analysis indicates that a single zone pebble bed core is self-stabilizing in terms of the pebble velocity profile and the effect of the temperature profile on the pebble flow is insignificant. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A diffusion synthetic acceleration scheme for rectangular geometries based on bilinear discontinuous finite elements

A DSA technique to accelerate the iterative convergence of S{sub n} transport solves is derived for bilinear discontinuous (BLD) finite elements on rectangular grids. The diffusion synthetic acceleration equations are discretized using BLD elements by adapting the Modified Interior Penalty technique, introduced in [4] for triangular grids. The MIP-DSA equations are SPD and thus are solved using a preconditioned CG technique. Fourier analyses and implementation of the technique in a BLD S{sub n} transport code show that the technique is stable is effective. (authors)

97 MATHEMATICS AND COMPUTING↗

Vibration and stress analysis of soft-bonded shuttle insulation tiles. Modal analysis with compact widely space stringers

An efficient iterative procedure is described for the vibration and modal stress analysis of reusable surface insulation (RSI) of multi-tiled space shuttle panels. The method, which is quite general, is rapidly convergent and highly useful for this application. A user-oriented computer program based upon this procedure and titled RESIST (REusable Surface Insulation Stresses) has been prepared for the analysis of compact, widely spaced, stringer-stiffened panels. RESIST, which uses finite element methods, obtains three dimensional tile stresses in the isolator, arrestor (if any) and RSI materials. Two dimensional stresses are obtained in the tile coating and the stringer-stiffened primary structure plate. A special feature of the program is that all the usual detailed finite element grid data is generated internally from a minimum of input data. The program can accommodate tile idealizations with up to 850 nodes (2550 degrees-of-freedom) and primary structure idealizations with a maximum of 10,000 degrees-of-freedom. The primary structure vibration capability is achieved through the development of a new rapid eigenvalue program named ALARM (Automatic LArge Reduction of Matrices to tridiagonal form).

Ojalvo, I. U.↗

Blade design and analysis using a modified Euler solver

An iterative method for blade design based on Euler solver and described in an earlier paper is used to design compressor and turbine blades providing shock free transonic flows. The method shows a rapid convergence, and indicates how much the flow is sensitive to small modifications of the blade geometry, that the classical iterative use of analysis methods might not be able to define. The relationship between the required Mach number distribution and the resulting geometry is discussed. Examples show how geometrical constraints imposed upon the blade shape can be respected by using free geometrical parameters or by relaxing the required Mach number distribution. The same code is used both for the design of the required geometry and for the off-design calculations. Examples illustrate the difficulty of designing blade shapes with optimal performance also outside of the design point.

Leonard, O.↗

Prediction of Individual Social-Demographic Role Based on Travel Behavior Variability Using Long-Term GPS Data

With the development of and advances in smartphones and global positioning system (GPS) devices, travelers’ long-term travel behaviors are not impossible to obtain. This study investigates the pattern of individual travel behavior and its correlation with social-demographic features. For different social-demographic groups (e.g., full-time employees and students), the individual travel behavior may have specific temporal-spatial-mobile constraints. The study first extracts the home-based tours, including Home-to-Home and Home-to-Non-Home, from long-term raw GPS data. The travel behavior pattern is then delineated by home-based tour features, such as departure time, destination location entropy, travel time, and driving time ratio. The travel behavior variability describes the variances of travelers’ activity behavior features for an extended period. After that, the variability pattern of an individual’s travel behavior is used for estimating the individual’s social-demographic information, such as social-demographic role, by a supervised learning approach, support vector machine. In this study, a long-term (18-month) recorded GPS data set from Puget Sound Regional Council is used. The experiment’s result is very promising. The sensitivity analysis shows that as the number of tours thresholds increases, the variability of most travel behavior features converges, while the prediction performance may not change for the fixed test data.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗

Boundary-Layer Transition Prediction Through Loose Coupling of OVERFLOW and LASTRAC

Transition prediction based on linear stability theory is expected to more accurately reflect the causality of transition onset than phenomenological transition models based on RANS-like transport equations. To help achieve the CFD vision 2030 aim of building a CFD tool chain with automated prediction of boundary layer transition, a technique to loosely tie the NASA OVERFLOW CFD solver with the LASTRAC stability analysis tool is described. The coupled solver is then used to compute transition over selected over a flat plate in a freestream with sufficiently low levels of turbulence, NLF(1)-0416 airfoil, the 6:1 prolate spheroid at an angle of attack, and a NASA juncture flow model with symmetric wing configuration. The findings show that the loosely coupled approach can reliably predict the transition location accurately in scenarios that are dominated by a single transition mechanism involving Tollmien-Schlichting instabilities, crossflow instabilities, or separation bubble-induced transition, or include a mixture of selected mechanisms. The toolset presents here appears to be robust to the prescription of the initial transition location, and it can lead to a converged solution in four or five rounds of the mean flow calculation and stability analysis, with minimal input from the user.

boundary layer transition↗

Advanced Capability of Monte Carlo Transport and Depletion Calculation Program SuperMC for Fusion and Fission Applications

Since the advantages of high simulation accuracy and flexible geometry representation, the Monte Carlo methods have been broadly adopted in nuclear design and analysis of nuclear energy systems. However, great challenges exist for the current Monte Carlo methods including the calculation modeling of complex geometries, slow convergence of complex calculation, etc. Super Monte Carlo Simulation Program for Nuclear and Radiation Process (SuperMC), developed by FDS team, is a general purpose, intelligent and multi-functional program for the design and safety evaluation of nuclear systems, which is designed to perform the comprehensive neutronics calculation. The latest version of SuperMC can accomplish the transport calculation of n,{gamma} for criticality and shielding problems, depletion calculation of isotopes, and has been integrated with the functions of automatic modeling, visualization and cloud computing. The latest progress will be introduced in this paper. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Structural Analysis of Test Flight Vehicles for Application of Multifunctional Energy Storage System

Under the NASA Aeronautics Research Mission Directorate (ARMD) Convergent Aeronautical Solutions (CAS) project, NASA Glenn Research Center has been leading Multifunctional Structures for High Energy Lightweight Load-bearing Storage (M-SHELLS) research efforts. The technology of integrating load-carrying structures with electrical energy storage capacity has the potential to reduce the overall weight of future electric aircraft. The proposed project goals were to develop M-SHELLS in the form of honeycomb coupons and subcomponents, integrate them into the structure, and conduct low-risk flight-tests onboard a remotely piloted small aircraft. Experimental M-SHELLS energy-storing coupons were fabricated and tested in the laboratory for their electrical and mechanical properties. In this report, finite element model development and structural analyses of two small test aircraft candidates are presented. The finite element analysis of the initial two-spar wing is described for strain, deflection, and weight estimation. After a test aircraft Tempest was acquired, a load-deflection test of the wing was conducted. A finite element model of the Tempest was then developed based on the test aircraft dimensions and construction detail. The component weight analyses from the finite element model and test measurements were correlated. Structural analysis results with multifunctional energy storage panels in the fuselage of the test vehicle are presented. Although the flight test was cancelled because of programmatic reasons and time constraints, the structural analysis results indicate that the mid-fuselage floor composite panel could provide structural integrity with minimal weight penalty while supplying electrical energy. To explore potential future applications of the multifunctional structure, analyses of the NASA X-57 Maxwell electric aircraft and a NASA N+3 Technology Conventional Configuration (N3CC) fuselage are presented. Secondary aluminum structures in the fuselage sub-floor and cargo area were partially replaced with reinforced five-layer composite panels with M-SHELLS honeycomb core. The N3CC fuselage weight reduction associated with each design without risking structural integrity are described. The structural analysis and weight estimation with the application of composite M-SHELLS panels to the N3CC fuselage indicate a 3.2% reduction in the fuselage structural weight, prior to accounting for the additional weight of core material required to complete the energy storage functionality.

Mukhopadhyay, Vivek↗

Ensuring Solution Uniqueness in Three-Phase Power System State Estimation

This paper is concerned with the issue of potential non-unique solutions in three-phase state estimation. Theory of observability analysis for positive sequence power system state estimation is based on certain assumptions that avoid possibility of multiple solutions. Also, it is shown that observability of a positive sequence network remains independent of the network parameters or the operating state. When extending single-phase observability analysis directly to the three-phase case, this paper considers the possibility of converging to multiple solutions, i.e. solution non-uniqueness, even for cases where state estimator successfully converges. The study illustrates via numerical examples the likelihood of converging to entirely different solutions for certain network parameters. It also examines how the operating state, particularly under unbalanced loading, leads to solution non-uniqueness. The paper then describes an alternative approach to ensure a unique solution in three-phase state estimation. This method aims to accurately and uniquely estimate the state of any unbalanced three-phase system, irrespective of load imbalance, network configuration, existence of synchronous generators or transformers.

Power System State Estimation, Three-Phase, Distri↗