SEARCH · Engineering Papers

Results for “Convergence analysis”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Piston driven converging shock waves in a stiffened gas

We present that the problem of a one-dimensional (1D) cylindrically or spherically symmetric shock wave converging into an inviscid, ideal gas was first investigated by Guderley[Starke kugelige und zylinrische verdichtungsstosse in der nahe des kugelmitterpunktes bzw. Der zylinderachse,” Luftfahrtforschung 19, 302 (1942)]. In the time since, many authors have discussed the practical notion of how Guderley-like flows might be generated. One candidate is a constant velocity, converging “cylindrical or spherical piston,” giving rise to a converging shock wave in the spirit of its classical, planar counterpart. A limitation of pre-existing analyses along these lines is the restriction to flows in materials described by an ideal gas equation of state (EOS) constitutive law. This choice is of course necessary for the direct comparison with the classical Guderley solution, which also features an ideal gas EOS. However, the ideal gas EOS is limited in its utility in describing a wide variety of physical phenomena and, in particular, the shock compression of solid materials. This work is thus intended to provide an extension of previous work to a nonideal EOS. The stiff gas EOS is chosen as a logical starting point due to not only its close resemblance to the ideal gas law but also its relevance to the shock compression of various liquid and solid materials. Using this choice of EOS, the solution of a 1D planar piston problem is constructed and subsequently used as the lowest order term in a quasi-self-similar series expansion intended to capture both curvilinear and nonideal EOS effects. The solution associated with this procedure provides correction terms to the 1D planar solution so that the expected accelerating shock trajectory and nontrivially varying state variable profiles can be obtained. This solution is further examined in the limit as the converging shock wave approaches the 1D curvilinear origin. Lastly, given the stiff gas EOS is not otherwise expected to admit a Guderley-like solution when coupled to the inviscid Euler equations, this work thus provides the semianalytical limiting behavior of a flow that cannot be otherwise captured using self-similar analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Stochastic Adaptive Droop Control in Frequency Regulation of Power Systems With Intermittent Generators

Modern power systems (MPSs), including microgrids (MGs), are increasingly incorporating multiple renewable energy sources (RESs) such as wind and solar power, as well as battery storage and controllable loads. While environmentally beneficial, these sources pose challenges for control and management due to their intermittent and stochastic nature, especially in maintaining frequency stability with multiple interconnected generators of varying capacities. Traditional droop control methods are effective in systems with generators that are dispatchable and have fixed generation capacities, but they fall short when applied to systems with RESs, where generation capacities are dynamic and affected by unpredictable environmental conditions. To address these challenges, this paper introduces a novel stochastic adaptive droop control (SADC) method for load frequency control (LFC). The proposed method adapts droop coefficients in real time, based on the measured stochastic data of power generation capacities, enabling more effective frequency regulation in systems with variable and intermittent power generation. Unlike traditional adaptive control methods, which assume constant or slowly-varying system parameters, this approach accounts for stochastic processes by modeling them as Markov chains, enabling robust performance under highly dynamic and unpredictable conditions. The key contributions of this work include the development of real-time droop coefficient adaptation algorithms, derivation of their stability and convergence properties, and the demonstration of the advantages of the method through simulations. Case studies highlight the improved performance of frequency regulation, particularly in addressing the impact of stochastic weather conditions and the benefits of reducing dependence on battery reserves in dealing with intermittency of RESs. Finally, this paper provides a comprehensive analysis of the theoretical foundations of the method, as well as practical implementation insights for future power systems with high penetration of RESs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1973

Stress, stability, and vibration of complex, branched shells of revolution.

A comprehensive computer program, designated BOSOR4, for analysis of the stress, stability, and vibration of segmented, ring-stiffened, branched shells of revolution and prismatic shells and panels is described. The program performs large-deflection axisymmetric stress analysis, small-deflection nonsymmetric stress analysis, modal vibration analysis with axisymmetric nonlinear prestress included, and buckling analysis with axisymmetric or nonsymmetric prestress. One of the main advantages of the code is the provision for realistic engineering details such as eccentric load paths, internal supports, arbitrary branching conditions, and a library of wall constructions. The program is based on the finite-difference energy method, which is very rapidly convergent with increasing numbers of mesh points.

