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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗
NASA NTRS1983

Development of a Multilevel Optimization Approach to the Design of Modern Engineering Systems

A general algorithm is proposed which carries out the design process iteratively, starting at the top of the hierarchy and proceeding downward. Each subproblem is optimized separately for fixed controls from higher level subproblems. An optimum sensitivity analysis is then performed which determines the sensitivity of the subproblem design to changes in higher level subproblem controls. The resulting sensitivity derivatives are used to construct constraints which force the controlling subproblems into chosing their own designs so as to improve the lower levels subproblem designs while satisfying their own constraints. The applicability of the proposed algorithm is demonstrated by devising a four-level hierarchy to perform the simultaneous aerodynamic and structural design of a high-performance sailplane wing for maximum cross-country speed. Finally, the concepts discussed are applied to the two-level minimum weight structural design of the sailplane wing. The numerical experiments show that discontinuities in the sensitivity derivatives may delay convergence, but that the algorithm is robust enough to overcome these discontinuities and produce low-weight feasible designs, regardless of whether the optimization is started from the feasible space or the infeasible one.

Barthelemy, J. F. M.↗
NASA NTRS1998

Effects of Buoyancy on Lean Premixed V-Flames Part I: Laminar and Turblent Flame Structure

Laser schlieren and planar laser-induced fluorescence techniques have been used to investigate laminar and turbulent v-flames in +g, -g, and micro g under flow conditions that span the regimes of momentum domination (Ri < 0. 1) and buoyancy domination (Ri > 0.1). Overall flame features shown by schlieren indicate that buoyancy dominates the entire flow field for conditions close to Ri = 1. With decreasing Ri, buoyancy effects are observed only in the far-field regions. Analyses of the mean flame angles demonstrate that laminar and turbulent flames do not have similar responses to buoyancy. Difference in the laminar +g and -g flame angles decrease with Ri (i.e., increasing Re) and converge to the microgravity flame angle at the momentum limit (Ri - 0). This is consistent with the notion that the effects of buoyancy diminish with increasing flow momentum. The +g and -g turbulent flame angles, however, do not converge at Ri = 0. As shown by OH-PLIF images, the inconsistency in +g and -g turbulent flame angles is associated with the differences in flame wrinkles. Turbulent flame wrinkles evolve more slowly in +g than in -g. The difference in flame wrinkle structures, however, cannot be explained in terms of buoyancy effects on flame instability mechanisms. It seems to be associated with the field effects of buoyancy that stretches the turbulent flame brushes in +g and compresses the flame brush in -g. Flame wrinkling offers a mechanism through which the flame responds to the field effects of buoyancy despite increasing flow momentum. These observations point to the need to include both upstream and downstream contributions in theoretical analysis of flame turbulence interactions.

Cheng, Robert K.↗
NASA NTRS1973

Comparison of finite-difference schemes for analysis of shells of revolution

Several finite difference schemes are applied to the stress and free vibration analysis of homogeneous isotropic and layered orthotropic shells of revolution. The study is based on a form of the Sanders-Budiansky first-approximation linear shell theory modified such that the effects of shear deformation and rotary inertia are included. A Fourier approach is used in which all the shell stress resultants and displacements are expanded in a Fourier series in the circumferential direction, and the governing equations reduce to ordinary differential equations in the meridional direction. While primary attention is given to finite difference schemes used in conjunction with first order differential equation formulation, comparison is made with finite difference schemes used with other formulations. These finite difference discretization models are compared with respect to simplicity of application, convergence characteristics, and computational efficiency. Numerical studies are presented for the effects of variations in shell geometry and lamination parameters on the accuracy and convergence of the solutions obtained by the different finite difference schemes. On the basis of the present study it is shown that the mixed finite difference scheme based on the first order differential equation formulation and two interlacing grids for the different fundamental unknowns combines a number of advantages over other finite difference schemes previously reported in the literature.

Noor, A. K.↗
NASA NTRS1984

Data analysis and model simulation of the initiation of convection on April 24, 1982

