Search NASASearch

SEARCH · Search NASA

Results for “structural model decomposition”

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 163 records · Page 9

Interpretable and flexible non-intrusive reduced-order models using reproducing kernel Hilbert spaces

This paper develops an interpretable, non-intrusive reduced-order modeling technique using regularized kernel interpolation. Existing non-intrusive approaches approximate the dynamics of a reduced-order model (ROM) by solving a data-driven least-squares regression problem for low-dimensional matrix operators. Our approach instead leverages regularized kernel interpolation, which yields an optimal approximation of the ROM dynamics from a user-defined reproducing kernel Hilbert space. We show that our kernel-based approach can produce interpretable ROMs whose structure mirrors full-order model structure by embedding judiciously chosen feature maps into the kernel. The approach is flexible and allows a combination of informed structure through feature maps and closure terms via more general nonlinear terms in the kernel. We also derive a computable a posteriori error bound that combines standard error estimates for intrusive projection-based ROMs and kernel interpolants. In conclusion, the approach is demonstrated in several numerical experiments that include comparisons to operator inference using both proper orthogonal decomposition and quadratic manifold dimension reduction.

Data-driven model reduction

GLAD-M35: a joint P and S global tomographic model with uncertainty quantification

We present our third and final generation joint P and S global adjoint tomography (GLAD) model, GLAD-M35, and quantify its uncertainty based on a low-rank approximation of the inverse Hessian. Starting from our second-generation model, GLAD-M25, we added 680 new earthquakes to the database for a total of 2160 events. New P-wave categories are included to compensate for the imbalance between P- and S-wave measurements, and we enhanced the window selection algorithm to include more major-arc phases, providing better constraints on the structure of the deep mantle and more than doubling the number of measurement windows to 40 million. Two stages of a Broyden–Fletcher–Goldfarb–Shanno (BFGS) quasi-Newton inversion were performed, each comprising five iterations. With this BFGS update history, we determine the model’s standard deviation and resolution length through randomized singular value decomposition.

58 GEOSCIENCES

Decomposition-aggregation stability analysis of the spinning Skylab

Stability of an 11-th order linear model of the spinning Skylab is determined by the decomposition-aggregation method based upon the comparison principle and vector Liapunov functions. To reduce the inherent conservativeness of the method an optimization problem is formulated and resolved producing the optimum comparison system. The system provides the best estimate of the stability region of the important structural parameter - asymmetry in the boom settings.

Cuk, S. M.

Measured and calculated acoustic attenuation rates of tuned resonator arrays for two surface impedance distribution models with flow

An experiment was performed to validate two analytical models for predicting low frequency attenuation of duct liner configurations built from an array of seven resonators that could be individually tuned via adjustable cavity depths. These analytical models had previously been developed for high frequency aero-engine inlet duct liner design. In the low frequency application, the liner surface impedance distribution is unavoidably spatially varying by virtue of available fabrication techniques. The characteristic length of this spatial variation may be a significant fraction of the acoustic wavelength. Comparison of measured and predicted attenuation rates and transmission losses for both modal decomposition and finite element propagation models were in good to excellent agreement for a test frequency range that included the first and second cavity resonance frequencies. This was true for either of two surface impedance distribution modeling procedures used to simplify the impedance boundary conditions. In the presence of mean flow, measurements revealed a fine scale structure of acoustic hot spots in the attenuation and phase profiles. These details were accurately predicted by the finite element model. Since no impedance changes due to mean flow were assumed, it is concluded that this fine scale structure was due to convective effects of the mean flow interacting with the surface impedance nonuniformities.

Parrott, Tony L.

A hybrid formulation of a component mode synthesis method

Component mode synthesis is a substructuring technique frequently employed in structural dynamics. In this method, a given structure is subdivided into components or substructures, each of which is analyzed independently for natural frequencies and for mode shapes. The substructure mode shapes are then assembled to give displacement shapes or load patterns of the original structure. An analytical justification of the basic concept is presented using spectral decompositions, and a variant substructuring approach where intersubstructure continuity is enforced in a weak form is derived. This leads to a hybrid formulation of the basic method which is particularly suitable for assembling heterogeneous substructures and analyzing nonconforming and incompatible finite element substructure models. For problems where both the basic and hybrid methods are applicable, the hybrid variant can be computationally more advantageous.

Farhat, Charbel

Task planning with uncertainty for robotic systems

