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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A velocity/stress mixed stabilized nodal finite element for elastodynamics: Analysis and computations with strongly and weakly enforced boundary conditions

A new nodal mixed finite element is proposed for the simulation of linear elastodynamics and wave propagation problems in time domain. Our method is based on equal-order interpolation discrete spaces for both the velocity (or displacement) and stress (or strain) tensor variables. The mixed form is derived using either the velocity/stress or velocity/strain pair of unknowns, the latter being instrumental in extensions of the method to nonlinear mechanics. The proposed approach works equally well on hexahedral or tetrahedral grids and, for this reason, it is suitable for time-domain engineering applications in complex geometry. The peculiarity of the proposed approach is the use of the rate form of the stress update equation, which yields a set of governing equations with the structure of a non-dissipative space/time Friedrichs’ system. We complement standard traction boundary conditions for the stress with strongly and weakly enforced boundary conditions for the velocity (or displacement). Weakly enforced boundary conditions are particularly suitable when considering complex geometrical shapes, because they do not require dedicated data structures for the imposition of the boundary degrees of freedom, but, rather, they utilize the structure of the variational formulation. Here we also show how the framework of weakly enforced boundary conditions can be used to develop variational forms for multi-domain simulations of heterogeneous media. A complete analysis including stability and convergence proofs is included, in the case of a space–time variational approach. A series of computational tests are used to demonstrate and verify the performance of the proposed approach.

42 ENGINEERING↗
NASA NTRS1995

An all-at-once reduced Hessian SQP scheme for aerodynamic design optimization

This paper introduces a computational scheme for solving a class of aerodynamic design problems that can be posed as nonlinear equality constrained optimizations. The scheme treats the flow and design variables as independent variables, and solves the constrained optimization problem via reduced Hessian successive quadratic programming. It updates the design and flow variables simultaneously at each iteration and allows flow variables to be infeasible before convergence. The solution of an adjoint flow equation is never needed. In addition, a range space basis is chosen so that in a certain sense the 'cross term' ignored in reduced Hessian SQP methods is minimized. Numerical results for a nozzle design using the quasi-one-dimensional Euler equations show that this scheme is computationally efficient and robust. The computational cost of a typical nozzle design is only a fraction more than that of the corresponding analysis flow calculation. Superlinear convergence is also observed, which agrees with the theoretical properties of this scheme. All optimal solutions are obtained by starting far away from the final solution.

Feng, Dan↗
NASA NTRS2006

High Fidelity System Simulation of Multiple Components in Support of the UEET Program

The High Fidelity System Simulation effort has addressed various important objectives to enable additional capability within the NPSS framework. The scope emphasized High Pressure Turbine and High Pressure Compressor components. Initial effort was directed at developing and validating intermediate fidelity NPSS model using PD geometry and extended to high-fidelity NPSS model by overlaying detailed geometry to validate CFD against rig data. Both "feedforward" and feedback" approaches of analysis zooming was employed to enable system simulation capability in NPSS. These approaches have certain benefits and applicability in terms of specific applications "feedback" zooming allows the flow-up of information from high-fidelity analysis to be used to update the NPSS model results by forcing the NPSS solver to converge to high-fidelity analysis predictions. This apporach is effective in improving the accuracy of the NPSS model; however, it can only be used in circumstances where there is a clear physics-based strategy to flow up the high-fidelity analysis results to update the NPSS system model. "Feed-forward" zooming approach is more broadly useful in terms of enabling detailed analysis at early stages of design for a specified set of critical operating points and using these analysis results to drive design decisions early in the development process.

Plybon, Ronald C.↗
NASA NTRS2001

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗
DOE OSTI2024

