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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

An evaluation of higher-order modal methods for calculating transient structural response

The present study evaluates a higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1990

Theoretical study of the C-H bond dissociation energy of C2H

A theoretical study of the convergence of the C-H bond dissociation energy D(0) in C2H with respect to both the one- and n-particle spaces is presented. The calculated C-H bond energies of C2H2 and C2H4, which are in excellent agreement with experiment, are used for calibration. The best estimate for D(0) of 112.4 + or - 2.0 kcal/mol is slightly below the recent experimental value of 116.3 + or - 2.6 kcal/mol, but substantially above a previous theoretical estimate of 102 kcal/mol. The remaining discrepancy with experiment may reflect primarily the uncertainty in the experimental D(0) value of C2 required in the analysis.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2026

Planning for Seattle at the Convergence of Resilience and Energy

Grid reliability metrics obscure important temporal, spatial, and categorical considerations for increasing energy resilience. Systemwide or feeder-level outage metrics do not identify which kinds of services are affected by outages, where, and for how long. These outage metrics indicate the impacts of outages but cannot measure the consequences to customers that could result from those outages. The consequences of power outages for surrounding community members are caused by disruptions to electricity-dependent critical services, rather than to electricity itself. Power outages can decrease a community's access to healthcare, fuel, safe indoor temperatures, and provisions like food and water. This project developed critical service access, a new consequence-focused resilience metric that quantifies the relative access to critical services provided to households by distribution infrastructure during normal conditions and major disruptions. We use a spatially granular grid analysis that facilitates targeted resilience interventions; dividing feeders into isolatable sections and combining those sections with the critical service access metric allows us to identify where energy improvements like solar-plus-storage microgrids could create the most benefits for community members by protecting access to food, fuel, health, shelter, and public safety services. We identify locations in a South Seattle study area that could be high priorities for resilience investment and summarize their potential neighborhood-scale resilience benefits. This analysis was complemented by direct feedback from study area residents collect through a survey and focus groups. Results can help utilities understand how and where long power outages can create real consequences for customers, set strategic targets based on that understanding, and measure the potential benefits of energy resilience upgrades for more informed decisions.

14 SOLAR ENERGY↗
NASA NTRS1989

Stability analysis of multiple-robot control systems

In a space telerobotic service scenario, cooperative motion and force control of multiple robot arms are of fundamental importance. Three paradigms to study this problem are proposed. They are distinguished by the set of variables used for control design. They are joint torques, arm tip force vectors, and an accelerated generalized coordinate set. Control issues related to each case are discussed. The latter two choices require complete model information, which presents practical modeling, computational, and robustness problems. Therefore, focus is on the joint torque control case to develop relatively model independent motion and internal force control laws. The rigid body assumption allows the motion and force control problems to be independently addressed. By using an energy motivated Lyapunov function, a simple proportional derivative plus gravity compensation type of motion control law is always shown to be stabilizing. The asymptotic convergence of the tracing error to zero requires the use of a generalized coordinate with the contact constraints taken into account. If a non-generalized coordinate is used, only convergence to a steady state manifold can be concluded. For the force control, both feedforward and feedback schemes are analyzed. The feedback control, if proper care has been taken, exhibits better robustness and transient performance.

Wen, John T.↗
NASA NTRS1993

Analysis of the anisotropic viscoplastic-damage response of composite laminates - Continuum basis and computational algorithms

The mathematical structure underlying the rate equations of a recently-developed constitutive model for the coupled viscoplastic-damage response of anisotropic composites is critically examined. In this regard, a number of tensor projection operators have been identified, and their properties were exploited to enable the development of a general computational framework for their numerical implementation using the Euler fully-implicit integration method. In particular, this facilitated (i) the derivation of explicit expressions of the (consistent) material tangent stiffnesses that are valid for both three-dimensional as well as subspace (e.g. plane stress) formulations, (ii) the implications of the symmetry or unsymmetry properties of these tangent operators from a thermodynamic standpoint, and (iii) the development of an effective time-step control strategy to ensure accuracy and convergence of the solution. In addition, the special limiting case of inviscid elastoplasticity is treated. The results of several numerical simulations are given to demonstrate the effectiveness of the schemes developed.