Bushnell, D.↗

NASA NTRS1976

A study of efficiency in low resistivity silicon solar cells

A general device-analysis program has been utilized to study the efficiency of silicon solar cells. The analysis is applied to specific geometries of both n(+)-p and n(+)-p-p(+) solar cells, and involves a numerical solution of the basic transport and continuity equations. This approach allows solutions that are free of typical limiting assumptions involved in solving the device equations as well as solutions relating to lifetime, mobility variations, and diffused-region profiles. The analysis includes available empirical information on diffusion length, mobility, and lifetime as a function of doping, as well as a Gaussian profile for the diffused region. Results are presented which illustrate the limitations of efficiency as a function of doping. It is found that the maximum efficiencies for both types of cell converge at lower resistivities to around 16% with air-mass-zero radiation and a single-layer absorbing-SiO antireflecting film. It is also found that the minority-carrier lifetime, both in the n(+) surface and p-type bulk regions, presents serious limitations to conversion efficiency, particularly in the low-resistivity cells.

Dunbar, P. M.↗

NASA NTRS1986

A new approximate LU factorization scheme for the Reynolds-averaged Navier-Stokes equations

A new approximate LU factorization scheme is developed to solve the steady state Reynolds-averaged Navier-Stokes (NS) equations. Central differencing is used for both implicit and explicit operator and special care is taken to obtain well-conditioned factors on the implicit side. The scheme is then analyzed and optimized according to a simple linear analysis. It is unconditionally stable for the model hyperbolic equation in both two- and three-dimensions. However, the requirement for well-conditioned factors has essentially limited the effective time step the scheme can achieve. Supersonic and transonic three-dimensinal flows past a hemisphere cylinder are computed to demonstrate the convergence characteristics of the scheme. A good convergence rate is achieved for the inviscid case. Finally, an explicit eigenvector annihilation procedure is adopted successfully to remove the stiffness caused by the fine grid spacing for viscous flows.

Pan, D.↗

NASA NTRS1991

Two-dimensional Navier-Stokes heat transfer analysis for rough turbine blades

A quasi-three-dimensional thin-layer Navier-Stokes analysis was used to predict heat transfer to rough surfaces. Comparisons are made between predicted and experimental heat transfer for turbine blades and flat plates of known roughness. The effect of surface roughness on heat transfer was modeled using a mixing length approach. The effect of near-wall grid spacing and convergence criteria on the accuracy of the heat transfer predictions are examined. An eddy viscosity mixing length model having an inner and outer layer was used. A discussion of the appropriate model for the crossover between the inner and outer layers is included. The analytic results are compared with experimental data for both flat plates and turbine blade geometries. Comparisons between predicted and experimental heat transfer showed that a modeling roughness effects using a modified mixing length approach results in good predictions of the trends in heat transfer due to roughness.

Boyle, R. J.↗

NASA NTRS1991

Two-dimensional Navier-Stokes heat transfer analysis for rough turbine blades

Boyle, R. J.↗

NASA NTRS1996

The Quest for Pi

This article gives a brief history of the analysis and computation of the mathematical constant Pi=3.14159 ..., including a number of the formulas that have been used to compute Pi through the ages. Recent developments in this area are then discussed in some detail, including the recent computation of Pi to over six billion decimal digits using high-order convergent algorithms, and a newly discovered scheme that permits arbitrary individual hexadecimal digits of Pi to be computed.

Bailey, David H.↗

NASA NTRS2002

Eigenvalues of Rectangular Waveguide Using FEM With Hybrid Elements

A finite element analysis using hybrid triangular-rectangular elements is developed to estimate eigenvalues of a rectangular waveguide. Use of rectangular vector-edge finite elements in the vicinity of the PEC boundary and triangular elements in the interior region more accurately models the physical nature of the electromagnetic field, and consequently quicken the convergence.

Deshpande, Manohar D.↗

NASA NTRS1997

An Improved Treatment of External Boundary for Three-Dimensional Flow Computations

We present an innovative numerical approach for setting highly accurate nonlocal boundary conditions at the external computational boundaries when calculating three-dimensional compressible viscous flows over finite bodies. The approach is based on application of the difference potentials method by V. S. Ryaben'kii and extends our previous technique developed for the two-dimensional case. The new boundary conditions methodology has been successfully combined with the NASA-developed code TLNS3D and used for the analysis of wing-shaped configurations in subsonic and transonic flow regimes. As demonstrated by the computational experiments, the improved external boundary conditions allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable speedup of convergence of the multigrid iterations.

