SEARCH · Engineering Papers

Results for “vector fitting”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2023

WeakIdent: Weak formulation for identifying differential equation using narrow-fit and trimming

Data-driven identification of differential equations is an interesting but challenging problem, especially when the given data are corrupted by noise. When the governing differential equation is a linear combination of various differential terms, the identification problem can be formulated as solving a linear system, with the feature matrix consisting of linear and nonlinear terms multiplied by a coefficient vector. This product is equal to the time derivative term, and thus generates dynamical behaviors. The goal is to identify the correct terms that form the equation to capture the dynamics of the given data. We propose a general and robust framework to recover differential equations using a weak formulation with two new mechanisms, narrow-fit and trimming, for both ordinary and partial differential equations (ODEs and PDEs). The weak formulation facilitates an efficient and robust way to handle noise, and two new mechanisms, narrow-fit and trimming, improve the coefficient support and value recoveries respectively. For each sparsity level, Subspace Pursuit is utilized to find an initial set of support from the large dictionary. Then, we focus on highly dynamic regions (rows of the feature matrix), and error normalize the feature matrix in the narrow-fit step. The support is further updated via trimming the terms that contribute the least. Finally, the support set of features with the smallest Cross-Validation error is chosen as the result. A comprehensive set of numerical experiments are presented for both systems of ODEs and PDEs with various noise levels. The proposed method gives a robust recovery of the coefficients, and a significant denoising effect which can handle up to 100% noise-to-signal ratio for some equations. We compare the proposed method with several state-of-the-art algorithms for the recovery of differential equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Using Downwelling Far- and Thermal-Infrared Hyperspectral Radiance for Cloud Phase Classification in the Antarctic

The cloud phase is one of the most important parameters of clouds. In this paper, we propose a method for cloud phase classification that synergistically utilizes the far- and thermal-infrared bands based on the Atmospheric Emitted Radiance Interferometer (AERI) at the Atmospheric Radiation Measurement West Antarctic Radiation Experiment (AWARE) observatory in 2016. The possible features in the far- and thermal-infrared bands are analyzed based on the differences in the simulated cloud brightness temperature (BT) spectra with different cloud phases. Using the support vector machine (SVM) algorithm, four features are determined to identify the cloud phase, which include the BT at 900 cm -1 , the slope of the fitted function of BT in the 900–1000 cm -1 interval, the BT difference (BTD) between 512 cm -1 and 726 cm -1 , and the BTD between 550 cm -1 and 726 cm -1 . Here, the performance of the proposed method is evaluated with Shupe’s and Turner’s method. The monthly average accuracy of the proposed method, the method without the two far-infrared features, and Turner’s method are about 76%, 36%, and 49%, respectively, which infer the good performance of the proposed method and also indicate that the far-infrared band features can effectively enhance cloud phase classification. It is notable that, compared to Shupe’s method, the accuracy for the proposed method is only 61% during the Antarctic summer, which results from the definitions of cloud phase and radiative effect. In addition, the accuracy is only 44% for Turner’s method in seasons with a low frequency of mixed clouds due to the significant effect of water vapor.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1975

An initial feasibility stage for Stoer's constrained least squares algorithm

A procedure is described for computing an initial feasible vector, x sub 0, for Stoer's algorithm for solving the linear least squares problem subject to linear equality and inequality constraints. The procedure described fits well with Stoer's algorithm since much of the computation performed to determine x sub 0 accomplishes initializing transformations of the problem data, which would otherwise be done in Stoer's algorithm after being given an x sub 0.

Lawson, C. L.↗

NASA NTRS1979

An experimental study of a jet issuing from a lifting wing

An experimental program was conducted to determine the behavior of a round turbulent jet issuing from a lifting two-dimensional wing in crossflow. The jet was located at 65% wing chord on an NACA 0021 airfoil fitted with a 30% chord NACA 4415 flap. The flowfield associated with the jet was surveyed extensively with directional pressure probes to determine local velocity vectors and pressures for three different values of lift coefficient at jet effective velocity ratios (square root of the ratio of the jet dynamic pressure to the freestream dynamic pressure) of 4, 6, and 8. Data describing the jet centerline and the path of the contrarotating vortices accompanying the deflected jet are presented and compared with similar data for a round jet issuing from a large flat plate. The spacing and strength of the vortices are calculated using a simple vortex model previously proposed for the flat plate case. The results show that the penetration of the jet and the vortices increases significantly with increasing lift for the range of test parameters covered in the study. The calculated vortex spacing and strength also show an increase with lift.

