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 577 records · Page 32

Refinement of Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Computational simulations of the flow within a streamline-traced, external-compression supersonic inlet for Mach 1.664 without and with vortex generators were performed to refine the characterization of the inlet performance as measured by the total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. The refinement of the simulations concerned two aspects: 1) refinement of the grid for the simulations to evaluate and reduce uncertainty, and 2) refinement of the modeling and design of the vortex generators. The vortex generators studied were rectangular vane-type vortex generators arranged in co-rotating arrays within the subsonic diffuser. The vortex generator geometric factors of interest included the height, circumferential spacing, and angle-of-incidence. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vortex generators were simulated using either a vortex generator model or with grids generated about each vortex generator. Statistical methods were used to compute confidence intervals and grid convergence indices to establish the uncertainties of the analyses with respect to grid refinement. Design-of-experiments methods were applied to quantify the effects of the geometric factors of the vortex generators. The analyses of the computed results illustrate the complexities of quantifying the uncertainties of the inlet performance and the implications of the uncertainties for the design of a vortex generator array for the STEX inlet.

Slater, John W.↗

On the nature of the interaction between structures and actuators in vibration suppression

This lecture presents a summary of work on the analysis of the interaction between a structure, an actuator used to suppress the vibration of the structure, and the control law implemented by the actuator. Two control applications are considered. First, a proof-mass actuator with experimentally verified dynamics, capable of being used in a space structure configuration is examined. This is connected to a cantilevered beam modeled as a single degree of freedom system. Secondly, an electric motor, again with experimentally verified dynamics, is used to slew a beam modeled by a partial differential equation. Both the experimental and numerical configuration predict the presence of potential instabilities in system performance if proper consideration is not given to interactions between the control law, the structure and the actuator. In addition, an understanding of the interaction is shown to greatly effect performance. In particular it is shown that judicious choice of actuator parameters greatly improves closed loop performance.

Inman, Daniel J.↗

Operator inference with roll outs for learning reduced models from scarce and low-quality data

Data-driven modeling has become a key building block in computational science and engineering. However, data that are available in science and engineering are typically scarce, often polluted with noise and affected by measurement errors and other perturbations, which makes learning the dynamics of systems challenging. Here, in this work, we propose to combine data-driven modeling via operator inference with the dynamic training via roll outs of neural ordinary differential equations. Operator inference with roll outs inherits interpretability, scalability, and structure preservation of traditional operator inference while leveraging the dynamic training via roll outs over multiple time steps to increase stability and robustness for learning from low-quality and noisy data. Numerical experiments with data describing shallow water waves and surface quasi-geostrophic dynamics demonstrate that operator inference with roll outs provides predictive models from training trajectories even if data are sampled sparsely in time and polluted with noise of up to 10%.

97 MATHEMATICS AND COMPUTING↗

Bumper 3 Update for IADC Protection Manual

The Bumper code has been the standard in use by NASA and contractors to perform meteoroid/debris risk assessments since 1990. It has undergone extensive revisions and updates [NASA JSC HITF website; Christiansen et al., 1992, 1997]. NASA Johnson Space Center (JSC) has applied BUMPER to risk assessments for Space Station, Shuttle, Mir, Extravehicular Mobility Units (EMU) space suits, and other spacecraft (e.g., LDEF, Iridium, TDRS, and Hubble Space Telescope). Bumper continues to be updated with changes in the ballistic limit equations describing failure threshold of various spacecraft components, as well as changes in the meteoroid and debris environment models. Significant efforts are expended to validate Bumper and benchmark it to other meteoroid/debris risk assessment codes. Bumper 3 is a refactored version of Bumper II. The structure of the code was extensively modified to improve maintenance, performance and flexibility. The architecture was changed to separate the frequently updated ballistic limit equations from the relatively stable common core functions of the program. These updates allow NASA to produce specific editions of the Bumper 3 that are tailored for specific customer requirements. The core consists of common code necessary to process the Micrometeoroid and Orbital Debris (MMOD) environment models, assess shadowing and calculate MMOD risk. The library of target response subroutines includes a board range of different types of MMOD shield ballistic limit equations as well as equations describing damage to various spacecraft subsystems or hardware (thermal protection materials, windows, radiators, solar arrays, cables, etc.). The core and library of ballistic response subroutines are maintained under configuration control. A change in the core will affect all editions of the code, whereas a change in one or more of the response subroutines will affect all editions of the code that contain the particular response subroutines which are modified. Note that the Bumper II program is no longer maintained or distributed by NASA.

