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,045 records · Page 58

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗

Environmental Conditions Affecting Global Mesoscale Convective System Occurrence

Abstract The ERA5 environments of mesoscale convective systems (MCSs), tracked from satellite observations, are assessed over a 20-yr period. The use of a large set of MCS tracks allows us to robustly test the sensitivity of the results to factors such as region, latitude, and diurnal cycle. We aim to provide novel information on environments of observed MCSs for assessments of global atmospheric models and to improve their ability to simulate MCSs. Statistical analysis of all tracked MCSs is performed in two complementary ways. First, we investigate the environments when an MCS has occurred at different spatial scales before and after MCS formation. Several environmental variables are found to show marked changes before MCS initiation, particularly over land. The vertically integrated moisture flux convergence shows a robust signal across different regions and when considering MCS initiation diurnal cycle. We also found spatial scale dependence of the environments between 200 and 500 km, providing new evidence of a natural length scale for use with MCS parameterization. In the second analysis, the likelihood of MCS occurrence for given environmental conditions is evaluated, by considering all environments and determining the probability of being in an MCS core or shield region. These are compared to analogous non-MCS environments, allowing discrimination between conditions suitable for MCS and non-MCS occurrence. Three environmental variables are found to be useful predictors of MCS occurrence: total column water vapor, midlevel relative humidity, and total column moisture flux convergence. Such relations could be used as trigger conditions for the parameterization of MCSs, thereby strengthening the dependence of the MCS scheme on the environment. Significance Statement Large storm systems called mesoscale convective systems form across Earth. These are collections of thunderstorms, with associated high-level clouds that produce substantial, lighter rainfall and modulate Earth’s energy balance. They produce hazardous weather conditions, such as floods and high winds, and are responsible for a high percentage of rainfall in many regions globally. We investigate the environmental conditions under which they form, so that we can understand the spatial extent of the environment which is important for their formation, and also where and when the effects of these storms might be felt. The novel information generated here should help improve the representation of these storms in weather and climate models, improving the prediction of rainfall, thunderclouds, and high-level clouds.

54 ENVIRONMENTAL SCIENCES↗

A space–angle DGFEM approach for the Boltzmann radiation transport equation with local angular refinement

In this paper a new method for the discretization of the radiation transport equation is presented, based on a discontinuous Galerkin method in space and angle that allows for local refinement in angle where any spatial element can support its own angular discretization. To cope with the discontinuous spatial nature of the solution, a generalized Riemann procedure is required to distinguish between incoming and outgoing contributions of the numerical fluxes. A new consistent framework is introduced that is based on the solution of a generalized eigenvalue problem. The resulting numerical fluxes for the various possible cases where neighboring elements have an equal, higher or lower level of refinement in angle are derived based on tensor algebra and the resulting expressions have a very clear physical interpretation. The choice of discontinuous trial functions not only has the advantage of easing local refinement, it also facilitates the use of efficient sweep-based solvers due to decoupling of unknowns on a large scale thereby approaching the efficiency of discrete ordinates methods with local angular resolution. The approach is illustrated by a series of numerical experiments. Results show high orders of convergence for the scalar flux on angular refinement. The generalized Riemann upwinding procedure leads to stable and consistent solutions. Further the sweep-based solver performs well when used as a preconditioner for a Krylov method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cooling rate and size effects on the medium-range structure of multicomponent oxide glasses simulated by molecular dynamics