Mcmahon, H. M.↗

NASA NTRS1979

A fast, conservative algorithm for solving the transonic full-potential equation

A fast, fully implicit approximate factorization (AF) algorithm designed to solve the conservative transonic full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., 'rotated differencing'), and thereby greatly enhances the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line overrelaxation algorithms.

Holst, T. L.↗

NASA NTRS2014

Conventional and Bimodal Nuclear Thermal Rocket (NTR) Artificial Gravity Mars Transfer Vehicle Concepts

A variety of countermeasures have been developed to address the debilitating physiological effects of "zero-gravity" (0-g) experienced by cosmonauts and astronauts during their approximately 0.5-1.2 year long stays in LEO (Low Earth Orbit). Longer interplanetary flights, combined with possible prolonged stays in Mars orbit, could subject crewmembers to up to approximately 2.5 years of weightlessness. In view of known and recently diagnosed problems associated with 0-g, an artificial gravity spacecraft offers many advantages and may indeed be an enabling technology for human flights to Mars. A number of important human factors must be taken into account in selecting the rotation radius, rotation rate, and orientation of the habitation module or modules. These factors include the gravity gradient effect, radial and tangential Coriolis forces, along with cross-coupled acceleration effects. Artificial gravity (AG) Mars transfer vehicle (MTV) concepts are presented that utilize both conventional NTR, as well as, enhanced "bimodal" nuclear thermal rocket (BNTR) propulsion. The NTR is a proven technology that generates high thrust and has a specific impulse (I (sub sp)) capability of approximately 900 s - twice that of today's best chemical rockets. The AG/MTV concepts using conventional NTP carry twin cylindrical "ISS-type" habitation modules with their long axes oriented either perpendicular or parallel to the longitudinal spin axis of the MTV and utilize photovoltaic arrays (PVAs) for spacecraft power. The twin habitat modules are connected to a central operations hub located at the front of the MTV via two pressurized tunnels that provide the rotation radius for the habitat modules. For the BNTR AG/MTV option, each engine has its own "closed" secondary helium-xenon gas loop and Brayton rotating unit that can generate tens of kilowatts (kW (sub e)) of spacecraft electrical power during the mission coast phase eliminating the need for large PVAs. A single inflatable "TransHab-type" habitation module is also used with multiple vertical floors oriented radial to the MTV spin axis. The BNTR MTV's geometry - long and linear - is naturally compatible with AG operation. By rotating the vehicle about its center-of-mass and perpendicular to its flight vector at approximately 3.0 - 5.2 rpm, a centrifugal force and AG environment corresponding to approximately 0.38 - 1.0 g can be established to help maintain crew fitness out to Mars and back. Vehicles using NTP/BNTP can more readily accommodate the heavier payload mass and increased RCS propellant loading associated with AG operation, and can travel faster to and from Mars thereby reducing the crew's exposure to galactic cosmic radiation and solar flares. Mission scenario descriptions, key vehicle features and operational characteristics for each propulsion options are presented using the lift capability and payload volumes estimated for the SLS-1A and HLV.

Nuclear Propulsion↗

DOE OSTI2022

Rare Higgs Processes at CMS and Precision Timing Detector Studies for HL-LHC CMS Upgrade

