Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse modeling”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation, surface-to-surface radiation exchange, and flowfield coupling. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
DOE OSTI2024

Quasilinear modelling of collisional trapped electron modes

The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.

fusion plasma↗
DOE OSTI2020

Nonadiabatic Energy Fluctuations of Scale-Invariant Quantum Systems in a Time-Dependent Trap

We consider the nonadiabatic energy fluctuations of a many-body system in a time-dependent harmonic trap. In the presence of scale-invariance, the dynamics becomes self-similar and the nondiabatic energy fluctuations can be found in terms of the initial expectation values of the second moments of the Hamiltonian, square position, and squeezing operators. Nonadiabatic features are expressed in terms of the scaling factor governing the size of the atomic cloud, which can be extracted from time-of-flight images. We apply this exact relation to a number of examples: the single-particle harmonic oscillator, the one-dimensional Calogero-Sutherland model, describing bosons with inverse-square interactions that includes the non-interacting Bose gas and the Tonks-Girdardeau gas as limiting cases, and the unitary Fermi gas. We illustrate these results for various expansion protocols involving sudden quenches of the trap frequency, linear ramps and shortcuts to adiabaticity. Our results pave the way to the experimental study of nonadiabatic energy fluctuations in driven quantum fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Machine Learning (ML) Classifier to Assist Metadata Creation

The Atmospheric Radiation Measurement (ARM) Data Center is responsible for the timely collection, archival, and curation of science data products. These products are freely available through an online data repository. Metadata creation is paramount for scientific users to find and access over seven petabytes of atmospheric science data. The hierarchical metadata structure allows users to search for information at both broad and narrow levels. This project aims to leverage 30 years’ worth of manually created metadata to enable machine predictions of broad-term classifications from narrow-term descriptions. These classification predictions would assist metadata coordinators with their term selections. This paper discusses the cleaning and preprocessing of the training data, the pipeline developed to determine the best model for this task, and the creation of an API metadata classifier for ARM measurement metadata. Our results show that the Linear Support Vector Classification (LinearSVC) algorithm, along with the Term Frequency – Inverse Document Frequency (TF-IDF) vectorizer, is well-suited for our multi-class classification task. Lengthier input training data led to better results, and artificial balancing was unnecessary for this particular use case. This predictive classifier enhances efficiency in metadata creation, as well as supports greater consistency and accuracy in metadata tagging.

Collier, Hannah [ORNL] (ORCID:0000000341284292)↗
NASA NTRS1977

The inversion of stratospheric aerosol and ozone vertical profiles from spacecraft solar extinction measurements

The inversions of multi-channel solar extinction measurements have been analyzed for the 0.35-1.0 micron wavelength region to retrieve stratospheric aerosol and ozone vertical profiles using both the constrained linear inversion scheme and the iterative scheme. The inversions of the multi-wavelength solar extinction data obtained from spacecraft have been analyzed based on the inversion of computer simulated data using various atmospheric models with differing amounts of aerosol and ozone in the stratosphere. The sensitivities of the inversion schemes to different experimental errors are discussed in terms of accuracy and resolution of the retrieved profiles.

Chu, W. P.↗
NASA NTRS1999

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q model) or a point dipole (the P model). The Q model has four unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P model has six unknown parameters and describes many intracloud (IC) discharges. In this paper we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a dme-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer-simulated sources and on natural lightning at the NASA Kennedy Space Center and U.S. Air Force Eastern Range.

Koshak, William J.↗
NASA NTRS2013

Shape and Stress Sensing of Multilayered Composite and Sandwich Structures Using an Inverse Finite Element Method

The marked increase in the use of composite and sandwich material systems in aerospace, civil, and marine structures leads to the need for integrated Structural Health Management systems. A key capability to enable such systems is the real-time reconstruction of structural deformations, stresses, and failure criteria that are inferred from in-situ, discrete-location strain measurements. This technology is commonly referred to as shape- and stress-sensing. Presented herein is a computationally efficient shape- and stress-sensing methodology that is ideally suited for applications to laminated composite and sandwich structures. The new approach employs the inverse Finite Element Method (iFEM) as a general framework and the Refined Zigzag Theory (RZT) as the underlying plate theory. A three-node inverse plate finite element is formulated. The element formulation enables robust and efficient modeling of plate structures instrumented with strain sensors that have arbitrary positions. The methodology leads to a set of linear algebraic equations that are solved efficiently for the unknown nodal displacements. These displacements are then used at the finite element level to compute full-field strains, stresses, and failure criteria that are in turn used to assess structural integrity. Numerical results for multilayered, highly heterogeneous laminates demonstrate the unique capability of this new formulation for shape- and stress-sensing.

Cerracchio, Priscilla↗
NASA NTRS1981

Radio occultation of Jupiter's ring - Bounds on optical depth and particle size and a comparison with infrared and optical results

