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 127 records · Page 7

Synthesis of Correct Digital Controller Models from Specifications by Model Transformation (21-0320)

The design of high consequence controllers (in weapons systems, autonomy, etc.) that do what they are supposed to do is a significant challenge. Testing simply does not come close to meeting the requirements for assurance. Today circuit designers at Sandia (and elsewhere) typically capture the core behavior of their components using state models in tools such as STATEFLOW. They then check that their models meet certain requirements (e.g. “The system bus must not deadlock” or “both traffic lights at an intersection must not be green at the same time”) using tools called model checkers. If the model checker returns “yes” then the property is guaranteed to be satisfied by the model. However, there are several drawbacks to this industry practice: (1) there is a lot of detail to get right, this is particularly challenging when there are multiple components requiring complex coordination (2) any errors returned by the model checker have to be traced back through the design and fixed, necessitating rework, (3) there are severe scalability problems with this approach, particularly when dealing with concurrency. All this places high demands on the designers who now face not only an accelerated schedule but also controllers of increasing complexity. This report describes a new and fundamentally different approach to the construction of safety-critical digital controllers. Instead of directly constructing a complete model and then trying to verify it, the designer can start with an initial abstract (think “sketch”) model plus the requirements, from which a correct concrete model is automatically synthesized. There is no need for post-hoc verification of required functional properties. Having tool to carry this out will significantly impact the nation’s ability to ensure the safety of high-consequence digital systems. The approach has been implemented in a prototype tool, along with a suite of examples, including ones that reflect actual problems faced by designers. Our approach operates on a variant of Statecharts developed at Sandia called Qspecs. Statecharts are a widely used formalism for developing concurrent reactive systems, supporting scalability through allowing state models containing composite states, which are the serial or parallel composition of substates which can themselves contain statecharts. Statecharts enable an incremental style of development, in which states are progressively refined to incorporate greater detail in an incremental model of software development. Our approach formulates a set of constraints from the structure of the models and the requirements and propagates these constraints to a fixpoint. The solution to the constraints is an inductive invariant along with guards on the transitions. We also show how our approach extends to implementation refinement, decomposition, composition, and elaboration. We currently handle safety requirements written in LTL (Linear Temporal Logic)

42 ENGINEERING

Reuse at the Software Productivity Consortium

The Software Productivity Consortium is sponsored by 14 aerospace companies as a developer of software engineering methods and tools. Software reuse and prototyping are currently the major emphasis areas. The Methodology and Measurement Project in the Software Technology Exploration Division has developed some concepts for reuse which they intend to develop into a synthesis process. They have identified two approaches to software reuse: opportunistic and systematic. The assumptions underlying the systematic approach, phrased as hypotheses, are the following: the redevelopment hypothesis, i.e., software developers solve the same problems repeatedly; the oracle hypothesis, i.e., developers are able to predict variations from one redevelopment to others; and the organizational hypothesis, i.e., software must be organized according to behavior and structure to take advantage of the predictions that the developers make. The conceptual basis for reuse includes: program families, information hiding, abstract interfaces, uses and information hiding hierarchies, and process structure. The primary reusable software characteristics are black-box descriptions, structural descriptions, and composition and decomposition based on program families. Automated support can be provided for systematic reuse, and the Consortium is developing a prototype reuse library and guidebook. The software synthesis process that the Consortium is aiming toward includes modeling, refinement, prototyping, reuse, assessment, and new construction.

Weiss, David M.

The Insoluble Carbonaceous Material of CM Chondrites as Possible Source of Discrete Organics During the Asteroidal Aqueous Phase

The larger portion of the organic carbon in carbonaceous chondrites (CC) is present as a complex and heterogeneous macromolecular material that is insoluble in acids and most solvents (IOM). So far, it has been analyzed only as a whole by microscopy (TEM) and spectroscopy (IR, NMR, EPR), which have offered and overview of its chemical nature, bonding, and functional group composition. Chemical or pyrolytic decomposition has also been used in combination with GC-MS to identify individual compounds released by these processes. Their value in the recognition of the original IOM structure resides in the ability to properly interpret the decomposition pathways for any given process. We report here a preliminary study of IOM from the Murray meteorite that combines both the analytical approaches described above, under conditions that would realistically model the IOM hydrothermal exposure in the meteorite parent body. The aim is to document the possible release of water and solvent soluble organics, determine possible changes in NMR spectral features, and ascertain, by extension, the effect of this loss on the frame of the IOM residue. Additional information is included in the original extended abstract.

