Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Manifold model”

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

Grassmannian diffusion maps based surrogate modeling via geometric harmonics

Abstract A novel surrogate model based on the Grassmannian diffusion maps (GDMaps) and utilizing geometric harmonics (GH) is developed for predicting the response of complex physical phenomena. The method utilizes GDMaps to obtain a low‐dimensional representation of the underlying behavior of physical/mathematical systems with respect to uncertain input parameters. Using this representation, GH, an out‐of‐sample extension technique, is employed to create a global map from the input parameter space to a Grassmannian diffusion manifold. GH is further employed to locally map points on the diffusion manifold onto the tangent space of a Grassmann manifold. The exponential map is then used to project the points in the tangent space onto the Grassmann manifold, where reconstruction of the full solution is performed. The performance of the proposed surrogate model is verified with three examples. The first problem is a toy example used to illustrate the technique. In the second example, errors associated with the various mappings are assessed by studying response predictions of the electric potential of a dielectric cylinder in a homogeneous electric field. The last example applies the method for uncertainty prediction in the strain field evolution in a model amorphous material using the shear transformation zone theory of plasticity.

42 ENGINEERING↗

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗

Experimental aerodynamic evaluation of the SSME two-duct, hot gas manifold

An investigation was performed on an improved Space Shuttle Main Engine (SSME) hot gas manifold (HGM) design using a full scale HGM model in cold-flow, air-blowdown tests. The new HGM design replaces the current three transfer ducts with two enlarged elliptical transfer ducts, faired duct inlets, and enlarged fuel bowl. Extensive experimental results have been obtained to verify previous results showing that the two-duct design improved the HGM internal flow path by eliminating separated flow regions, reducing high local velocities, and providing more uniform flow.

Kim, B. Y.↗

Time domain identification of element matrices

Accurate mathematical models are clearly an important part of the design, modification, control, and damage assessment of large space structures. A critical part of model determination of any type of structure is the use of system identification (SI), a process for using measured excitation and response data to improve the form of a mathematical model or the values of the parameters in it. The reasons for using SI are manifold, but most of them involve uncertainties in either the form of or the parameters in the mathematical model. Some of the more common uncertainties are the following: the nature of the damping, the characteristics of sloppy or sticky joints, inelastic material properties, and parameters in a reduced Degree-of-Freedom (DOF) model. This paper contains a description of the SI algorithm followed by two illustrations: one in which all parameters in a shear building model are determined using pull back and quick release data, and another in which an equivalent reduced DOF model is obtained for a two story frame.

Matzen, Vernon C.↗

Optimization on Manifolds via Graph Gaussian Processes

This paper integrates manifold learning techniques within a Gaussian process upper confidence bound algorithm to optimize an objective function on a manifold. Our approach is motivated by applications where a full representation of the manifold is not available and querying the objective is expensive. We rely on a point cloud of manifold samples to define a graph Gaussian process surrogate model for the objective. Query points are sequentially chosen using the posterior distribution of the surrogate model given all previous queries. We establish regret bounds in terms of the number of queries and the size of the point cloud. Several numerical examples complement the theory and illustrate the performance of our method.

Bayesian optimization↗

Lattice Clifford fractons and their Chern-Simons-like theory

We use Dirac matrix representations of the Clifford algebra to build fracton models on the lattice and their effective Chern-Simons-like theory. As an example, we build lattice fractons in odd D spatial dimensions and their (D+1) spacetime dimensional effective theory. The model possesses an anti-symmetric K matrix resembling that of hierarchical quantum Hall states. The gauge charges are conserved in sub-dimensional manifolds which ensures the fractonic behavior. The construction extends to any lattice fracton model built from commuting projectors and with tensor products of spin-1/2 degrees of freedom at the sites.

Fontana, Weslei↗

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗

Revealing the hidden structure of disordered materials by parameterizing their local structural manifold