A comparison of Voyager 1 data of the Jovian ring at radio, infrared and optical wavelengths suggests a population density that either falls more rapidly than the inverse square of the linear size, or is sharply bounded in maximum particle size. A fragmentation power law is used to estimate a minimum particle radius of 1-2 microns, and a specific model is developed which consists of a power law distribution of lossless or slightly absorbent Mie scatterers with a refractive index near 1.62, a power law index of -3.5, and a minimum particle size of 1.5 microns. The maximum particle radius is not critical, provided that it is greater than 4 microns, and the power law polydispersion of lossless or low absorbent Mie scatterers is consistent with previously reported observations as long as the parameter size minimum is 20 and the power index is about -3.5.

Tyler, G. L.↗
DOE OSTI2024

Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P∝ k m , where k is the wavenumber. The hydrodynamic condition is the same as in the $\theta$-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The $\theta$-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ~1.86 shock and gamma-law = 5/3 fluids with Atwood number A = 0.5. The bubble amplitude h B is found to grow in two stages. Initially, h B undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, h B grows as a power law $~t^{\theta}$ when $\langle{k}\rangle{|h_B|}$ > O(1), where $\langle{k}\rangle$ is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Inversion of Crater Morphometric Data to Gain Insight on the Cratering Process

In recent years, morphometric data for Venus and several outer planet satellites have been collected, so we now have observational data of complex Craters formed in a large range of target properties. We present general inversion techniques that can utilize the morphometric data to quantitatively test various models of complex crater formation. The morphometric data we use in this paper are depth of a complex crater, the diameter at which the depth-diameter ratio changes, and onset diameters for central peaks, terraces, and peak rings. We tested the roles of impactor velocities and hydrostatic pressure vs. crustal strength, and we tested the specific models of acoustic fluidization (Melosh, 1982) and nonproportional growth (Schultz, 1988). Neither the acoustic fluidization model nor the nonproportional growth in their published formulations are able to successfully reproduce the data. No dependence on impactor velocity is evident from our inversions. Most of the morphometric data is consistent with a linear dependence on the ratio of crustal strength to hydrostatic pressure on a planet, or the factor c/pg.

Herrick, Robert R.↗
NASA NTRS2009

Efficient Method for Analyzing Hyperspectral Remote Sensing Data

The current and next generation satellite infrared sensors such as AIRS, IASI, CrIS, and CLARREO all have one thing in common: large number of spectral channels. In order to fully utilize the information content of these sensors, a large number of Radiative Transfer (RT) calculations through the inhomogeneous atmosphere are needed. It is also computationally intensive to invert atmospheric and surface properties using all the channels. Usually, only subsets of channels are used to perform physical inversions for atmospheric profiles. We will explore ways to speed up RT calculations and inversions and test the algorithms using data from the Infrared Atmospheric Sounding Interferometer (IASI) satellite instrument. We will describe a retrieval algorithm using a Principal Component-based Radiative Transfer Model (PCRTM) for generating atmospheric temperature/moisture/ozone profiles and surface properties. The retrieval algorithm performs both the radiative transfer calculations and inversions in the Principal Component (PC) domain. The inversion algorithm is based on a non-linear Levenberg- Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the ultraspectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels. We will also use data collected during the Joint Airborne IASI Validation Experiment (JAIVEx) field campaign to validate the algorithm and IASI retrievals.

Liu, Xu↗
DOE OSTI2022

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRSDate not supplied

Adjusting a Full Body Model to Mitigate Inverse Kinematics Artifacts in Opensim

BACKGROUND: In support of Vibration Isolation and Stabilization (VIS) system development for Human Health Countermeasures (HHC) exercise systems in space, such as the European Enhanced Exploration Exercise Device (E4D) [1], dynamic quantities required to model the response of a proposed VIS while considering the effect of VIS motion on the forces between the human and VIS platform were obtained using motion capture data [2]. On occasion, large-amplitude oscillatory spikes were found in the subject’s linear and angular momentum derivatives, affecting analyses that depend on forces and moments derived from motion capture. The purpose of this investigation was to identify causes of these artifacts and techniques for their resolution. METHODS AND RESULTS: To obtain the required human dynamic quantities to drive the VIS simulation, motion capture data was collected containing recorded trajectories of passive retroreflective markers on body landmarks of an exercising subject. Since the full body Rajagopal model [3] was originally used to enhance gait analysis, upper body joints did not require large Ranges of Motion (ROM). We thus modified the Rajagopal model [4, 5] to allow it to be used for upper body intensive exercises like those common to the E4D. OpenSim Inverse Kinematics (IK) [6] was performed using these scaled subject models to generate the joint angles throughout the exercise while minimizing marker error. At times, the arms were observed to ‘snap’ from one configuration to another, causing spike artifacts. Following IK, a custom OpenSim plugin [7] was used to determine the required dynamic quantities including the rates of change of the linear and angular momenta of the human. Since motion capture is recorded at a larger time step than required by the VIS simulation, the human center of mass location was fit with splines and a second derivative taken to obtain the momentum derivative, allowing the VIS simulation to maintain conservation of momentum when appropriate. In a few cases during this stage, artifacts much larger than the expected noise of the second derivatives were introduced. Investigation of cases containing artifacts revealed several modifications that could be made to the OpenSim model to improve IK results. Since OpenSim models use Euler angles and rotation sequences, ‘gimbal lock’ would be encountered in the arms when raised 90 degrees to the side (e.g., T-pose, some hang clean exercise, etc.). This was resolved by reorienting the horizontal axes at the shoulder joint by 45 degrees, placing ‘gimbal lock’ outside common arm ROM, with the arm ROMs adjusted following this change. Elbow and wrist ROMs could also be adjusted to allow realistic motion while at the same time limiting the likelihood of unrealistic orientations. On occasion, the arms flipped backwards when raised above the head. This was prevented by using medial elbow markers in scaling and IK. When medial markers were not available, the acromial joint location in the unscaled model was shifted before model scaling to better align the arm with available markers. Lastly, artifacts which became apparent after taking the second derivatives of spline-fit data were found to occur in cases when the pelvis rotation limit prevented the full range of motion of an exercise. These issues were resolved by unclamping the pelvis rotation limit. Through this investigation, an understanding of conditions leading to IK artifacts was acquired allowing the automation of artifact detection. These artifact detection and mitigation techniques can be applied toward modeling of upper body motions in aerospace and other fields for improved IK results.

