Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test vector”

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 1,081 records · Page 60

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗

GRB Progenitor Classification from Gamma-Ray Burst Prompt and Afterglow Observations

Using an established classification technique, we leverage standard observations and analyses to predict the progenitors of gamma-ray bursts (GRBs). This technique, utilizing support vector machine (SVM) statistics, provides a more nuanced prediction than the previous two-component Gaussian mixture in duration of the prompt gamma-ray emission. Based on further covariance testing from Fermi-GBM, Swift-BAT, and Swift-XRT data, we find that our classification based only on prompt emission properties gives perspective on the recent evidence that mergers and collapsars exist in both “long” and “short” GRB populations.

P Nuessle↗

Using support vector machines to improve elemental ion identification in macromolecular crystal structures

A method to automatically identify possible elemental ions in X-ray crystal structures has been extended to use support vector machine (SVM) classifiers trained on selected structures in the PDB, with significantly improved sensitivity over manually encoded heuristics. In the process of macromolecular model building, crystallographers must examine electron density for isolated atoms and differentiate sites containing structured solvent molecules from those containing elemental ions. This task requires specific knowledge of metal-binding chemistry and scattering properties and is prone to error. A method has previously been described to identify ions based on manually chosen criteria for a number of elements. Here, the use of support vector machines (SVMs) to automatically classify isolated atoms as either solvent or one of various ions is described. Two data sets of protein crystal structures, one containing manually curated structures deposited with anomalous diffraction data and another with automatically filtered, high-resolution structures, were constructed. On the manually curated data set, an SVM classifier was able to distinguish calcium from manganese, zinc, iron and nickel, as well as all five of these ions from water molecules, with a high degree of accuracy. Additionally, SVMs trained on the automatically curated set of high-resolution structures were able to successfully classify most common elemental ions in an independent validation test set. This method is readily extensible to other elemental ions and can also be used in conjunction with previous methods based on a priori expectations of the chemical environment and X-ray scattering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiple Shaker Placement for Ground Vibration Test of X-59 Aircraft Using Topology Optimization

A multiple shaker placement methodology is developed and tested using a topology optimization technique. Current multiple shaker placement methodology requires optimum accelerometer placement and optimum single-shaker placement techniques. The proposed methodology is tested using a finite element model of the X-59 Low Boom Flight Demonstrator aircraft. The effective independence and the driving point acceleration transfer function (DPATF) methods are used for the accelerometer placement study. In this study, four shakers are used to excite each mode more effectively during the ground vibration test; all the modes of interest thus are separated into four groups. Each shaker takes care of a separate group of modes. Grouping the modes of interest is performed utilizing topology optimization. The number of modes for each group therefore will be automatically decided during grouping. For each group of modes, perform the following two steps to determine optimal location of four shakers: 1) At each accelerometer location, compare the magnitude of DPATF values at natural frequencies, select the minimum value, and make a vector with these minimum values of the DPATF magnitudes for each group; and 2) Select the degrees of freedom corresponding to the maximum value of this vector. The objective function value is the maximum value of the vector with minimum value of the magnitude of the superposed acceleration transfer function. This objective function value is maximized by changing the modes for each group. Forty accelerometers are enough to have good correlation between mode shapes obtained from the reduced order model and the simulated ground vibration test.

Optimization Analysis↗

EGS Collab Experiment 1: In-situ observation of pre-, co- and post-seismic shear slip preceding hydraulic fracturing

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and a transient afterslip phase. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment. This submission includes the pressure and deformation data recorded by the SIMFIP probe during the first injection at the 164 ft (50 m) notch of borehole E1-I. The injection was performed on on 05/22/2018 as part of Experiment 1 of the EGS Collab project. This data accompanies a manuscript submitted to GRL, linked in this submission.

15 GEOTHERMAL ENERGY↗

Aeroelastic Model Structure Computation for Envelope Expansion

Structure detection is a procedure for selecting a subset of candidate terms, from a full model description, that best describes the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modelling may be of critical importance in the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion which may save significant development time and costs. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of nonlinear aeroelastic systems. The LASSO minimises the residual sum of squares by the addition of an l(sub 1) penalty term on the parameter vector of the traditional 2 minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudolinear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. Applicability of this technique for model structure computation for the F/A-18 Active Aeroelastic Wing using flight test data is shown for several flight conditions (Mach numbers) by identifying a parsimonious system description with a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

