Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equivalence Method”

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,117 records · Page 62

Response of joint dominated space structures

An efficient linearization method is presented for calculating the transient response of nonlinear systems due to initial disturbances. The method is an extension of the describing function approach in which the steady state response of the system is calculated by representing the nonlinear element, typically joints in the case of space structures, by impedances which are functions of the amplitude of response. Thus, the problem of solving the differential equation for the steady state response becomes one of solving a set of nonlinear algebraic equations involving the steady state amplitudes and phases of the system. It is shown that for the transient case the steady state impedances can be averaged over the range of responses in order to provide equivalent values of stiffness and damping that, for a given set of initial displacements, may be treated as being constant for purposes of calculating system response. Single degree of freedom system are used to demonstrate the method and to develop an approach for optimizing the joint's characteristics so as to minimize transient response times. The use of this method for response estimation and optimization in multiple degree of freedom systems is investigated.

Source record↗

Development of a fractional-step method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A fractional step method is developed for solving the time-dependent three-dimensional incompressible Navier-Stokes equations in generalized coordinate systems. The primitive variable formulation uses the pressure, defined at the center of the computational cell, and the volume fluxes across the faces of the cells as the dependent variables, instead of the Cartesian components of the velocity. This choice is equivalent to using the contravariant velocity components in a staggered grid multiplied by the volume of the computational cell. The governing equations are discretized by finite volumes using a staggered mesh system. The solution of the continuity equation is decoupled from the momentum equations by a fractional step method which enforces mass conservation by solving a Poisson equation. This procedure, combined with the consistent approximations of the geometric quantities, is done to satisfy the discretized mass conservation equation to machine accuracy, as well as to gain the favorable convergence properties of the Poisson solver. The momentum equations are solved by an approximate factorization method, and a novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two- and three-dimensional laminar test cases are computed and compared with other numerical and experimental results to validate the solution method. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗

Initial Development of a Metric to Describe the Level of Safety Associated with Piloting an Aircraft with Synthetic Vision Systems (SVS) Displays

Synthetic Vision Systems (SVS) displays provide pilots with a continuous view of terrain combined with integrated guidance symbology in an effort to increase situation awareness (SA) and decrease workload during operations in Instrument Meteorological Conditions (IMC). It is hypothesized that SVS displays can replicate the safety and operational flexibility of flight in Visual Meteorological Conditions (VMC), regardless of actual out-the-window (OTW) visibility or time of day. Significant progress has been made towards evolving SVS displays as well as demonstrating their ability to increase SA compared to conventional avionics in a variety of conditions. While a substantial amount of data has been accumulated demonstrating the capabilities of SVS displays, the ability of SVS to replicate the safety and operational flexibility of VMC flight performance in all visibility conditions is unknown to any specific degree. In order to more fully quantify the relationship of flight operations in IMC with SVS displays to conventional operations conducted in VMC, a fundamental comparison to current day general aviation (GA) flight instruments was warranted. Such a comparison could begin to establish the extent to which SVS display concepts are capable of maintaining an "equivalent level of safety" with the round dials they could one day replace, for both current and future operations. A combination of subjective and objective data measures were used to quantify the relationship between selected components of safety that are associated with flying an approach. Four information display methods ranging from a "round dials" baseline through a fully integrated SVS package that includes terrain, pathway based guidance, and a strategic navigation display, were investigated in this high fidelity simulation experiment. In addition, a broad spectrum of pilots, representative of the GA population, were employed for testing in an attempt to enable greater application of the results and determine if "equivalent levels of safety" are achievable through the incorporation of SVS technology regardless of a pilot's flight experience.

Bartolone, Anthony P.↗

Lithium-Ion Battery Diagnostics Using Electrochemical Impedance via Machine-Learning