C A Bell↗
NASA NTRS1989

Nonlinear and distributed parameter models of the mini-mast truss

Large spacecraft such as Space Station Freedom employ large trusses in their construction. The structural dynamics of such trusses often exhibit nonlinear behavior and little damping which can impact significantly the performance of control systems. The Mini-MAST truss was constructed to research such structural dynamics and control systems. The Mini-MAST truss is an object of study for the guest investigator program as part of NASA's controls-structures interaction program. The Mini-MAST truss is deployable and about 65 ft long. Although the bending characteristics of the Mini-MAST truss are essentially linear, the angular deflection under torsional loading has exhibited significant hysteresis and nonlinear stiffness. It is the purpose to develop nonlinear and distributed parameter models of the truss and to compare the model dynamics with actual measurements. Distributed parameter models have the advantage of requiring fewer model parameters. A tangent function is used to describe the nonlinear stiffness in torsion, partly because of the convenience of its easily expressed inverse. Hysteretic slip elements are introduced and extended to a continuum to account for the observed hysteresis in torsion. The contribution of slipping to the structural damping is analyzed and found to be strongly dependent on the applied loads. Because of the many factors which affect the damping and stiffness in a truss, it is risky to assume linearity.

Taylor, Lawrence W., Jr.↗
DOE OSTI2022

Clump survival and migration in VDI galaxies: an analytical model versus simulations and observations

ABSTRACT We address the nature of the giant clumps in high-z galaxies that undergo violent disc instability, distinguishing between long-lived and short-lived clumps. We study the evolution of long-lived clumps during migration through the disc via an analytical model tested by simulations and confront theory with CANDELS-HST observations. The clump ‘bathtub’ model, which considers gas and stellar gain and loss, involves four parameters: the accretion efficiency α, the star formation rate (SFR) efficiency ϵd, and the outflow mass-loading factors for gas and stars, η and ηs. The corresponding time-scales are comparable to the migration time, two-three orbital times. The accretion-rate dependence on clump mass, gas, and stars, allows an analytical solution involving exponential growing and decaying modes. For the fiducial parameter values there is a main evolution phase where the SFR and gas mass are constant and the stellar mass is rising linearly with time. This makes the inverse specific SFR an observable proxy for clump age. When η or ϵd are high, or α is low, the decaying mode induces a decline of SFR and gas mass till the migration ends. Later, the masses and SFR approach an hypothetical exponential growth with a constant specific SFR. The model matches simulations with different, moderate feedbacks, both in isolated and cosmological settings. The observed clumps agree with our predictions, indicating that the massive clumps are long-lived and migrating. A challenge is to model feedback that is non-disruptive in massive clumps but suppresses SFR to match the galactic stellar-to-halo mass ratio.

Dekel, Avishai↗
NASA NTRS2000

X-38 Experimental Controls Laws

X-38 Experimental Control Laws X-38 is a NASA JSC/DFRC experimental flight test program developing a series of prototypes for an International Space Station (ISS) Crew Return Vehicle, often called an ISS "lifeboat." X- 38 Vehicle 132 Free Flight 3, currently scheduled for the end of this month, will be the first flight test of a modem FCS architecture called Multi-Application Control-Honeywell (MACH), originally developed by the Honeywell Technology Center. MACH wraps classical P&I outer attitude loops around a modem dynamic inversion attitude rate loop. The dynamic inversion process requires that the flight computer have an onboard aircraft model of expected vehicle dynamics based upon the aerodynamic database. Dynamic inversion is computationally intensive, so some timing modifications were made to implement MACH on the slower flight computers of the subsonic test vehicles. In addition to linear stability margin analyses and high fidelity 6-DOF simulation, hardware-in-the-loop testing is used to verify the implementation of MACH and its robustness to aerodynamic and environmental uncertainties and disturbances.

Munday, Steve↗