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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Angular filter refractometry analysis using simulated annealing [An improved method for characterizing plasma density profiles using angular filter refractometry]

Here, a new method of analysis for angular filter refractometry images was developed to characterize laser-produced, long-scale-length plasmas using an annealing algorithm to iterative converge upon a solution. Angular filter refractometry (AFR) is a novel technique used to characterize the density pro files of laser-produced, long-scale-length plasmas. A synthetic AFR image is constructed by a user-defined density profile described by eight parameters, and the algorithm systematically alters the parameters until the comparison is optimized. The optimization and statistical uncertainty calculation is based on a minimization of the $\chi$2 test statistic. The algorithm was successfully applied to experimental data of plasma expanding from a flat, laser-irradiated target, resulting in average uncertainty in the density profile of 5-10% in the region of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Rigorous justification for the space–split sensitivity algorithm to compute linear response in Anosov systems

Ruelle gave a formula for linear response of transitive Anosov diffeomorphisms. Recently, practically computable realizations of Ruelle's formula have emerged that potentially enable sensitivity analysis of certain high-dimensional chaotic numerical simulations encountered in the applied sciences. In this paper, we provide full mathematical justification for the convergence of one such efficient computation, the space–split sensitivity, or S3, algorithm. In S3, Ruelle's formula is computed as a sum of two terms obtained by decomposing the perturbation vector field into a coboundary and a remainder that is parallel to the unstable direction. Such a decomposition results in a splitting of Ruelle's formula that is amenable to efficient computation. Here, we prove the existence of the S3 decomposition and the convergence of the computations of both resulting components of Ruelle's formula.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Fast Active-Set Thresholding Method for Nonnegative Least Squares

Nonnegative Least Squares (NNLS) is a fundamental constrained optimization problem encountered in many applications such as image deblurring, signal processing, nonnegative matrix factorization, magnetic microscopy, and hyperspectral imaging. Active-set based methods are a common class of algorithms for solving NNLS which identify the optimal variable set of the NNLS solution. They do so by iteratively solving a series of unconstrained least squares problems, identifying which variables violate the nonnegativity constraints, and then swapping variables in/out of consideration until the optimal set of variables is found. Several variations improving upon this method exist in the literature. In this work, we propose an active-set swap heuristic which further improves upon existing active-set based methods for NNLS. Our optimizations are based upon adding multiple variables to the passive set within a threshold of the smallest gradient value and removing variables within a similar threshold of the closest boundary constraint. We leverage these optimizations to yield a Fast Active-Set Thresholding NNLS (FAST-NNLS) algorithm which significantly outperforms the existing state-of-the-art NNLS algorithms for a wide range of problems. Rigorous convergence guarantees are proven for the proposed method. We demonstrate the effectiveness of our proposed method on multiple synthetic datasets and two realworld text analysis applications. In doing so, we present the most comprehensive NNLS solver comparison in the literature to date.

Cobb, Benjamin [Georgia Institute of Technology]↗
DOE OSTI2022

The Role of Ocean Circulation in Southern Ocean Heat Uptake, Transport, and Storage Response to Quadrupled CO 2

In response to quadrupled CO 2 , the Southern Ocean primarily uptakes excess heat around 60°S, which is then redistributed by the northward ocean heat transport (OHT) and mostly stored in the ocean or released back to the atmosphere around 45°S. However, the relative roles of mean ocean circulation and ocean circulation change in the uptake and redistribution of heat in the Southern Ocean remain controversial. Here, in this study, a set of climate model experiments embedded with a novel partial coupling technique are used to separate the roles of mean ocean circulation (passive component) and ocean circulation change (active component). For the ocean heat uptake (OHU) response, the mean ocean circulation and ocean circulation change are of equal importance. The OHT response south of 50°S is mainly determined by mean ocean circulation, while the ocean circulation change generates an anomalous southward OHT north of 50°S. A heat budget analysis finds that the divergence of passive OHT acts to balance the passive surface heat gain to the south of ~50°S, while the convergence of active OHT acts to balance the active surface heat loss to the north of ~50°S. Intriguingly, all the increase in ocean heat storage (OHS) is attributable to the passive component, with the ocean circulation change playing almost no role. In the Southern Ocean, both the active and the passive ocean heat transports are overcompensated by the reverse atmospheric heat transport via the Bjerknes compensation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Vibrations of cantilevered shallow cylindrical shells of rectangular planform

