Search NASA⌕ Search

SEARCH · Search NASA

Results for “structural equation 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Sensitivity Analysis of Tunable Equation of State Material Model In Pulsed Mercury Target Simulation

A pulsed neutron spallation target is subjected to very short but intense loads from repeated proton pulses. Approximately 60% of the energy from each proton pulse is deposited into the mercury target material and the stainless-steel target structure, leading to a high-pressure region in both the stationary target structure and the flowing mercury. The high-pressure region propagates and leads to fluid-structure interaction. The resultant loading on the target structure containing liquid mercury is difficult to predict, although various simulation approaches and material models for the mercury have been tried. To date, the best match of simulation to experimental data is obtained by using an equation of state (EOS) material model with a specified tensile cutoff pressure, which simulates the cavitation threshold. The inclusion of a threshold to represent cavitation is key to the successful predictions of stress waves triggered by the high-energy pulse striking the mercury and vessel. However, recent measurements of target structure strain show that significant discrepancies remain between the measured and simulated strain values in the EOS mercury model. These differences grow when noncondensable helium gas is intentionally injected into the flowing mercury to reduce the loading on the structure. An EOS-based proportional–integral–derivative (PID) mercury model has been proposed to reduce the gap between the measured and simulated vessel strain responses for targets with gas injection. The conceptual and numerical description and initial investigation of the PID model are presented in previous work. Further studies of this PID model — including the sensitivity of the structure’s strain response to model parameters (the tensile cutoff, PID parameters Kp, Ki, and Kd) — are reported in this article. Results show the strain response is more sensitive to changes in the tensile cutoff value than to changes in the model parameters Kp, Ki, and Kd. These results will aid in future work where the model parameters will be optimized to match simulation data to strain measurements.

Lin, Lianshan↗
DOE OSTI2020

Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

In this work, we present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system’s governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh–Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗
DOE OSTI2025

Symmetry structure of a Riccati equation appearing in penetration mechanics

In the design of projectile penetration experiments a matter of considerable interest is scaling: that is, the potential relevance of small-scale experiments to their full-scale counterparts, in a manner analogous to that most often encountered in the context of fluid mechanics. From the theoretical standpoint, phenomena associated with scaling and scalability can be assessed using the well-established tools of dimensional analysis and the Buckingham-Pi Theorem. However, the familiar precepts of dimensional analysis are themselves a specific manifestation of the broader group invariance properties or symmetries of a mathematical model. Here, this work explores these notions – that is, dimensional analysis, the conditions for realizing complete similarity, and any additional symmetry structures – in the context of a Riccati differential equation appearing in the context of penetration mechanics. The aim of the investigation is twofold: 1) to complement existing empirical considerations with a concrete theoretical basis, and 2) to provide a deeper theoretical understanding of the projectile penetration model and its many implications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Applications of the Dulmage–Mendelsohn decomposition for debugging nonlinear optimization problems

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

42 ENGINEERING↗
NASA NTRS2020

Global Aerodynamic Modeling Using Automated Local Model Networks in Real Time

A novel method is presented for automated real-time global aerodynamic modeling using local model networks, known as Smoothed Partitioning with Localized Trees in Real Time (SPLITR), as part of NASA’s Learn-to-Fly technology development initiative. The global nonlinear aerodynamics are partitioned into several local regions known as cells, with the dimension, location, and timing of each partition automatically selected based on a residual characterization procedure, under the constraints of real-time operation. Regression trees represent the successive partitioning of the global flight envelope and describe the evolution of the cell structure. Recursive equation-error least-squares parameter estimation in the time domain is used to estimate a model that represents the local aerodynamics in each region, so that it can be updated independently with non-contiguous data in the range of each cell over time. A weighted superposition of these piecewise local models across the flight envelope forms a global nonlinear model that also accurately captures the local aerodynamics. The SPLITR approach is demonstrated using both simulation and flight data, and the results are analyzed in terms of model predictive capabilities as well as interpretability. The results show that SPLITR can be used to automatically partition complex nonlinear aerodynamic behavior, produce an accurate model, and provide valuable physical insight into the local and global aerodynamics.