Abstract Durable interest in developing a framework for the detailed structure of glassy materials has produced numerous structural descriptors that trade off between general applicability and interpretability. However, none approach the combination of simplicity and wide-ranging predictive power of the lattice-grain-defect framework for crystalline materials. Working from the hypothesis that the local atomic environments of a glassy material are constrained by enthalpy minimization to a low-dimensional manifold in atomic coordinate space, we develop a generalized distance function, the Gaussian Integral Inner Product (GIIP) distance, in connection with agglomerative clustering and diffusion maps, to parameterize that manifold. Applying this approach to a two-dimensional model crystal and a three-dimensional binary model metallic glass results in parameters interpretable as coordination number, composition, volumetric strain, and local symmetry. In particular, we show that a more slowly quenched glass has a higher degree of local tetrahedral symmetry at the expense of cyclic symmetry. While these descriptors require post-hoc interpretation, they minimize bias rooted in crystalline materials science and illuminate a range of structural trends that might otherwise be missed.

36 MATERIALS SCIENCE↗

Saturn - A study of the 3 nu sub 3 methane band.

Photometric spectra are given for various manifolds of the R-branch of this methane band of Saturn. The spectra were obtained by measurements, primarily in the winter of 1970, along the planet's central meridian. Analysis of the manifold shapes in terms of a reflecting-layer model proved largely self-consistent. It implies a rotational temperature of about 140 K and a C/H ratio in Saturn's atmosphere which is close to that on the sun.

Trafton, L.↗

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗

The structure of random discrete spacetime

The usual picture of spacetime consists of a continuous manifold, together with a metric of Lorentzian signature which imposes a causal structure on the spacetime. A model, first suggested by Bombelli et al., is considered in which spacetime consists of a discrete set of points taken at random from a manifold, with only the causal structure on this set remaining. This structure constitutes a partially ordered set (or poset). Working from the poset alone, it is shown how to construct a metric on the space which closely approximates the metric on the original spacetime manifold, how to define the effective dimension of the spacetime, and how such quantities may depend on the scale of measurement. Possible desirable features of the model are discussed.

Brightwell, Graham↗

Structure of random discrete spacetime

The usual picture of spacetime consists of a continuous manifold, together with a metric of Lorentzian signature which imposes a causal structure on the spacetime. A model, first suggested by Bombelli et al., is considered in which spacetime consists of a discrete set of points taken at random from a manifold, with only the causal structure on this set remaining. This structure constitutes a partially ordered set (or poset). Working from the poset alone, it is shown how to construct a metric on the space which closely approximates the metric on the original spacetime manifold, how to define the effective dimension of the spacetime, and how such quantities may depend on the scale of measurement. Possible desirable features of the model are discussed.

Brightwell, Graham↗

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗

Boundary Layer Control for Hypersonic Airbreathing Vehicles

Active and passive methods for tripping hypersonic boundary layers have been examined in NASA Langley Research Center wind tunnels using a Hyper-X model. This investigation assessed several concepts for forcing transition, including passive discrete roughness elements and active mass addition (or blowing), in the 20-Inch Mach 6 Air and the 31-Inch Mach 10 Air Tunnels. Heat transfer distributions obtained via phosphor thermography, shock system details, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. The comparisons between the active and passive methods for boundary layer control were conducted at test conditions that nearly match the Hyper-X nominal Mach 7 flight test-point of an angle-of-attack of 2-deg and length Reynolds number of 5.6 million. For passive roughness, the primary parametric variation was a range of trip heights within the calculated boundary layer thickness for several trip concepts. The passive roughness study resulted in a swept ramp configuration, scaled to be roughly 0.6 of the calculated boundary layer thickness, being selected for the Mach 7 flight vehicle. For the active blowing study, the manifold pressure was systematically varied (while monitoring the mass flow) for each configuration to determine the jet penetration height, with schlieren, and transition movement, with the phosphor system, for comparison to the passive results. All the blowing concepts tested, which included various rows of sonic orifices (holes), two- and three-dimensional slots, and random porosity, provided transition onset near the trip location with manifold stagnation pressures on the order of 40 times the model surface static pressure, which is adequate to ensure sonic jets. The present results indicate that the jet penetration height for blowing was roughly half the height required with passive roughness elements for an equivalent amount of transition movement.

