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,459 records · Page 81

MR Cygni revisited

New analysis tools and additional unanalyzed observations justify a reanalysis of MR Cygni. The reanalysis applied successively more restrictive physical models, each with an optimization program. The final model assigned separate first and second order limb darkening coefficients, from model atmospheres, to individual grid points. Proper operation of the optimization procedure was tested on simulated observational data, produced by light synthesis with assigned system parameters, and modulated by simulated observational error. The iterative solution converged to a weakly-determined mass ratio of 0.75. Assuming the B3 primary component is on the main sequence, the HR diagram location of the secondary from the light ratio (ordinate) and adjusted T sub eff (abscissa) was calculated. The derived mass ratio, together with a main-sequence mass for the B3 component, implies a main-sequence secondary spectral type of B4. The photometrically-determined secondary radii agree with this spectral type, in marginal disagreement with the B7 type from the HR diagram analysis. The individual masses, derived from the radial velocity curve of the primary component, the photometrically-determined i, and alternative values of derived mass ratio are seriously discrepant with main sequence objects. The imputed physical status of the system is in disagreement with representations that have appeared in the literature.

Linnell, Albert P.↗

Difference-Equation/Flow-Graph Circuit Analysis

Numerical technique enables rapid, approximate analyses of electronic circuits containing linear and nonlinear elements. Practiced in variety of computer languages on large and small computers; for circuits simple enough, programmable hand calculators used. Although some combinations of circuit elements make numerical solutions diverge, enables quick identification of divergence and correction of circuit models to make solutions converge.

Mcvey, I. M.↗

A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds

Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

97 MATHEMATICS AND COMPUTING↗

Phase-field simulations of gas density within bubbles under irradiation

Phase-field simulations are used to study the evolution of gas density within irradiation-induced bubbles. In our simulations, the dpa rate, gas production rate, and defect diffusivities are systematically varied to understand their effect on bubble nucleation rates, bubble densities, and the distribution of gas concentration within bubbles and in the solid regions. We find that gas densities within bubbles fluctuate drastically in the early nucleation stages, when growth rates are highest, but converge to steady-state values during the later coarsening stages. The steady-state gas densities within bubbles correspond with the ratio of total accumulated vacancy content divided by the total accumulated gas content, in agreement with a thermodynamic analysis concerning free-energy minimization.

36 MATERIALS SCIENCE↗

High-accuracy emulators for observables in ΛCDM, N eff, Σ m ν, and w cosmologies

ABSTRACT We use the emulation framework CosmoPower to construct and publicly release neural network emulators of cosmological observables, including the cosmic microwave background (CMB) temperature and polarization power spectra, matter power spectrum, distance-redshift relation, baryon acoustic oscillation (BAO) and redshift-space distortion (RSD) observables, and derived parameters. We train our emulators on Einstein–Boltzmann calculations obtained with high-precision numerical convergence settings, for a wide range of cosmological models including ΛCDM, wCDM, ΛCDM + Neff, and ΛCDM + Σmν. Our CMB emulators are accurate to better than 0.5 per cent out to ℓ = 104, which is sufficient for Stage-IV data analysis, and our P(k) emulators reach the same accuracy level out to $k=50 \, \, \mathrm{Mpc}^{-1}$, which is sufficient for Stage-III data analysis. We release the emulators via an online repository (CosmoPower Organisation), which will be continually updated with additional extended cosmological models. Our emulators accelerate cosmological data analysis by orders of magnitude, enabling cosmological parameter extraction analyses, using current survey data, to be performed on a laptop. We validate our emulators by comparing them to class and camb and by reproducing cosmological parameter constraints derived from Planck TT, TE, EE, and CMB lensing data, as well as from the Atacama Cosmology Telescope Data Release 4 CMB data, Dark Energy Survey Year-1 galaxy lensing and clustering data, and Baryon Oscillation Spectroscopic Survey Data Release 12 BAO and RSD data.

Astronomy & Astrophysics↗

Generation of anisotropy in turbulent flows subjected to rapid distortion