Diagnosing battery states such as health, state-of-charge, or temperature is crucial for ensuring the safety and reliability of electrochemical energy storage systems. While some states, such as temperature, may be measured using cheap sensors, accurate diagnosis of battery health metrics usually requires time-consuming performance measurements, making them infeasible for use in real-world operation. These health metrics can be measured during lab-testing and then estimated on-line using predictive life models or via state observer algorithms such as Kalman filters, but these predictive methods should be supplemented by actual measurement of battery health whenever possible to ensure reliability. Rapid measurement of battery health may be done by various types of fast diagnostic techniques such as electrochemical impedance spectroscopy (EIS), which can be performed in only a few minutes and require only a fraction of the energy and power needed for a full charge and discharge measurement. But there is a substantial challenge for estimating battery health using EIS data, as EIS is sensitive to cell temperature, state-of-charge, current, and resting time in addition to health. Thus, utilizing EIS data to predict battery capacity requires correcting for all these additional variables, a task that is extremely difficult to handle analytically. This talk utilizes machine-learning methods to estimate the effectiveness of battery capacity prediction from EIS data, leveraging a data set of hundreds of EIS measurements recorded at varying temperature and state-of-charge throughout a 500-day aging study of 32 commercial, large-format NMC-Graphite lithium-ion batteries. Using EIS as input to machine-learning models is complicated by the nonlinear response of impedance to battery health, temperature, and state-of-charge, as well as the collinearity between the impedance response at neighboring frequencies, which can easily lead to overfit models. To train robust models, features from EIS data need to be extracted from the data or some subset of critical frequencies selected. Many approaches for extracting and selecting features from EIS data from electrochemical analysis and machine-learning fields were identified for analysis: using the entire raw spectra; selection of one, two, or many frequencies from the entire spectra; selecting interesting points from the EIS measurement using domain knowledge; fitting EIS with an equivalent-circuit model; calculating statistics on the raw impedance values; and reducing the dimensionality of the data using unsupervised linear (principal component analysis) and non-linear (uniform manifold approximation and projection) methods. These approaches were rigorously compared using a machine-learning pipeline approach, training linear, Gaussian process, and random forest regression models and quantifying performance using cross-validation as well as a held-out test set. An artificial neural network model trained on the raw spectra was also tested. Promising pipelines were fine-tuned via Bayesian hyperparameter optimization using cross-validation loss and training with class-specific weights to counter data set imbalance. The most reliable method for utilizing impedance in this work was the selection of two optimal frequencies through an exhaustive search, resulting in about 2% mean absolute error on test data for both Gaussian process and random forest model architectures. Interrogation of a variety of models reveals critical frequencies of 100 Hz and 103 Hz for this data set, though the optimal set of frequencies is not necessarily intuitive, i.e., the best performing models are not simply those that use impedance at frequencies that have the highest correlation to the relative discharge capacity. The best performing model is an ensemble model, which is able to predict battery capacity with 1.9% mean absolute error for unseen cells using impedance recorded at a variety of temperatures and states-of-charge.

battery↗

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling of oxygen-methane reacting coaxial jets using x-ray fluorescence to measure mixture fraction

Stoichiometric mixing length L s of reacting coaxial jet flames is a critical scaling parameter for liquid rocket engine combustors. Previous studies have shown that L s for shear coaxial flames can be scaled like their non-reacting counterparts using a non-dimensional momentum flux ratio J. In addition, stoichiometric mixing lengths of reacting and nonreacting coaxial jets collapse upon a single line by altering J using an effective outer flow gas density. This effective density is calculated from a modified version of the equivalence principle, originally developed by Tacina and Dahm [1, 2] and accounts for the effects of heat release on mixing. However, previous studies also required a second nonphysical scaling constant S c for the reacting jets, which is not predicted by the equivalence principle [3]. It was originally hypothesized that S c is attributed to the limitation of hydroxyl (OH) planar laser-induced fluorescence, which only infers L s . Direct quantitative measurement of conserved scalar fields using conventional optical diagnostics is difficult due to the lack of a tracer that easily fluoresces, survives high temperature oxygen flames, and is not dominated by quenching effects. To measure a conserved scalar field, this work implements x-ray fluorescence of Kr and Ar tracers to obtain quantitative mixture fraction fields. From these mixture fraction fields, stoichiometric mixing lengths for two CH 4 /O 2 flames are calculated and scaled against nonreacting coaxial mixing lengths using the equivalence principle. By directly measuring the stoichiometric mixing length, it is established that the additional constant is a byproduct of the OH measurement technique and the equivalence principle fully captures the scaling. Finally, comparison with high-fidelity simulation of the flame further supports this conclusion. In addition to further strengthening this scaling method, this work represents the first use of x-ray fluorescence to make quantitative conserved scalar measurements in turbulent flames.

