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 523 records · Page 29

Reconstructing quintessence

We present an Effective Field Theory based reconstruction of quintessence models of dark energy directly from cosmological data. We show that current cosmological data possess enough constraining power to test several quintessence model properties for redshifts z ϵ [0, 1.5] with no assumptions about the behavior of the scalar field potential. We use measurements of the cosmic microwave background, supernovae distances, and the clustering and lensing of galaxies to constrain the evolution of the dark energy equation of state, Swampland Conjectures, the shape of the scalar field reconstructed potential, and the structure of its phase space. The standard cosmological model still remains favored by data and, within quintessence models, deviations from its expansion history are bounded to be below the 10% level at 95% confidence at any redshift below z = 1.5.

79 ASTRONOMY AND ASTROPHYSICS↗

Domain-aware Control-oriented Neural Models for Autonomous Underwater Vehicles

Conventional physics-based modeling is a time-consuming bottleneck in control design for complex nonlinear systems like autonomous underwater vehicles (AUVs). In contrast, purely data-driven models, require a large number of observations and lack operational guarantees for safety-critical systems. Data-driven models leveraging available partially characterized dynamics have potential to provide reliable systems models in a typical data-limited scenario for high value complex systems, thereby avoiding months of expensive expert modeling time. In this work we explore this middle-ground between expert-modeled and pure data-driven modeling. We present control-oriented parametric models with varying levels of domain-awareness that exploit known system structure and prior physics knowledge to create constrained deep neural dynamical system models. We employ universal differential equations to construct data-driven blackbox and graybox representations of the AUV dynamics. In addition, we explore a hybrid formulation that explicitly models the residual error related to imperfect graybox models. We compare the prediction performance of the learned models for different distributions of initial conditions and control inputs to assess their suitability for control.

Shaw Cortez, Wenceslao E.↗

Neutron Star Equation of State in Light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. Here, we construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Enabling Real-time Scattering Data Analysis with Scalable Optimization [Slides]

Diffraction experiments produce datasets with rich multidimensional physics information such as microstructure, equations of state, crystal structure, elastoplastic properties, and other key inputs to LANL mission-essential multiphysics models. This information is typically extracted through a process called Rietveld refinement, which involves selecting appropriate models of the instrument, crystal structure, and microstructure, identifying suitable starting values, and then fitting often hundreds of model parameters using a sequence of empirical parameter turnon/off sequences within a non-global gradient-based optimization. Extensive user expertise is required to properly setup a refinement, identify appropriate models, and select initial parameter values close to truth, such that the refinement will yield parameter values that are optimally predictive. This is a very tedious manual process performed far after the beamline campaign has ended. As facilities have become capable of generating larger volumes of data, the limitation in throughput due to Rietveld refinement has led to a dramatic increase in unanalyzed data as opposed to an intended increase in new science. In our FY22 TED, we demonstrated an integrated toolset providing near real-time automated Rietveld analysis. If this toolset can be optimized to provide automated Rietveld analysis in real-time, this could alleviate the bottleneck in unanalyzed diffraction data, aid in decision-making during experiments, and increase efficiency of the facility.

74 ATOMIC AND MOLECULAR PHYSICS↗

Review of magnetic field observations

Recent observations of magnetic fields in the magnetosphere are reviewed, and critical experiments and data are identified for theoretical analysis and interpretation. Quantitative studies of the solar wind interaction with the earth's magnetic field, regional measurements near the earth's equator at R = 2-8 R sub E, the polar cusp region of the geomagnetosphere, and structural models of the neutral sheet region in the geomagnetic tail are considered.

Ness, N. F.↗

An analysis of the Venus measurements