A computational tool for the anisotropic time-evolution of the spectral velocity correlation tensor is presented. We operate in the linear, rapid distortion limit of the mean-field-coupled equations. Each term of the equations is written in the form of an expansion to arbitrary order in the basis of irreducible representations of the SO(3) symmetry group. The computational algorithm for this calculation solves a system of coupled equations for the scalar weights of each generated anisotropic mode. The analysis demonstrates that rapid distortion rapidly but systematically generates higher-order anisotropic modes. To maintain a tractable computation, the maximum number of rotational modes to be used in a given calculation is specified a priori . The computed Reynolds stress converges to the theoretical result derived by Batchelor and Proudman [Quart. J. Mech. Appl. Math. 7, 83 (1954)] if a sufficiently large maximum number of rotational modes is utilized; more modes are required to recover the solution at later times. In conclusion, the emergence and evolution of the underlying multidimensional space of functions is presented here using a 64-mode calculation. Alternative implications for modeling strategies are discussed.

42 ENGINEERING↗

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculations of the pressure distribution on axisymmetric boattails including effects of viscous interactions and exhaust jets in subsonic flow

A method of calculating the pressure distributions on boattails is proposed. This method accounts for viscous effects including the presence of a separated region for base flows by combining an inviscid analysis with a boundary layer analysis in an iterative calculation. Details of the reversed flow region are not considered. Some preliminary results have been obtained for boattails at subsonic free stream Mach number with turbulent boundary layers separating at the boattail base. In some cases convergence could not be obtained using the present computer program. It is possible, in principle, to extend this method to the calculation of boattail flows with pressure gradient induced separation on the boattail.

Rom, J.↗

Transient airload computer analysis for simulating wind induced impulsive noise conditions of a hovering helicopter rotor

A numerical analysis was developed to determine the airloads on helicopter rotors operating under near-hovering flight conditions capable of producing impulsive noise. A computer program was written in which the solutions for the rotor tip vortex geometry, inflow, aeroelastic response, and airloads are solved in a coupled manner at sequential time steps, with or without the influence of an imposed steady ambient wind or transient gust. The program was developed for future applications in which predicted airloads would be incorporated in an acoustics analysis to attempt to predict and analyze impulsive noise (blade slap). The analysis was applied to a hovering full-scale rotor for which impulsive noise was recorded in the presence of ambient wind. The predicted tip vortex coordinates are in reasonable agreement with the test data, and the blade airload solutions converged to a periodic behavior for an imposed steady ambient wind conditions.

Hall, G. F.↗

A theoretical and experimental investigation of shock-associated noise in supersonic jets

This paper examines the fundamental mechanism of broadband shock noise generation in an improperly expanded supersonic jet. The study includes circular convergent-divergent as well as circular convergent nozzles. The main source of shock noise is determined to be the transient interaction between the shock front and the convected vorticity within the jet plume. The discussion of the noise-generation mechanism is based on detailed numerical analysis, theoretical aerodynamic and acoustic modeling, refined measurements of the jet mean flow, shock-cell structure, turbulence, and noise. These results provide a broad-based generalization for the Harper-Bourne and Fisher analysis and prediction method.

Pao, S. P.↗

Numerical solution of transonic wing flow fields

A fast, fully implicit, approximate factorization algorithm designed to solve the conservative full-potential equation is used to compute lifting-wing flow fields with embedded transonic flow. The computer program (TWING - transonic wing analysis) uses an elliptic-solver numerical grid-generation routine and is capable of analyzing arbitrary wing geometries with sweep, twist, taper, and section variation. Computed flowfield results for several wing geometries are presented; they demonstrate substantial improvement in convergence speed relative to the FL028 computer code.

Holst, T. L.↗

Finite element analysis of wrinkling membranes

The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.

Miller, R. K.↗

Finite element analysis of partly wrinkled membranes

The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.

Miller, R. K.↗

High degree interpolation polynomial in Newton form

Polynomial interpolation is an essential subject in numerical analysis. Dealing with a real interval, it is well known that even if f(x) is an analytic function, interpolating at equally spaced points can diverge. On the other hand, interpolating at the zeroes of the corresponding Chebyshev polynomial will converge. Using the Newton formula, this result of convergence is true only on the theoretical level. It is shown that the algorithm which computes the divided differences is numerically stable only if: (1) the interpolating points are arranged in a different order, and (2) the size of the interval is 4.

Tal-Ezer, Hillel↗

Computing Stress, Stability, And Vibration Of Shells