This thesis describes the search for two rare Higgs processes. The first analysis describes the CMS Run 2 search for $H$ $\rightarrow$ $\mu$$\mu$ decays, with 137.3 fb$^{-1}$ of data at $\sqrt{s}$ = 13 TeV. The analysis targeted four different Higgs production modes: the gluon fusion (ggH), the vector boson fusion (VBF), the Higgs-strahlung process (VH), and the production in association with a pair of top quarks (ttH). Each category used a dedicated machine learning based classifier to separate the signal from the background processes. A combined fit from all these categories saw a slight excess in the data corresponding to 3.0 standard deviations at $M$$_{H}$ = 125.38 GeV, and gave the first evidence for the Higgs boson decay to second-generation fermions. The best-fit signal strength and the corresponding 68% CL interval was found to be +0.17?????? = 1.19 $_{-0.39}^{+0.41}$ (stat)$_{-0.16}^{+0.17}$(syst) at $M$$_{H}$ = 125.38 GeV. The second analysis describes the CMS Run 2 search for 𝐻𝐻 → 𝑏𝑏𝑏𝑏 with highly boosted Higgs bosons. This analysis used a dedicated jet identification algorithm based on graph neural networks (ParticleNet) to identify boosted H→ bb jets. This search targeted the gluon fusion and the vector boson fusion HH production modes, and put constraints on the allowed values of the various Higgs couplings as: 𝜅𝜆 ∈ [−9.9, 16.9] when 𝜅𝑉 = 1, 𝜅2𝑉 = 1; 𝜅𝑉 ∈ [−1.17, −0.79] ∪ [0.81, 1.18] when 𝜅𝜆 = 1, 𝜅2𝑉 = 1; 𝜅2𝑉 ∈ [0.62, 1.41] when 𝜅𝜆 = 1, 𝜅𝑉 = 1. A scenario with 𝜅2𝑉 = 0 was excluded with a significance of 6.3 standard deviations for the first time, when other H couplings are fixed to their SM values. The combined observed (expected) 95% upper limit on the HH production cross section was found to be 9.9 (5.1) × SM. Finally, this thesis also discusses the planned MIP Timing Detector (MTD) upgrade for CMS at the HL-LHC. The MTD will be a time-of-flight (TOF) detector, designed to provide a precision timing information for charged particles using SiPMs + LYSO scintillating crystals, with a time resolution of ∼30 ps. This thesis describes several R&D tests that have been performed for characterizing the sensor properties (time resolution, light yield, etc.) and optimizing the sensor design geometry. This thesis also contains a description of mock test setups for cooling the sensors, since it is known to be an effective way of mitigating the increased dark current rates in the sensors due to radiation damage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Novel Geometric Operations for Linear Programming

This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Application of multivariate autoregressive spectrum estimation to ULF waves

The estimation of the power spectrum of a time series by fitting a finite autoregressive model to the data has recently found widespread application in the physical sciences. The extension of this method to the analysis of vector time series is presented here through its application to ULF waves observed in the magnetosphere by the ATS 6 synchronous satellite. Autoregressive spectral estimates of the power and cross-power spectra of these waves are computed with computer programs developed by the author and are compared with the corresponding Blackman-Tukey spectral estimates. The resulting spectral density matrices are then analyzed to determine the direction of propagation and polarization of the observed waves.

Ioannidis, G. A.↗

DOE OSTI2022

Simulation of the Impact of Point Defects and Edge Dislocations on X-Ray Diffraction in Hexagonal (Ni,Co) 1+2 x Ti 1– x O 3 Thin Films

In this work, a computer code for simulating high-resolution X-ray diffraction (XRD) data from disordered crystals with arbitrary spatial composition and local lattice parameters is developed. Simulated patterns are compared with the experimental data collected on a single phase, highly crystalline (Ni 0.42 Co 0.58 ) 2.22 Ti 0.39 O 3 solid-solution thin film exhibiting a large number of subsidiary minima. As a case study, the edge dislocations in hexagonal (Ni,Co) 3 O 3 thin films with Burgers vector parallel and perpendicular to the a-axis and c-axis, respectively, are modeled. No peak profiles are assumed and thus issues related to profile fitting are avoided. Both macroscopic features, such as film thickness, and atomic-scale structure, such as dislocations, are simultaneously modeled. Simulations are run on a desktop machine. Commonly applied database-based phase identification routines can result in wrong phase identification, unless intensities are properly modeled. Modeling tools of diffractometer manufacturers are compared with the present approach. An example of how simulated reciprocal lattice patterns can be used to choose relevant measurement geometries in defected thin films is given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1991

Vector magnetic fields in sunspots. I - Stokes profile analysis using the Marshall Space Flight Center magnetograph