A set of molecular dynamics simulations were performed to investigate the effect of cooling rate and system size on the medium-range structure of melt-derived multicomponent silicate glasses, represented by the quaternary 45S5 Bioglass composition. Given the significant impact of the glass degradation on applications of these materials in biomedicine and nuclear waste disposal, bulk structural features which directly affect the glass dissolution process are of particular interest. Connectivity of the silicate matrix, ion clustering and nanosegregation, distribution of ring and chain structural patterns represent critical features in this context, which can be directly extracted from the models. A key issue is represented by the effect of the computational approach on the corresponding glass models, especially in light of recent indications questioning the suitability of conventional MD approaches (that is, involving melt-and-quench of systems containing ∼10{sup 3} atoms at cooling rates of 5-10 K/ps) when applied to model these glasses. The analysis presented here compares MD models obtained with conventional and nonconventional cooling rates and system sizes, highlighting the trend and range of convergence of specific structural features in the medium range. The present results show that time-consuming computational approaches involving much lower cooling rates and/or significantly larger system sizes are in most cases not necessary in order to obtain a reliable description of the medium-range structure of multicomponent glasses. We identify the convergence range for specific properties and use them to discuss models of several glass compositions for which a possible influence of cooling-rate or size effects had been previously hypothesized. The trends highlighted here represent an important reference to obtain reliable models of multicomponent glasses and extract converged medium-range structural features which affect the glass degradation and thus their application in different fields. In addition, as a first application of the present findings, the fully converged structure of the 45S5 glass was further analyzed to shed new light on several dissolution-related features whose interpretation has been rather controversial in the past.

36 MATERIALS SCIENCE↗

Calculation of Thomson scattering spectral fits for interpenetrating flows

Collective mode optical Thomson scattering has been used to investigate the interactions of radially convergent ablation flows in Tungsten wire arrays. These experiments were carried out at the Magpie pulsed power facility at Imperial College, London. Analysis of the scattered spectra has provided direct evidence of ablation stream interpenetration on the array axis, and has also revealed a previously unobserved axial deflection of the ablation streams towards the anode as they approach the axis. It is has been suggested that this deflection is caused by the presence of a static magnetic field, advected with the ablation streams, stagnated and accrued around the axis. Analysis of the Thomson scattering spectra involved the calculation and fitting of the multi-component, non-relativistic, Maxwellian spectral density function S (k, ω). The method used to calculate the fits of the data are discussed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HALF-BROTHERS IN THE SCHULHOF FAMILY?

Disruptive collisions in the asteroid belt produced groups of fragments known as the asteroid families. The studies of identified asteroid families help us to better understand issues related to impact physics, space weathering, asteroid interior, and collisional evolution of the main belt. Here, we analyze a family near the central main belt asteroid (2384) Schulhof. We show that the previously found group of objects around (81337) 2000 GP36 is actually a sub-cluster in the larger Schulhof family. Using backward integrations we demonstrate that the orbits of sub-cluster asteroids converge to that of (2384) Schulhof at 780 {+-} 100 kyr ago, suggesting that the breakup event happened very recently. Interestingly, a similar analysis of the two newly discovered members of the Schulhof family may indicate a second event {approx}< 100 kyr ago (e.g., secondary collision, fission, satellite instability). If confirmed, the formation history of the Schulhof family would suggest that small asteroids may have very colorful lifetimes. Additional astrometric observations of the two new member asteroids will be needed to improve their present orbit and better constrain their past histories.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhancing DESI DR1 full-shape analyses using HOD-informed priors

We present an analysis of DESI Data Release 1 (DR1) that incorporates Halo Occupation Distribution (HOD)-informed priors into Full-Shape (FS) modeling of the power spectrum based on cosmological perturbation theory (PT). By leveraging physical insights from the galaxy-halo connection, these HOD-informed priors on nuisance parameters substantially mitigate projection effects in extended cosmological models that allow for dynamical dark energy. The resulting credible intervals now encompass the posterior maximum from the baseline analysis using gaussian priors, eliminating a significant posterior shift observed in baseline studies. In the ΛCDM framework, a combined DESI DR1 FS information and constraints from the DESI DR1 baryon acoustic oscillations (BAO) — including Big Bang Nucleosynthesis (BBN) constraints and a weak prior on the scalar spectral index — yields Ω m = 0.2994 ± 0.0090 and σ 8 = 0.836$^{+0.024}_{-0.027}$, representing improvements of approximately 4% and 23% over the baseline analysis, respectively. For the w 0 w a CDM model, our results from various data combinations are highly consistent, with all configurations converging to a region with w 0 > -1 and w a < 0. This convergence not only suggests intriguing hints of dynamical dark energy but also underscores the robustness of our HOD-informed prior approach in delivering reliable cosmological constraints.

59 BASIC BIOLOGICAL SCIENCES↗

Residential energy demand, emissions, and expenditures at regional and income-decile level for alternative futures