33 ADVANCED PROPULSION SYSTEMS↗

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗

Allocating Reserves in Active Distribution Systems for Tertiary Frequency Regulation

This paper proposes a cooperative game theory-based approach for reserve optimization to enable distributed energy resources (DERs) participate in tertiary frequency regulation. Tertiary frequency regulation schemes ensure that reserve requirements of primary and secondary frequency regulation are fulfilled with a minimum cost. While the available reserve from a single distribution system may not suffice tertiary frequency regulation, stacked reserve from several distribution systems can enable them participate in tertiary frequency regulation at scale. In this paper, a two-stage strategy is proposed to effectively and precisely allocate spinning reserve requirement from each DER in distribution systems. In the first stage, two types of characteristic functions are computed: worthiness index (WI) and power loss reduction (PLR). In the second stage, the equivalent Shapley values are computed based on the characteristic functions, which are used to determine distribution factors for reserve allocation among DERs. The effectiveness of the proposed method for allocating reserves among DERs is demonstrated through several case studies on modified versions of the IEEE 13-node and 33-node distribution systems.

Gautam, Mukesh↗

Sulfur Pellets Responses to a Bare and Steel Reflected Pulse of the Oak Ridge National Laboratory Health Physics Research Reactor

The experiments analyzed in this report were conducted at the Health Physics Research Reactor (HPRR), also known as the $\textit{Fast Burst Reactor}$. The reactor was designed and built at Oak Ridge National Laboratory (ORNL) in 1961. The HPRR was an unmoderated, unshielded fast reactor that used highly enriched uranium and molybdenum alloy as fuel. The reactor was initially sent to the Nevada Test Site in 1962, where it was used to evaluate radiation doses received as a result of the Hiroshima and Nagasaki bombings during World War II. A few years later, the reactor was sent back to ORNL to be part of the Dosimetry Application Research (DOSAR) facility shown in Figure 1, which included a reactor building shown on the left (west) of the picture and a control and laboratory building in the upper right corner (northeast). The critical assembly was used for numerous technical studies, including systems calibration, dosimetry, radiobiology of plants and animals, testing of radiation alarms, as well as teaching and training in radiation dosimetry and nuclear engineering. Between 1963 and 1987, the HPRR was operated for thousands of hours, achieving criticality close to 10,000 times and motivating many publications. The HPRR was decommissioned in 1987. The goal of this effort was to use historical data from operation of the HPRR to create a criticality accident alarm system (CAAS) benchmark to be included in the $\textit{International Handbook of Evaluated Criticality Safety Benchmark Experiments}$ (ICSBEP Handbook). A thorough inspection was performed of all available documentation and information available. The most promising experiments that were selected for evaluation were those described in the 1987 ORNL report entitled $\textit{Health Physics Research Reactor Reference Dosimetry}$, ORNL-6240. The report includes reference dosimetry results of the shielded and unshielded configurations of the HPRR after burst operations. Because of changes to the reactor positioning and storage systems that were made in 1985, the previous dosimetry reports became obsolete, and the newly designed experiments were needed to create the HPRR’s adjusted dosimetry data. The various results reported in ORNL-6240 include reference doses and dose equivalents from different conventions at different distances and elevations as determined using the detected neutron fluence and conversion factors. The HPRR neutron fluence was obtained through different methods, including sulfur pellet analysis and threshold detector unit data. Information about the HPRR spectrum was also obtained through Bonner sphere measurements. This benchmark is focused on a part of the measured sulfur fluences reported in Appendix H of ORNL-6240. Standard commercial sulfur pellets were placed at different distances from the HPRR centerline during burst operation and were activated due to the 32 S(n,p) 32 P reaction. The resulting 32 P activity was then measured and the information about the corresponding sulfur fluence and/or neutron dose could be extracted. Many of those measurements have 7 been performed with the HPRR in its bare configuration or with different shields (combinations of Lucite, concrete, steel). All the necessary, precise information about material and/or dimensions of the different shields was not found, so it was decided to focus only on the unshielded and steel-shielded configurations to minimize the benchmark uncertainty. A total of 31 cases (24 unshielded and 7 shielded cases at different positions) of sulfur fluence were selected before evaluation to develop the benchmark.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Optimization of a traveling wave superconducting radiofrequency cavity for upgrading the International Linear Collider