In a practical robotic system, it is important to represent and plan sequences of operations and to be able to choose an efficient sequence from them for a specific task. During the generation and execution of task plans, different kinds of uncertainty may occur and erroneous states need to be handled to ensure the efficiency and reliability of the system. An approach to task representation, planning, and error recovery for robotic systems is demonstrated. Our approach to task planning is based on an AND/OR net representation, which is then mapped to a Petri net representation of all feasible geometric states and associated feasibility criteria for net transitions. Task decomposition of robotic assembly plans based on this representation is performed on the Petri net for robotic assembly tasks, and the inheritance of properties of liveness, safeness, and reversibility at all levels of decomposition are explored. This approach provides a framework for robust execution of tasks through the properties of traceability and viability. Uncertainty in robotic systems are modeled by local fuzzy variables, fuzzy marking variables, and global fuzzy variables which are incorporated in fuzzy Petri nets. Analysis of properties and reasoning about uncertainty are investigated using fuzzy reasoning structures built into the net. Two applications of fuzzy Petri nets, robot task sequence planning and sensor-based error recovery, are explored. In the first application, the search space for feasible and complete task sequences with correct precedence relationships is reduced via the use of global fuzzy variables in reasoning about subgoals. In the second application, sensory verification operations are modeled by mutually exclusive transitions to reason about local and global fuzzy variables on-line and automatically select a retry or an alternative error recovery sequence when errors occur. Task sequencing and task execution with error recovery capability for one and multiple soft components in robotic systems are investigated.

Cao, Tiehua

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.

Application of the Schwinger multichannel formulation to electron-impact excitation of the b 3Sigma(+) state of CO

The Chi 1Sigma - b 3Sigma(+) transition of CO is described using the Schwinger multichannel method in the two-state approximation, in an energy range from threshold to 20 eV. The resonance structure is analyzed by performing a partial-wave decomposition and the resonance positions are established based on the typical discontinuous behavior in the K-matrix elements as well as the pi radian change in the partial-wave eigenphases. The resonance behavior is related to the concept of core-excited resonances and the 'grandparent' model of resonances. The results are related to the formation of the negative ion by carrying out bound-state calculations on the (5 sigma)(3 s sigma)-squared core-excited Rydberg state of CO(-), and the position is found to agree well with the low-energy resonance positions.

Weatherford, Charles A.

PROACTIVE Focus Area 4: Structured Decision Metrics Analysis (Final Report)

A structured decision metrics analysis was developed as one of the tasks under PROACTIVE’s Focus Area 4 (FA4) and used structured decision metrics for processes, items, and facilities (PIF) to determine the efficacy of an M&V system in meeting treaty goals and technical objectives. The method starts with the functional decomposition of a “treaty” with the M&V system overlaid on it. A Functional Decomposition Rubric (FDR) for each PIF is used to determine a score ranging from 0 (weak) to 4 (strong) for assurance, security, and burden. The individual scores are combined to provide an overall Verification System Score (VSS) which assesses the overall verification system’s suitability given goals; scores for each step of the process are also given. The outcome of the evaluation is to identify needs or modifications to the enterprise model, verification system, or proposed testbed capabilities. VeriScore was used to evaluate FA4’s Spiral 0 exercise; those results are included in this report. The VeriScore and FDR framework can also be used in a non-treaty environment, as any goal/objective/method structure will work, thus expanding its applicability beyond traditional arms control structures.

99 GENERAL AND MISCELLANEOUS

Parallel processing in finite element structural analysis

A brief review is made of the fundamental concepts and basic issues of parallel processing. Discussion focuses on parallel numerical algorithms, performance evaluation of machines and algorithms, and parallelism in finite element computations. A computational strategy is proposed for maximizing the degree of parallelism at different levels of the finite element analysis process including: 1) formulation level (through the use of mixed finite element models); 2) analysis level (through additive decomposition of the different arrays in the governing equations into the contributions to a symmetrized response plus correction terms); 3) numerical algorithm level (through the use of operator splitting techniques and application of iterative processes); and 4) implementation level (through the effective combination of vectorization, multitasking and microtasking, whenever available).

Noor, Ahmed K.

Mapping Thermal Conductivity at the Atomic Scale: A Step toward the Thermal Design of Materials

We describe a spatial decomposition of the thermal conductivity, termed site-projected thermal conductivity, which quantifies the thermal conduction activity at each atomic site—a critical parameter for the thermal design of materials. The method is based on the Green–Kubo formula and the harmonic approximation and requires the force-constant and dynamical matrices, as well as a relaxed structural model. Throughout the study, it uses high quality models previously tested and compared to many experiments. It discusses the method and underlying approximations for amorphous silicon, carries the detailed analysis for amorphous silicon, and then examines an amorphous-crystal silicon interface and representative carbon materials. This study identifies the sites and local structures that reduce heat transport, and quantify these (estimate the spatial range) over which these “thermal defects” are effective. It identifies filamentary structures in the amorphous silicon network which impact heat transport and electronic structure (the Urbach edge) and electronic transport.