Cosmology from cross-correlation of ACT-DR4 CMB lensing and DES-Y3 cosmic shear

Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy weak lensing measured by the Dark Energy Survey (DES) Y3 data. Our baseline analysis uses the CMB convergence map derived from ACT-DR4 and Planck data, where most of the contamination due to the thermal Sunyaev Zel’dovich effect is removed, thus avoiding important systematics in the cross-correlation. In our modelling, we consider the nuisance parameters of the photometric uncertainty, multiplicative shear bias and intrinsic alignment of galaxies. The resulting cross-power spectrum has a signal-to-noise ratio = 7.1 and passes a set of null tests. We use it to infer the amplitude of the fluctuations in the matter distribution (S_8 ≡ σ_8(Ω_m/0.3)^0.5 = 0.782 ± 0.059) with informative but well-motivated priors on the nuisance parameters. We also investigate the validity of these priors by significantly relaxing them and checking the consistency of the resulting posteriors, finding them consistent, albeit only with relatively weak constraints. This cross-correlation measurement will improve significantly with the new ACT-DR6 lensing map and form a key component of the joint 6×2pt analysis between DES and ACT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Runge–Kutta method for solving fuzzy differential equations under generalized differentiability

In this paper, we interpret a fuzzy differential equation by using the strongly generalized differentiability concept. Utilizing the Generalized Characterization Theorem, we investigate the problem of finding a numerical approximation of solutions. The Runge–Kutta approximation method is implemented and its error analysis, which guarantees pointwise convergence, is given. The method applicability is illustrated by solving a linear first-order fuzzy differential equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Structures of drug-specific monoclonal antibodies bound to opioids and nicotine reveal a common mode of binding

Opioid-related fatal overdoses have reached epidemic proportions. Because existing treatments for opioid use disorders offer limited long-term protection, accelerating the development of newer approaches is critical. Monoclonal antibodies (mAbs) are an emerging treatment strategy that targets and sequesters selected opioids in the bloodstream, reducing drug distribution across the blood-brain barrier, thus preventing or reversing opioid toxicity. We previously identified a series of murine mAbs with high affinity and selectivity for oxycodone, morphine, fentanyl, and nicotine. To determine their binding mechanism, we used X-ray crystallography to solve the structures of mAbs bound to their respective targets, to 2.2 Å resolution or higher. Structural analysis showed a critical convergent hydrogen bonding mode that is dependent on a glutamic acid residue in the mAbs’ heavy chain and a tertiary amine of the ligand. Further, characterizing drug-mAb complexes represents a significant step toward rational antibody engineering and future manufacturing activities to support clinical evaluation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Symmetry determination on Pb-free piezoceramic 0.5Ba(Zr{sub 0.2}Ti{sub 0.8})O{sub 3}-0.5(Ba{sub 0.7}Ca{sub 0.3})TiO{sub 3} using convergent beam electron diffraction method

(1−x)(Ba(Zr{sub 0.2}Ti{sub 0.8})O{sub 3}-x(Ba{sub 0.7}Ca{sub 0.3})TiO{sub 3} (BZT-xBCT) Pb-free piezoceramic has been reported showing ultrahigh piezoelectric performance in its morphotropic phase boundary (MPB) region. However, the crystal structure characteristic for the MPB composition of BZT-xBCT is still under debate—between single orthorhombic phase and tetragonal + rhombohedral two phase mixture. In the present study, we perform the local symmetry determination on the MPB composition x = 0.5 using convergent beam electron diffraction analysis (CBED). Our CBED results from multiple zone axes suggest that there are two coexisting phases with the point group symmetries of 4 mm (tetragonal) and 3 m (rhombohedral) respectively, which agree with two phase mixture model. The strong piezoelectricity can thus be understood by considering the polarization rotation between tetragonal and rhombohedral phases by external field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

IMPROVING THE PRECISION OF TIME-DELAY COSMOGRAPHY WITH OBSERVATIONS OF GALAXIES ALONG THE LINE OF SIGHT

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line-of-sight convergence {kappa}{sub ext}. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with {kappa}{sub ext} obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of {kappa}{sub ext} from {sigma}{sub {kappa}} {approx} 0.06 to {approx}0.04 for very overdense LOSs like that of the system B1608+656. For more common LOSs, {sigma}{sub {kappa}} is reduced to {approx}<0.03, corresponding to an uncertainty of {approx}< 3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain {kappa}{sub ext} almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable {kappa}{sub ext} estimator gives {sigma}{sub {kappa}} = 0.047 down from 0.065 using only galaxy counts. Although deeper multiband observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Extratropical Cloud Feedback Constrained by Cloud Sources and Sinks in Cyclones

Constraining cloud feedback in global climate models (GCMs) using observations is important for establishing accurate predictions of future climate. Uncertainty in shortwave cloud feedback (SW FB ) dominates uncertainty in total cloud feedback. Recent studies show a shift toward more positive extratropical SW FB in the latest generations of GCMs leading to the emergence of very high equilibrium climate sensitivity (ECS). In this study, we use precipitation efficiency and albedo susceptibility to constrain liquid water path (LWP) response to warming and SW FB in the Southern Ocean (SO; 50°–80°S). We analyze precipitation in extratropical cyclones (ECs) to learn about extratropical condensed water sink processes, combined with observations of clouds and moisture convergence, and use the analysis to better understand and constrain SW FB . We utilize a perturbed parameter ensemble (PPE) hosted in the Community Atmosphere Model, version 6 (CAM6), to provide a constraint on SW FB based on observations from Clouds and the Earth’s Radiant Energy System (CERES) and Multisensor Advanced Climatology of LWP (MAC-LWP). We apply Gaussian process regression to emulate the model response to all parameters perturbed in the PPE. Confronting the emulator output with observations provides a new estimated response of Earth to global warming. Furthermore, our new estimates of SO LWP reduce the PPE range by 66%–72%, which results in a shortwave cloud radiative effect estimated range that is 27%–34% less than the PPE range. Observations suggest a more positive SO SW FB than the Community Earth System Model, version 2 (CESM2), and consequently do not reject the high climate sensitivity GCMs emerging from the Coupled Model Intercomparison Project phase 6 (CMIP6).

Atmosphere↗
DOE OSTI2022

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗
DOE OSTI2016

DREAMING OF ATMOSPHERES

Here, we introduce the RobERt (Robotic Exoplanet Recognition) algorithm for the classification of exoplanetary emission spectra. Spectral retrieval of exoplanetary atmospheres frequently requires the preselection of molecular/atomic opacities to be defined by the user. In the era of open-source, automated, and self-sufficient retrieval algorithms, manual input should be avoided. User dependent input could, in worst-case scenarios, lead to incomplete models and biases in the retrieval. The RobERt algorithm is based on deep-belief neural (DBN) networks trained to accurately recognize molecular signatures for a wide range of planets, atmospheric thermal profiles, and compositions. Reconstructions of the learned features, also referred to as the “dreams” of the network, indicate good convergence and an accurate representation of molecular features in the DBN. Using these deep neural networks, we work toward retrieval algorithms that themselves understand the nature of the observed spectra, are able to learn from current and past data, and make sensible qualitative preselections of atmospheric opacities to be used for the quantitative stage of the retrieval process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

FIRST OBSERVATIONAL SIGNATURE OF ROTATIONAL DECELERATION IN A MASSIVE, INTERMEDIATE-AGE STAR CLUSTER IN THE MAGELLANIC CLOUDS

While the extended main-sequence turnoffs (eMSTOs) found in almost all 1–2 Gyr old star clusters in the Magellanic Clouds are often explained by postulating extended star formation histories (SFHs), the tight subgiant branches (SGBs) seen in some clusters challenge this popular scenario. Puzzlingly, the SGB of the eMSTO cluster NGC 419 is significantly broader at bluer than at redder colors. We carefully assess and confirm the reality of this observational trend. If we would assume that the widths of the features in color–magnitude space were entirely owing to a range in stellar ages, the SFHs of the eMSTO stars and the blue SGB region would be significantly more prolonged than that of the red part of the SGB. This cannot be explained by assuming an internal age spread. We show that rotational deceleration of a population of rapidly rotating stars, a currently hotly debated alternative scenario, naturally explains the observed trend along the SGB. Our analysis shows that a “converging” SGB could be produced if the cluster is mostly composed of rapidly rotating stars that slow down over time owing to the conservation of angular momentum during their evolutionary expansion from main-sequence turnoff stars to red giants.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Optimization of a corrugated stiffened composite panel under uniaxial compression

An approach of structural optimization has been used to optimize the weight of a simply supported, corrugated hat stiffened composite panel under uniaxial compression. The approach consists of the employment of nonlinear mathematical programming techniques to reach an optimum solution. Some simplifying assumptions are made in the stress analysis to obtain faster convergence to an optimum solution. With these simplifying assumptions the number of unknown design parameters is reduced to twelve.

Agarwal, B. L.↗