Tsynkov, Semyon V.↗

NASA NTRS2004

Magnetic-Nozzle Studies for Fusion Propulsion Applications: Gigawatt Plasma Source Operation and Magnetic Nozzle Analysis

This project has been a multiyear effort to assess the feasibility of a key process inherent to virtually all fusion propulsion concepts: the expansion of a fusion-grade plasma through a diverging magnetic field. Current fusion energy research touches on this process only indirectly through studies of plasma divertors designed to remove the fusion products from a reactor. This project was aimed at directly addressing propulsion system issues, without the expense of constructing a fusion reactor. Instead, the program designed, constructed, and operated a facility suitable for simulating fusion reactor grade edge plasmas, and to examine their expansion in an expanding magnetic nozzle. The approach was to create and accelerate a dense (up to l0(exp 20)/m) plasma, stagnate it in a converging magnetic field to convert kinetic energy to thermal energy, and examine the subsequent expansion of the hot (100's eV) plasma in a subsequent magnetic nozzle. Throughout the project, there has been a parallel effort between theoretical and numerical design and modelling of the experiment and the experiment itself. In particular, the MACH2 code was used to design and predict the performance of the magnetoplasmadynamic (MPD) plasma accelerator, and to design and predict the design and expected behavior for the magnetic field coils that could be added later. Progress to date includes the theoretical accelerator design and construction, development of the power and vacuum systems to accommodate the powers and mass flow rates of interest to out research, operation of the accelerator and comparison to theoretical predictions, and computational analysis of future magnetic field coils and the expected performance of an integrated source-nozzle experiment.

Gilland, James H.↗

DOE OSTI2015

Pathways and transformations of dissolved methane and dissolved inorganic carbon in Arctic tundra watersheds: Evidence from analysis of stable isotopes

Arctic soils contain a large pool of terrestrial C and are of interest due to their potential for releasing significant carbon dioxide (CO 2 ) and methane (CH 4 ) to the atmosphere. Due to substantial landscape heterogeneity, predicting ecosystem-scale CH 4 and CO 2 production is challenging. This study assessed dissolved inorganic carbon (DIC = Σ (total) dissolved CO 2 ) and CH 4 in watershed drainages in Barrow, Alaska as critical convergent zones of regional geochemistry, substrates, and nutrients. In July and September of 2013, surface waters and saturated subsurface pore waters were collected from 17 drainages. Based on simultaneous DIC and CH 4 cycling, we synthesized isotopic and geochemical methods to develop a subsurface CH 4 and DIC balance by estimating mechanisms of CH 4 and DIC production and transport pathways and oxidation of subsurface CH 4 . We observed a shift from acetoclastic (July) toward hydrogenotropic (September) methanogenesis at sites located toward the end of major freshwater drainages, adjacent to salty estuarine waters, suggesting an interesting landscape-scale effect on CH 4 production mechanism. The majority of subsurface CH 4 was transported upward by plant-mediated transport and ebullition, predominantly bypassing the potential for CH 4 oxidation. Thus, surprisingly, CH 4 oxidation only consumed approximately 2.51± 0.82% (July) and 0.79 ± 0.79% (September) of CH 4 produced at the frost table, contributing to <0.1% of DIC production. DIC was primarily produced from respiration, with iron and organic matter serving as likely e- acceptors. Furthermore, this work highlights the importance of spatial and temporal variability of CH 4 production at the watershed scale and suggests broad scale investigations are required to build better regional or pan-Arctic representations of CH 4 and CO 2 production.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Conforming window functions for meshfree methods

Window functions provide a base for the construction of approximation functions in many meshfree methods. They control the smoothness and extent of the approximation functions and are commonly defined using Euclidean distances which helps eliminate the need for a meshed discretization, simplifying model development for some classes of problems. However, for problems with complicated geometries such as nonconvex or multi-body domains, poor solution accuracy and convergence can occur unless the extents of the window functions, and thus approximation functions, are carefully controlled, often a time consuming or intractable task. Here in this paper, we present a method to provide more control in window function design, allowing efficient and systematic handling of complex geometries. “Conforming” window functions are constructed using Bernstein–Bézier splines defined on local triangulations with constraints imposed to control smoothness. Graph distances are used in conjunction with Euclidean metrics to provide adequate information for shaping the window functions. The conforming window functions are demonstrated using the Reproducing Kernel Particle Method showing improved accuracy and convergence rates for problems with challenging geometries. Conforming window functions are also demonstrated as a means to simplify the imposition of essential boundary conditions.

, rapid design-to-analysis↗

NASA NTRS1971

Analysis of transient, linear wave propagation in shells by the finite difference method

The applicability of the finite difference method to propagation problems in shells, and the response of a cylindrical shell with cutouts to both longitudinal and radial transient excitations are investigated. It is found that the only inherent limitation of the finite difference method is its inability to reproduce accurately response discontinuities. The short wave length limitations of thin shell theory create significant convergence difficulties may often be overcome through proper selection of finite difference mesh dimensions and temporal or spatial smoothing of the excitation. Cutouts produce moderate changes in early and intermediate time response of a cylindrical shell to axisymmetric pulse loads applied at one end. The cutouts may facilitate the undesirable late-time transfer of load-injected extensional energy into nonaxisymmetric flexural response.

Geers, T. L.↗

NASA NTRS2004

Rayleigh-Sommerfield Diffraction vs Fresnel-Kirchhoff, Fourier Propagation and Poisson's Spot

The boundary conditions imposed on the diffraction problem in order to obtain the Fresnel-Kirchhoff (FK) solution are well-known to be mathematically inconsistent and to be violated by the solution when the observation point is close to the diffracting screen 1-3. These problems are absent in the Rayleigh-Sommerfeld (RS) solution. The difference between RS and FK is in the inclination factor and is usually immaterial because the inclination factor is approximated by unity. But when this approximation is not valid, FK can lead to unacceptable answers. Calculating the on-axis intensity of Poisson s spot provides a critical test, a test passed by RS and failed by FK. FK fails because (a) convergence of the integral depends on how it is evaluated and (b) when the convergence problem is xed, the predicted amplitude at points near the obscuring disk is not consistent with the assumed boundary conditions.

RS (RAYLEIGH SOMMERFELD)↗

NASA NTRSDate not supplied

ANALYSIS OF THE MSL/MEDLI ENTRY DATA WITH COUPLED CFD AND MATERIAL RESPONSE.

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing 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). 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 Open-FOAM (PATO) 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, 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.

