Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

An Impact Crater Origin for the InSight Landing Site at Homestead Hollow, Mars: Implications for Near Surface Stratigraphy, Surface Processes, and Erosion Rates

The InSight mission to Mars landed within Homestead hollow on an Early Amazonian lavaplain. The hollow is a 27‐m‐diameter, 0.3‐m‐deep quasi‐circular depression that shares morphologic and sedimentologic characteristics to degraded impact craters. Unlike the intercrater plains outside of the hollow, the interior lacks cobbles and is dominated by loose sand, granules, and pebbles. Fresher craters near the landing site exhibit meter‐scale bedforms in their ejecta and on their floors due to sediment trapping.The sedimentology of the interior fill of Homestead hollow suggests similar trapping. The hollow falls along amorphologic continuum that requires low rates of rim degradation and fill. Crater degradation rates (rim erosion plus filling) in the landing site decline nonlinearly through time from 10−2to 10−4m/Myr as craters evolve to a hollow‐like form. Rim erosion rates are lower initially, at 10−3m/Myr, but converge with degradation rates to 10−4m/Myr. This implies that while filling plays an important role soon after crater formation, it is limited in later stages. Crater statistics indicate that the bulk of the fill occurred in the first~50 Myr for Homestead hollow. The estimated maximum age of the hollow is ~400 to 700 Myr. This requires near‐zero fill aggradation and long‐term soil stability for the bulk of the crater's history. Fill stability manifests in Homestead hollow as a ~5‐to 10‐cm‐thick duricrust, formed by exchanges of atmospheric water vapor with soil. The estimated degradation in the hollow requires ~2 to 3 m of sedimentary fill beneath the lander.

Maria Elaine Banks↗
NASA NTRS1995

Locomotion training of legged robots using hybrid machine learning techniques

In this study artificial neural networks and fuzzy logic are used to control the jumping behavior of a three-link uniped robot. The biped locomotion control problem is an increment of the uniped locomotion control. Study of legged locomotion dynamics indicates that a hierarchical controller is required to control the behavior of a legged robot. A structured control strategy is suggested which includes navigator, motion planner, biped coordinator and uniped controllers. A three-link uniped robot simulation is developed to be used as the plant. Neurocontrollers were trained both online and offline. In the case of on-line training, a reinforcement learning technique was used to train the neurocontroller to make the robot jump to a specified height. After several hundred iterations of training, the plant output achieved an accuracy of 7.4%. However, when jump distance and body angular momentum were also included in the control objectives, training time became impractically long. In the case of off-line training, a three-layered backpropagation (BP) network was first used with three inputs, three outputs and 15 to 40 hidden nodes. Pre-generated data were presented to the network with a learning rate as low as 0.003 in order to reach convergence. The low learning rate required for convergence resulted in a very slow training process which took weeks to learn 460 examples. After training, performance of the neurocontroller was rather poor. Consequently, the BP network was replaced by a Cerebeller Model Articulation Controller (CMAC) network. Subsequent experiments described in this document show that the CMAC network is more suitable to the solution of uniped locomotion control problems in terms of both learning efficiency and performance. A new approach is introduced in this report, viz., a self-organizing multiagent cerebeller model for fuzzy-neural control of uniped locomotion is suggested to improve training efficiency. This is currently being evaluated for a possible patent by NASA, Johnson Space Center. An alternative modular approach is also developed which uses separate controllers for each stage of the running stride. A self-organizing fuzzy-neural controller controls the height, distance and angular momentum of the stride. A CMAC-based controller controls the movement of the leg from the time the foot leaves the ground to the time of landing. Because the leg joints are controlled at each time step during flight, movement is smooth and obstacles can be avoided. Initial results indicate that this approach can yield fast, accurate results.

Simon, William E.↗
DOE OSTI2022

High-repetition-rate, multi-MeV deuteron acceleration from converging heavy water microjets at laser intensities of 10 21 W/cm 2

Herein we demonstrate high repetition-rate deuteron acceleration by irradiating a continuously flowing, ambient temperature liquid heavy water jet with the high-intensity ALEPH laser. The laser delivered up to 5.5 J (120 TW, 1.2 × 10 21 W/cm 2 ) laser energy on target at 0.5 Hz. A high repetition-rate Thomson parabola spectrometer measured the deuteron beam energy spectra on each shot for 60 sequential shots (two minutes). Peak fluxes of 5 x 10 10 deuterons/sr/pulse, corresponding to an average flux of 1.5 x 10 12 deuterons/sr/min, were demonstrated with deuteron energies reaching up to 4.4 MeV. High shot-to-shot stability is observed up to 40%–50% of the maximum deuteron energy. These deuteron beams are suited for fast neutron production through deuteron breakup in a converter yielding energies similar to deuteron–deuteron (D–D, 2.45 MeV) fusion reactions of importance for material damage studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRS1976