Adaptive Learning for Reliability Analysis using Support Vector Machines

A novel algorithm is presented for adaptive learning of an unknown function that separates two regions of a domain.In the context of reliability analysis these two regions represent the failure domain, where a set of constraints or requirements are violated, and a safe domain where they are satisfied. The Limit State Function (LSF) separates these two regions. Evaluating the constraints for a given parameter point requires the evaluation of a computational model that may well be expensive. For this reason we wish to construct a meta-model that can estimate the LSFas accurately as possible, using only a limited amount of training data. This work presents an adaptive strategy employing a Support Vector Machine (SVM) as a meta-model to provide a semi-algebraic approximation of the LSF.We describe an optimization process that is used to select informative parameter points to add to training data at each iteration to improve the accuracy of this approximation. A formulation is introduced for bounding the predictions of the meta-model; in this way we seek to incorporate this aspect of Gaussian Process Models (GPMs) within anSVM meta-model. Finally, we apply our algorithm to two benchmark test cases, demonstrating performance that is comparable with, if not superior, to a standard technique for reliability analysis that employs GPMs

Adaptive learning↗

Seismic moment tensor classification using elliptical distribution functions on the hypersphere

Discrimination of underground explosions from naturally occurring earthquakes and other anthropogenic sources is one of the fundamental challenges of nuclear explosion monitoring. In an operational setting, the number of events that can be thoroughly investigated by analysts is limited by available resources. The capability to rapidly screen out events that can be robustly identified as not being explosions is, therefore, of great potential benefit. Nevertheless, possible mis-classification of explosions as earthquakes currently limits the use of screening methods for verification of test-ban treaties. Moment tensors provide a physics-based classification tool for the characterization of different seismic sources and have enabled the advent of new techniques for discriminating between earthquakes and explosions. Following normalization and projection of their six-degree vectors onto the hypersphere, existing screening approaches use spherically symmetric metrics to determine whether any new moment tensor may have been an explosion. Here, we show that populations of moment tensors for both earthquakes and explosions are anisotropically distributed on the hypersphere. Distributions possessing elliptical symmetry, such as the scaled von Mises–Fisher distribution, therefore provide a better description of these populations than the existing spherically symmetric models. We describe a method that uses these elliptical distributions in combination with a Bayesian classifier to achieve successful classification rates of 99 per cent for explosions and 98 per cent for earthquakes using existing catalogues of events from the western United States. The 1983 May 5 Crowdie underground nuclear test and 2018 July 20 DAG-1 deep-borehole chemical explosion are the only two explosions out of 140 that are incorrectly classified. Application of the method to the 2006–2017 nuclear tests in the Democratic People’s Republic of Korea yields 100 per cent identification rates and we provide a simple routine MTid for general usage. The approach provides a means to rapidly assess the likelihood of an event being an explosion and can be built into monitoring workflows that rely on simultaneously assessing multiple different discrimination metrics.

58 GEOSCIENCES↗

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validating time-distance helioseismology with realistic quiet-sun simulations

Linear time-distance helioseismic inversions are carried out for vector flow velocities using travel times measured from two ∼100{sup 2} Mm{sup 2} × 20 Mm realistic magnetohydrodynamic quiet-Sun simulations of about 20 hr. The goal is to test current seismic methods on these state-of-the-art simulations. Using recent three-dimensional inversion schemes, we find that inverted horizontal flow maps correlate well with the simulations in the upper ∼3 Mm of the domains for several filtering schemes, including phase-speed, ridge, and combined phase-speed and ridge measurements. In several cases, however, the velocity amplitudes from the inversions severely underestimate those of the simulations, possibly indicating nonlinearity of the forward problem. We also find that, while near-surface inversions of the vertical velocities are best using phase-speed filters, in almost all other example cases these flows are irretrievable due to noise, suggesting a need for statistical averaging to obtain better inferences.

79 ASTRONOMY AND ASTROPHYSICS↗

Birefringence tests of gravity with multimessenger binaries