Saleeb, A. F.↗
DOE OSTI2024

Extreme Precipitation Over the Southern Slope of the Tibetan Plateau and the Associated Atmospheric Circulation Anomalies

The southern slope of the Tibetan Plateau (SSTP) is one of the rainiest regions in the world where geological hazards caused by extreme precipitation often occur. This study investigates the characteristics and mechanisms of extreme precipitation over SSTP from June to September during 2001–2020 using Global Precipitation Measurement satellite observation. The extreme precipitation days are defined as the days with top 5% of regional-mean daily precipitation over SSTP in this period, which has an average precipitation of 27.2 mm/d. Averaging over the extreme precipitation days, precipitation peaks at an altitude of about 300 m, coinciding with the climatological maximum precipitation, but with a much larger value of 37.2 mm/d than the climatology of 11.9 mm/d. Composite analysis of circulations on extreme days reveals significant circulation anomalies in both the lower and upper troposphere. Specifically, the lower-tropospheric circulations are characterized by significant westerly anomalies over northern India, and the upper-tropospheric circulations are characterized by northerly anomalies over the central Tibetan Plateau, which are statistically independent. The lower-tropospheric westerly anomalies blowing toward SSTP are blocked by the topography, favoring extreme precipitation over SSTP. The upper-tropospheric northerly anomalies, on the other hand, correspond to anomalous northeasterlies north of SSTP in the middle troposphere and southeasterlies to the south in the lower troposphere, and the convergence of these circulation anomalies favors extreme precipitation over SSTP. Lastly, the lower-tropospheric westerly and upper-tropospheric northerly anomalies, respectively, correspond to less precipitation over the South Asian monsoon region and the Tibetan Plateau.

Extreme precipitation, Tibetan Plateau↗
DOE OSTI2026

Wavelet flow for extragalactic foreground simulations

Extragalactic foregrounds in cosmic microwave background (CMB) observations are both a source of cosmological and astrophysical information and a nuisance to the CMB. Effective field-level modeling that captures their non-Gaussian statistical distributions is increasingly important for optimal information extraction, particularly given the low-noise observations from current and upcoming experiments. Here, we explore the use of Wavelet Flow (WF) models to tackle the novel task of modeling the field-level probability distributions of multi-component CMB secondaries and foregrounds. Specifically, we jointly train correlated CMB lensing convergence (κ) and cosmic infrared background (CIB) maps with a WF model and obtain a network that statistically recovers the input to high accuracy — the trained network generates samples of κ and CIB fields whose average power spectra are within a few percent of the inputs across all scales, and whose Minkowski functionals are similarly accurate compared to the inputs. Leveraging the multiscale architecture of these models, we fine-tune both the model parameters and the priors at each scale independently, optimizing performance across different resolutions. These results demonstrate that WF models can accurately simulate correlated components of CMB secondaries, supporting improved analysis of cosmological data. Our code and trained models can be found on this GitHub repo.

cosmological simulations↗
DOE OSTI2023

Synthetic light-cone catalogues of modern redshift and weak lensing surveys with abacussummit

ABSTRACT The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy, and dark matter. However, the development of high-fidelity theoretical models is a necessary stepping stone. Here, we present public high-resolution weak lensing maps on the light-cone, generated using the N-body simulation suite abacussummit, and accompanying weak lensing mock catalogues, tuned to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument. Available in this release are maps of the cosmic shear, deflection angle, and convergence fields at source redshifts ranging from z = 0.15 to 2.45 as well as cosmic microwave background convergence maps for each of the 25 base-resolution simulations ($L_{\rm box} = 2000\, h^{-1}\, {\rm Mpc}$ and Npart = 69123) as well as for the two huge simulations ($L_{\rm box} = 7500\, h^{-1}\, {\rm Mpc}$ and Npart = 86403) at the fiducial abacussummit cosmology. The pixel resolution of each map is 0.21 arcmin, corresponding to a healpix Nside of 16 384. The sky coverage of the base simulations is an octant until z ≈ 0.8 (decreasing to about 1800 deg2 at z ≈ 2.4), whereas the huge simulations offer full-sky coverage until z ≈ 2.2. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper Suprime-Cam data sets. The mock catalogues are validated against theoretical predictions for various clustering and lensing statistics, such as correlation multipoles, galaxy–shear, and shear–shear, showing excellent agreement. All products can be downloaded via a Globus endpoint (see Data Availability section).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Vergence-dependent adaptation of the vestibulo-ocular reflex

The gain of the vestibulo-ocular reflex (VOR) normally depends on the distance between the subject and the visual target, but it remains uncertain whether vergence angle can be linked to changes in VOR gain through a process of context-dependent adaptation. In this study, we examined this question with an adaptation paradigm that modified the normal relationship between vergence angle and retinal image motion. Subjects were rotated sinusoidally while they viewed an optokinetic (OKN) stimulus through either diverging or converging prisms. In three subjects the diverging prisms were worn while the OKN stimulus moved out of phase with the head, and the converging prisms were worn when the OKN stimulus moved in-phase with the head. The relationship between the vergence angle and OKN stimulus was reversed in the fourth subject. After 2 h of training, the VOR gain at the two vergence angles changed significantly in all of the subjects, evidenced by the two different VOR gains that could be immediately accessed by switching between the diverged and converged conditions. The results demonstrate that subjects can learn to use vergence angle as the contextual cue that retrieves adaptive changes in the angular VOR.

NASA Discipline Developmental Biology↗
NASA NTRS1988

Computational aspects of sensitivity calculations in transient structural analysis

A key step in the application of formal automated design techniques to structures under transient loading is the calculation of sensitivities of response quantities to the design parameters. This paper considers structures with general forms of damping acted on by general transient loading and addresses issues of computational errors and computational efficiency. The equations of motion are reduced using the traditional basis of vibration modes and then integrated using a highly accurate, explicit integration technique. A critical point constraint formulation is used to place constraints on the magnitude of each response quantity as a function of time. Three different techniques for calculating sensitivities of the critical point constraints are presented. The first two are based on the straightforward application of the forward and central difference operators, respectively. The third is based on explicit differentiation of the equations of motion. Condition errors, finite difference truncation errors, and modal convergence errors for the three techniques are compared by applying them to a simple five-span-beam problem. Sensitivity results are presented for two different transient loading conditions and for both damped and undamped cases.

Greene, William H.↗
NASA NTRS1987

An evaluation of higher-order model methods for calculating transient structural response

A higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method is evaluated. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force-derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1987

A general method for dynamic analysis of structures overview

The presented research deals with the development of a dynamic analysis method for structural systems. The modeling approach is essentially a finite element method in the sense that the structure is divided into n elements. An element is defined as any structural unit whose degree of freedom (dofs) can be categorized as either interface or non-interface dofs. An element could be a fundamental unit such as a rod, a beam, a plate etc., or it could be an entire structural component. Furthermore, the parameters for the element could be distributed or lumped. The choice of elements is totally arbitrary and is a matter of user convenience. In particular, issues of accuracy and convergence do not enter on the level of example that bookkeeping is reduced to a minimum. Each element is modeled using a set of interface constraint modes (ICM) combined with a set of interface restrained normal models (IRNM). The next step is the solution of the system eigenvalue problem. The procedure calls for the sequential solution of a number of small eigenvalue problems based on a truncation principle for IRNM. In addition, the form of these eigenvalue problems is very simple such that an escalator type of eigenvalue problem solver can be used which is extremely cost-effective and fast.

Engels, Remi C.↗
NASA NTRS2007

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

Mesoscale weather conditions can significantly affect the space launch and landing operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). During the summer months, land-sea interactions that occur across KSC and CCAFS lead to the formation of a sea breeze, which can then spawn deep convection. These convective processes often last 60 minutes or less and pose a significant challenge to the forecasters at the National Weather Service (NWS) Spaceflight Meteorology Group (SMG). The main challenge is that a "GO" forecast for thunderstorms and precipitation at the Shuttle Landing Facility is required at the 90 minute deorbit decision for End Of Mission (EOM) and at the 30 minute Return To Launch Site (RTLS) decision. Convective initiation, timing, and mode also present a forecast challenge for the NWS in Melbourne, FL (MLB). The NWS MLB issues such tactical forecast information as Terminal Aerodrome Forecasts (TAF5), Spot Forecasts for fire weather and hazardous materials incident support, and severe/hazardous weather Watches, Warnings, and Advisories. Lastly, these forecasting challenges can also affect the 45th Weather Squadron (45 WS), which provides comprehensive weather forecasts for shuttle launch, as well as ground operations, at KSC and CCAFS. The need for accurate mesoscale model forecasts to aid in their decision making is crucial. This study specifically addresses the skill of different model configurations in forecasting warm season convective initiation. Numerous factors influence the development of convection over the Florida peninsula. These factors include sea breezes, river and lake breezes, the prevailing low-level flow, and convergent flow due to convex coastlines that enhance the sea breeze. The interaction of these processes produces the warm season convective patterns seen over the Florida peninsula. However, warm season convection remains one of the most poorly forecast meteorological parameters. To determine which configuration options are best to address this specific forecast concern, the Weather Research and Forecasting (WRF) model, which has two dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM) was employed. In addition to the two dynamical cores, there are also two options for a "hot-start" initialization of the WRF model - the Local Analysis and Prediction System (LAPS; McGinley 1995) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS; Brewster 1996). Both LAPS and ADAS are 3- dimensional weather analysis systems that integrate multiple meteorological data sources into one consistent analysis over the user's domain of interest. This allows mesoscale models to benefit from the addition of highresolution data sources. Having a series of initialization options and WRF cores, as well as many options within each core, provides SMG and MLB with considerable flexibility as well as challenges. It is the goal of this study to assess the different configurations available and to determine which configuration will best predict warm season convective initiation.

Watson, Leela R.↗
NASA NTRS1990

Heat Pipe Vapor Dynamics

The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.

Issacci, Farrokh↗
NASA NTRS1996

Stereo viewing 3-component, planar PIV utilizing fuzzy inference

An all electronic 3-D Digital Particle Image Velocimetry (DPIV) system has been developed for use in high velocity (supersonic) flows. Two high resolution CCD cameras mounted in a stereo viewing configuration are used to determine the out-of-plane velocity component from the difference of the in-plane velocity measurements. Double exposure image frames are acquired and Fuzzy inference techniques are used to maximize the validity of the velocity estimates obtained from the auto-correlation analysis. The CCD cameras are tilted relative to their respective lens axes to satisfy Scheimpflug's condition. Tilting the camera film plane ensures that the entire image plane is in focus. Perspective distortion still results, but can be corrected by proper calibration of the optical system. A calibration fixture is used to determine the experimental setup parameters and to assess the accuracy to which the z-plane displacements can be estimated. The details of the calibration fixture and procedure are discussed in the text. A pair of pulsed Nd:YAG lasers operating at 532 nm are used to illuminate the seeded flow from a convergent nozzle operated in an underexpanded condition. The light sheet was oriented perpendicular to the nozzle flow, yielding planar cross-sections of the 3-component velocity field at several axial stations. The key features of the supersonic jet are readily observed in the cross-plane vector plots.

Wernet, Mark P.↗
DOE OSTI2022

Natural gas maximal load delivery for multi-contingency analysis

An increasing dependence on natural gas has amplified existing vulnerabilities to the power grid, including disruptions to gas transmission networks from natural and man-made disasters. To address the operational challenges arising from these disruptions, we, in this study, consider the problem of estimating the steady-state operating capacity of a damaged gas pipeline network while ensuring the maximal delivery of load. Specifically, we formulate the mixed-integer nonconvex maximal load delivery (MLD) problem, which proves difficult to solve on large-scale networks. To address this challenge, we present a relaxation of the MLD problem and use it to determine bounds on the transport capacity of a gas pipeline system. A rigorous computational evaluation over network models ranging in size from 11 to 4,197 junctions shows that the relaxation-based method is suitable for analyzing the impacts of multi-contingency network disruptions, often converging to the optimal solution of the relaxation in less than ten seconds.

03 NATURAL GAS↗
DOE OSTI2020

The many-body expansion for aqueous systems revisited: 1. Water-water interactions

We revisit the Many-Body Expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the Basis Set Superposition Error (BSSE) correction, on the various terms for selected sizes of water clusters up to n = 21. The analysis is performed at the second order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality for the (H2O)n, n = 7, 10, 13, 16 and 21 clusters, for which we report either the complete MBE (n = 7, 10) or the one through the 6-body (n = 13) and the 5-body terms (n = 16, 21). Our results suggest that any sizeable contributions to the total cluster binding energy arising from the 5-body and larger terms are solely an artifact of the size of the basis set. Indeed, all terms above the 4-body converge monotonically to practically zero at the Complete Basis Set (CBS) limit. The magnitudes of the 5-body and above terms are accurately reproduced even with the smaller basis set of the series (aug-cc-pVDZ) once the BSSE correction is taken into account. The same level of theory (MP2/aug-cc-pVDZ, BSSEcorrected) also accurately reproduces the MP2/CBS values of the 3- and 4-body terms. The contribution of electron correlation to the 3- and 4-body terms is quite small, so that neglecting the correlation contribution in all terms above the 3-body results in an error of ??(0.1%). The BSSE correction to the largest 2-body term in the MBE was accurately estimated from the function ?? [1 + erf(-?? ??)], which is proportional to the common (overlapping) area between two Gaussian distributions whose centers are separated by R, with the constants ??, ?? fitted to the calculated BSSE corrections for the individual 2-body terms of the clusters with each basis set. Our results unquestionably demonstrate that the MBE for water-water interactions converges to the 4-body term, since any finite terms above the 4-body are artifacts of the size of the basis set and can be safely truncated at the 4-body term when either a very large basis set is used or BSSE corrections are taken into account even with the aug-cc-pVDZ basis set. We expect these findings to have important consequences in the pursuit of accurate many-body molecular dynamics simulations for aqueous systems.

Heindel, Joseph↗
NASA NTRS1988

Weak rotating flow disturbances in a centrifugal compressor with a vaneless diffuser

A theory is presented to predict the occurrence of weak rotating waves in a centrifugal compression system with a vaneless diffuser. As in a previous study of axial systems, an undisturbed performance characteristic is assumed known. Following an inviscid analysis of the diffuser flow, conditions for a neutral rotating disturbance are found. The solution is shown to have two branches; one with fast rotation, the other with very slow rotation. The slow branch includes a dense set of resonant solutions. The resonance is a feature of the diffuser flow, and therefore such disturbances must be expected at the various resonant flow coefficients regardless of the compressor characteristic. Slow solutions seem limited to flow coefficients less than about 0.3, where third and fourth harmonics appear. Fast waves seem limited to a first harmonic. These fast and slow waves are described, and effects of diffuser-wall convergence, backward blade angles, and partial recovery of exit velocity head are assessed.

Moore, F. K.↗