Christiansen, Eric L.↗

A multi-ion generalized transport model of the polar wind

The higher-order generalizations of the equations of standard hydrodynamics, known collectively as generalized transport theories, have been used since the early 1980s to describe the terrestrial polar wind. Inherent in the structure of generalized transport theories is the ability to describe not only interparticle collisions but also certain non-Maxwellian processes, such as heat flow and viscous stress, that are characteristic of any plasma flow that is not collision dominated. Because the polar wind exhibits a transition from collision-dominated to collisionless flow, generalized transport theories possess advantages for polar wind modeling not shared by either collision-dominated models (such as standard hydrodynamics) or collisionless models (such as those based on solving the collisionless Boltzmann equation). In general, previous polar wind models have used generalized transport equations to describe electrons and only one species of ion (H(+)). If other ion species were included in the models at all, it was in a simplified or semiempirical manner. The model described in this paper is the first polar wind model that uses a generalized transport theory (bi-Maxwellian-based 16-moment theory) to describe all of the species, both major and minor, in the polar wind plasma. In the model, electrons and three ion species (H(+), He(+), O(+)) are assumed to be major and several ion species are assumed to be minor (NO(+), Fe(+), O(++)). For all species, a complete 16-moment transport formulation is used, so that profiles of density, drift velocity, parallel and perpendicular temperatures, and the field-aligned parallel and perpendicular energy flows are obtained. In the results presented here, emphasis is placed on describing those constituents of the polar wind that have received little attention in past studies. In particular, characteristic solutions are presented for supersonic H(+) outflow and for both supersonic and subsonic outflows of the major ion He(+). Solutions are also presented for various minor ions, both atomic and molecular and both singly and multiply charged.

Demars, H. G.↗

The structure of the planets Jupiter and Saturn

Planetary models for Jupiter and Saturn are computed using a fourth-order theory and an equation of state for the molecular hydrogen and helium planetary envelopes taken from the Monte Carlo calculations of Slattery and Hubbard (1976). Models for Jupiter are found that have a small amount of heavy elements either mixed with hydrogen and helium throughout the interior of the planet or concentrated in a small dense core. Saturn is modeled with a solar-composition hydrogen and helium envelope and a small dense core. It is concluded that the molecular equation of state linked with suitable interior equations of state can produce Jovian models which satisfy the observational data. The planetary models show that the enrichment of heavy elements (relative to solar composition) is approximately 3 times for Jupiter and 10 times for Saturn.

Slattery, W. L.↗

The equilibrium structure of thin magnetic flux tubes. I

A model atmosphere within a thin magnetic flux tube that is embedded in an arbitrarily stratified medium is presently constructed by solving the radiative transfer equation in the two-stream approximation for gray opacity, under the assumption that the magnetic field is sufficiently strong to warrant the neglect of both thermal conduction and convective diffusion; energy inside the flux tube therefore being transported solely by radiation. The structure of the internal atmosphere is determined on the basis of the hydrostatic and radiative equilibrium conditions of the tube embedded in an external atmosphere. The gas temperature along the axis of the tube is computed, and the geometry of the flux tube is determined on the basis of magnetic flux conservation and total pressure equilibrium.

Ferrari, A.↗