Mars Science Laboratory↗

DOE OSTI2025

Effects of threshold resummation for large-x PDF in large momentum effective theory

Parton distribution functions (PDFs) at large x are challenging to extract from experimental data, yet they are essential for understanding hadron structure and searching for new physics beyond the Standard Model. Within the framework of the large momentum P z expansion of lattice quasi-PDFs, we investigate large x PDFs, where the matching coefficient is factorized into the hard kernel, related to the active quark momentum xP z , and the threshold soft function, associated with the spectator momentum (1 − x)P z . The renormalization group equation of the soft function enables the resummation of the threshold double logarithms α k ln 2k (1 − x), which is crucial for a reliable and controllable calculation of large x PDFs. Our analysis with pion valence PDFs indicates that perturbative matching breaks down when the spectator momentum (1 − x)P z approaches Λ QCD , but remains valid when both xP z and (1 − x)P z are much larger than Λ QCD . Additionally, we incorporate leading renormalon resummation within the threshold framework, demonstrating good perturbative convergence in the region where both spectator and active quark momenta are perturbative scales.

factorization↗

DOE OSTI2018

Development of the Theory of Branched Continued Fractions in 1996–2016

We perform the analysis of investigations in the theory of branched continued fractions carried out for the last 20 years in the directions developed under the general supervision of Prof. V. Ya. Skorobohat’ko (18.07.1927–04.07.1996) including, in particular, the interpolation of functions of several variables by branched continued fractions, the determination of efficient criteria of convergence and computational stability of these fractions, the correspondence between multiple power series and functional branched continued fractions, the investigation of various classes of functional fractions, and the application of branched continued fractions.

97 MATHEMATICS AND COMPUTING↗