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 775 records · Page 43

All sources NASA NTRS DOE OSTI

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 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 OSTI2022

Kingdom-Wide Analysis of Fungal Protein-Coding and tRNA Genes Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Large Scale Bilevel Optimization for N-K SCOPF Using Adversarial Robustness

Ensuring a secure dispatch against multiple simultaneous outages has long been desired to maintain grid security in the presence of severe events, such as extreme weather phenomena. Traditionally denoted as N-k security constrained optimal power flow (N-k SCOPF), this problem is intractable to solve due to its size being combinatorial in the number of simultaneous outages and due to the non-convex nature of the AC network constraints. This hinders the use of N-k SCOPF for operating realistic-scale systems. In this paper, we introduce a methodology to scalably solve an AC-feasible dispatch that improves security over k simultaneous outages. Our methodology poses N-k SCOPF as a bilevel optimization problem and solves it using an adversarial robustness approach. We develop new efficient methods to solve each level of the bilevel optimization by employing knowledge of the physics of the underlying system. This yields significant improvements in speed and convergence that enable us to address the N-k SCOPF problem at scale. We demonstrate the effectiveness of our method by conducting a comprehensive analysis of an N-3 SCOPF for a 500-bus network. Furthermore, we emphasize the ability of our physics-driven techniques to handle larger systems by successfully scaling up to 12,000 buses.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Unitarity cuts of the worldsheet

We compute the imaginary parts of genus-one string scattering amplitudes. Following Witten’s i\varepsilon i ε prescription for the integration contour on the moduli space of worldsheets, we give a general algorithm for computing unitarity cuts of the annulus, Möbius strip, and torus topologies exactly in \alpha' α ′ . With the help of tropical analysis, we show how the intricate pattern of thresholds (normal and anomalous) opening up arises from the worldsheet computation. The result is a manifestly-convergent representation of the imaginary parts of amplitudes, which has the analytic form expected from Cutkosky rules in field theory, but bypasses the need for performing laborious sums over the intermediate states. We use this representation to study various physical aspects of string amplitudes, including their behavior in the (s,t) ( s , t ) plane, exponential suppression, decay widths of massive strings, total cross section, and low-energy expansions. We find that planar annulus amplitudes exhibit a version of low-spin dominance: at any finite energy, only a finite number of low partial-wave spins give an appreciable contribution to the imaginary part.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Effects of asymmetry on the dynamic stability of aircraft

The oblique wing concept for transonic aircraft was proposed to reduce drag. The dynamic stability of the aircraft was investigated by analytically determining the stability derivatives at angles of skew ranging from 0 and 45 deg and using these stability derivatives in a linear analysis of the coupled aircraft behavior. The stability derivatives were obtained using a lifting line aerodynamic theory and found to give reasonable agreement with derivatives developed in a previous study for the same aircraft. In the dynamic analysis, no instability or large changes occurred in the root locations for skew angles varying from 0 to 45 deg with the exception of roll convergence. The damping in roll, however, decreased by an order of magnitude. Rolling was a prominent feature of all the oscillatory mode shapes at high skew angles.

Fantino, R. E.↗

NASA NTRS1990

Edge detection applied to SST fields

An algorithm designed to detect fronts automatically in satellite-derived sea-surface temperature (SST) fields is presented. The algorithm is operated at different levels to detect and differentiate between false and true edges. For purposes of comparison, the algorithm is applied to a test set of 98 SST images to detect the northern edge of the Gulf Stream. The algorithm successfully detected valid temperature fronts and ignored false edges, and also produced statistics about the temperature fronts that are useful in the subsequent analysis of these fronts. It is assumed that the algorithm performs equally well on other SST fronts such as those associated with rings, the subtropical convergence, or the shelf/slope fronts.

Cayula, Jean-Francois↗

NASA NTRS1994

Dynamics of Numerics and CFD

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Yee, H. C.↗

NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗

NASA NTRS1985

Analysis of experimental shaft seal data for high-performance turbomachines, as for Space Shuttle main engines

High-pressure, high-temperature seal flow (leakage) data for nonrotating and rotating Raleigh-step and convergent-tapered-bore seals were characterized in terms of a normalized flow coefficient. The data for normalized Rayleigh-steip and nonrotating tapered-bore seals were in reasonable agreement with theory, but data for the rotating tapered-bore seals were not. The tapered-bore-seal operational clearances estimated from the flow data were significantly larger than calculated. Although clearances are influenced by wear from conical to cylindrical geometry and errors in clearance corrections, the problem was isolated to the shaft temperature - rotational speed clearance correction. The geometric changes support the use of some conical convergence in any seal. Under these conditions rotation reduced the normalized flow coefficiently by nearly 10 percent.

Hendricks, R. C.↗

NASA NTRS1987

Analysis of experimental shaft seal data for high-performance turbomachines - As for Space Shuttle main engines

High-pressure, high-temperature seal flow (leakage) data for nonrotating and rotating Raleigh-step and convergent-tapered-bore seals were characterized in terms of a normalized flow coefficient. The data for normalized Rayleigh-step and nonrotating tapered-bore seals were in reasonable agreement with theory, but data for the rotating tapered-bore seals were not. The tapered-bore-seal operational clearances estimated from the flow data were significantly larger than calculated. Although clearances are influenced by wear from conical to cylindrical geometry and errors in clearance corrections, the problem was isolated to the shaft temperature - rotational speed clearance correction. The geometric changes support the use of some conical convergence in any seal. Under these conditions rotation reduced the normalized flow coefficiently by nearly 10 percent.

Hendricks, R. C.↗

NASA NTRS2004

Numerical Investigation of Dual-Mode Scramjet Combustor with Large Upstream Interaction