The numerical evaluation of the maximum-likelihood estimate of a subset of mixture proportions

Necessary and sufficient conditions are given for a maximum likelihood estimate of a subset of mixture proportions. From these conditions, likelihood equations are derived satisfied by the maximum-likelihood estimate and a successive-approximations procedure is discussed as suggested by equations for numerically evaluating the maximum-likelihood estimate. It is shown that, with probability one for large samples, this procedure converges locally to the maximum-likelihood estimate whenever a certain step-size lies between zero and two. Furthermore, optimal rates of local convergence are obtained for a step-size which is bounded below by a number between one and two.

Peters, B. C., Jr.↗
NASA NTRS1990

First GPS baseline results from the North Andes

The CASA Uno GPS experiment (January-February 1988) has provided the first epoch baseline measurements for the study of plate motions and crustal deformation in and around the North Andes. Two dimensional horizontal baseline repeatabilities are as good as 5 parts in 10 to the 8th for short baselines (100-1000 km), and better than 3 parts in 10 to the 8th for long baselines (greater than 1000 km). Vertical repeatabilities are typically 4-6 cm, with a weak dependence on baseline length. The expected rate of plate convergence across the Colombia Trench is 6-8cm/yr, which should be detectable by the repeat experiment planned for 1991. Expected deformation rates within the North Andes are of the order of 1 cm/yr, which may be detectable with the 1991 experiment.

Kellogg, James N.↗
DOE OSTI2024

Variance-Reduced Accelerated First-Order Methods: Central Limit Theorems and Confidence Statements

In this paper, we consider a strongly convex stochastic optimization problem and propose three classes of variable sample-size stochastic first-order methods: (i) the standard stochastic gradient descent method, (ii) its accelerated variant, and (iii) the stochastic heavy-ball method. In each scheme, the exact gradients are approximated by averaging across an increasing batch size of sampled gradients. We prove that when the sample size increases at a geometric rate, the generated estimates converge in mean to the optimal solution at an analogous geometric rate for schemes (i)–(iii). Based on this result, we provide central limit statements, whereby it is shown that the rescaled estimation errors converge in distribution to a normal distribution with the associated covariance matrix dependent on the Hessian matrix, the covariance of the gradient noise, and the step length. If the sample size increases at a polynomial rate, we show that the estimation errors decay at a corresponding polynomial rate and establish the associated central limit theorems (CLTs). Under certain conditions, we discuss how both the algorithms and the associated limit theorems may be extended to constrained and nonsmooth regimes. As a result, we provide an avenue to construct confidence regions for the optimal solution based on the established CLTs and test the theoretical findings on a stochastic parameter estimation problem.

Lei, Jinlong↗
NASA NTRS1985

New stability theorems for averaging and their application to the convergence analysis of adaptive identification and control schemes

New stability theorems are developed for the convergence analysis of a class of one and two-time scale time-varying nonlinear systems using averaging theory. These theorems are applied to a class of continuous time adaptive identifiers and model reference adaptive controllers to obtain estimates of the parameter rate of exponential convergence.

Fu, L.-C.↗
NASA NTRS1984

Mesoscale budgets of heat and moisture in a convective system over the central United States

Mesoscale heat and moisture budget calculations from the April 10-11 AVE-SESAME central U.S. region are calculated in order to diagnose the effect of midlatitude organized convection on its environment. When averaged over an area greater than the observational scale, the noise levels of 5 C/day for the heat budget and 2 g/kg per day for the moisture budget are considerably reduced. The credibility of such area-averaged budget results is indicated by the general agreement obtained between the observed rainfall rate and the vertically integrated heat sources and moisture sinks. The vertical convective heating that is diagnosed from this organized midlatitude convection differs significantly from tropical structures diagnosed on much larger temporal and spatial scales. A time lag of several hours between the moisture convergence and rainfall rate is noted.

Kuo, Y.-H.↗
NASA NTRS1998

Ultra-Fast Fracture Strength of Advanced Ceramics at Elevated Temperatures

An attempt was made to determine elevated-temperature, 'ultra'-fast fracture strengths of one alumina, two silicon nitrides and one silicon carbide by using constant stress-rate ('dynamic fatigue') testing with a series of 'ultra'-fast test rates. Of the materials tested, the alumina exhibited a convergence of strength at stress rates below 3.3 x 10(exp 4) MPa/s. The strength approached approximately the room-temperature inert strength. By contrast, the silicon nitrides and silicon carbide did not reveal a strength approach, but exhibited elevated-temperature strengths 10 and 20% lower than their respective room-temperature strengths. Although the analytical results imply that the elevated-temperature 'inert' strength of a ceramic material can be obtained by using sufficiently high stress rates, the experimental testing rates were only sufficient to demonstrate convergence for the alumina.