Berry, Scott A.↗

Viscous flow effects on hydrogen leaks from cracks in the Orbiter Challenger main engines

An analytical model was developed to provide additional insight and understanding of the factors that influence the simulation and prediction of leak rates from small cracks in pressurized containers. Specifically, the analysis was aimed at developing an analytical model capable of predicting the hydrogen leak rates from a crack in the combustion chamber coolant discharge manifold on main engine 1 of the Orbiter Challenger that was discovered during flight readiness firings 1 and 2. This model was based on viscous pipe flow analyses and calibrated for the crack geometry by using helium leak-rate data obtained from both low- and high-pressure tests used to simulate the flight readiness firing test conditions. In addition, this model includes the effects of crack width changes caused by different working stresses associated with the different test conditions. Because of the combination of the small crack dimensions and the wide range of pressures used for the test conditions, either laminar or turbulent viscous effects dominated the flows at all test conditions. This model was used to illustrate the sensitivity of the predicted leak rates to considerations of test conditions, viscous flow effects, and geometric features of the crack. In addition, the model was certified by comparing the hydrogen leak-rate prediction for the flight readiness firing test condition to the actual measured leak rate. The prediction was within 9 percent of the measured value.

Goodrich, W. D.↗

A Comparison of Active and Passive Methods for Control of Hypersonic Boundary Layers on Airbreathing Configurations

Active and passive methods for control of hypersonic boundary layers have been experimentally examined in NASA Langley Research Center wind tunnels on a Hyper-X model. Several configurations for forcing transition using passive discrete roughness elements and active mass addition, or blowing, methods were compared in two hypersonic facilities, the 20-Inch Mach 6 Air and the 31-Inch Mach 10 Air tunnels. Heat transfer distributions, obtained via phosphor thermography, shock system details, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. The comparisons between the active and passive methods for boundary layer control were conducted at test conditions that nearly match the nominal Mach 7 flight trajectory of an angle-of-attack of 2-deg and length Reynolds number of 5.6 million. For the passive roughness examination, the primary parametric variation was a range of trip heights within the calculated boundary layer thickness for several trip concepts. The prior passive roughness study resulted in a swept ramp configuration being selected for the Mach 7 flight vehicle that was scaled to be roughly 0.6 of the calculated boundary layer thickness. For the active jet blowing study, the blowing manifold pressure was systematically varied for each configuration, while monitoring the mass flow, to determine the jet penetration height with schlieren and transition movement with the phosphor system for comparison to the passive results. All the blowing concepts tested were adequate for providing transition onset near the trip location with manifold stagnation pressures on the order of 40 times the model static pressure or higher.

Berry, Scott A.↗

Preliminary SP-100/Stirling heat exchanger designs

Analytic modeling of several heat exchanger concepts to couple the SP-100 nuclear reactor primary lithium loop and the Space Stirling Power Convertor (SSPC) was performed. Four 25 kWe SSPC's are used to produce the required 100 kW of electrical power. This design work focused on the interface between a single SSPC and the primary lithium loop. Manifolding to separate and collect the four channel flow was not modeled. This work modeled two separate types of heat exchanger interfaces (conductive coupling and radiative coupling) to explore their relative advantages and disadvantages. The minimum mass design of the conductively coupled concepts was 18 kg or 0.73 kg/kWe for a single 25 kWe convertor. The minimum mass radiatively coupled concept was 41 kg or 1.64 kg/kWe. The direct conduction heat exchanger provides a lighter weight system because of its ability to operate the Stirling convertor evaporator at higher heat fluxes than those attainable by the radiatively coupled systems. Additionally the conductively coupled concepts had relatively small volumes and provide potentially simpler assembly. Their disadvantages were the tight tolerances and material joining problems associated with this refractory to superalloy interface. The advantages of the radiatively coupled designs were the minimal material interface problems.

Schmitz, Paul↗