Yabuta, H.

Distribution of SLS Integrated Load Uncertainty to Surface Pressures and Sectional Loads

Aerodynamic loads that are important to launch vehicle programs such as NASA’s Space Launch System (SLS) include both integrated loads such as the force & moment on the entire vehicle and distributed loads. In this work two examples of distributed loads are considered: the pressure field on the surface of the vehicle and sectional loads, which are one-dimensional distributions along the axis of the launch vehicle. In some modern flight programs, the integrated loads, such as lift and drag, used to design the guidance and control laws for the vehicle come from wind tunnel testing, while distributed loads are produced using Computational Fluid Dynamics(CFD). The first task that is addressed in this paper, then, is to provide a formal method to adjust the distributed loads so that integrating them matches the prescribed integrated load. In addition, the integrated loads in a launch vehicle typically include an uncertainty estimate. The second task is to distribute this prescribed integrated uncertainty to each point in a distributed load. Both tasks are addressed using the same technique, which is to create distributed load profiles that isolate adjustments to one integrated load while leaving the others unaffected. These adjustments are informed by Proper Orthogonal Decomposition (POD) of the entire CFD-based distributed load database. Once applied, the adjusted distributed loads can be used to evaluate any scalar quantity of interest that might be needed by downstream users such as structural analysis or trajectory modelers.

SLS

The deep structure of lunar basins - Implications for basin formation and modification

Models for the crustal structure in the vicinity of nine impact basins, from an inversion of gravity and topographic data from the lunar nearside are presented. The models display a low-density nonmare crustal layer and a mare basalt layer, both of variable thicknesses. Assuming that topography in mare areas is isostatically compensated before the emplacement of mare basalts and that compensation of mare basalt units may be neglected, a decomposition of the gravity anomaly into contributions from Moho relief and mare fill is permitted. Minimum values for mare basalt thicknesses are obtained but because mare basalts and mantle material are similar in density, the thicknesses of the nonmare crust are estimated. An important constraint is the crustal thickness inferred from the Apollo 12 and 14 landing sites from seismic observations. The crustal thickness model indicates that the crust is thinner beneath each of the major nearside basins than in surrounding areas. New bounds on the volume of material ejected from each basin are derived. The geological implications of structural differences among basins for the processes of basin formation and modification are evaluated as functions of time on the moon.

Bratt, S. R.

Smagglce: Surface Modeling and Grid Generation for Iced Airfoils: Phase 1 Results

SmaggIce (Surface Modeling and Grid Generation for Iced Airfoils) is a software toolkit used in the process of aerodynamic performance prediction of iced airfoils with grid-based Computational Fluid Dynamics (CFD). It includes tools for data probing, boundary smoothing, domain decomposition, and structured grid generation and refinement. SmaggIce provides the underlying computations to perform these functions, a GUI (Graphical User Interface) to control and interact with those functions, and graphical displays of results, it is being developed at NASA Glenn Research Center. This paper discusses the overall design of SmaggIce as well as what has been implemented in Phase 1. Phase 1 results provide two types of software tools: interactive ice shape probing and interactive ice shape control. The ice shape probing tools will provide aircraft icing engineers and scientists with an interactive means to measure the physical characteristics of ice shapes. On the other hand, the ice shape control features of SmaggIce will allow engineers to examine input geometry data, correct or modify any deficiencies in the geometry, and perform controlled systematic smoothing to a level that will make the CFD process manageable.

Vickerman, Mary B.

Goated: goal-oriented tensor decompositions in python