Update on the Status of the FLUKA Monte Carlo Transport Code

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. Here we review the progresses achieved in the last year on the physics models. From the point of view of hadronic physics, most of the effort is still in the field of nucleus--nucleus interactions. The currently available version of FLUKA already includes the internal capability to simulate inelastic nuclear interactions beginning with lab kinetic energies of 100 MeV/A up the the highest accessible energies by means of the DPMJET-II.5 event generator to handle the interactions for greater than 5 GeV/A and rQMD for energies below that. The new developments concern, at high energy, the embedding of the DPMJET-III generator, which represent a major change with respect to the DPMJET-II structure. This will also allow to achieve a better consistency between the nucleus-nucleus section with the original FLUKA model for hadron-nucleus collisions. Work is also in progress to implement a third event generator model based on the Master Boltzmann Equation approach, in order to extend the energy capability from 100 MeV/A down to the threshold for these reactions. In addition to these extended physics capabilities, structural changes to the programs input and scoring capabilities are continually being upgraded. In particular we want to mention the upgrades in the geometry packages, now capable of reaching higher levels of abstraction. Work is also proceeding to provide direct import into ROOT of the FLUKA output files for analysis and to deploy a user-friendly GUI input interface.

Pinsky, L.↗

Effects of equatorially trapped ions on refilling of the plasmasphere

The kinetic aspects of shock formation in response to the equatorial heating of ions are reported by means of small-scale kinetic simulations of countersteaming plasma flows along a spatially varying magnetic field having the features of the dipolar magnetic field lines. Attention is given to the trapping of the field-aligned flowing ions in response to the simulated equatorial heating of ions and the associated self-consistent structures in the electric potentials and the f-plasma flow pattern. Trapping is found to lead to the formation of electrostatic shocks, which propagate away from the 'equator' as seen in the hydrodynamic model of Singh and Torr (1990). The shocks form near the mirror points of the 'equatorially' heated ions. The second stage of the refilling with equatorially trapped ions is reported. A variety of structures in density and temperature distributions as the refilling proceeds slowly is found.

Singh, Nagendra↗

Simulation of Gas Mixtures in RF Discharges for Nitride Deposition

Tungsten and titanium nitride films have long been grown by chemical vapor deposition (CVD) methods. However, there has been recent interest in low temperature growth using plasma enhanced CVD. For the present work, we focus on the radio frequency (BE) discharge characteristics of gas mixtures used in nitride deposition (for example, WF6 and ammonia). Because the radial variations for a standard 200 mm, parallel plate reactor are limited to a small zone near the edges of the electrodes, a 1-D (one-dimensional) analysis is considered. This model consists of a self-consistent, 3-D (three-dimensional) moment fluid simulation that solves the continuity, momentum, and energy equations for neutral and charged species. The results in terms of plasma structure, radical concentrations, and local deposition rate will be presented. We will also compare the 1-D results with those obtained from a 2-D hybrid plasma equipment model (HPEM) developed at the University of Illinois.

Hwang, Helen H.↗

Exploring the Neutron Substructure with Advanced Polarized Helium-3 Targets (Or: How I Learned to Stop Worrying and Love Spectroscopy)