Income and its distribution profile are important determinants of residential energy demand and carry direct implications for human well-being and climate. We explore the sensitivity of residential energy systems to income growth and distribution across shared socioeconomic pathway-representative concentration pathways scenarios using a global, integrated, multisector dynamics model, Global Change Analysis Model, which tracks national/regional household energy services and fuel choice by income decile. Nation/region energy use patterns across deciles tend to converge over time with aggregate income growth, as higher-income consumers approach satiation levels in floorspace and energy services. However, in some regions, existing within-region inequalities in energy consumption persist over time due to slow income growth in lower income groups. Due to continued differences in fuel types, lower income groups will have higher exposure to household air pollution, despite lower contributions to greenhouse gas emissions. We also find that the share of income dedicated to energy is higher for lower deciles, with strong regional differences.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Cosmology from weak lensing, galaxy clustering, CMB lensing, and tSZ – I. 10 × 2pt modelling methodology

ABSTRACT The overlap of galaxy surveys and cosmic microwave background (CMB) experiments presents an ideal opportunity for joint cosmological data set analyses. In this paper we develop a halo model-based method for the first joint analysis combining these two experiments using 10 correlated two-point functions (10 × 2pt) derived from galaxy position, galaxy shear, CMB lensing convergence, and Compton-y fields. We explore this method using the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) as examples. We find such LSS × CMB joint analyses lead to significant improvement in Figure-of-Merit of Ωm and S8 over the constraints from using LSS-only probes within Λ cold dark matter (ΛCDM) model. We identify that the shear–y and y–y correlations are the most valuable additions when thermal Sunyaev–Zel’dolvich (tSZ) is included. We further identify the dominant sources of halo model uncertainties in the small-scale modelling, and investigate the impact of halo self-calibration due to the inclusion of small-scale tSZ information.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct determination of one-dimensional interphase structures using normalized crystal truncation rod analysis

Here, a one-dimensional non-iterative direct method was employed for normalized crystal truncation rod analysis. The non-iterative approach, utilizing the Kramers–Kronig relation, avoids the ambiguities due to an improper initial model or incomplete convergence in the conventional iterative methods. The validity and limitations of the present method are demonstrated through both numerical simulations and experiments with Pt(111) in a 0.1 M CsF aqueous solution. The present method is compared with conventional iterative phase-retrieval methods.

36 MATERIALS SCIENCE↗

Final Scientific Report for "The Interhemispheric Pattern in 20th Century and Future Abrupt Change in Regional Tropical Rainfall"

This is the final scientific report for grant DOE-FG02-08ER64588, "The Interhemispheric Pattern in 20th Century and Future Abrupt Change in Regional Tropical Rainfall."The project investigates the role of the interhemispheric pattern in surface temperature – i.e. the contrast between the northern and southern temperature changes – in driving rapid changes to tropical rainfall changes over the 20th century and future climates. Previous observational and modeling studies have shown that the tropical rainband – the Intertropical Convergence Zone (ITCZ) over marine regions, and the summer monsoonal rainfall over land – are sensitive to the interhemispheric thermal contrast; but that the link between the two has not been applied to interpreting long-term tropical rainfall changes over the 20th century and future.The specific goals of the project were to i) develop dynamical mechanisms to explain the link between the interhemispheric pattern to abrupt changes of West African and Asian monsoonal rainfall; ii) Undertake a formal detection and attribution study on the interhemispheric pattern in 20th century climate; and iii) assess the likelihood of changes to this pattern in the future. In line with these goals, our project has produced the following significant results: 1.We have developed a case that suggests that the well-known abrupt weakening of the West African monsoon in the late 1960s was part of a wider co-ordinated weakening of the West African and Asian monsoons, and driven from an abrupt cooling in the high latitude North Atlantic sea surface temperature at the same time. Our modeling work suggests that the high-latitude North Atlantic cooling is effective in driving monsoonal weakening, through driving a cooling of the Northern hemisphere that is amplified by positive radiative feedbacks. 2.We have shown that anthropogenic sulfate aerosols may have partially contributed to driving a progressively southward displacement of the Atlantic Intertropical Convergence Zone (ITCZ) over the course of the 20th century prior to the 1980s. This is based on our detection and attribution analysis of 20th century simulations done by international modeling groups as part of the Coupled Model Intercomparison Project phase 3 (CMIP3). We repeated the same analysis with the current CMIP5 multimodel simulations, with essentially similar results. 3.Future projections of the global interhemispheric thermal gradient suggest a pronounced trend that well exceeds the 20th century range of behavior. The major cause of this trend is due to anthropogenic greenhouse gas emissions, acting in such a way as to warm the North more than the South. This result is based on our analysis of the CMIP3 and 5 simulations of future scenarios. The underlying suggestion is that tropical rainfall may concentrate more northwards in the future climate, though further research is required to more firmly establish that result.Taken together, our results shows the important role of the interhemispheric thermal gradient in determining tropical rainfall changes in the 20th century and future. Our analysis specifically highlights high-latitude North Atlantic sea surface temperature, and anthropogenic sulfate aerosols, as important drivers of the interhemispheric gradient over the 20th century; and anthropogenic greenhouse gases in the 21st. The PI has written a review paper in order to promote the awareness of the interhemispheric gradient amongst the climate science community.Our project was instrumental in developing the career of a postdoctoral scholar, as well as contributing to the research training of three Ph.D. candidates.