Dual-mode scramjet combustor configuration with significant upstream interaction is investigated numerically, The possibility of scaling the domain to accelerate the convergence and reduce the computational time is explored. The supersonic combustor configuration was selected to provide an understanding of key features of upstream interaction and to identify physical and numerical issues relating to modeling of dual-mode configurations. The numerical analysis was performed with vitiated air at freestream Math number of 2.5 using hydrogen as the sonic injectant. Results are presented for two-dimensional models and a three-dimensional jet-to-jet symmetric geometry. Comparisons are made with experimental results. Two-dimensional and three-dimensional results show substantial oblique shock train reaching upstream of the fuel injectors. Flow characteristics slow numerical convergence, while the upstream interaction slowly increases with further iterations. As the flow field develops, the symmetric assumption breaks down. A large separation zone develops and extends further upstream of the step. This asymmetric flow structure is not seen in the experimental data. Results obtained using a sub-scale domain (both two-dimensional and three-dimensional) qualitatively recover the flow physics obtained from full-scale simulations. All results show that numerical modeling using a scaled geometry provides good agreement with full-scale numerical results and experimental results for this configuration. This study supports the argument that numerical scaling is useful in simulating dual-mode scramjet combustor flowfields and could provide an excellent convergence acceleration technique for dual-mode simulations.

Mohieldin, T. O.↗

DOE OSTI2025

Nontrivial critical behavior at magnetic transitions: A case study of Sm 7 ⁢Pd 3

We present a comprehensive analysis of the critical behavior of Sm 7 ⁢Pd 3 in the vicinity of its second-order magnetoelastic transition at 𝑇 c =173 K. The critical exponents (CEs) 𝛽 and 𝛾, determined using both the standard convergence procedure and the average normalized slope (ANS) method, diverge at 𝑇 c –-a characteristic typically associated with first-order transitions. Notably, none of the established universality classes satisfactorily describe the critical behavior of Sm 7 ⁢Pd 3 , and we discuss the possible origins of this deviation in the context of the strong spin-lattice coupling intrinsic to the sample. We emphasize the importance of accurately selecting the critical temperature and magnetic field ranges to ensure robust critical behavior analysis, and propose a quantitative approach to assess the reliability of the extracted CEs. Additionally, we demonstrate that in the ANS method, the critical exponents 𝛽 and 𝛾 should be calculated separately using data for 𝑇 ⩽ 𝑇 c and 𝑇 ⩾ 𝑇 c , respectively. In conclusion, our findings underscore the need for a revised theoretical framework to accurately describe second-order magnetoelastic transitions.

Ferrimagnets↗

DOE OSTI2020

A Mixed Discontinuous Galerkin Method for Linear Elasticity with Strongly Imposed Symmetry

In this paper, we study a mixed discontinuous Galerkin (MDG) method to solve linear elasticity problem with arbitrary order discontinuous finite element spaces in d -dimension ( \(d=2,3\) ). This method uses polynomials of degree \(k+1\) for the stress and of degree k for the displacement ( \(k\ge 0\) ). The mixed DG scheme is proved to be well-posed under proper norms. Specifically, we prove that, for any \(k \ge 0\) , the \(H(\mathrm{div})\) -like error estimate for the stress and \(L^2\) error estimate for the displacement are optimal. We further establish the optimal \(L^2\) error estimate for the stress provided that the \({\mathcal {P}}_{k+2}-{\mathcal {P}}_{k+1}^{-1}\) Stokes pair is stable and \(k \ge d\) . Finally, we also provide numerical results of MDG showing that the orders of convergence are actually sharp.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

On the accuracy of limiters and convergence to steady state solutions

This paper addresses the practical problem of obtaining convergence to steady state solutions when limiters are used in conjunction with upwind schemes on unstructured grids. The base scheme forms a gradient and limits it by imposing monotonicity conditions in the reconstruction stage. It is shown by analysis in one dimension that such an approach leads to various schemes meeting TVD requirements in one dimension. It is further shown that these formally second order accurate schemes are less than second order accurate in practice because of the action of the limiter function in smooth regions of the solution. Modifications are proposed to the limiter that restore the second order accuracy. In multiple dimensions these schemes produce steady state solutions that are monotone and devoid of oscillations. However, convergence stalls after a few orders of reduction in the residual. With the modified limiter, on the other hand, it is shown that converged steady state solutions can be obtained.

Venkatakrishnan, V.↗

DOE OSTI2020

Phase retrieval for short wavelength orbital angular momentum beams using knife-edge diffraction

Laser beams operating in high order transverse modes, especially orbital angular momentum (OAM) modes, have been used increasingly in many frontier research areas. The characterization of these laser modes and the related component analysis can be carried out using the retrieved optical phase. However, conventional iterative phase retrieval methods, when applied to pure OAM modes, may fail because their convergence depends on the choice of the initial phase. Here, we present an improved phase retrieval method which can incorporate additional phase information obtained using knife-edge diffraction. Using no additional optics or few optics available in a wide spectral range, this method is particularly suitable for short wavelength lasers such as free-electron lasers (FELs) in the VUV, EUV, and X-ray regimes. Finally, this method has been tested using simulated beams, and successfully applied to experimentally generated beams from a HeNe laser and an FEL.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1973

A convergence theory for a class of nonlinear programming problems.

A recent convergence theory of Elkin concerning methods for unconstrained minimization is extended to a certain class of nonlinear programming problems. As in Elkin's original approach, the analysis of a variety of step-length algorithms is treated entirely separately from that of several direction algorithms. This allows for their combination into many different methods for solving the constrained problem. These include some of the methods of Rosen and Zoutendijk. We also extend the results of Topkis and Veinott to nonconvex sets and drop their requirement of the uniform feasibility of a subsequence of the search directions.

Rauch, S. W.↗