A cantilevered, shallow shell of circular cylindrical curvature and rectangular planform exhibits free vibration behavior which differs considerably from that of a cantilevered beam or of a flat plate. Some numerical results can be found for the problem in the previously published literature, mainly obtained by using various finite element methods. The present paper is the first definitive study of the problem, presenting accurate non-dimensional frequency parameters for wide ranges of aspect ratio, shallowness ratio and thickness ratio. The analysis is based upon shallow shell theory. Numerical results are obtained by using the Ritz method, with algebraic polynomial trial functions for the displacements. Convergence is investigated, with attention being given both to the number of terms taken for each co-ordinate direction and for each of the three components of displacement. Accuracy of the results is also established by comparison with finite element results for shallow shells and with other accurate flat plate solutions.

Leissa, A. W.↗
NASA NTRS1993

A thermal NO(x) prediction model - Scalar computation module for CFD codes with fluid and kinetic effects

A thermal NO(x) prediction model is developed to interface with a CFD, k-epsilon based code. A converged solution from the CFD code is the input to the postprocessing model for prediction of thermal NO(x). The model uses a decoupled analysis to estimate the equilibrium level of (NO(x))e which is the constant rate limit. This value is used to estimate the flame (NO(x)) and in turn predict the rate of formation at each node using a two-step Zeldovich mechanism. The rate is fixed on the NO(x) production rate plot by estimating the time to reach equilibrium by a differential analysis based on the reaction: O + N2 = NO + N. The rate is integrated in the nonequilibrium time space based on the residence time at each node in the computational domain. The sum of all nodal predictions yields the total NO(x) level.

Mcbeath, Giorgio↗
DOE OSTI2018

Asymptotically safe f(R)-gravity coupled to matter I: The polynomial case

We use the functional renormalization group equation for the effective average action to study the non-Gaussian renormalization group fixed points (NGFPs) arising within the framework of f(R)-gravity minimally coupled to an arbitrary number of scalar, Dirac, and vector fields. Based on this setting we provide comprehensible estimates whose gravity–matter systems give rise to NGFPs suitable for rendering the theory asymptotically safe. The analysis employs an exponential split of the metric fluctuations and retains a 7-parameter family of coarse-graining operators allowing the inclusion of non-trivial endomorphisms in the regularization procedure. For vanishing endomorphisms, it is established that gravity coupled to the matter content of the standard model of particle physics and many beyond the standard model extensions exhibit NGFPs whose properties are strikingly similar to the case of pure gravity: there are two UV-relevant directions and the position and critical exponents converge rapidly when higher powers of the scalar curvature are included. Conversely, none of the phenomenologically interesting gravity–matter systems exhibits a stable NGFP when a Type II coarse graining operator is employed. Our analysis resolves this tension by demonstrating that the NGFPs seen in the two settings belong to different universality classes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Structural, magnetic, and oxygen storage properties of hexagonal Dy{sub 1−x}Y{sub x}MnO{sub 3+δ}