Extensions to General Relativity (GR) allow the polarization of gravitational waves (GW) from astrophysical sources to suffer from amplitude and velocity birefringence, which respectively induce changes in the ellipticity and orientation of the polarization tensor. Here, we introduce a multi-messenger approach to test this polarization behavior of GWs during their cosmological propagation using binary sources, for which the initial polarization is determined by the inclination and orientation angles of the orbital angular momentum vector with respect to the line of sight. In particular, we use spatially resolved radio imaging of the jet from a binary neutron star (BNS) merger to constrain the orientation angle and hence the emitted polarization orientation of the GW signal at the site of the merger, and compare to that observed on Earth by GW detectors. For GW170817, using past measurements of the inclination angle, we constrain the deviation from GR due to amplitude birefringence to κ A = $-0.12^{+0.60}_ {-0.61}$, while the velocity birefringence parameter κ V remains unconstrained. The inability to constrain κ V is due to the low amplitude of GW170817 in the Virgo detector, and measurements of the polarization orientation require information from a combination of multiple detectors with different alignments. For this reason, we also mock future BNS mergers with resolved afterglow proper motion and project that κ V could be constrained to a precision of 5 rad (corresponding to an angular shift of the GW polarization of δ φV ≈ 0.2 rad for a BNS at 100 Mpc) by a future network of third-generation ground-based GW detectors such as Cosmic Explorer and the radio High Sensitivity Array. Crucially, this velocity birefringence effect cannot be constrained with dark binary mergers as it requires polarization information at the emission time, which can be provided only by electromagnetic emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Extreme data compression for the CMB

We apply the Karhunen-Loéve methods to cosmic microwave background (CMB) data sets, and show that we can recover the input cosmology and obtain the marginalized likelihoods in Λ cold dark matter cosmologies in under a minute, much faster than Markov chain Monte Carlo methods. This is achieved by forming a linear combination of the power spectra at each multipole l, and solving a system of simultaneous equations such that the Fisher matrix is locally unchanged. Instead of carrying out a full likelihood evaluation over the whole parameter space, we need evaluate the likelihood only for the parameter of interest, with the data compression effectively marginalizing over all other parameters. The weighting vectors contain insight about the physical effects of the parameters on the CMB anisotropy power spectrum C l . The shape and amplitude of these vectors give an intuitive feel for the physics of the CMB, the sensitivity of the observed spectrum to cosmological parameters, and the relative sensitivity of different experiments to cosmological parameters. We test this method on exact theory C l as well as on a Wilkinson Microwave Anisotropy Probe (WMAP)-like CMB data set generated from a random realization of a fiducial cosmology, comparing the compression results to those from a full likelihood analysis using CosmoMC. Furthermore, after showing that the method works, we apply it to the temperature power spectrum from the WMAP seven-year data release, and discuss the successes and limitations of our method as applied to a real data set.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Prophylactic treatment with PEGylated bovine IFNλ3 effectively bridges the gap in vaccine-induced immunity against FMD in cattle

Foot-and-mouth disease (FMD) is a vesicular disease of cloven-hoofed animals with devastating economic implications. The current FMD vaccine, routinely used in enzootic countries, requires at least 7 days to induce protection. However, FMD vaccination is typically not recommended for use in non-enzootic areas, underscoring the need to develop new fast-acting therapies for FMD control during outbreaks. Interferons (IFNs) are among the immune system’s first line of defense against viral infections. Bovine type III IFN delivered by a replication defective adenovirus (Ad) vector has effectively blocked FMD in cattle. However, the limited duration of protection—usually only 1–3 days post-treatment (dpt)—diminishes its utility as a field therapeutic. Here, we test whether polyethylene glycosylation (PEGylation) of recombinant bovine IFNλ3 (PEGboIFNλ3) can extend the duration of IFN-induced prevention of FMDV infection in both vaccinated and unvaccinated cattle. We treated groups of heifers with PEGboIFNλ3 alone or in combination with an adenovirus-based FMD O1Manisa vaccine (Adt-O1M) at either 3 or 5 days prior to challenge with homologous wild type FMDV. We found that pre-treatment with PEGboIFNλ3 was highly effective at preventing clinical FMD when administered at either time point, with or without co-administration of Adt-O1M vaccine. PEGboIFNλ3 protein was detectable systemically for >10 days and antiviral activity for 4 days following administration. Furthermore, in combination with Adt-O1M vaccine, we observed a strong induction of FMDV-specific IFNγ+ T cell response, demonstrating its adjuvanticity when co-administered with a vaccine. Our results demonstrate the promise of this modified IFN as a pre-exposure prophylactic therapy for use in emergency outbreak scenarios.

