Engineering Papers⌕ Search

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.

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
DOE OSTI2024

Neural network based emulation of galaxy power spectrum covariances: A reanalysis of BOSS DR12 data

We train neural networks to quickly generate redshift-space galaxy power spectrum covariances from a given parameter set (cosmology and galaxy bias). This covariance emulator utilizes a combination of traditional fully connected network layers and transformer architecture to accurately predict covariance matrices for the high redshift, north galactic cap sample of the BOSS DR12 galaxy catalog. We run simulated likelihood analyses with emulated and brute-force computed covariances, and we quantify the network’s performance via two different metrics: (1) difference in Χ 2 and (2) likelihood contours for simulated BOSS DR 12 analyses. We find that the emulator returns excellent results over a large parameter range. We then use our emulator to perform a reanalysis of the BOSS HighZ NGC galaxy power spectrum, and find that varying covariance with cosmology along with the model vector produces Ω m = $0.27⁢6$$^{+0.013}_{–0.015}$, H 0 = 70.2 ± 1.9 km/s/Mpc, and σ 8 = $0.67⁢4$$^{+0.058}_{–0.077}$. These constraints represent an average 0.46⁢σ shift in best-fit values and a 5% increase in constraining power compared to fixing the covariance matrix (Ω m = 0.293 ± 0.017, H 0 = 70.3 ± 2.0 km/s/Mpc, σ 8 = $0.70⁢2$$^{+0.063}_{–0.075}$). As a result, this work demonstrates that emulators for more complex cosmological quantities than second-order statistics can be trained over a wide parameter range at sufficiently high accuracy to be implemented in realistic likelihood analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Recent update of the RPLUS2D/3D codes

The development of the RPLUS2D/3D codes is summarized. These codes utilize LU algorithms to solve chemical non-equilibrium flows in a body-fitted coordinate system. The motivation behind the development of these codes is the need to numerically predict chemical non-equilibrium flows for the National AeroSpace Plane Program. Recent improvements include vectorization method, blocking algorithms for geometric flexibility, out-of-core storage for large-size problems, and an LU-SW/UP combination for CPU-time efficiency and solution quality.

Tsai, Y.-L. Peter↗
NASA NTRS2011

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole s accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamical model, we find a. > 0.92 (3 ). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle, and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk s low luminosity.

Gou, Lijun↗
NASA NTRS1978

Roll tracking effects of G-vector tilt and various types of motion washout

In a dogfight scenario, the task was to follow the target's roll angle while suppressing gust disturbances. All subjects adopted the same behavioral strategies in following the target while suppressing the gusts, and the MFP-fitted math model response was generally within one data symbol width. The results include the following: (1) comparisons of full roll motion (both with and without the spurious gravity tilt cue) with the static case. These motion cues help suppress disturbances with little net effect on the visual performance. Tilt cues were clearly used by the pilots but gave only small improvement in tracking errors. (2) The optimum washout (in terms of performance close to real world, similar behavioral parameters, significant motion attenuation (60 percent), and acceptable motion fidelity) was the combined attenuation and first-order washout. (3) Various trends in parameters across the motion conditions were apparent, and are discussed with respect to a comprehensive model for predicting adaptation to various roll motion cues.

Jex, H. R.↗
NASA NTRS1984

Modification to shock fitting program

A modified form of the Lepping - Argentiero single spacecraft, shock normal determination procedure is presented. The modified method incorporates a simple predictor-corrector algorithm which allows a faster convergence rate and the use of average values of the parameters for the starting vector.

Acuna, M. H.↗
NASA NTRS1979

Numerical solution of the Navier-Stokes equations for laminar, transonic flows

An implicit finite difference solution of the Navier-Stokes equations yielded time histories of the transonic laminar flow development about a circular cylinder and NACA-0018 airfoil. Reynolds numbers ranged from those corresponding to purely laminar flow to those corresponding to significant turbulence in the boundary layer. Body thermal conditions of an adiabatic wall and a specified body temperature were considered. Versatility in treating arbitrary bodies was incorporated by using numerically generated, body-fitted coordinate transformations. Solution of the simultaneous difference equations for the dependent variables was obtained using an accelerated Gauss-Seidel iterative scheme. Computational results are presented in the form of velocity vector fields, Mach number contours, aerodynamic coefficients, heat transfer rates at the body surface, and body temperature distributions. Truncation analyses of first and second derivative difference approximations resulted in general criteria for numerically generated coordinates so that flow near a body is more accurately represented.

Turner, L., III↗
NASA NTRS1986

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗
NASA NTRS1992

Interactive Spectral Analysis and Computation (ISAAC)

Isaac is a task in the NSO external package for IRAF. A descendant of a FORTRAN program written to analyze data from a Fourier transform spectrometer, the current implementation has been generalized sufficiently to make it useful for general spectral analysis and other one dimensional data analysis tasks. The user interface for Isaac is implemented as an interpreted mini-language containing a powerful, programmable vector calculator. Built-in commands provide much of the functionality needed to produce accurate line lists from input spectra. These built-in functions include automated spectral line finding, least squares fitting of Voigt profiles to spectral lines including equality constraints, various filters including an optimal filter construction tool, continuum fitting, and various I/O functions.