A combined modelling and data analysis is used to examine the initiation and development of a hail-producing mesoscale convective system which developed over the Texas panhandle near Amarillo, Texas, on April 24, 1982, and propagated east into Oklahoma. This event occurred during the NASA AVE-VAS IV experiment, as a dense upper air mesoscale network was collecting data at three-hour intervals. The results show that both synoptic and mesoscale processes separately and in concert help to force the convection. The destabilization of the atmosphere due to cold advection is important in removing any cap to deep convective development. The upper level moistening due to flow from the northwest also helps to increase the convective instability of the atmosphere. The mesoscale convergence due to cloud shading appears dominant over that of topography, but the latter cannot be ignored, especially in its sustained effect. The developing mesoscale circulation also appears to tap a surface supply of moisture feeding the convection.

Mcnider, R. T.↗
NASA NTRS2015

Space Radiation Transport Code Development: 3DHZETRN

The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.

Wilson, John W.↗
NASA NTRS2021

Adaptive Aerostructures for Revolutionary Civil Supersonic Transportation

To enable commercially viable civil supersonic transport (SST) aircraft, innovative solutions must be developed to meet noise and efficiency requirements for overland flight. This research effort consists of a multidisciplinary team of academic and industrial experts exploring for the first time the potential of small real-time geometric outer mold line (OML) reconfigurations to minimize sonic boom signatures and aircraft drag in response to changing ambient conditions, thereby enabling noise-compliant overland supersonic flight. The team utilizes recent advances in supersonic computational fluid dynamic (CFD) methods, new noise prediction tools, and new design approaches to consider embedded highly energy-dense shape memory alloy (SMA) actuators for local shape modifications to an SST aircraft leading to optimal low boom signature and low drag in different environments. This university-led program will provide strategic leadership toward technology convergence that advances NASA's Aerospace Research Mission Directorate's (ARMD) research objectives with regard to Thrust 2: “Innovation in Commercial Supersonic Aircraft” by exploring for the first time enabling low-boom operation across a range of flight conditions via structural adaptivity, and will promote education of the next generation of engineers. The overall research strategy is to pursue three critical areas: the design of configurations for reducing boom, SMA material development and modeling, and technology feasibility demonstration in a relevant environment. Initially, the team will identify potential applications where structure or geometry adaptivity provides a benefit in noise or drag across the entire flight envelope. For selected applications/structural locations, required OML geometry changes will be determined based on analysis of sonic boom ground signature and drag reduction using new design tools, trade studies, and atmospheric sensing techniques. Designs will be developed and evaluated against requirements on loading, stroke length, and operational temperature. New alloy formulations will be developed tailored for both autonomous and controlled actuation modes. As the SMA material development matures, integrated system-level factors will be investigated. Optimized designs for small-scale distributed adaptivity applications of maximum benefit will then be matured and tested, moving toward demonstration of the innovative technology approaches at a TRL 4-5 and showing that sonic booms can be reduced by reconfiguration on demand.

Koushik Datta↗
DOE OSTI2021

National direct-drive program on OMEGA and the National Ignition Facility

A major advantage of the laser direct-drive (DD) approach to ignition is the increased fraction of laser drive energy coupled to the hot spot and relaxed hot-spot requirements for the peak pressure and convergence ratios relative to the indirect-drive approach at equivalent laser energy. With the goal of a successful ignition demonstration using DD, the recently established national strategy has several elements and involves multiple national and international institutions. These elements include the experimental demonstration on OMEGA cryogenic implosions of hot-spot conditions relevant for ignition at MJ-scale energies available at the National Ignition Facility (NIF) and developing an understanding of laser-plasma interactions and laser coupling using DD experiments on the NIF. DD designs require reaching central stagnation pressures in excess of 100 Gbar. The current experiments on OMEGA have achieved inferred peak pressures of 56 Gbar (Regan et al 2016 Phys. Rev. Lett. 117 025001). Extensive analysis of the cryogenic target experiments and two- and three-dimensional simulations suggest that power balance, target offset, and target quality are the main limiting factors in target performance. In addition, cross-beam energy transfer (CBET) has been identified as the main mechanism reducing laser coupling. Reaching the goal of demonstrating hydrodynamic equivalence on OMEGA includes improving laser power balance, target position, and target quality at shot time. CBET must also be significantly reduced and several strategies have been identified to address this issue. (2016-01-01)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local bilinear computation of Jacobi sets

Abstract We propose a novel method for the computation of Jacobi sets in 2D domains. The Jacobi set is a topological descriptor based on Morse theory that captures gradient alignments among multiple scalar fields, which is useful for multi-field visualization. Previous Jacobi set computations use piecewise linear approximations on triangulations that result in discretization artifacts like zig-zag patterns. In this paper, we utilize a local bilinear method to obtain a more precise approximation of Jacobi sets by preserving the topology and improving the geometry. Consequently, zig-zag patterns on edges are avoided, resulting in a smoother Jacobi set representation. Our experiments show a better convergence with increasing resolution compared to the piecewise linear method. We utilize this advantage with an efficient local subdivision scheme. Finally, our approach is evaluated qualitatively and quantitatively in comparison with previous methods for different mesh resolutions and across a number of synthetic and real-world examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Sea level variation

Published values for the long-term, global mean sea level rise determined from tide gauge records range from about one to three mm per year. The scatter of the estimates appears to arise largely from the use of data from gauges located at convergent tectonic plate boundaries where changes of land elevation give fictitious sea level trends, and the effects of large interdecadal and longer sea level variations on short (less than 50+ years) or sappy records. In addition, virtually all gauges undergo subsidence or uplift due to isostatic rebound from the last deglaciation at a rate comparable to or greater than the secular rise of sea level. Modeling rebound by the ICE-3G model of Tushingham and Peltier (1990) and avoiding tide gauge records in areas of converging tectonic plates produces a highly consistent set of long sea level records. A global set of 21 such stations in nine oceanic regions with an average record length of 76 years during the period 1880-1980 yields the global sea level rise value 1.8 mm/year +/- 0.1. Greenhouse warming scenarios commonly forecast an additional acceleration of global sea level in the next 5 or 6+ decades in the range 0.1-0.2 mm/yr2. Because of the large power at low frequencies in the sea level spectrum, very long tide gauge records (75 years minimum) have been examined for past apparent sea level acceleration. For the 80-year period 1905-1985, 23 essentially complete tide gauge records in 10 geographic groups are available for analysis. These yielded the apparent global acceleration -0.011 (+/- 0.012) mm/yr2. A larger, less uniform set of 37 records in the same 10 groups with 92 years average length covering the 141 years from 1850-1991 gave 0.001 (+/- 0.008) mm/yr2. Thus there is no evidence for an apparent acceleration in the past 100+ years that is significant either statistically, or in comparison to values associated with global warming. Unfortunately, the large interdecadal fluctuations of sea level severely affect estimates of global sea level acceleration for time spans of less than about 50 years. This means that tide gauges alone cannot serve as a reliable leading indicator of climate change in less than many decades. This time required can be significantly reduced if the interdecadal fluctuations of sea level can be understood in terms of their forcing mechanisms, and then removed from the tide gauge records.

Douglas, Bruce C.↗
NASA NTRS2020

Robust Conceptual Design of Transonic Airfoils

This paper describes an integrated, multi-fidelity analysis and heuristic design approach that can be used to derive initial airfoil designs for transonic flight. If successful, the final result is a geometry that can be expected to produce reasonable aerodynamic performance when used with higher order analysis methods. A key aspect of the methodology is the use of a “sonic-plateau” pressure distribution as the target distribution for inverse design. The sonic- plateau distribution is easily parameterized and has the advantage of automatically resulting in a smooth airfoil shape without any discontinuities built into the surface due to the presence of a shock in the target pressure distribution. Inverse design is performed on each airfoil using a parametrically defined pressure distribution at a reduced lift coefficient and Mach number from the design condition. The methodology is demonstrated by designing an airfoil at 38% of the wing semispan for a 737-200-like aircraft. The demonstration problem shows that the methodology is able to achieve rapid and robust convergence to the solution. The calculated designed airfoil was found to be sufficiently higher than Mach number, and the maximum thickness was close to the targeted value.

Olson, Erik D.↗
DOE OSTI2022

Toward designing effective exascale scientific computing workflows: experiences and best practices

Many fields within scientific computing have embraced advances in big-data analysis and machine learning, which often requires the deployment of large, distributed and complicated workflows that may combine training neural networks, performing simulations, running inference, and performing database queries and data analysis in asynchronous, parallel and pipelined execution frameworks. Such a shift has brought into focus the need for scalable, efficient workflow management solutions with reproducibility, error and provenance handling, traceability, and checkpoint-restart capabilities, among other needs. Here, we discuss challenges and best-practices for deploying exascale-generation computational science workflows on resources at the Oak Ridge Leadership Computing Facility (OLCF). We present our experiences with large-scale deployment of distributed workflows on the Summit supercomputer, including for bioinformatics and computational biophysics, materials science, and deep learning model optimization. We also present problems and solutions created by working within a Python-centric software base on traditional HPC systems, and discuss steps that will be required before the convergence of HPC, AI, and data science can be fully realized. Our results point to a wealth of exciting new possibilities for harnessing this convergence to tackle new scientific challenges.

Coletti, Mark↗
DOE OSTI2021

MADLens, a python package for fast and differentiable non-Gaussian lensing simulations

Here, we present MADLens a python package for producing non-Gaussian lensing convergence maps at arbitrary source redshifts with unprecedented precision. MADLens is designed to achieve high accuracy while keeping computational costs as low as possible. A MADLens simulation with only particles produces convergence maps whose power agrees with theoretical lensing power spectra up to within the accuracy limits of HaloFit. This is made possible by a combination of a highly parallelizable particle-mesh algorithm, a sub-evolution scheme in the lensing projection, and a machine-learning inspired sharpening step. Further, MADLens is fully differentiable with respect to the initial conditions of the underlying particle-mesh simulations and a number of cosmological parameters. These properties allow MADLens to be used as a forward model in Bayesian inference algorithms that require optimization or derivative-aided sampling. Another use case for MADLens is the production of large, high resolution simulation sets as they are required for training novel deep-learning-based lensing analysis tools. We make the MADLens package publicly available under a Creative Commons License

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Comparative genomics provides insights into the cold adaptation of endophytic fungi associated with Deschampsia antarctica

Endophytic fungi from Deschampsia antarctica , the southernmost flowering plant, provide insights into the cold adaptation mechanisms of plant-associated fungi in extreme environments. This study presents the genome sequences and comparative analysis of eight fungal isolates from D. antarctica leaves. These Antarctic fungal isolates were analyzed alongside 121 plant-associated fungal genomes to uncover signatures of adaptation and endophytic specialization. Antarctic endophytes show striking patterns, including reduced genome size (∼26.3 Mb on average), streamlined gene content (∼8844 genes), and notably small secretomes (∼288 proteins). Despite this reduced gene repertoire, they maintain a robust set of genes encoding carbohydrate-active enzymes (CAZymes) but lack those for lignin and bacterial cell wall degradation, indicating a symbiotic lifestyle that avoids host damage and predation. One isolate, Alternaria sp. UNIPAMPA017 stood out, with 26% of its genome occupied by transposable elements. Lifestyle, rather than phylogeny, was the main driver of CAZyme and secretome profiles, underscoring ecological convergence. Compared to endophytes from Arabidopsis and Populus, D. antarctica endophytes harbor fewer pectin-degrading enzymes, reflecting their adaptation to the cell wall structure of their monocot host. Together, these fungi reveal a pattern of genomic reduction and functional fine-tuning, hallmarks of life adapted to persist in cold, nutrient-scarce niches.

Ascomycota↗
NASA NTRS1982

Doppler-radar observation of the evolution of downdrafts in convective clouds

A detailed analysis of the 20 July 1977 thunderstorm complex which formed and evolve over the South Park region in Central Colorado is presented. The storm was extensively analyzed using multiple Doppler radar and surface mesonet data, developed within an environment having very weak wind shear. The storm owed its intensification to the strength of the downdraft, which was nearly coincident with the region where the cloud had grown. The noteworthy features of this storm were its motion to the right of the cloud-level winds, its multicellular nature and discrete propagation, its north-south orientation, and its relatively large storm size and high reflectivity factor (55 dBZ). This scenario accounts for the observed mesoscale and cloud-scale event. A line of convergence was generated at the interface between the easterly upslope winds and westerly winds. During stage II, the convergence line subsequently propagated down the slopes of the Mosquito Range, and was the main forcing mechanism for the development of updraft on the west flank of the storm. The formation of downdraft on the eastern side of updraft blacked surface inflow, and created a detectable gust front. As the original downdraft intensified, the accumulation of evaporatively-chilled air caused the intensification of the mesohigh, which likely destroyed the earlier convergence line and created a stronger convergence line to the east, which forced up-lifting of the moist, westerly inflow and caused the formation of updraft to the east. An organized downdraft circulation, apparently maintained by precipitation drag and evaporational cooling, was responsible in sustaining a well-defined gust front. The storm attained its highest intensity as a consequence of merging with a neighboring cloud. The interaction of downdrafts or gust fronts from two intense cells appeared to be the primary mechanism of this merging process as suggested by Simpson et al. (1980). The merging process coincided with more rain than occurred in unmerged echoes.