biotherapeutics↗

Casting the Coronal Magnetic Field Reconstruction Tools in 3D Using the MHD Bifrost Model

Quantifying the coronal magnetic field remains a central problem in solar physics. Nowadays, the coronal magnetic field is often modeled using nonlinear force-free field (NLFFF) reconstructions, whose accuracy has not yet been comprehensively assessed. Here we perform a detailed casting of the NLFFF reconstruction tools, such as π -disambiguation, photospheric field preprocessing, and volume reconstruction methods, using a 3D snapshot of the publicly available full-fledged radiative MHD model. Specifically, from the MHD model, we know the magnetic field vector in the entire 3D domain, which enables us to perform a “voxel-by-voxel” comparison of the restored and the true magnetic fields in the 3D model volume. Our tests show that the available π -disambiguation methods often fail in the quiet-Sun areas dominated by small-scale magnetic elements, while they work well in the active region (AR) photosphere and (even better) chromosphere. The preprocessing of the photospheric magnetic field, although it does produce a more force-free boundary condition, also results in some effective “elevation” of the magnetic field components. This “elevation” height is different for the longitudinal and transverse components, which results in a systematic error in absolute heights in the reconstructed magnetic data cube. The extrapolations performed starting from the actual AR photospheric magnetogram are free from this systematic error, while other metrics are comparable with those for extrapolations from the preprocessed magnetograms. This finding favors the use of extrapolations from the original photospheric magnetogram without preprocessing. Our tests further suggest that extrapolations from a force-free chromospheric boundary produce measurably better results than those from a photospheric boundary.

79 ASTRONOMY AND ASTROPHYSICS↗

Vector plotting as an indication of the approach to flutter

A binary flexure-torsion analysis was made to check theoretically a method for predicting flutter which depends on plotting vectorially the amplitudes of response relative to the exciting force and extracting the relevant damping rate. The results of this calculation are given in graphs both of the vector plots themselves and of the estimated damping rate against forward speed. The estimated damping rates are compared with calculated values. The method has the advantage that in a flight flutter test damping can be estimated from continuous excitation records: the method is an extension of the Kennedy and Pancu technique used in ground resonance testing.

Broadbent, E. G.↗

Efficient reanalysis of structures by a direct modification method

A procedure for the local stiffness modifications of large structures is described. It enables structural modifications without an a priori definition of the changes in the original structure and without loss of efficiency due to multiple loading conditions. The solution procedure, implemented in NASTRAN, involved the decomposed stiffness matrix and the displacement vectors of the original structure. It solves the modified structure exactly, irrespective of the magnitude of the stiffness changes. In order to investigate the efficiency of the present procedure and to test its applicability within a design environment, several real and large structures were solved. The results of the efficiency studies indicate that the break-even point of the procedure varies between 8% and 60% stiffness modifications, depending upon the structure's characteristics and the options employed.

Raibstein, A. I.↗

Characterization of the non axial thrust generated by large solid propellant rocket motors in three axis stabilized ascent

Nonaxial thrusts produced by solid rocket motors during three-axis stabilized attitude control have been determined from ascent experience on twenty three Burner II, Burner IIA and Block 5D-1 upper stage vehicles. A data base representing four different rocket motor designs (three spherical and one extended spherical) totaling twenty five three-axis stabilized firings is generated. Solid rocket motor time-varying resultant and lateral side force vector magnitudes, directions and total impulses, and roll torque couple magnitudes, directions, and total impulses are tabulated in the appendix. Population means and three sigma deviations are plotted. Existing applicable ground test side force and roll torque magnitudes and total impulses are evaluated and compared to the above experience data base. Within the spherical motor population, the selected AEDC ground test data consistently underestimated experienced motor side forces, roll torques and total impulses. Within the extended spherical motor population, the selected AEDC test data predicted experienced motor side forces, roll torques, and total impulses, with surprising accuracy considering the very small size of the test and experience populations.

Kosmann, W. J.↗