The standing wave TESLA niobium-based superconducting radio frequency structure is limited to an accelerating gradient of about 50 MV/m by the critical rf magnetic field. To break through this barrier, we explore the option of niobium-based traveling wave (TW) structures. Optimization of TW structures was done considering experimentally known limiting electric and magnetic fields. It is shown that a TW structure can have an accelerating gradient above 70 MeV/m that is about 1.5 times higher than contemporary standing wave structures with the same critical magnetic field. The other benefit of TW structures shown is R/Q about 2 times higher than the TESLA structure that reduces the dynamic heat load by a factor of 2. A method is proposed how to make TW structures multipactor-free. Some design proposals are offered to facilitate fabrication. Further increase of the real-estate gradient (equivalent to 80 MV/m active gradient) is also possible by increasing the length of the accelerating structure because of higher group velocity and cell-to-cell coupling. Realization of this work opens paths to International Linear Collider energy upgrades beyond 1 to 3 TeV in competition with CLIC. The paper will discuss corresponding opportunities and challenges.

43 PARTICLE ACCELERATORS↗

Sparse Data Machine Learning Integration with Theory, Experiment and Uncertainty Quantification: Process-Structure-Property-Performance of Friction Deformation Processing

Computer vision and deep learning tools that advance the ability to establish processing-structure-property-performance (PSPP) relations are presented. The Bayesian binning method for image segmentation enables quantitative analysis of microstructural features in an automated way, while the analysis of shapes and relative orientation of these features reveals local deformation maps indicative of both, material flow and residual stresses due to materials processing. The deep learning method leads to the previous knowledge agnostic mapping of empirically observed microstructural zones in friction stir welding (FSW) process and synthetic microstructure generation capability that is statistically equivalent to experimentally collected data.

97 MATHEMATICS AND COMPUTING↗

Quasicoordinate equations for flexible spacecraft

Equations of motion equivalent to Lagrange's equations for generalized coordinates can be formulated in quasi-coordinate form. An analysis of the resulting equations is conducted and some alternative methods are considered. It is found that Lagrange's quasi-coordinate equations produce precisely the same equations of motion for any material system as would energy from the combination of Lagrange's generalized coordinate equations and Euler's rotational equation in vector form.

Likins, P. W.↗

Penalty function finite element analysis of steady viscous incompressible flow in rotating coordinates

Finite element methods for incompressible viscous flow in turbomachines have not been presented in the literature previously. This paper develops a penalty function primitive variable method including Coriolis and centrifugal force terms for steady flow in a rotating coordinate system. Simplex elements are used with the result of solution times comparable to equivalent finite different solutions. Example cases considered are Couette flow, Poiseuile flow, flow over a step and flow in a rotating channel. Both laminar and turbulent flows are discussed. The accuracy of computed solutions compares well with theoretical solutions and experimenal measurements.

Rosen, M. C.↗

Asymptotic high frequency analysis of the electromagnetic backscattering from an inlet model consisting of piecewise linearly tapered sections