The Marshall Space Flight Center (MSFC) vector magnetograph is a tunable filter magnetograph with a bandpass of 125 mA. Results are presented of the inversion of Stokes polarization profiles observed with the MSFC vector magnetograph centered on a sunspot to recover the vector magnetic field parameters and thermodynamic parameters of the spectral line forming region using the Fe I 5250.2 A spectral line using a nonlinear least-squares fitting technique. As a preliminary investigation, it is also shown that the recovered thermodynamic parameters could be better understood if the fitted parameters like Doppler width, opacity ratio, and damping constant were broken down into more basic quantities like temperature, microturbulent velocity, or density parameter.

Balasubramaniam, K. S.↗

NASA NTRS2010

General Pressurization Model in Simscape

System integration is an essential part of the engineering design process. The Ares I Upper Stage (US) is a complex system which is made up of thousands of components assembled into subsystems including a J2-X engine, liquid hydrogen (LH2) and liquid oxygen (LO2) tanks, avionics, thrust vector control, motors, etc. System integration is the task of connecting together all of the subsystems into one large system. To ensure that all the components will "fit together" as well as safety and, quality, integration analysis is required. Integration analysis verifies that, as an integrated system, the system will behave as designed. Models that represent the actual subsystems are built for more comprehensive analysis. Matlab has been an instrument widely use by engineers to construct mathematical models of systems. Simulink, one of the tools offered by Matlab, provides multi-domain graphical environment to simulate and design time-varying systems. Simulink is a powerful tool to analyze the dynamic behavior of systems over time. Furthermore, Simscape, a tool provided by Simulink, allows users to model physical (such as mechanical, thermal and hydraulic) systems using physical networks. Using Simscape, a model representing an inflow of gas to a pressurized tank was created where the temperature and pressure of the tank are measured over time to show the behavior of the gas. By further incorporation of Simscape into model building, the full potential of this software can be discovered and it hopefully can become a more utilized tool.

Servin, Mario↗

NASA NTRS2016

Conventional and Bimodal Nuclear Thermal Rocket (NTR) Artificial Gravity Mars Transfer Vehicle Concepts

A variety of countermeasures have been developed to address the debilitating physiological effects of zero-gravity (0-g) experienced by cosmonauts and astronauts during their approximately 0.5 to 1.2 year long stays in low Earth orbit (LEO). Longer interplanetary flights, combined with possible prolonged stays in Mars orbit, could subject crewmembers to up to approximately 2.5 years of weightlessness. In view of known and recently diagnosed problems associated with 0-g, an artificial gravity (AG) spacecraft offers many advantages and may indeed be an enabling technology for human flights to Mars. A number of important human factors must be taken into account in selecting the rotation radius, rotation rate, and orientation of the habitation module or modules. These factors include the gravity gradient effect, radial and tangential Coriolis forces, along with cross-coupled acceleration effects. Artificial gravity Mars transfer vehicle (MTV) concepts are presented that utilize both conventional NTR, as well as, enhanced bimodal nuclear thermal rocket (BNTR) propulsion. The NTR is a proven technology that generates high thrust and has a specific impulse (Isp) capability of approximately 900 s-twice that of today's best chemical rockets. The AG/MTV concepts using conventional Nuclear Thermal Propulsion (NTP) carry twin cylindrical International Space Station (ISS)- type habitation modules with their long axes oriented either perpendicular or parallel to the longitudinal spin axis of the MTV and utilize photovoltaic arrays (PVAs) for spacecraft power. The twin habitat modules are connected to a central operations hub located at the front of the MTV via two pressurized tunnels that provide the rotation radius for the habitat modules. For the BNTR AG/MTV option, each engine has its own closed secondary helium(He)-xenon (Xe) gas loop and Brayton Rotating Unit (BRU) that can generate 10s of kilowatts (kWe) of spacecraft electrical power during the mission coast phase eliminating the need for large PVAs. A single inflatable TransHab-type habitation module is also used with multiple vertical floors oriented radial to the MTV spin axis. The BNTR MTV's geometry-long and linear-is naturally compatible with AG operation. By rotating the vehicle about its center-of-mass (CM) and perpendicular to its flight vector at approximately 3.0 to 5.2 rpm, a centrifugal force and AG environment corresponding to approximately 0.38 to 1.0 g can be established to help maintain crew fitness out to Mars and back. Vehicles using NTP/ Bimodal Nuclear Thermal Propulsion (BNTP) can more readily accommodate the heavier payload mass and increased RCS propellant loading associated with AG operation, and can travel faster to and from Mars thereby reducing the crew's exposure to galactic cosmic radiation and solar flares. Mission scenario descriptions, key vehicle features and operational characteristics for each propulsion option are presented using the lift capability and payload volumes estimated for the Space Launch System (SLS)-1B and followon Heavy Lift Vehicle (HLV).