Plots of the Mariner 5 and Venera 4, 5, and 6 pressure vs temperature show that the Venera profiles are essentially congruent with the Mariner 5 day and night profiles, but are displaced 28 percent higher in pressure. Data suggest a variability in the atmospheric water vapor content in the region from 30 to 50 km. The Venera 7 measurements are interpreted as evidence for updrafts, down-drafts, horizontal wind layers, and nonaqueous precipitation. The previously observed band of retrograde winds which circle the equator with an average speed of 110 m/s is found to extend downward to the one atmosphere level at the equatorial morning terminator. The possibility of a low altitude equator-to-pole circulation with warm gas rising at the poles is inferred. Venera 7 temperature data used with radar topography and microwaves interferometer measurements suggest that the variation of surface temperature with altitude in a band about the equator is less than 5 K/km. The available data are used to calculate a model of the structure of the Venus atmosphere for the first 75 km above the equatorial region.

Ainsworth, J. E.↗

Physicochemical and dynamical processes in cometary ionospheres. I - The basic flow profile

The paper deals with the basic aspects of solar wind interaction with the ionosphere of a large comet and with detailed hydrodynamic models developed to describe the flow exterior to the tangential discontinuity surface that separates the outflowing purely cometary ions from the inflowing contaminated solar wind. The equation of the basic flow model and the results of a preliminary analysis of these equations are given, and the shape and structure of the shock layer are analyzed. The shock layer structure associated with fan-shaped comet emanating from the cometary nucleus in the general direction to the sun is identified.

Houpis, H. L. F.↗

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.↗

Stable feedback control of distributed parameter systems: Time and frequency domain conditions

Large space structures, or any mechanically flexible structures, are inherently distributed parameter systems (DPS) whose dynamics are modeled by partial, rather than ordinary, differential equations. Such DPS are described by operator equations on an infinite-dimensional Hilbert (or Banach) space. However, any feedback controller for such a DPS must be a finite-dimensional (and discrete-time) system in order to be implemented with on-line digital computers and a finite (small) number of actuators and sensors. There are many ways to synthesize such controllers; the Galerkin or finite-element approach is emphasized.

Balas, M. J.↗

A modal model for SCOLE structural dynamics

Information on a modal model for the Spacecraft Control Laboratory (SCOLE) is given in viewgraph form. A partial differential equation model covering roll bending, pitch bending, torsion, shear forces, and bending moments is given.

Joshi, Suresh M.↗

Environmental effects on long term behavior of composite laminates

Model equations are presented for approximate methods simulating the long-term behavior of composite materials and structures in hot/humid service environments. These equations allow laminate property upgradings with time, and can account for the effects of service environments on creep response. These methodologies are illustrated for various individual and coupled temperature/moisture, longitudinal/transverse, and composite material type cases. Creep deformation is noted to rise dramatically for cases of matrix-borne, but not of fiber-borne, loading in hot, humid environments; the coupled influence of temperature and moisture is greater than a mere combination of their individual influences.

Singhal, S. N.↗

A three dimensional multigrid Reynolds-averaged Navier-Stokes solver for unstructured meshes

A three-dimensional unstructured mesh Reynolds averaged Navier-Stokes solver is described. Turbulence is simulated using a single field-equation model. Computational overheads are minimized through the use of a single edge-based data-structure, and efficient multigrid solution technique, and the use of multi-tasking on shared memory multi-processors. The accuracy and efficiency of the code are evaluated by computing two-dimensional flows in three dimensions and comparing with results from a previously validated two-dimensional code which employs the same solution algorithm. The feasibility of computing three-dimensional flows on grids of several million points in less than two hours of wall clock time is demonstrated.

Mavriplis, D. J.↗

Helioseismic Constraints on New Solar Models from the MoSEC Code