SAND2026-20464O Goated performs goal-oriented tensor decompositions in Python, enabling efficient compression of multi-dimensional simulation data. It extends common tensor decomposition methods by incorporating domain-specific knowledge, such as conservation laws in physics, through a penalty term in the optimization process. This approach improves data compression and modeling accuracy across various applications, including physics simulations, by using specialized algorithms and structure-aware subroutines to accelerate solver performance. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

SciDAC

Closed-form solutions of performability

Methods which yield closed form performability solutions for continuous valued variables are developed. The models are similar to those employed in performance modeling (i.e., Markovian queueing models) but are extended so as to account for variations in structure due to faults. In particular, the modeling of a degradable buffer/multiprocessor system is considered whose performance Y is the (normalized) average throughput rate realized during a bounded interval of time. To avoid known difficulties associated with exact transient solutions, an approximate decomposition of the model is employed permitting certain submodels to be solved in equilibrium. These solutions are then incorporated in a model with fewer transient states and by solving the latter, a closed form solution of the system's performability is obtained. In conclusion, some applications of this solution are discussed and illustrated, including an example of design optimization.

Meyer, J. F.

Toward a More Robust Pruning Procedure for MLP Networks

Choosing a proper neural network architecture is a problem of great practical importance. Smaller models mean not only simpler designs but also lower variance for parameter estimation and network prediction. The widespread utilization of neural networks in modeling highlights an issue in human factors. The procedure of building neural models should find an appropriate level of model complexity in a more or less automatic fashion to make it less prone to human subjectivity. In this paper we present a Singular Value Decomposition based node elimination technique and enhanced implementation of the Optimal Brain Surgeon algorithm. Combining both methods creates a powerful pruning engine that can be used for tuning feedforward connectionist models. The performance of the proposed method is demonstrated by adjusting the structure of a multi-input multi-output model used to calibrate a six-component wind tunnel strain gage.

Stepniewski, Slawomir W.

Bulge Growth and Quenching Since Z=2.5 in Candels/3D-HST

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift informationprovided by 3D-HST, we investigate the relation between structure and stellar populations fora mass-selected sample of 6764 galaxies above 1010 M, spanning the redshift range 0.5 z 2.5.For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e.,bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar massmaps reconstructed from resolved stellar population modeling. We confirm that the increased bulgeprominence among quiescent galaxies, as reported previously based on rest-optical observations, remainsin place when considering the distributions of stellar mass. Moreover, we observe an increaseof the typical Sersic index and bulge-to-total ratio (with median BT reaching 40-50) among starforminggalaxies above 1011 M. Given that quenching for these most massive systems is likely tobe imminent, our findings suggest that significant bulge growth precedes a departure from the starformingmain sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge andtotal mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy thanthe total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model bySomerville et al. In the latter, bulges and black holes grow hand in hand through merging andordisk instabilities, and AGN-feedback shuts off star formation. Further observations will be requiredto pin down star formation quenching mechanisms, but our results imply they must be internal to thegalaxies and closely associated with bulge growth.

Quenching

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Stabilization of large-scale systems - A spinning flexible spacecraft

Stabilization of high-order constant systems by multilevel control is proposed in the framework of the decomposition-aggregation method for stability analysis of large-scale systems. An 11-order linear model of a spinning flexible spacecraft is stabilized by the method using both the passive and the active stabilization devices. The method can take advantage of the special structural features of the model to provide the best estimate of the stability region for an important parameter which couples the wobble and the spin motion of the spacecraft.

Siljak, D. D.

Energy Finite Element Analysis for Computing the High Frequency Vibration of the Aluminum Testbed Cylinder and Correlating the Results to Test Data