Nuclear Propulsion↗

DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗

DOE OSTI2023

Anisotropic magnetism of the Shastry-Sutherland lattice material BaNd 2 PtO 5

For this work, single crystals were grown and characterized to investigate the physical properties and magnetic ground state of BaNd 2 PtO 5 , a candidate to host physics of the Shastry-Sutherland model. Analysis of single crystal x-ray diffraction data yields an updated crystal structure, similar to the prior report but now in space group P4/mbm. Magnetization and specific heat measurements reveal an antiferromagnetic transition at T N =1.9K and a large magnetic anisotropy with in-plane magnetization much larger than out-of-plane magnetization. Single crystal neutron diffraction at zero field reveals a propagation vector of ($\frac{1}{2}$$\frac{1}{2}$$\frac{1}{2}$) for the magnetic ground state as compared to the ($\frac{1}{2}$$\frac{1}{2}$0) wave vector observed in the related body-centered material BaNd 2 ZnO 5 . The ground state is found to be a fully compensated antiferromagnet with diffraction data well fitted within the magnetic space group P S –1 (BNS setting #2.7). The ordered moment is mostly in the basal plane with nearest neighbors forming ferromagnetic dimers, however it also has a finite out-of-plane component unlike the related easy-plane materials BaNd 2 ZnO 5 and BaNd 2 ZnS 5 . Field-induced transitions are observed below T N when the field is applied within the basal plane, and in-plane anisotropy of the associated critical fields is observed. The magnetization is strongly impacted by misalignment of the field away from high symmetry directions. These results suggest complex magnetic structures may arise in the field-induced states, and they highlight the need for extreme care when studying this and related Shastry-Sutherland materials.

36 MATERIALS SCIENCE↗

NASA NTRS1973

Measurement of auroral Birkeland currents and energetic particle fluxes.

A rocket-borne experiment containing a vector magnetometer and a set of charged-particle detectors was launched from Poker Flat, Alaska, at 2217 LT on February 13, 1971, over a single auroral arc. Particle data in the energy range 0.5-20 kev, obtained during part of the flight, show a peak incident flux of 40 millions electrons per sq cm sec ster kev in the energy range 1-2 kev. The pitch angle distributions of these precipitating electrons were found to be fairly isotropic from 0 to 60 deg and to decrease from 60 to 90 deg. The results of the vector magnetometer indicate the existence of a system of Birkeland currents with magnitudes of 5 microamperes/sq m in the vicinity of the auroral arc. Two different models that fit the data are discussed.

Cloutier, P. A.↗

NASA NTRS1980

Equivalent source modeling of the main field using Magsat data

A graphic software package was developed to plot the dipole positions for a particular model on a world map. An arrow is drawn at each dipole in the direction of the horizontal magnetization vector, with length proportional to the horizontal magnitude. A contouring package represents the radical component of the magnetization vector. Based on quiet Magsat data, several dipole models were derived. These are being verified against each other and against MGST (6/80) data set as to goodness of fit and stability of solution. The spatial power spectra computed from spherical harmonic expansions of the dipole models is being analyzed relative to crustal and core content. The power spectra for one of the dipole models and a spherical harmonic model based on Magsat data through degree and order 23 is presented.

Source record↗

NASA NTRS1982

Three-dimensional inviscid analysis of the scramjet inlet flow field

A computer code has been developed to analyze the inviscid flow field in a supersonic combustion ramjet (scramjet) inlet. The code uses the three-dimensional Euler equations in full conservation form to describe the inlet flow. An algebraic numerical coordinate transformation is used to generate a set of boundary-fitted curvilinear coordinates. The governing equations are solved by a time-asymptotic, unsplit, two-step, finite-difference method. This method is highly efficient on the vector processing computers for which the current code is written. Detailed results are presented for two scramjet inlet configurations over a range of Mach numbers. The calculated results are compared with the available experimental and theoretical results.

Kumar, A.↗