As we seek to understand the smallest, physical aspects of our universe, we cannot simply rely on our senses to probe the world around us as we did in the past. The smallest physical elements of our universe behave in strange, probabilistic ways and are completely invisible to the naked eye/ear/etc. So, we design clever experiments (such as scattering experiments) to probe these minute realms. Then, just as with the larger, observable world, we devise models and equations to describe what we think is happening. Due to the nature of the physical universe at the quantum scale and with the aid of symmetries such as Lorentz invariance, we can write down equations that describe the scattering, but the expressions contain functions, which we call ?form factors? and ?structure functions?, that we cannot compute from first principles. We can, however, formulate models that make predictions for these functions. By comparing our predictions with the observed data, we can gain insight into the validity of our models and thus a better physical understanding of what is happening at these minuscule scales. Studying the constituents inside of the nucleus of an atom adds another layer of difficulty if we can?t remove those components from the nucleus. This is the case with the neutron. When not bound in the nucleus with protons and other neutrons, the neutron will decay into a proton after about 15 minutes. So, we?re forced to study the neutron while it is still bound in the nucleus of an atom such as helium-3 (3He). For the last 1,000 years (rounding up), our group has developed high quality, polarized 3He targets made of an aluminosilicate glass. These targets are made in order to perform experiments at Jefferson Lab (JLab), experiments which let us determine the form factors and structure functions of the neutron by scattering polarized electrons from polarized neutrons (or rather polarized 3He). The specific experiments reported on in this thesis push the bounds of our understanding of the internal structure of the neutron. Good science is often about pushing experimental techniques to a new level. Toward that goal we study our polarized 3He targets both to advance the technology and to choose the best ones for our experiments. We do this using a process called nuclear magnetic resonance (NMR) to gauge the maximum polarization of a target and how fast the polarization decays with time. While these tests primarily provide us information that make analysis of our experimental scattering data possible, they also let us determine whether or not a target-cell is useful or even, dare I say, of spectacular quality. Our latest targets utilize a novel convection design allowing 3He to be polarized and quickly moved in front of the electron-beam, making it possible to use larger targets with higher electron-beam currents than ever before. This means more electrons scatter and we get more data. And by studying our targets in detail prior to using them in our experiments, we have found techniques to take effects which could have been detrimental to target quality and turn them to our advantage! It?s a real case of making lemonade out of lemons. We also use laser spectroscopy to study the absorption lines of alkali-metals in the target (potassium and rubidium, specifically). We add these alkali-metals to our target to facilitate polarizing the 3He. We can use the measurement of these pressure broadened absorption lines to determine the 3He density inside of the target with great precision. Historically, we understood the width of these lines would be dependent on the temperature of the target. Specifically, if I raise the temperature, the width should get bigger. I found that was not the case, which was very confusing at first, though very exciting now that I realize the data are self-consistent and suggestive of unexpected behavior. This thesis details the development of high quality, glass, polarized 3He targets for the 2020 An 1 /dn 2 and 2023 Gn E experiments, which utilized the first 3He convection targets and broke records in target quality. This thesis also covers the initial development of metal windows for the next-generation of 3He target-cells. Finally, this thesis documents the temperature dependence of the width of potassium (K) and rubidium (Rb) absorption lines as measured with laser spectroscopy.

Jantzi, Christopher↗

Phase Transitions and Thermal Equation of State of Fe‐9wt.%Si Applied to the Moon and Mercury

Abstract Accurate knowledge of the phase transitions and thermoelastic properties of candidate iron alloys, such as Fe‐Si alloys, is essential for understanding the nature and dynamics of planetary cores. The phase diagrams of some Fe‐Si alloys between 1 atm and 16 GPa have been back‐extrapolated from higher pressures, but the resulting phase diagram of Fe 83.6 Si 16.4 (9 wt.% Si) is inconsistent with temperature‐induced changes in its electrical resistivity between 6 and 8 GPa. This study reports in situ synchrotron X‐ray diffraction (XRD) measurements on pre‐melted and powder Fe 83.6 Si 16.4 samples from ambient conditions to 60 GPa and 900 K using an externally heated diamond‐anvil cell. Upon compression at 300 K, the bcc phase persisted up to ∼38 GPa. The hcp phase appeared near 8 GPa in the pre‐melted sample, and near 17 GPa in the powder sample. The appearance of the hcp phase in the pre‐melted sample reconciles the reported changes in electrical resistivity of a similar sample, thus resolving the low‐pressure region of the phase diagram. The resulting high‐temperature Birch‐Murnaghan equation of state (EoS) and thermal EoS based on the Mie‐Gruneisen‐Debye model of the bcc and hcp structures are consistent with, and complement the literature data at higher pressures. The calculated densities based on the thermal EoS of Fe‐9wt.%Si indicate that both bcc and hcp phases agree with the reported core density estimates for the Moon and Mercury.

Berrada, Meryem↗

Fortran mimetic abstraction language (Formal) v0.1.