36 MATERIALS SCIENCE

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS

Hierarchical analysis of spatial pattern and processes of Douglas-fir forests

There has been an increased interest in the quantification of pattern in ecological systems over the past years. This interest is motivated by the desire to construct valid models which extend across many scales. Spatial methods must quantify pattern, discriminate types of pattern, and relate hierarchical phenomena across scales. Wavelet analysis is introduced as a method to identify spatial structure in ecological transect data. The main advantage of the wavelet transform over other methods is its ability to preserve and display hierarchical information while allowing for pattern decomposition. Two applications of wavelet analysis are illustrated, as a means to: (1) quantify known spatial patterns in Douglas-fir forests at several scales, and (2) construct spatially-explicit hypotheses regarding pattern generating mechanisms. Application of the wavelet variance, derived from the wavelet transform, is developed for forest ecosystem analysis to obtain additional insight into spatially-explicit data. Specifically, the resolution capabilities of the wavelet variance are compared to the semi-variogram and Fourier power spectra for the description of spatial data using a set of one-dimensional stationary and non-stationary processes. The wavelet cross-covariance function is derived from the wavelet transform and introduced as a alternative method for the analysis of multivariate spatial data of understory vegetation and canopy in Douglas-fir forests of the western Cascades of Oregon.

Bradshaw, G. A.

Using SpF to Achieve Petascale for Legacy Pseudospectral Applications

Pseudospectral (PS) methods possess a number of characteristics (e.g., efficiency, accuracy, natural boundary conditions) that are extremely desirable for dynamo models. Unfortunately, dynamo models based upon PS methods face a number of daunting challenges, which include exposing additional parallelism, leveraging hardware accelerators, exploiting hybrid parallelism, and improving the scalability of global memory transposes. Although these issues are a concern for most models, solutions for PS methods tend to require far more pervasive changes to underlying data and control structures. Further, improvements in performance in one model are difficult to transfer to other models, resulting in significant duplication of effort across the research community. We have developed an extensible software framework for pseudospectral methods called SpF that is intended to enable extreme scalability and optimal performance. Highlevel abstractions provided by SpF unburden applications of the responsibility of managing domain decomposition and load balance while reducing the changes in code required to adapt to new computing architectures. The key design concept in SpF is that each phase of the numerical calculation is partitioned into disjoint numerical kernels that can be performed entirely inprocessor. The granularity of domain decomposition provided by SpF is only constrained by the datalocality requirements of these kernels. SpF builds on top of optimized vendor libraries for common numerical operations such as transforms, matrix solvers, etc., but can also be configured to use open source alternatives for portability. SpF includes several alternative schemes for global data redistribution and is expected to serve as an ideal testbed for further research into optimal approaches for different network architectures. In this presentation, we will describe our experience in porting legacy pseudospectral models, MoSST and DYNAMO, to use SpF as well as present preliminary performance results provided by the improved scalability.

DYNAMO

Search for UnderUtilized Airspace for Extensible Traffic Management Operations Based on Air Traffic Patterns

This paper presents a new method to facilitate integrating new vehicle traffic operations, operating with existing air traffic operations. The extensible traffic management (xTM) concept assumes that the new vehicles can operate in a dedicated Cooperative Area (CA) with minimal interaction with conventional air traffic and requiring minimal air traffic supervision. Our method assumes that a new xTM CA can be created when an underutilized airspace with little or no traffic can be identified. Our approach involves modeling airspace as a tree data structure and iteratively subdividing it into smaller cells, with underutilized airspace defined as any cells without flight tracks. The benefits of our approach include its applicability to all xTM scenarios, the ability to handle both 2D and 3D space using a unique tree data structure, and computational efficiency for key functions such as space decomposition, labeling of connected cells, and searching of cells containing a given point. By automatically searching for underutilized airspace based on operating air traffic patterns, we can optimize airspace utilization and improve air traffic management. Our proposed approach can quantitatively determine when and where to allow xTM operations in the National Airspace System.

extensible traffic management

Interannual Variability and Trends of Extratropical Ozone: Northern Hemisphere - Part 1