Motallebi, N.↗
NASA NTRS1989

Development of parallel algorithms for electrical power management in space applications

The application of parallel techniques for electrical power system analysis is discussed. The Newton-Raphson method of load flow analysis was used along with the decomposition-coordination technique to perform load flow analysis. The decomposition-coordination technique enables tasks to be performed in parallel by partitioning the electrical power system into independent local problems. Each independent local problem represents a portion of the total electrical power system on which a loan flow analysis can be performed. The load flow analysis is performed on these partitioned elements by using the Newton-Raphson load flow method. These independent local problems will produce results for voltage and power which can then be passed to the coordinator portion of the solution procedure. The coordinator problem uses the results of the local problems to determine if any correction is needed on the local problems. The coordinator problem is also solved by an iterative method much like the local problem. The iterative method for the coordination problem will also be the Newton-Raphson method. Therefore, each iteration at the coordination level will result in new values for the local problems. The local problems will have to be solved again along with the coordinator problem until some convergence conditions are met.

Berry, Frederick C.↗
NASA NTRS1999

Design of a Low Speed Fan Stage for Noise Suppression

This report describes the design of a low tip speed, moderate pressure rise fan stage for demonstration of noise reduction concepts. The fan rotor is a fixed-pitch configuration delivering a design pressure ratio of 1.378 at a specific flow of 43.1 lbm/sec/sq ft. Four exit stator configurations were provided to demonstrate the effectiveness of circumferential and axial sweep in reducing rotor-stator interaction tone noise. The fan stage design was combined with an axisymmetric inlet, conical convergent nozzle, and nacelle to form a powered fan-nacelle subscale model. This model has a 22-inch cylindrical flow path and employs a rotor with a 0.30 hub-to-tip radius ratio. The design is fully compatible with an existing NASA force balance and rig drive system. The stage aerodynamic and structural design is described in detail. Three-dimensional (3-D) computational fluid dynamics (CFD) tools were used to define optimum airfoil sections for both the rotor and stators. A fan noise predictive system developed by Pratt & Whitney under contract to NASA was used to determine the acoustic characteristics of the various stator configurations. Parameters varied included rotor-to-stator spacing and vane leading edge sweep. The structural analysis of the rotor and stator are described herein. An integral blade and disk configuration was selected for the rotor. Analysis confirmed adequate low cycle fatigue life, vibratory endurance strength, and aeroelastic suitability. A unique load carrying stator arrangement was selected to minimize generation of tonal noise due to sources other than rotor-stator interaction. Analysis of all static structural components demonstrated adequate strength, fatigue life, and vibratory characteristics.

Dalton, W. N.↗
DOE OSTI2022

An Evaluation of Tropical Cyclone Rainfall Structures in the HighResMIP Simulations against Satellite Observations

This study evaluates tropical cyclone (TC) rainfall structures in the CMIP6 HighResMIP global climate model (GCM) simulations against satellite rainfall retrievals. Here, we specifically focus on TCs within the deep tropics (25°S–25°N). Analysis of TC rain rate composites indicates that in comparison to the satellite observations at the same intensity, many HighResMIP simulations tend to overproduce rain rates around TCs, in terms of both maximum rain rate magnitude and area-averaged rain rates. In addition, as model horizontal resolution increases, the magnitude of the peak rain rate appears to increase. However, the area-averaged rain rates decrease with increasing horizontal resolution, partly due to the TC eyewall being located closer to the TC center, thus occupying a smaller area and contributing less to the area-averaged rain rates. The effect of ocean coupling is to lower the TC rain rates, bringing them closer to the satellite observations, due to reduced horizontal moisture flux convergence and surface latent heat flux beneath TCs. Examination of horizontal rain rate distributions indicates that vertical wind shear–induced rainfall asymmetries in HighResMIP-simulated TCs are qualitatively consistent with the observations. In addition, a positive relationship is observed between the area-averaged inner-core rainfall and TC intensification likelihoods across the HighResMIP simulations, as GCM simulations producing stronger TCs more frequently have the greater rainfall close to the center, in agreement with previous theoretical and GCM simulation results.

54 ENVIRONMENTAL SCIENCES↗