Choi, Sung R.↗
DOE OSTI2022

Pseudodiagonalization Method for Accelerating Nonlinear Subspace Diagonalization in Density Functional Theory

In density functional theory, each self-consistent field (SCF) nonlinear step updates the discretized Kohn-Sham orbitals by solving a linear eigenvalue problem. The concept of pseudodiagonalization is to solve this linear eigenvalue problem approximately, and specifically utilizing a method involving a small number of Jacobi rotations that takes advantage of the good initial guess to the solution given by the approximation to the orbitals from the previous SCF iteration. The approximate solution to the linear eigenvalue problem can be very rapid, particularly for those steps near SCF convergence. Here, we adapt pseudodiagonalization to finite-temperature and metallic systems, where partially-occupied orbitals must be individually resolved with some accuracy. We apply pseudodiagonalization to the subspace eigenvalue problem that arises in Chebyshev-filtered subspace iteration. In tests on metallic and other systems for a range of temperatures, we show that pseudodiagonalization achieves similar rates of SCF convergence to exact diagonalization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Precipitation Efficiency in the Tropical Deep Convective Regime

Precipitation efficiency in the tropical deep convective regime is analyzed based on a 2-D cloud resolving simulation. The cloud resolving model is forced by the large-scale vertical velocity and zonal wind and large-scale horizontal advections derived from TOGA COARE for a 20-day period. Precipitation efficiency may be defined as a ratio of surface rain rate to sum of surface evaporation and moisture convergence (LSPE) or a ratio of surface rain rate to sum of condensation and deposition rates of supersaturated vapor (CMPE). Moisture budget shows that the atmosphere is moistened (dryed) when the LSPE is less (more) than 100 %. The LSPE could be larger than 100 % for strong convection. This indicates that the drying processes should be included in cumulus parameterization to avoid moisture bias. Statistical analysis shows that the sum of the condensation and deposition rates is bout 80 % of the sum of the surface evaporation rate and moisture convergence, which ads to proportional relation between the two efficiencies when both efficiencies are less han 100 %. The CMPE increases with increasing mass-weighted mean temperature and creasing surface rain rate. This suggests that precipitation is more efficient for warm environment and strong convection. Approximate balance of rates among the condensation, deposition, rain, and the raindrop evaporation is used to derive an analytical solution of the CMPE.

Li, Xiaofan↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Self-regulating behavior of hybrid membrane systems as demonstrated in an element-scale forward osmosis-reverse osmosis hybrid system

Hybrid membrane systems can be difficult to design due to the requisite flow rate matching between up- and downstream unit operations. In this work, we use a forward osmosis-reverse osmosis (FO-RO) hybrid system to demonstrate how some membrane systems can exhibit self-regulating behavior due to osmotic coupling. This can reduce the need for complex control systems for flow balancing. We show this behavior using a module-scale test bed that can mimic the behavior of larger scale operations. The system shows permeate flow rate near-convergence between the FO and RO modules after startup or when perturbed by a change in RO module pressure. The behavior of this hybrid system demonstrates that some membrane operations can exploit osmotic interdependence, rather than expensive control systems, to achieve steady state operation.

Debottlenecking↗
DOE OSTI2024

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

An Analytical Investigation of Wing-Jet Interaction

The aerodynamic interaction between the wing and an inviscid jet with Mach number nonuniformity is formulated by using a two-vortex-sheet model for the jet. One of the vortex sheets accounts for the induced jet flow and the other the induced outer flow. No additional source distribution is needed for the jet at an angle of attack. The above problem is solved by satisfying the jet and wing tangency and the jet pressure-continuity conditions and using a quasi vortex lattice method for computing the induced flow field. The latter method is derived through theoretical consideration by properly accounting for singularities present in the equations and possesses the same simplicity and generality as the conventional vortex lattice method but has a better rate of numerical convergence. The resulting system of algebraic equations is solved by Purcell's vector method. The numerical formulation is first applied to the wing-slipstream interaction problem. Results for one centered-jet configuration are compared with those predicted by some existing theories.

Lan, E. C.↗
NASA NTRS1988

Karhunen-Loeve expansion of Burgers' model of turbulence

The properties of the Karhunen-Loeve expansion of a strongly inhomogeneous random process are examined with emphasis on applications to turbulent flow fields. The ability of the KL expansion to represent functions that have both slow and rapid variations in a relatively small number of expansion terms is tested on a one-dimensional model based on the forced Burgers' equation. The rate of the convergence of the expansion is evaluated, and its dependence on the Reynolds number is determined. It is shown that the KL eigenfunctions possess wall boundary layers attached to outer structures that are independent of the Reynolds number (at high Reynolds numbers). It is also shown that the spectrum of eigenvalues is broad at large Reynolds numbers, requiring many terms to represent higher-order derivatives of the function.

Chambers, D. H.↗