The authors apply principal component analysis (PCA) to the extratropical total column ozone from the combined merged ozone data product and the European Centre for Medium-Range Weather Forecasts assimilated ozone from January 1979 to August 2002. The interannual variability (IAV) of extratropical O-3 in the Northern Hemisphere (NH) is characterized by four main modes. Attributable to dominant dynamical effects, these four modes account for nearly 60% of the total ozone variance in the NH. The patterns of variability are distinctly different from those derived for total O-3 in the tropics. To relate the derived patterns of O-3 to atmospheric dynamics, similar decompositions are performed for the 30 100-Wa geopotential thickness. The results reveal intimate connections between the IAV of total ozone and the atmospheric circulation. The first two leading modes are nearly zonally symmetric and represent the connections to the annular modes and the quasi-biennial oscillation. The other two modes exhibit in-quadrature, wavenumber-1 structures that, when combined, describe the displacement of the polar vortices in response to planetary waves. In the NH, the extrema of these combined modes have preferred locations that suggest fixed topographical and land-sea thermal forcing of the involved planetary waves. Similar spatial patterns and trends in extratropical column ozone are simulated by the Goddard Earth Observation System chemistryclimate model (GEOS-CCM). The decreasing O-3 trend is captured in the first mode. The largest trend occurs at the North Pole, with values similar to-1 Dobson Unit (DU) yr(-1). There is almost no trend in tropical O-3. The trends derived from PCA are confirmed using a completely independent method, empirical mode decomposition, for zonally averaged O-3 data. The O-3 trend is also captured by mode 1 in the GEOS-CCM, but the decrease is substantially larger than that in the real atmosphere.

Yung, Yuk L.

SpF: Enabling Petascale Performance for Pseudospectral Dynamo Models

Pseudospectral (PS) methods possess a number of characteristics (e.g., efficiency, accuracy, natural boundary conditions) that are extremely desirable for dynamo models. Unfortunately, dynamo models based upon PS methods face a number of daunting challenges, which include exposing additional parallelism, leveraging hardware accelerators, exploiting hybrid parallelism, and improving the scalability of global memory transposes. Although these issues are a concern for most models, solutions for PS methods tend to require far more pervasive changes to underlying data and control structures. Further, improvements in performance in one model are difficult to transfer to other models, resulting in significant duplication of effort across the research community.We have developed an extensible software framework for pseudospectral methods called SpF that is intended to enable extreme scalability and optimal performance. High-level abstractions provided by SpF unburden applications of the responsibility of managing domain decomposition and load balance while reducing the changes in code required to adapt to new computing architectures. The key design concept in SpF is that each phase of the numerical calculation is partitioned into disjoint numerical kernels that can be performed entirely in-processor. The granularity of domain-decomposition provided by SpF is only constrained by the data-locality requirements of these kernels. SpF builds on top of optimized vendor libraries for common numerical operations such as transforms, matrix solvers, etc., but can also be configured to use open source alternatives for portability. SpF includes several alternative schemes for global data redistribution and is expected to serve as an ideal testbed for further research into optimal approaches for different network architectures.In this presentation, we will describe the basic architecture of SpF as well as preliminary performance data and experience with adapting legacy dynamo codes. We will conclude with a discussion of planned extensions to SpF that will provide pseudospectral applications with additional flexibility with regard to time integration, linear solvers, and discretization in the radial direction.

Pseudospectral (PS)

Myco-Architecture off Planet: Growing Surface Structures at Destination

Our work focused on filling major select key technical knowledge gaps, and technical aspects to be addressed in a Phase II proposal were identified such as the secretion of glues or plastics by the fungi or bacteria to form biocomposites. Specifically, we raised the technical TRL by assessing the growth of in-house mycelial-producing strains on potential food substrates and analyzed the advantages and disadvantages to their use off planet. Growth of the mycelia on sawdust and powdered nutrients including dried, powdered algae, was used as a baseline. We measured growth as a function of temperature, tested the relevant material properties of the mycelia products produced, and measured the mass of input and output volume for production per volume of material. Mycotecture was produced in a bag simulating the proposed mission implementation. Post-production such as heat treatment (as is done with terrestrial mycotecture)was assessed.The material properties of the dried, frozen and baked mycelial outputs included density, strength,thermal insulation and flame retardation. Tensile/compression testing on biocomposite material specimens was performed to obtain three-dimensional elastic constants, as well as to examine deformation and fracture behavior. These mechanical tests were performed on both dormant and activated samples to understand the evolution of the structural material. In addition, experimental data obtained from the mechanical testing was used to build a failure prediction model that accounts for material anisotropy. This aided in future structural design through a quantitative understanding of the mechanical limits of the material. To examine the texture, defects, fracture surfaces of the material, high-resolution microstructural imaging was utilized before and after the mechanical tests. This microstructural analysis informed us about the macro-structural behavior and influenced structural design. Thermal decomposition analysis was also performed tounderstand the thermal limits.

Biocomposites