BOSOR4 computer program developed as comprehensive program for analysis of stress, stability, and vibration of complex, branched shells of revolution made of elastic materials. Used to analyze prismatic shells and panels. 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 main advantage, provisions for such realistic engineering details as eccentric load paths, internal supports, arbitrary branching conditions, and "library" of wall constructions. Based on finite-difference energy method and offers very rapid convergence with increasing numbers of mesh points. Written in FORTRAN 77.

Bushnell, David↗

Research highlights: June 1990 - May 1991

Linear instability calculations at MSFC have suggested that the Geophysical Fluid Flow Cell (GFFC) should exhibit classic baroclinic instability at accessible parameter settings. Interest was in the mechanisms of transition to temporal chaos and the evolution of spatio-temporal chaos. In order to understand more about such transitions, high resolution numerical experiments for the physically simplest model of two layer baroclinic instability were conducted. This model has the advantage that the numerical code is exponentially convergent and can be efficiently run for very long times, enabling the study of chaotic attractors without the often devastating effects of low-order trunction found in many previous studies. Numerical algorithms for implementing an empirical orthogonal function (EOF) analysis of the high resolution numerical results were completed. Under conditions of rapid rotation and relatively low differential heating, convection in a spherical shell takes place as columnar banana cells wrapped around the annular gap, but with axes oriented along the axis of rotation; these were clearly evident in the GFFC experiments. The results of recent numerical simulations of columnar convection and future research plans are presented.

Source record↗

Advanced Subsonic Technology (AST) Area of Interest (AOI) 6: Develop and Validate Aeroelastic Codes for Turbomachinery

AlliedSignal Engines, in cooperation with NASA GRC (National Aeronautics and Space Administration Glenn Research Center), completed an evaluation of recently-developed aeroelastic computer codes using test cases from the AlliedSignal Engines fan blisk and turbine databases. Test data included strain gage, performance, and steady-state pressure information obtained for conditions where synchronous or flutter vibratory conditions were found to occur. Aeroelastic codes evaluated included quasi 3-D UNSFLO (MIT Developed/AE Modified, Quasi 3-D Aeroelastic Computer Code), 2-D FREPS (NASA-Developed Forced Response Prediction System Aeroelastic Computer Code), and 3-D TURBO-AE (NASA/Mississippi State University Developed 3-D Aeroelastic Computer Code). Unsteady pressure predictions for the turbine test case were used to evaluate the forced response prediction capabilities of each of the three aeroelastic codes. Additionally, one of the fan flutter cases was evaluated using TURBO-AE. The UNSFLO and FREPS evaluation predictions showed good agreement with the experimental test data trends, but quantitative improvements are needed. UNSFLO over-predicted turbine blade response reductions, while FREPS under-predicted them. The inviscid TURBO-AE turbine analysis predicted no discernible blade response reduction, indicating the necessity of including viscous effects for this test case. For the TURBO-AE fan blisk test case, significant effort was expended getting the viscous version of the code to give converged steady flow solutions for the transonic flow conditions. Once converged, the steady solutions provided an excellent match with test data and the calibrated DAWES (AlliedSignal 3-D Viscous Steady Flow CFD Solver). However, efforts expended establishing quality steady-state solutions prevented exercising the unsteady portion of the TURBO-AE code during the present program. AlliedSignal recommends that unsteady pressure measurement data be obtained for both test cases examined for use in aeroelastic code validation.

Gardner, Kevin D.↗

The Influence of Environmental State on Lightning and Convective Parameter Distributions

A merged and highly reduced database of TRMM level 1 (precipitation radar, microwave imager, lightning) and NCEP reanalysis (basic state, radiative and surface flux) data has been assembled for three years of the TRMM mission. This allows direct examination of the dependence of convective spectra (as observed through radar reflectivity, microwave brightness temperature and lightning flash rate) on environmental basic states and anomalies. Such analysis may be more physically justified and instructive than traditional geographic and/or seasonal binning. The dependence of convective spectra on several environmental forcing parameters is presented, including surface Bowen ratio (sensible heat to total turbulent flux), net atmospheric radiative flux convergence and net atmospheric enthalpy flux convergence. The latter are basic drivers of net moisture convergence in simple quasi-equilibrium models of tropical atmospheric convection.

Boccippio, D.↗