Single-phase polycrystalline samples of hexagonal Dy{sub 1−x}Y{sub x}MnO{sub 3+δ} were synthesized over the full solubility range under reducing oxygen partial pressure conditions. Rietveld refinements using neutron powder diffraction data (NPD) confirmed the persistence of the parent hexagonal P6{sub 3}cm (Space group #185) structure for all the stoichiometric (δ=0) samples. The (Mn–O) bond lengths derived from NPD data were found to be consistent with non-degenerate Mn{sup +3} high-spin state e“{sup 2}e′{sup 2} in trigonal–bipyramidal coordination. Detailed thermogravimetric measurements show large reversible oxygen storage/release capabilities (within the range of oxygen content 3.0–3.3) at unusually low temperatures of 200–400 °C in air or oxygen atmospheres. A phase transition to a new hexagonal phase (Hex1) was observed to take place around δ∼0.29 for the oxygen-loaded Dy-rich samples (x<0.5). Analysis of x-ray and neutron diffraction data using Rietveld refinements demonstrate that the superstructure is formed by tripling the c-axis of the original P6{sub 3}cm unit cell. Modeling of the Hex1 superstructure converged to a structural symmetry consistent with the noncentrosymmetric space group R3c (#161). In magnetization measurements, a ferrimagnetic order for DyMnO{sub 3.0} was observed below 8 K. The dilution of Dy magnetism with nonmagnetic yttrium results in reducing both the transition temperature and the effective paramagnetic moment. - Graphical abstract: Superstructure of the (R3c) Dy{sub 0.7}Y{sub 0.3}MnO{sub 3.29} in the lattice range c=0.27611–0.42829 (left panel) along with its low-temperature oxygen absorption/desorption capability in pure O{sub 2} (right panel). - Highlights: • Polycrystalline samples of Dy{sub 1−x}Y{sub x}MnO{sub 3+δ} were synthesized in the hexagonal phase. • They exhibit large oxygen storage/release capabilities at low oxidation temperatures. • A phase transition to the R3c superstructure was observed around δ=0.29. • Both the Tc and effective paramagnetic moment were reduced with the Dy substitution.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SCALE Code System

The SCALE Code System is a widely-used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.SCALE 6.2 provides many new capabilities and significant improvements of existing features.New capabilities include:• ENDF/B-VII.1 nuclear data libraries CE and MG with enhanced group structures,• Neutron covariance data based on ENDF/B-VII.1 and supplemented with ORNL data,• Covariance data for fission product yields and decay constants,• Stochastic uncertainty and correlation quantification for any SCALE sequence with Sampler,• Parallel calculations with KENO,• Problem-dependent temperature corrections for CE calculations,• CE shielding and criticality accident alarm system analysis with MAVRIC,• CE depletion with TRITON (T5-DEPL/T6-DEPL),• CE sensitivity/uncertainty analysis with TSUNAMI-3D,• Simplified and efficient LWR lattice physics with Polaris,• Large scale detailed spent fuel characterization with ORIGAMI and ORIGAMI Automator,• Advanced fission source convergence acceleration capabilities with Sourcerer,• Nuclear data library generation with AMPX, and• Integrated user interface with Fulcrum.Enhanced capabilities include:• Accurate and efficient CE Monte Carlo methods for eigenvalue and fixed source calculations,• Improved MG resonance self-shielding methodologies and data,• Resonance self-shielding with modernized and efficient XSProc integrated into most sequences,• Accelerated calculations with TRITON/NEWT (generally 4x faster than SCALE 6.1),• Spent fuel characterization with 1470 new reactor-specific libraries for ORIGEN,• Modernization of ORIGEN (Chebyshev Rational Approximation Method [CRAM] solver, API for high-performance depletion, new keyword input format)• Extension of the maximum mixture number to values well beyond the previous limit of 2147 to ~2 billion,• Nuclear data formats enabling the use of more than 999 energy groups,• Updated standard composition library to provide more accurate use of natural abundances, andvi• Numerous other enhancements for improved usability and stability.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2014

Conjugate heat and mass transfer in the lattice Boltzmann equation method

An interface treatment for conjugate heat and mass transfer in the lattice Boltzmann equation method is proposed based on our previously proposed second-order accurate Dirichlet and Neumann boundary schemes. The continuity of temperature (concentration) and its flux at the interface for heat (mass) transfer is intrinsically satisfied without iterative computations, and the interfacial temperature (concentration) and their fluxes are conveniently obtained from the microscopic distribution functions without finite-difference calculations. The present treatment takes into account the local geometry of the interface so that it can be directly applied to curved interface problems such as conjugate heat and mass transfer in porous media. For straight interfaces or curved interfaces with no tangential gradient, the coupling between the interfacial fluxes along the discrete lattice velocity directions is eliminated and thus the proposed interface schemes can be greatly simplified. Several numerical tests are conducted to verify the applicability and accuracy of the proposed conjugate interface treatment, including (i) steady convection-diffusion in a channel containing two different fluids, (ii) unsteady convection-diffusion in the channel, (iii) steady heat conduction inside a circular domain with two different solid materials, and (iv) unsteady mass transfer from a spherical droplet in an extensional creeping flow. The accuracy and order of convergence of the simulated interior temperature (concentration) field, the interfacial temperature (concentration), and heat (mass) flux are examined in detail and compared with those obtained from the "half-lattice division" treatment in the literature. The present analysis and numerical results show that the half-lattice division scheme is second-order accurate only when the interface is fixed at the center of the lattice links, while the present treatment preserves second-order accuracy for arbitrary link fractions. For curved interfaces, the present treatment yields second-order accurate interior and interfacial temperatures (concentrations) and first-order accurate interfacial heat (mass) flux. An increase of order of convergence by one degree is obtained for each of these three quantities compared with the half-lattice division scheme. The surface-averaged Sherwood numbers computed in test (iv) agree well with published results.

Li, LK↗
NASA NTRS1988

Steady and unsteady transonic small disturbance analysis of realistic aircraft configurations

A transonic unsteady aerodynamic and aeroelasticity code CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) has been developed for application to realistic aircraft configurations. The code uses a time-accurate approximate factorization (AF) algorithm for solution of the unsteady transonic small-disturbance potential equation. The paper gives an overview of the CAP-TSD code development effort and reports on recent algorithm modifications. The algorithm modifications include: an Engquist-Osher (E-O) type-dependent switch to treat regions of supersonic flow, extension of the E-O switch for second-order spatial accuracy, nonisentropic effects to treat strong-shock cases, nonreflecting far field boundary conditions for unsteady applications, and several modifications to accelerate convergence to steady state. The modifications have significantly improved the stability of the AF algorithm and hence the reliability of the CAP-TSD code in general. Calculations are also presented from a flutter analysis of a 45-deg sweptback wing which agree well with experimental data. The paper present descriptions of the CAP-TSD code and algorithm details along with results.