Evolutionary solar models are computed using a new stellar evolution code, MOSEC (Modular Stellar Evolution Code). This code has been designed with carefully controlled truncation errors in order to achieve a precision which reflects the increasingly accurate determination of solar interior structure by helioseismology. A series of models is constructed to investigate the effects of the choice of equation of state (OPAL or MHD-E, the latter being a version of the MHD equation of state recalculated by the author), the inclusion of helium and heavy-element settling and diffusion, and the inclusion of a simple model of mixing associated with the solar tachocline. The neutrino flux predictions are discussed, while the sound speed of the computed models is compared to that of the sun via the latest inversion of SOI-NMI p-mode frequency data. The comparison between models calculated with the OPAL and MHD-E equations of state is particularly interesting because the MHD-E equation of state includes relativistic effects for the electrons, whereas neither MHD nor OPAL do. This has a significant effect on the sound speed of the computed model, worsening the agreement with the solar sound speed. Using the OPAL equation of state and including the settling and diffusion of helium and heavy elements produces agreement in sound speed with the helioseismic results to within about +.-0.2%; the inclusion of mixing slightly improves the agreement.

Elliott, J. R.↗

Development and Application of Agglomerated Multigrid Methods for Complex Geometries

We report progress in the development of agglomerated multigrid techniques for fully un- structured grids in three dimensions, building upon two previous studies focused on efficiently solving a model diffusion equation. We demonstrate a robust fully-coarsened agglomerated multigrid technique for 3D complex geometries, incorporating the following key developments: consistent and stable coarse-grid discretizations, a hierarchical agglomeration scheme, and line-agglomeration/relaxation using prismatic-cell discretizations in the highly-stretched grid regions. A signi cant speed-up in computer time is demonstrated for a model diffusion problem, the Euler equations, and the Reynolds-averaged Navier-Stokes equations for 3D realistic complex geometries.

Nishikawa, Hiroaki↗

G-DYN Multibody Dynamics Engine

G-DYN is a multi-body dynamic simulation software engine that automatically assembles and integrates equations of motion for arbitrarily connected multibody dynamic systems. The algorithm behind G-DYN is based on a primal-dual formulation of the dynamics that captures the position and velocity vectors (primal variables) of each body and the interaction forces (dual variables) between bodies, which are particularly useful for control and estimation analysis and synthesis. It also takes full advantage of the spare matrix structure resulting from the system dynamics to numerically integrate the equations of motion efficiently. Furthermore, the dynamic model for each body can easily be replaced without re-deriving the overall equations of motion, and the assembly of the equations of motion is done automatically. G-DYN proved an essential software tool in the simulation of spacecraft systems used for small celestial body surface sampling, specifically in simulating touch-and-go (TAG) maneuvers of a robotic sampling system from a comet and asteroid. It is used extensively in validating mission concepts for small body sample return, such as Comet Odyssey and Galahad New Frontiers proposals.

Acikmese, Behcet↗

Turbulent Axial Odometer Model

A Galilean-invariant field equation is proposed and tested on standard turbulence model test cases. The field equation provides an additional non-dimensional outer scale which allows the turbulence models to reproduce the axial normal stress increase with Re seen in high Reynolds numbers experiments. The field equation provides a Reynolds number for every point based on the length of turbulent flow upstream of that point in the domain. This outer scale equation can be considered an odometer that gives a length scale conjectured to be related to the large stream wise structures that are seen in turbulent flow and that require run length to develop. A new RANS model using this additional scale is able to match the Reynolds number variation of the normal stresses seen at high Reynolds number. Furthermore, the good attached flow prediction capabilities of current RANS models appears to be attained. Using this scale equation, the entire Reynolds-stress state appears to be predicted correctly, over a large run length Reynolds number range such as experienced in aircraft design.

Olsen, Michael E.↗

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A revised version of the transfer matrix method to analyze one-dimensional structures

A new and general method to analyze both free and forced vibration characteristics of one-dimensional structures is discussed in this paper. This scheme links for the first time the classical transfer matrix method with the recently developed integrating matrix technique to integrate systems of differential equations. Two alternative approaches to the problem are presented. The first is based upon the lumped parameter model to account for the inertia properties of the structure. The second releases that constraint allowing a more precise description of the physical system. The free vibration of a straight uniform beam under different support conditions is analyzed to test the accuracy of the two models. Finally some results for the free vibration of a 12th order system representing a curved, rotating beam prove that the present method is conveniently extended to more complicated structural dynamics problems.

Nitzsche, F.↗