The Energy Finite Element Analysis (EFEA) is a finite element based computational method for high frequency vibration and acoustic analysis. The EFEA solves with finite elements governing differential equations for energy variables. These equations are developed from wave equations. Recently, an EFEA method for computing high frequency vibration of structures either in vacuum or in contact with a dense fluid has been presented. The presence of fluid loading has been considered through added mass and radiation damping. The EFEA developments were validated by comparing EFEA results to solutions obtained by very dense conventional finite element models and solutions from classical techniques such as statistical energy analysis (SEA) and the modal decomposition method for bodies of revolution. EFEA results have also been compared favorably with test data for the vibration and the radiated noise generated by a large scale submersible vehicle. The primary variable in EFEA is defined as the time averaged over a period and space averaged over a wavelength energy density. A joint matrix computed from the power transmission coefficients is utilized for coupling the energy density variables across any discontinuities, such as change of plate thickness, plate/stiffener junctions etc. When considering the high frequency vibration of a periodically stiffened plate or cylinder, the flexural wavelength is smaller than the interval length between two periodic stiffeners, therefore the stiffener stiffness can not be smeared by computing an equivalent rigidity for the plate or cylinder. The periodic stiffeners must be regarded as coupling components between periodic units. In this paper, Periodic Structure (PS) theory is utilized for computing the coupling joint matrix and for accounting for the periodicity characteristics.

Vlahopoulos, Nickolas

New Synthetic Route to 4,6‐Diamino‐5,7‐dinitro‐benzo‐furazan, Important Decomposition Product of 1,3,5‐Triamino‐2,4,6‐trinitrobenzene

The important 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) decomposition byproduct, 4,6-diamino-5,7-dinitro-benzo-furazan (F 1 ), was prepared in a four-step reaction sequence starting with commercially available materials through a key intermediate, 4,6-dichloro-5,7-dinitro-benzo-furazan. The introduction of azido to 4,6-dichloro-benzo-furazan, followed by nitration and amination gave F 1 in good yield. This new method showed a greatly improved yield and ease of purification than previous methods (the total yield was 30%–40% in seven steps). F 1 is a stable yellow solid with a melting point of 308.7°C and a thermal decomposing temperature of 310.3°C (differential scanning calorimetry at 10°C/min heating rate). The solid-state structure was solved by single-crystal X-ray analyses at low temperatures (100 K). The compound crystallizes in a P 2 1 c space group with the formula C 6 H 4 N 6 O 5 •H 2 O with 12 molecules of F 1 /H 2 O in the unit cell. This compound has eluded complete characterization for almost 50 years. With this information, more advanced decomposition models can be created to further improve critical safety margins in handling the energetic material, TATB.

decomposition

The multicomponent structure of bulges

The morphology of disk galaxies is usually described by two major components, the centrally concentrated spheroidal component, called the bulge, and an oblate disk. The ratio of there contribution to the total luminosity - the bulge-to-disk ratio - is one of the parameters characterizing the Hubble sequence. Following de Vaucouleurs (1948), for most galaxies the radial distribution of the outer spheroid is fairly well described by the r exp 1/4 law I(r)=I(sub 0) exp(-(alpha)r), whereas the radial luminosity distribution of the disk follows an exponential law: I(r)=I(sub 0) exp(-alpha(r exp 1/4)) (Freeman 1970), with r the radial distance from the center. I(sub 0) and alpha are characteristic constants for each individual galaxy. Parameters for the structural properties of these components give important constraints for models of the formation and evolution of galaxies. Therefore we have tried to decompose disk and bulge components from high S/N CCD observations of a sample of edge-on disk galaxies. A common procedure for the decomposition is to model one component in a region where it dominates and subtract it from the combined light distribution. This technique was successfully carried out e.g. by van der Kruit & Earl (1981, 1982) and Wakamatsu & Hamabe (1984, 1989). Here we present two more examples of bulge-dominated edge-on SO galaxies, namely ESO 506-G33 and NGC 7123, which show an additional small and concentrated central component besides disk and 'bulge'.

Dettmar, Ralf-Juergen

Multifluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. The interstellar pickup protons are treated in the model as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. The simulation assumes that the background magnetic field on the Sun is either a dipole (aligned or tilted with respect to the solar rotation axis) or one that is deduced from solar magnetograms.

Goldstein, Melvyn L.

Magnetofluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. In this report, the interstellar pickup protons are treated as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. Where possible, the model is compared with Voyager data. Initial results from generalization to a three-fluid model is described elsewhere in this session.

Goldstein, Melvyn L.

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING