SEARCH · Engineering Papers

Results for “convergence analysis”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Systematic Assessment of Retrieval Methods for Canopy Far-Red Solar-Induced Chlorophyll Fluorescence Using High-Frequency Automated Field Spectroscopy

Remote sensing of solar-induced chlorophyll fluorescence (SIF) offers potential to infer photosynthesis across scales and biomes. Many retrieval methods have been developed to estimate top-of-canopy SIF using ground-based spectroscopy. However, inconsistencies among methods may confound interpretation of SIF dynamics, eco-physiological/environmental drivers, and its relationship with photosynthesis. Using high temporal- and spectral resolution ground-based spectroscopy, we aimed to (1) evaluate performance of SIF retrieval methods under diverse sky conditions using continuous field measurements; (2) assess method sensitivity to fluctuating light, reflectance, and fluorescence emission spectra; and (3) inform users for optimal ground-based SIF retrieval. Analysis included field measurements from bi-hemispherical and hemispherical-conical systems and synthetic upwelling radiance constructed from measured downwelling radiance, simulated reflectance, and simulated fluorescence for benchmarking. Fraunhofer-based differential optical absorption spectroscopy (DOAS) and singular vector decomposition (SVD) retrievals exhibit convergent SIF-PAR relationships and diurnal consistency across different sky conditions, while O2A-based spectral fitting method (SFM), SVD, and modified Fraunhofer line discrimination (3FLD) exhibit divergent SIF-PAR relationships across sky conditions. We find that such behavior holds across system configurations, though hemispherical-conical systems diverge less across sky conditions. O 2 A retrieval accuracy, influenced by atmospheric distortion, improves with a narrower fitting window and when training SVD with temporally local spectra. This may impact SIF-photosynthesis relationships interpreted by previous studies using O 2 A-based retrievals with standard (759–767.76 nm) fitting windows. Fraunhofer-based retrievals resist atmospheric impacts but are noisier and more sensitive to assumed SIF spectral shape than O 2 A-based retrievals. We recommend SVD or SFM using reduced fitting window (759.5–761.5 nm) for robust far-red SIF retrievals across sky conditions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Assessing correlated truncation errors in modern nucleon-nucleon potentials

We test the BUQEYE model of correlated effective field theory (EFT) truncation errors on Reinert, Krebs, and Epelbaum's semilocal momentum-space implementation of the chiral EFT (𝜒⁢EFT ) expansion of the nucleon-nucleon (NN) potential. This Bayesian model hypothesizes that dimensionless coefficient functions extracted from the order-by-order corrections to NN observables can be treated as draws from a Gaussian process (GP). We combine a variety of graphical and statistical diagnostics to assess when predicted observables have a 𝜒⁢EFT convergence pattern consistent with the hypothesized GP statistical model. Our conclusions are that, first, the BUQEYE model is generally applicable to the potential investigated here, which enables statistically principled estimates of the impact of higher EFT orders on observables. Second, parameters defining the extracted coefficients such as the expansion parameter 𝑄 must be well chosen for the coefficients to exhibit a regular convergence pattern—a property we exploit to obtain posterior distributions for such quantities. Third, the assumption of GP stationarity across lab energy and scattering angle is not generally met; this necessitates adjustments in future work. We provide a workflow and interpretive guide for our analysis framework, and show what can be inferred about probability distributions for 𝑄, the EFT breakdown scale Λ 𝑏 , the scale associated with soft physics in the 𝜒⁢EFT potential 𝑚 eff , and the GP hyperparameters. All our results can be reproduced using a publicly available Jupyter notebook, which can be straightforwardly modified to analyze other 𝜒⁢EFT NN potentials.

Bayesian methods↗

DOE OSTI2025

Identification of key steps in the evolution of anaerobic methanotrophy in Candidatus Methanovorans (ANME-3) archaea

Despite their large environmental impact and multiple independent emergences, the processes leading to the evolution of anaerobic methanotrophic archaea (ANME) remain unclear. This work uses comparative metagenomics of a recently evolved but understudied ANME group, “Candidatus Methanovorans” (ANME-3), to identify evolutionary processes and innovations at work in ANME, which may be obscured in earlier evolved lineages. We identified horizontal transfer of hdrA homologs and convergent evolution in carbon and energy metabolic genes as potential early steps in Methanovorans evolution. We also identified the erosion of genes required for methylotrophic methanogenesis along with horizontal acquisition of multiheme cytochromes and other loci uniquely associated with ANME. The assembly and comparative analysis of multiple Methanovorans genomes offers important functional context for understanding the niche-defining metabolic differences between methane-oxidizing ANME and their methanogen relatives. Furthermore, this work illustrates the multiple evolutionary modes at play in the transition to a globally important metabolic niche.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

UHT-CAMANCHE: Ultra-High Temperature Ceramic Additively Manufactured Compact Heat Exchangers

The conceptual basis for this project is the convergence of advanced ultra-high temperature ceramic materials and additive manufacturing technologies to produce compact ceramic heat exchangers with complex internal flow path geometries. Task areas were broadly divided into materials and manufacturing development, heat exchanger design, component testing, and techno-economic analysis. Technical challenges included the design and commissioning of new test facilities, improving feature resolution and deposition rate of ceramic additive manufacturing techniques, establishing process-structure-property relationships in additively manufactured ultra-high temperature ceramics, and assessing high temperature materials compatibility in CO 2 environments. The primary candidate material evaluated in this work is a composite comprising zirconium diboride (ZrB2) with 30 vol. % silicon carbide (SiC) which was selected based on its desirable combination of high temperature mechanical properties, high thermal and electrical conductivities, and oxidation resistance. High solids loaded ZrB2-SiC pastes suitable for extrusion-based additive manufacturing were developed for the first time as part of this work. Materials compatibility studies indicate this material oxidizes in CO 2 to form a protective borosilicate scale which transforms to pure silica above 1000°C. Parts made by additive manufacturing displayed enlarged grain sizes produced by pressureless sintering as compared to hot-press sintering. Increases in microstructural coarseness have outsized effect on oxidation performance up to 1400°C due to incomplete oxidation of coarse large diameter SiC particles resulting in lower amounts of silica that apparently inhibit protective scale formation. Additive manufacturing as a forming technique did not appear to significantly affect thermal conductivity, hardness, or elastic modulus, though flexural strength was reduced by half or more as compared to traditionally hot-pressed materials. This effect was attributed to the presence of strength-limiting flaws (ca. 40 microns in size) originating from extrudate inhomogeneities that could potentially be eliminated with further process improvements. Attempts to attain economies of scale for production of multi-kilowatt scale heat exchangers by ceramic additive manufacturing proved difficult. Lack of automation and a modest extrusion rate while retaining fine feature resolution made the overall process labor intensive and limited experimental throughput. A number of full-scale components were taken through post-process heat treatments including drying, binder burnout, and sintering; however, none survived without significant flaws or cracks. Therefore, no operational data from a newly installed heat exchanger test loop were able to be obtained during the performance period. Continued research and development is recommended to improve economic feasibility of ceramic additive manufacturing by standardizing the use of advanced sensors, artificial intelligence, and automation tools to reduce associated labor costs and accelerate production rates. The materials and manufacturing techniques demonstrated in this work are likely to find applications in defense and energy applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1989

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗

NASA NTRS1991

Three-Dimensional Navier-Stokes Simulation of Space Shuttle Main Propulsion 17-inch Disconnect Valves

A steady incompressible three-dimensional viscous flow analysis has been conducted for the Space Shuttle external tank/orbiter propellant feed line disconnect flapper valves with upstream elbows. The Navier-Stokes code, INS3D, is modified to handle interior obstacles and a simple turbulence model. The flow solver is tested for stability and convergence in the presence of interior flappers. An under-relaxation scheme has been incorporated to improve the solution stability. Important flow characteristics such as secondary flows, recirculation, vortex and wake regions, and separated flows are observed. Computed values for forces, moments, and pressure drop are in satisfactory agreement with water flow test data covering a maximum tube Reynolds number of 3.5 million. The predicted hydrodynamical stability of the flappers correlates well with the measurements.

Kandula, M.↗

NASA NTRS2000

Advancement of Bi-Level Integrated System Synthesis (BLISS)

Bi-Level Integrated System Synthesis (BLISS) is a method for optimization of an engineering system, e.g., an aerospace vehicle. BLISS consists of optimizations at the subsystem (module) and system levels to divide the overall large optimization task into sets of smaller ones that can be executed concurrently. In the initial version of BLISS that was introduced and documented in previous publications, analysis in the modules was kept at the early conceptual design level. This paper reports on the next step in the BLISS development in which the fidelity of the aerodynamic drag and structural stress and displacement analyses were upgraded while the method's satisfactory convergence rate was retained.

Sobieszczanski-Sobieski, J.↗

NASA NTRS2000

Advancement of Bi-Level Integrated System Synthesis (BLISS)

Sobieszczanski-Sobieski, Jaroslaw↗

NASA NTRS2000

Advancement of Bi-Level Integrated System Synthesis (BLISS)

Sobieszczanski-Sobieski, Jaroslaw↗

DOE OSTI2020

PtychoShelves , a versatile high-level framework for high-performance analysis of ptychographic data

Over the past decade, ptychography has been proven to be a robust tool for non-destructive high-resolution quantitative electron, X-ray and optical microscopy. It allows for quantitative reconstruction of the specimen's transmissivity, as well as recovery of the illuminating wavefront. Additionally, various algorithms have been developed to account for systematic errors and improved convergence. With fast ptychographic microscopes and more advanced algorithms, both the complexity of the reconstruction task and the data volume increase significantly. PtychoShelves is a software package which combines high-level modularity for easy and fast changes to the data-processing pipeline, and high-performance computing on CPUs and GPUs.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2004

Subsonic Aircraft With Regression and Neural-Network Approximators Designed

At the NASA Glenn Research Center, NASA Langley Research Center's Flight Optimization System (FLOPS) and the design optimization testbed COMETBOARDS with regression and neural-network-analysis approximators have been coupled to obtain a preliminary aircraft design methodology. For a subsonic aircraft, the optimal design, that is the airframe-engine combination, is obtained by the simulation. The aircraft is powered by two high-bypass-ratio engines with a nominal thrust of about 35,000 lbf. It is to carry 150 passengers at a cruise speed of Mach 0.8 over a range of 3000 n mi and to operate on a 6000-ft runway. The aircraft design utilized a neural network and a regression-approximations-based analysis tool, along with a multioptimizer cascade algorithm that uses sequential linear programming, sequential quadratic programming, the method of feasible directions, and then sequential quadratic programming again. Optimal aircraft weight versus the number of design iterations is shown. The central processing unit (CPU) time to solution is given. It is shown that the regression-method-based analyzer exhibited a smoother convergence pattern than the FLOPS code. The optimum weight obtained by the approximation technique and the FLOPS code differed by 1.3 percent. Prediction by the approximation technique exhibited no error for the aircraft wing area and turbine entry temperature, whereas it was within 2 percent for most other parameters. Cascade strategy was required by FLOPS as well as the approximators. The regression method had a tendency to hug the data points, whereas the neural network exhibited a propensity to follow a mean path. The performance of the neural network and regression methods was considered adequate. It was at about the same level for small, standard, and large models with redundancy ratios (defined as the number of input-output pairs to the number of unknown coefficients) of 14, 28, and 57, respectively. In an SGI octane workstation (Silicon Graphics, Inc., Mountainview, CA), the regression training required a fraction of a CPU second, whereas neural network training was between 1 and 9 min, as given. For a single analysis cycle, the 3-sec CPU time required by the FLOPS code was reduced to milliseconds by the approximators. For design calculations, the time with the FLOPS code was 34 min. It was reduced to 2 sec with the regression method and to 4 min by the neural network technique. The performance of the regression and neural network methods was found to be satisfactory for the analysis and design optimization of the subsonic aircraft.

Patnaik, Surya N.↗

NASA NTRS2019

Structural Analysis and Performance-Based Validation of a Composite Wing Spar

Electric-motor powered aircraft possess the ability to operate with efficient energy delivery, but lack the operational range of internal combustion engine powered aircraft. This range limitation requires the use of high aspect ratio, thin-chord wings to minimize aerodynamic drag losses, which results in highly loaded composite spar structures. High aspect ratio wings are required to increase mission durations for a NASA-developed experimental multi-rotor electric powered aircraft denoted as the Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) or “X-57”. This paper examines the structural performance of the composite main wing spars to validate spar strength using ply-based laminate finite element methods. Geometric scaling of a main spar test-section was initially proposed for proof-testing but sacrificed stability. Ply-based structures modeling with local structural features was implemented as a risk-reduction methodology. Ply-based modeling was selected to augment the conventional “building block” approach to reduce risk, and leverage a “performance-based” approval processes encouraged in Federal Aviation Administration (FAA) design guidance. Therefore, ply-based laminate modeling of the full-scale main spar and forward spar shear-web attachments were subsequently undertaken to determine load path complexity with predicted flight loads. Ply-based modeling included stress concentrations and interlaminate behavior at interface locations that can be obscured in traditional finite element sizing models. Analysis of the wing spar laminate ply-based models compared with bearing test coupon performance was used to reduce future wing assembly proof-testing burden and facilitate performance-based flight hardware safety for the X-57 experimental aircraft.

Moore, James B.↗

NASA NTRS1975

Numerical analysis and parametric studies of the buckling of composite orthotropic compression and shear panels

A computer program is presented which was developed for the combined compression and shear of stiffened variable thickness orthotropic composite panels on discrete springs: boundary conditions are general and include elastic boundary restraints. Buckling solutions are obtained by using a newly developed trigonometric finite difference procedure which improves the solution convergence rate over conventional finite difference methods. The classical general shear buckling results which exist only for simply supported panels over a limited range of orthotropic properties, were extended to the complete range of these properties for simply supported panels and, in addition, to the complete range of orthotropic properties for clamped panels. The program was also applied to parametric studies which examine the effect of filament orientation upon the buckling of graphite-epoxy panels. These studies included an examination of the filament orientations which yield maximum shear or compressive buckling strength for panels having all four edges simply supported or clamped over a wide range of aspect ratios. Panels with such orientations had higher buckling loads than comparable, equal weight, thin skinned aluminum panels. Also included among the parameter studies were examinations of combined axial compression and shear buckling and examinations of panels with rotational elastic edge restraints.

Housner, J. M.↗

NASA NTRS1984

Perturbation analysis of the limit cycle of the free van der Pol equation

A power series expansion in the damping parameter, epsilon, of the limit cycle of the free van der Pol equation is constructed and analyzed. Coefficients in the expansion are computed in exact rational arithmetic using the symbolic manipulation system MACSYMA and using a FORTRAN program. The series is analyzed using Padeapproximants. The convergence of the series for the maximum amplitude of the limit cycle is limited by two pair of complex conjugate singularities in the complex epsilon-plane. A new expansion parameter is introduced which maps these singularities to infinity and leads to a new expansion for the amplitude which converges for all real values of epsilon. Amplitudes computed from this transformed series agree very well with reported numerical and asymptotic results. For the limit cycle itself, convergence of the series expansion is limited by three pair of complex conjugate branch point singularities. Two pair remain fixed throughout the cycle, and correspond to the singularities found in the maximum amplitude series, while the third pair moves in the epsilon-plane as a function of t from one of the fixed pairs to the other. The limit cycle series is transformed using a new expansion parameter, which leads to a new series that converges for larger values of epsilon.

Dadfar, M. B.↗

NASA NTRS2010

Orion Launch Abort Vehicle Separation Analysis Using OVERFLOW

This slide presentation reviews the use of OVERFLOW, a flow solver, to analyze the effect of separation for a launch abort vehicle (i.e., Orion capsule) if required. Included in the presentation are views of the geometry, and the Overset grids, listing of the assumptions, the general run strategy, inputs into the Overflow solver, the required computational resources, the results of the convergence study. Charts and graphics are presented to show the results.

Booth, Tom↗

DOE OSTI2025

Benchmarking third-order cluster perturbation theory for electronically excited states

In this study, we investigate the reliability of cluster perturbation (CP) theory applied to the calculation of electronically excited states through a comprehensive benchmark. In CP theory, perturbative corrections are added to the properties of a parent excitation space, which converge toward the properties of a target excitation space. For the CPS(D-n) model, perturbative corrections through order n are added to the coupled cluster singles (CCS) excitation energies to target the coupled cluster singles and doubles (CCSD) excitation energies. Through a comparative analysis of excitation energy calculations across a diverse set of molecules and wavefunction methods, we present a comprehensive evaluation of the accuracy of the third-order CPS(D) model, CPS(D-3), in calculating excitation energies. Further, our findings demonstrate that CPS(D-3) is a reliable alternative to established methods, particularly CCSD, while systematically overestimating the excitation energies compared to high-level coupled cluster methods such as CC3. These results highlight the strengths and limitations of CPS(D-3), as well as the promising directions for its future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗

DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