Batina, John T.↗
NASA NTRS2017

New Flutter Analysis Technique for Time-Domain Computational Aeroelasticity

A new time-domain approach for computing flutter speed is presented. Based on the time-history result of aeroelastic simulation, the unknown unsteady aerodynamics model is estimated using a system identification technique. The full aeroelastic model is generated via coupling the estimated unsteady aerodynamic model with the known linear structure model. The critical dynamic pressure is computed and used in the subsequent simulation until the convergence of the critical dynamic pressure is achieved. The proposed method is applied to a benchmark cantilevered rectangular wing.

HWB aircraft↗
DOE OSTI2025

Spike-Free Adaptive Sliding Mode Control: Application to Permanent Magnet Synchronous Motors

A new methodology for adaptive sliding mode control (ASMC) has been widely used to improve the control performance in various systems. This method exhibits several advantages, including low sliding mode control (SMC) chattering, no knowledge of the system disturbance bound, and no overestimation of the control gain. Despite its advantages, this method can be hampered by the spike phenomenon, slow control gain convergence, and difficulty in achieving optimal performance under varying disturbances. Consequently, this article proposes a spike-free ASMC method with a disturbance observer (DOB) to address these problems. Previous ASMC methods have been analyzed via simulations to verify the aforementioned problems. Here, this analysis highlights the need for disturbance compensation and improvements in the SMC gain adaptation law. Therefore, a DOB is designed to mitigate the spike phenomenon by compensating for disturbances. Subsequently, an SMC gain adaptation law based on disturbance error estimation is designed to eliminate the spike phenomenon completely. The proposed adaptation law makes the SMC gain to converge to a slightly higher value than the disturbance estimation error. Consequently, the proposed method not only eliminates the spike phenomenon, but also ensures optimal performance under varying disturbances. The performance of the proposed method is experimen tally validated through a comparative study.

42 ENGINEERING↗
DOE OSTI2021

Revisiting the ODE Method for Recursive Algorithms: Fast Convergence Using Quasi Stochastic Approximation

Several decades ago, Profs. Sean Meyn and Lei Guo were postdoctoral fellows at ANU, where they shared interest in recursive algorithms. It seems fitting to celebrate Lei Guo's 60th birthday with a review of the ODE Method and its recent evolution. The method has been regarded as a technique for algorithm analysis. It is argued that this viewpoint is backwards: The original stochastic approximation method was surely motivated by an ODE, and tools for analysis came much later (based on establishing robustness of Euler approximations). The paper presents a brief survey of recent research in machine learning that shows the power of algorithm design in continuous time, following by careful approximation to obtain a practical recursive algorithm. While these methods are usually presented in a stochastic setting, this is not a prerequisite. In fact, recent theory shows that rates of convergence can be dramatically accelerated by applying techniques inspired by quasi Monte-Carlo. Subject to conditions, the optimal rate of convergence can be obtained by applying the averaging technique of Polyak and Ruppert. The conditions are not universal, but theory suggests alternatives to achieve acceleration. The theory is illustrated with applications to gradient-free optimization, and policy gradient algorithms for reinforcement learning.

learning and adaptive systems in artificial intell↗
NASA NTRS1979

Origin of longitudinal grooving in Tiu Vallis, Mars - Isolation of responsible fluid-types