Electromagnetic backscattering from an open ended three dimensional inlet model is analyzed and computed patterns are compared with results of experimental measurements. The model is comprised of two sections. The first section consists of a linearly tapered waveguide with a rectangular opening at one end and the other end is connected to the second section which is a uniform rectangular waveguide with a planar perfectly conducting termination. The model is electrically large so that many propagating modes are excited. The method of analysis contains conventional aperture integration and modal techniques combined with high frequency techniques, which employ concepts such as modal rays, geometrical theory of diffraction and equivalent currents. For the cases considered, it is shown that only a few of the many propagating modes contribute appreciably to the backscattered field. These modes are selected according to their modal ray angle directions.

Altintas, A.↗

Q branches of the nu7 fundamental of ethane (C2H6) Integrated intensity measurements for atmospheric measurement applications

Laboratory spectra covering the nu7 band of ethane (C2H6) have been recorded, and measurements of integrated intensities of selected Q branches from these spectra are reported. The method by which the spectra were obtained is described, and a typical spectrum covering the PQ3 branch at 2976.8/cm is shown along with a plot of equivalent width vs. optical density for this branch. The values of the integrated intensities reported for each branch are the means of five different optical densities.

Rinsland, C. P.↗

Design of ground test suspension systems for verification of flexible space structures

A simple model demonstrates the frequency-increasing effects of a simple cable suspension on flexible test article/suspension systems. Two passive suspension designs, namely a negative spring mechanism and a rolling cart mechanism, are presented to alleviate the undesirable frequency-increasing effects. Analysis methods are provided for systems in which the augmentations are applied to both discrete and continuous representations of test articles. The damping analyses are based on friction equivalent viscous damping. Numerical examples are given for comparing the two augmentations with respect to minimizing frequency and damping increases.

Cooley, V. M.↗

Lumbar spine disc height and curvature responses to an axial load generated by a compression device compatible with magnetic resonance imaging

STUDY DESIGN: Axial load-dependent changes in the lumbar spine of supine healthy volunteers were examined using a compression device compatible with magnetic resonance imaging. OBJECTIVE: To test two hypotheses: Axial loading of 50% body weight from shoulder to feet in supine posture 1) simulates the upright lumbar spine alignment and 2) decreases disc height significantly. SUMMARY OF BACKGROUND DATA: Axial compression on the lumbar spine has significantly narrowed the lumbar dural sac in patients with sciatica, neurogenic claudication or both. METHODS: Using a device compatible with magnetic resonance imaging, the lumbar spine of eight young volunteers, ages 22 to 36 years, was axially compressed with a force equivalent to 50% of body weight, approximating the normal load on the lumbar spine in upright posture. Sagittal lumbar magnetic resonance imaging was performed to measure intervertebral angle and disc height before and during compression. RESULTS: Each intervertebral angle before and during compression was as follows: T12-L1 (-0.8 degrees +/- 2.5 degrees and -1.5 degrees +/- 2.6 degrees ), L1-L2 (0.7 degrees +/- 1.4 degrees and 3.3 degrees +/- 2.9 degrees ), L2-L3 (4.7 degrees +/- 3.5 degrees and 7.3 degrees +/- 6 degrees ), L3-L4 (7.9 degrees +/- 2.4 degrees and 11.1 degrees +/- 4.6 degrees ), L4-L5 (14.3 degrees +/- 3.3 degrees and 14.9 degrees +/- 1.7 degrees ), L5-S1 (25.8 degrees +/- 5.2 degrees and 20.8 degrees +/- 6 degrees ), and L1-S1 (53.4 degrees +/- 11.9 degrees and 57.3 degrees +/- 16.7 degrees ). Negative values reflect kyphosis, and positive values reflect lordosis. A significant difference between values before and during compression was obtained at L3-L4 and L5-S1. There was a significant decrease in disc height only at L4-L5 during compression. CONCLUSIONS: The axial force of 50% body weight in supine posture simulates the upright lumbar spine morphologically. No change in intervertebral angle occurred at L4-L5. However, disc height at L4-L5 decreased significantly during compression.

Non-NASA Center↗