Weinstein, Rose↗
NASA NTRS1990

Approximations of thermoelastic and viscoelastic control systems

Well-posed models and computational algorithms are developed and analyzed for control of a class of partial differential equations that describe the motions of thermo-viscoelastic structures. An abstract (state space) framework and a general well-posedness result are presented that can be applied to a large class of thermo-elastic and thermo-viscoelastic models. This state space framework is used in the development of a computational scheme to be used in the solution of a linear quadratic regulator (LQR) control problem. A detailed convergence proof is provided for the viscoelastic model and several numerical results are presented to illustrate the theory and to analyze problems for which the theory is incomplete.

Burns, J. A.↗
NASA NTRS2006

Structural Analysis and Testing of the Inflatable Re-entry Vehicle Experiment (IRVE)

The Inflatable Re-entry Vehicle Experiment (IRVE) is a 3.0 meter, 60 degree half-angle sphere cone, inflatable aeroshell experiment designed to demonstrate various aspects of inflatable technology during Earth re-entry. IRVE will be launched on a Terrier-Improved Orion sounding rocket from NASA s Wallops Flight Facility in the fall of 2006 to an altitude of approximately 164 kilometers and re-enter the Earth s atmosphere. The experiment will demonstrate exo-atmospheric inflation, inflatable structure leak performance throughout the flight regime, structural integrity under aerodynamic pressure and associated deceleration loads, thermal protection system performance, and aerodynamic stability. Structural integrity and dynamic response of the inflatable will be monitored with photogrammetric measurements of the leeward side of the aeroshell during flight. Aerodynamic stability and drag performance will be verified with on-board inertial measurements and radar tracking from multiple ground radar stations. In addition to demonstrating inflatable technology, IRVE will help validate structural, aerothermal, and trajectory modeling and analysis techniques for the inflatable aeroshell system. This paper discusses the structural analysis and testing of the IRVE inflatable structure. Equations are presented for calculating fabric loads in sphere cone aeroshells, and finite element results are presented which validate the equations. Fabric material properties and testing are discussed along with aeroshell fabrication techniques. Stiffness and dynamics tests conducted on a small-scale development unit and a full-scale prototype unit are presented along with correlated finite element models to predict the in-flight fundamental mod

Lindell, Michael C.↗
NASA NTRSDate not supplied

Effects of Rotation on Internal Structure and Dynamics of Main-Sequence Stars

Current state-of-the-art computer simulations allow us to build 3D dynamical and radiative models of stars from physical first principles with a high degree of realism. The radiative 3D dynamical stellar models obtained with the StellarBox code take into account the effects of turbulence, stellar abundances, a realistic equation of state, and radiative energy transport. In this talk, I will discuss the effects of rotation on the turbulent dynamics and surface structure for a 1.47Msun star for rotational periods of 1 and 14 days. The simulations are performed with the computational domain at various latitudes. The models reproduce stellar granulation, the subsurface shear layer, structural changes in convection, and the tachocline, which is the interface between the inner radiative zone and the outer convection zone and plays a crucial role in stellar variability. In particular, the model results reveal the formation of differential rotation and meridional circulation.

SMD↗
DOE OSTI2023

Scientific machine learning for modeling and simulating complex fluids

The formulation of rheological constitutive equations—models that relate internal stresses and deformations in complex fluids—is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these “digital fluid twins” to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation—a task that is not achievable using any previously developed data-driven rheological equation of state.

Science & Technology - Other Topics↗
DOE OSTI2024

Nonlinear Nano‐Imaging of Interlayer Coupling in 2D Graphene‐Semiconductor Heterostructures

The emergent electronic, spin, and other quantum properties of 2D heterostructures of graphene and transition metal dichalcogenides are controlled by the underlying interlayer coupling and associated charge and energy transfer dynamics. However, these processes are sensitive to interlayer distance and crystallographic orientation, which are in turn affected by defects, grain boundaries, or other nanoscale heterogeneities. This obfuscates the distinction between interlayer charge and energy transfer. Here, nanoscale imaging in coherent four‐wave mixing (FWM) and incoherent two‐photon photoluminescence (2PPL) is combined with a tip distance‐dependent coupled rate equation model to resolve the underlying intra‐ and inter‐layer dynamics while avoiding the influence of structural heterogeneities in mono‐ to multi‐layer graphene/WSe 2 heterostructures. With selective insertion of hBN spacer layers, it is shown that energy, as opposed to charge transfer, dominates the interlayer‐coupled optical response. From the distinct nano‐FWM and ‐2PPL tip‐sample distance‐dependent modification of interlayer and intralayer relaxation by tip‐induced enhancement and quenching, an interlayer energy transfer time of τ ET ≈ ( 0 . 35 − 0.15 + 0.65 ) $\tau _{\rm ET} \approx (0.35^{+0.65}_{-0.15})$ ps consistent with recent reports is derived. As a local probe technique, this approach highlights the ability to determine intrinsic sample properties even in the presence of large sample heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Aemulus ν: precision halo mass functions in wνCDM cosmologies

Precise and accurate predictions of the halo mass function for cluster mass scales in wνCDM cosmologies are crucial for extracting robust and unbiased cosmological information from upcoming galaxy cluster surveys. Here, we present a halo mass function emulator for cluster mass scales (≳ 1013 M ⊙/h) up to redshift z = 2 with comprehensive support for the parameter space of wνCDM cosmologies allowed by current data. Based on the Aemulus ν suite of simulations, the emulator marks a significant improvement in the precision of halo mass function predictions by incorporating both massive neutrinos and non-standard dark energy equation of state models. This allows for accurate modeling of the cosmology dependence in large-scale structure and galaxy cluster studies. We show that the emulator, designed using Gaussian Process Regression, has negligible theoretical uncertainties compared to dominant sources of error in future cluster abundance studies. Our emulator is publicly available (https://github.com/DelonShen/aemulusnu_hmf), providing the community with a crucial tool for upcoming cosmological surveys such as LSST and Euclid.

cluster counts↗
DOE OSTI2022

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. Here we calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Marshall system for aerospace simulation (MARSYAS)

System is simple flexible language which can be coded by users unfamiliar with computer programming. It is designed for engineers with little experience in simulation, who desire to simulate large physical systems. User has ability to mix differential equations with diagrams in his model. With few exceptions, there is no rigid statement-operator structure within given module.

Tranboth, H. H.↗
NASA NTRS1975

Summary of theoretical and experimental investigation of grating type, silicon photovoltaic cells

Theoretical and experimental aspects are summarized for single crystal, silicon photovoltaic devices made by forming a grating pattern of p/n junctions on the light receiving surface of the base crystal. Based on the general semiconductor equations, a mathematical description is presented for the photovoltaic properties of such grating-like structures in a two dimensional form. The resulting second order elliptical equation is solved by computer modeling to give solutions for various, reasonable, initial values of bulk resistivity, excess carrier concentration, and surface recombination velocity. The validity of the computer model is established by comparison with p/n devices produced by alloying an aluminum grating pattern into the surface of n-type silicon wafers. Current voltage characteristics and spectral response curves are presented for cells of this type constructed on wafers of different resistivities and orientations.

Chen, L. Y.↗
NASA NTRS1985

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Two coupled, nonlinear differential equations are proposed for the modeling of the elastic and rate (time) dependent inelastic behavior of structural metals in the absence of recovery and aging. The structure of the model is close to the unified theories but contains essential differences. It is shown that the model reproduces almost elastic regions upon initial loading and in the unloading regions of the hysteresis loop. Under loading, unloading and reloading in strain control the model simulated the experimentally observed sharp transition from nearly elastic to inelastic behavior. When a formulation akin to existing unified theories is adopted the almost elastic regions reduce the points and the transition upon reloading is very gradual. For different formulations the behavior under sudden in(de)creases of the strain rate by two orders of magnitude is simulated by numerical experiments and differences are noted. The model represents cyclically neutral behavior and contains three constants and two positive, decreasing functions. The determination of constants and functions from monotonic loading with strain rate changes and relaxation periods is described.

Krempl, E.↗