54 ENVIRONMENTAL SCIENCES↗

Final Report on Subcontract B605152. Multigrid Methods for Systems of PDEs

The project team has continued with work on developing aggressive coarsening techniques for AMG methods. Of particular interest is the idea to use aggressive coarsening with polynomial smoothing. Using local Fourier analysis the optimal values for the parameters involved in defining the polynomial smoothers are determined automatically in a way to achieve fast convergence of cycles with aggressive coarsening. Numerical tests have the sharpness of the theoretical results. The methods are highly parallelizable and efficient multigrid algorithms on structured and semistructured grids in two and three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗

Economic and Physical Linkages of the Information and Communication Technology (ICT) Service Industry to Key Industries of the Economy: An Ad Hoc Report

This report examines the information and communications technology (ICT) services industry in response to an inquiry by the Department of Energy’s (DOE’s) Office of Energy Policy and Systems Analysis. The report answers several key questions: •How has the reliance on ICT services evolved in recent years for key infrastructure services such as air travel, freight transport, electricity and natural gas distribution, financial services, and critical health care, and for the household sector? •What ICT industry trends explain continued strong linkage to and reliance upon ICT? •What is the ICT industry’s reliance on grid-supplied power, uninterruptible power supplies, emergency generators and back-up energy storage technologies? •What are the observed direct effects of ICT disruptions induced by electrical system failures in recent history and how resilient are the components of the ICT industry?

99 GENERAL AND MISCELLANEOUS↗

2021 GeoAI Workshop Report: The Trillion Pixel Challenge

The convergence of geospatial big data with advancements from artificial intelligence, cloud infrastructure, and high-performance computing continues to revolutionize mapping and analysis of Earth's surface in unprecedented detail. Rapid innovations in sensing technologies will soon collect geospatial data in even higher resolution and throughput. These developments offer the potential for breakthroughs in science, policy, and national security via end-to-end GeoAI systems that can provide fresh insights into how humans occupy and alter their environment over time. At the 2021 GeoAI Trillion Pixel workshop, international subject matter experts from government, academia, industry, and nonprofit organizations gathered virtually to discuss the Trillion Pixel GeoAI Challenge. The event focused on six major themes currently influencing scientific innovation and breakthroughs. Particular focus was paid to societal impacts. As an additional takeaway message, the gathering identified remaining application gaps and challenges that are in need of stronger community partnerships and collaborations.

58 GEOSCIENCES↗

Performance of a transpiration-regenerative cooled rocket thrust chamber

The analysis, design, fabrication, and testing of a liquid rocket engine thrust chamber which is gas transpiration cooled in the high heat flux convergent portion of the chamber and water jacket cooled (simulated regenerative) in the barrel and divergent sections of the chamber are described. The engine burns LOX-hydrogen propellants at a chamber pressure of 600 psia. Various transpiration coolant flow rates were tested with resultant local hot gas wall temperatures in the 800 F to 1400 F range. The feasibility of transpiration cooling with hydrogen and helium, and the use of photo-etched copper platelets for heat transfer and coolant metering was successfully demonstrated.

Valler, H. W.↗

A General-Purpose Optimization Engine for Multi-Disciplinary Design Applications

A general purpose optimization tool for multidisciplinary applications, which in the literature is known as COMETBOARDS, is being developed at NASA Lewis Research Center. The modular organization of COMETBOARDS includes several analyzers and state-of-the-art optimization algorithms along with their cascading strategy. The code structure allows quick integration of new analyzers and optimizers. The COMETBOARDS code reads input information from a number of data files, formulates a design as a set of multidisciplinary nonlinear programming problems, and then solves the resulting problems. COMETBOARDS can be used to solve a large problem which can be defined through multiple disciplines, each of which can be further broken down into several subproblems. Alternatively, a small portion of a large problem can be optimized in an effort to improve an existing system. Some of the other unique features of COMETBOARDS include design variable formulation, constraint formulation, subproblem coupling strategy, global scaling technique, analysis approximation, use of either sequential or parallel computational modes, and so forth. The special features and unique strengths of COMETBOARDS assist convergence and reduce the amount of CPU time used to solve the difficult optimization problems of aerospace industries. COMETBOARDS has been successfully used to solve a number of problems, including structural design of space station components, design of nozzle components of an air-breathing engine, configuration design of subsonic and supersonic aircraft, mixed flow turbofan engines, wave rotor topped engines, and so forth. This paper introduces the COMETBOARDS design tool and its versatility, which is illustrated by citing examples from structures, aircraft design, and air-breathing propulsion engine design.

Patnaik, Surya N.↗

Parallel Aircraft Trajectory Optimization with Analytic Derivatives

Trajectory optimization is an integral component for the design of aerospace vehicles, but emerging aircraft technologies have introduced new demands on trajectory analysis that current tools are not well suited to address. Designing aircraft with technologies such as hybrid electric propulsion and morphing wings requires consideration of the operational behavior as well as the physical design characteristics of the aircraft. The addition of operational variables can dramatically increase the number of design variables which motivates the use of gradient based optimization with analytic derivatives to solve the larger optimization problems. In this work we develop an aircraft trajectory analysis tool using a Legendre-Gauss-Lobatto based collocation scheme, providing analytic derivatives via the OpenMDAO multidisciplinary optimization framework. This collocation method uses an implicit time integration scheme that provides a high degree of sparsity and thus several potential options for parallelization. The performance of the new implementation was investigated via a series of single and multi-trajectory optimizations using a combination of parallel computing and constraint aggregation. The computational performance results show that in order to take full advantage of the sparsity in the problem it is vital to parallelize both the non-linear analysis evaluations and the derivative computations themselves. The constraint aggregation results showed a significant numerical challenge due to difficulty in achieving tight convergence tolerances. Overall, the results demonstrate the value of applying analytic derivatives to trajectory optimization problems and lay the foundation for future application of this collocation based method to the design of aircraft with where operational scheduling of technologies is key to achieving good performance.

aircraft↗

A thermodynamically consistent discontinuous Galerkin formulation for interface separation

Our paper describes the formulation of an interface damage model, based on the discontinuous Galerkin (DG) method, for the simulation of failure and crack propagation in laminated structures. The DG formulation avoids common difficulties associated with cohesive elements. Specifically, it does not introduce any artificial interfacial compliance and, in explicit dynamic analysis, it leads to a stable time increment size which is unaffected by the presence of stiff massless interfaces. This proposed method is implemented in a finite element setting. Convergence and accuracy are demonstrated in Mode I and mixed-mode delamination in both static and dynamic analyses. Significantly, numerical results obtained using the proposed interface model are found to be independent of the value of the penalty factor that characterizes the DG formulation. By contrast, numerical results obtained using a classical cohesive method are found to be dependent on the cohesive penalty stiffnesses. The proposed approach is shown to yield more accurate predictions pertaining to crack propagation under mixed-mode fracture because of the advantage. Furthermore, in explicit dynamic analysis, the stable time increment size calculated with the proposed method is found to be an order of magnitude larger than the maximum allowable value for classical cohesive elements.

42 ENGINEERING↗