Fluid instability analysis developed for the growth or decay of large, longitudinal roller vortex perturbations imposed at the base of fluids of exponentially stratified viscosities is presented. The origin of longitudinal grooving and fluting in a region of converging flow in Tiu Vallis, Mars, is analyzed, and the propensity of fluids, such as water, mudflow/debris flow, and lava, to develop longitudinal roll patterns as part of their flow field is tested. It is found that the grooving pattern in Tiu Vallis was probably formed by vortex motion in a mudflow or debris-laden slurry, but exact representation of sediment/debris concentration or sediment transport machanism is not yet possible. It is concluded that longitudinal vortex motion modelled in this analysis can develop in water, lava, or ice, but that the specific flow characteristics of these models which include the vortex motion are unreasonable for flow in Tiu Vallis.

Thompson, D. E.↗
DOE OSTI2017

A conservative nonlocal convection–diffusion model and asymptotically compatible finite difference discretization

Here in this paper, we first propose a nonlocal convection–diffusion model, in which the convection term is constructed in a special upwind manner so that mass conservation and maximum principle are maintained in any space dimension. The well-posedness of the proposed nonlocal model and its convergence to the classical local convection–diffusion model are established. A quadrature-based finite difference discretization is then developed to numerically solve the nonlocal problem and it is shown to be consistent and unconditionally stable. We further demonstrate that the numerical scheme is asymptotically compatible, that is, the approximate solutions converge to the exact solution of the corresponding local problem when δ → 0 and h → 0. Numerical experiments are also performed to complement the theoretical analysis.

42 ENGINEERING↗
DOE OSTI2024

Entire expressed peripheral blood transcriptome in pediatric severe malarial anemia

Abstract This study on severe malarial anemia (SMA: Hb < 6.0 g/dL), a leading global cause of childhood morbidity and mortality, compares the entire expressed whole blood host transcriptome between Kenyan children (3-48 mos.) with non-SMA (Hb ≥ 6.0 g/dL, n = 39) and SMA ( n = 18). Differential expression analyses reveal 1403 up-regulated and 279 down-regulated transcripts in SMA, signifying impairments in host inflammasome activation, cell death, and innate immune and cellular stress responses. Immune cell profiling shows decreased memory responses, antigen presentation, and immediate pathogen clearance, suggesting an immature/improperly regulated immune response in SMA. Module repertoire analysis of blood-specific gene signatures identifies up-regulation of erythroid genes, enhanced neutrophil activation, and impaired inflammatory responses in SMA. Enrichment analyses converge on disruptions in cellular homeostasis and regulatory pathways for the ubiquitin-proteasome system, autophagy, and heme metabolism. Pathway analyses highlight activation in response to hypoxic conditions [Hypoxia Inducible Factor (HIF)−1 target and Reactive Oxygen Species (ROS) signaling] as a central theme in SMA. These signaling pathways are also top-ranking in protein abundance measures and a Ugandan SMA cohort with available transcriptomic data. Targeted RNA-Seq validation shows strong concordance with our entire expressed transcriptome data. These findings identify key molecular themes in SMA pathogenesis, offering potential targets for new malaria therapies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Structure of the sensory domain of McpX from Sinorhizobium meliloti, the first known bacterial chemotactic sensor for quaternary ammonium compounds

The α-proteobacterium Sinorhizobium meliloti can live freely in the soil or engage in a symbiosis with its legume host. S. meliloti facilitates nitrogen fixation in root nodules, thus providing pivotal, utilizable nitrogen to the host. The organism has eight chemoreceptors, namely McpT to McpZ and IcpA that facilitate chemotaxis. McpX is the first known bacterial sensor of quaternary ammonium compounds (QACs) such as choline and betaines. Because QACs are exuded at chemotaxis-relevant concentrations by germinating alfalfa seeds, McpX has been proposed to contribute to host-specific chemotaxis. We have determined the crystal structure of the McpX periplasmic region (McpX PR ) in complex with the proline betaine at 2.7 Å resolution. In the crystal, the protein forms a symmetric dimer with one proline betaine molecule bound to each monomer of McpX PR within membrane-distal CACHE module. The ligand is bound through cation–π interactions with four aromatic amino acid residues. Mutational analysis in conjunction with binding studies revealed that a conserved aspartate residue is pivotal for ligand binding. We discovered that, in a striking example of convergent evolution, the ligand-binding site of McpX PR resembles that of a group of structurally unrelated betaine-binding proteins including ProX and OpuAC. Through this comparison and docking studies, we rationalized the specificity of McpX PR for this specific group of ligands. Collectively, our structural, biochemical, and molecular docking data have revealed the molecular determinants in McpX that are crucial for its rare ligand specificity for QACs.

Shrestha, Manisha↗