Lytle, D. M.↗
DOE OSTI2026

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗
NASA NTRS1984

Modification to shock fitting program

A modified form of the Lepping-Argentiero single spacecraft, shock normal determination procedure is presented. The modified method incorporates a simple predictor-corrector algorithm which allows a faster convergence rate and the use of average values of the parameters for the starting vector.

Acuna, M. H.↗
DOE OSTI2023

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗
NASA NTRS1985

Evaluation of data obtained from atmospheric laser Doppler velocimeter

The Doppler lidar velocimeter produces a variety of kinds of data. Besides the line of sight velocity components, there are the reflected amplitudes and the Doppler widths. Computer graphics software for displaying these data was produced. Different methods of presentation are needed for the various items. A picture was produced as pictures are often the best way to gain understanding. The individual lidar scans form a crosshatch pattern. Thus one-dimensional measurements fit together to form a two-dimensional whole. A pair of velocity measurements at a point combine to form a wind measurement with direction and magnitude. This gives a forest of wind vectors. The goal is to find a recognizable pattern to these trees. Often it is necessary to show only part of the information. That is, show only the wind direction not its magnitude or reduce the wind to streamlines of air flow. In other cases data are reduced to contour plots. Just enough contour lines are included to show the picture described.

Keck, T. S.↗
DOE OSTI2020

Has AMS-02 observed two-component dark matter?

There is convincing observational evidence for an increasing cosmic-ray positron-to-electron ratio at energies larger than ~10 GeV, at odds with expectations from secondary positron production. The most recent AMS-02 data exhibit an interesting spectral feature consisting of a bump at an energy around 300 GeV followed by a drop around ~800 GeV. A possible explanation for the most recent data is that the excess positron flux originates from decaying dark matter. Here, we show that models consisting of two dark matter particle species contributing equally to the global cosmological dark matter density provide good fits to the data. The favored models, with a best-fit with χ 2 / d . o . f ~ 2, consist of a first species weighing 750 GeV decaying with a lifetime τ χ ~ 10 26 s to τ lepton pairs (or to a pair of vector bosons subsequently decaying to a τ pair each), and a second species with a mass around 2.3 TeV decaying to μ lepton pairs. We provide a few possible concrete realizations for this scenario.

Physics↗
DOE OSTI2024

A New Vehicle-to-Vehicle Communication Technology Based on Binary Light Code

The proliferation of Connected and Autonomous Vehicles (CAVs) has necessitated the development of efficient, reliable communication methods between vehicles and their surroundings. This study presents a novel methodology, the VECTOR system, that addresses this need by converting dynamic vehicle data into a binary code and displaying it on an LED panel. The system involves a three-step process of data collection, polynomial fitting of velocity data, and encoding the polynomial parameters into binary form using a Cyclic Redundancy Check (CRC). The study also explores the practical application of this system by conducting experiments with modified Lincoln MKZ hybrids, demonstrating the feasibility and efficiency of this approach. The paper presents both the methodology and results of these experiments, further expanding upon the potential applications and implications for CAV technology. By comparing the final accuracy, we find this approach can achieve 80% of the original velocity data.

Ma, Ke↗
NASA NTRS1982

Zonal harmonic model of Saturn's magnetic field from Voyager 1 and 2 observations

An analysis of the magnetic field of Saturn is presented which takes into account both the Voyager 1 and 2 vector magnetic field observations. The analysis is based on the traditional spherical harmonic expansion of a scale potential to derive the magnetic field within 8 Saturn radii. A third-order zonal harmonic model fitted to Voyager 1 and 2 observations is found to be capable of predicting the magnetic field characteristics at one encounter based on those observed at another, unlike models including dipole and quadrupole terms only. The third-order model is noted to lead to significantly enhanced polar surface field intensities with respect to dipole models, and probably represents the axisymmetric part of a complex dynamo field.

Connerney, J. E. P.↗
NASA NTRS1982

Sensitivity of constrained linear inversions to the selection of the Lagrange multiplier

The influence of the choice of the Lagrange multiplier on constrained linear inversions is explored, with reference made to applications in inferring the columnar aerosol size distributions from spectral aerosol optical depth measurements. A range of the Lagrange multiplier is examined to find all positive solutions for the solution vector, which represents modifying factors to the assumed form of the size distribution. An iterative method is devised to constrain the calculations to consideration of only positive quantities and a requirement that the regression fit to data be consistent with measurement errors. The determination of the variances and covariances is formulated and applied to existing data sets for optical depth. Variances in the solution are found to be large for particle radii when the information content of the data is small.

King, M. D.↗
NASA NTRS1988

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗
NASA NTRS1988

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