The Fortran mimetic abstraction language ("Formal") is a domain-specific language (DSL) embedded in Fortran 202Y [1]. Formal provides novel software abstractions for simulating phenomena governed by the partial differential equations (PDEs) of vector and tensor calculus. Such equations model an extremely broad set of physical phenomena, ranging from atmospheric winds to light propagation. Formal's data structures and algorithms mimic in form and behavior continuous functions and operators. Formal supports these mathematical constructs using mimetic discretizations that define a discrete calculus satisfying various tensor calculus theorems, thereby ensuring high-fidelity representations of the physics being modeled. [2] Formal 0.1.0 also lays a foundation for the future use of Fortran 202Y type-safe templates to facilitate the formal verification of tensor contractions in computational physics and artificial intelligence [3]. [1] "Fortran 202Y" is Fortran standard committee's informal designation for the next Fortran revision, which will likely be "Fortran 2028". [2] Corbino, J. and Castillo, J. (2020) Journal of Computational and Applied Mathematics, https://doi.org/10.1016/j.cam.2019.06.042. [3] Haveraaen, M., Järvi, J., & Rouson, D. (2019). Reflecting on Generics for Fortran. https://j3-fortran.org/doc/year/19/19-188.pdf.

Rouson, Damian [Lawrence Berkeley National Laborat↗

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

A numerical investigation of severe thunderstorm gust fronts

A numerical model was developed to simulate the evolution and structure of severe thunderstorm gust fronts. The model is a non-hydrostatic, fine resolution, cross-sectional primitive equation model. Two-dimensional horizontal and vertical equations of motion, the continuity equation, and the thermodynamic energy equation were utilized. It was shown that two dominant factors influencing gust front configuration are surface friction and the solenoidal field coincident with the front. It is suggested that solenoidal accelerations oppose the deceleration of surface friction. After a downdraft is initiated in the model, these opposing tendencies soon reach a balance and the gust front achieves a quasi-steady configuration. Thus, the experiments indicate that surface friction does not induce a cycle of front formation and collapse. In addition, the effect of evaporative cooling in producing a vigorous downdraft was parameterized by a local cooling function. Greater cooling in the downdraft results in a more intense gust front that exhibits stronger wind maximums and greater shears. The ambient air stability was shown to be an important factor influencing the depth of the cold outflow.

Mitchell, K. E.↗

Fokker-Planck response of stochastic satellites

The present investigation is concerned with the effects of stochastic geometry and random environmental torques on the pointing accuracy of spinning and three-axis stabilized satellites. The study of pointing accuracies requires a knowledge of the rates of error growth over and above any criteria for the asymptotic stability of the satellites. For this reason the investigation is oriented toward the determination of the statistical properties of the responses of the satellites. The geometries of the satellites are considered stochastic so as to have a phenomenological model of the motions of the flexible structural elements of the satellites. A widely used method of solving stochastic equations is the Fokker-Planck approach where the equations are assumed to define a Markoff process and the transition probability densities of the responses are computed directly as a function of time. The Fokker-Planck formulation is used to analyze the response vector of a rigid satellite.

Huang, T. C.↗

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The Helicopter Equations for Stability and Loads (HESL) program is extended to derive the governing equations of motion for the blade, and a Lagrangian formulation is used to obtain the equations in generalized coordinates. The program generates the steady-state and linearized perturbation equations in symbolic form and then codes them into FORTRAN subroutines. The coefficients for each equation and for each mode are identified through a numerical program; the latter can also be used to obtain the harmonic balance equations. The governing multiblade equations are derived explicitly using HESL. These equations can accommodate any number of elastic blade modes. Stability results are presented for several hingeless rotor blade structural models, and the influence of dynamic inflow in forward flight with an elastic hingeless rotor is investigated.

Reddy, T. S. R.↗

Core plasma in the magnetosphere

The paper presents a review of new findings related to core plasma (which includes ions and electrons with energies less than 50 eV) identified in studies of magnetosphere during 1983-1986. These findings include the cleft ion fountain; the plasmaspheric heavy ion torus; core molecular ions; high core-plasma densities in the plasma-sheet boundary layer; intense transverse heating of core plasmas in the equatorial regions; the supersonic polar wind (H/+/ and He/+/); toroidal or ring heavy-ion distribution functions in the auroral region; conic-to-field-aligned ion signatures of inverted-V auroral events; and the large-scale electron density structure of the magnetosphere. The progress made on modeling core plasma includes simulation of thermal helium heating at the equator; kinetic models of the plasma transport from localized ionospheric regions; time-dependent hydrodynamic models of the polar wind and plasmasphere-ionosphere coupling; and kinetic models of plasmasphere refilling.

Horwitz, J. L.↗