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

Identification of Low Order Equivalent System Models From Flight Test Data

Identification of low order equivalent system dynamic models from flight test data was studied. Inputs were pilot control deflections, and outputs were aircraft responses, so the models characterized the total aircraft response including bare airframe and flight control system. Theoretical investigations were conducted and related to results found in the literature. Low order equivalent system modeling techniques using output error and equation error parameter estimation in the frequency domain were developed and validated on simulation data. It was found that some common difficulties encountered in identifying closed loop low order equivalent system models from flight test data could be overcome using the developed techniques. Implications for data requirements and experiment design were discussed. The developed methods were demonstrated using realistic simulation cases, then applied to closed loop flight test data from the NASA F-18 High Alpha Research Vehicle.

Morelli, Eugene A.↗

Aladdin: Exploration and Sample Return from the Moons of Mars

Aladdin is a remote sensing and sample return mission focused on the two small moons of Mars, Phobos and Deimos. Understanding the moons of Mars will help us to understand the early history of Mars itself. Aladdin's primary objective is to acquire well documented, representative samples from both moons and return them to Earth for detailed analyses. Samples arrive at Earth within three years of launch. Aladdin addresses several of NASA's highest priority science objectives: the origin and evolution of the Martian system (one of two silicate planets with satellites) and the composition and nature of small bodies (the building blocks of the solar system). The Aladdin mission has been selected as a finalist in both the 1997 and 1999 Discovery competitions based on the high quality of science it would accomplish. The equivalent of Aladdin's Phase A development has been successfully completed, yielding a high degree of technical maturity. Aladdin uses an innovative flyby sample acquisition method, which has been validated experimentally and does not require soft landing or anchoring. An initial phasing orbit at Mars reduces mission propulsion requirements, enabling Aladdin to use proven, low-risk chemical propulsion with good mass margin. This phasing orbit is followed by a five month elliptical mission during which there are redundant opportunities for acquisition of samples and characterization of their geologic context using remote sensing. The Aladdin mission is a partnership between Brown University, the Johns Hopkins University Applied Physics Laboratory, Lockheed Martin Astronautics, and NASA Johnson Space Center.

Pieters, C.↗

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mixed-Fidelity Approach for Design of Low-Boom Supersonic Aircraft

This paper documents a mixed-fidelity approach for the design of low-boom supersonic aircraft with a focus on fuselage shaping.A low-boom configuration that is based on low-fidelity analysis is used as the baseline. The fuselage shape is modified iteratively to obtain a configuration with an equivalent-area distribution derived from computational fluid dynamics analysis that attempts to match a predetermined low-boom target area distribution and also yields a low-boom ground signature. The ground signature of the final configuration is calculated by using a state-of-the-art computational-fluid-dynamics-based boom analysis method that generates accurate midfield pressure distributions for propagation to the ground with ray tracing. The ground signature that is propagated from a midfield pressure distribution has a shaped ramp front, which is similar to the ground signature that is propagated from the computational fluid dynamics equivalent-area distribution. This result supports the validity of low-boom supersonic configuration design by matching a low-boom equivalent-area target, which is easier to accomplish than matching a low-boom midfield pressure target.

Li, Wu↗

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

Wing Shape Sensing from Measured Strain

A new two-step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. In the first step, a measured strain is fitted using a piecewise least-squares curve fitting method together with the cubic spline technique. These fitted strains are integrated twice to obtain deflection data along the fibers. In the second step, computed deflection along the fibers are combined with a finite element model of the structure in order to interpolate and extrapolate the deflection and slope of the entire structure through the use of the System Equivalent Reduction and Expansion Process. The theory is first validated on a computational model, a cantilevered rectangular plate wing. The theory is then applied to test data from a cantilevered swept-plate wing model. Computed results are compared with finite element results, results using another strainbased method, and photogrammetry data. For the computational model under an aeroelastic load, maximum deflection errors in the fore and aft, lateral, and vertical directions are -3.2%, 0.28%, and 0.09%, respectively; and maximum slope errors in roll and pitch directions are 0.28% and -3.2%, respectively. For the experimental model, deflection results at the tip are shown to be accurate to within 3.8% of the photogrammetry data and are accurate to within 2.2% in most cases. In general, excellent matching between target and computed values are accomplished in this study. Future refinement of this theory will allow it to monitor the deflection and health of an entire aircraft in real time, allowing for aerodynamic load computation, active flexible motion control, and active induced drag reduction.

shape sensing↗

On the accuracy of compressibility transformations

This study highlights the importance of satisfying the eddy viscosity equivalence below the logarithmic layer, to deriving accurate compressibility transformations. First, we analyze the ability of known transformations to satisfy the eddy viscosity equivalence and show that the accuracy of these transformations is strongly dependent on this ability. Second, in a step-by-step manner, we devise new transformations that satisfy this hypothesis. An approach based on curve fitting of the incompressible Direct Numerical Simulation data for eddy viscosity profiles below the logarithmic layer provides an extremely accurate transformation, which motivates self-contained methods, making use of mixing length formulas in the inner region. It is shown that the accuracy of existing transformations can be significantly improved by applying these ideas, below the logarithmic layer. Motivated by the effectiveness of the formulations derived from eddy viscosity equivalence, we introduce a new integral transformation based on Reynolds number equivalence between compressible and incompressible flows. This approach is based on defining a new compressible velocity scale, which affects the accuracy of transformations. Several choices for the velocity scale are tested, and in each attempt, it is shown that the eddy viscosity equivalence plays a very important role for the accuracy of compressibility transformations.

42 ENGINEERING↗

Safe Tracking Control of an Uncertain Euler-Lagrange System with Full-State Constraints using Barrier Functions

This paper presents a novel, safe tracking control design method that learns the parameters of an uncertain Euler-Lagrange (EL) system online using adaptive learning laws. A barrier function (BF) is first used to transform the full-state constrained EL-dynamics into an equivalent unconstrained dynamics. An adaptive tracking controller is then developed along with the parameter update law in the transformed state space such that the states remain bounded for all time within a prescribed bound. A stability analysis is developed that considers the EL-dynamics’ uncertainty, yielding a semi-globally uniformly ultimately bounded (SGUUB) tracking error and the parameter estimation error. The controller design is validated in simulations using a two-link planar manipulator. The results show the proposed method’s ability to track the reference trajectory while remaining inside each of the predefined state bounds.

Robots↗

Radiation Measured with Different Dosimeters for ISS-Expedition 18-19/ULF2 on Board International Space Station during Solar Minimum

Radiation field of particles in low Earth orbit (LEO) is mainly composed of galactic cosmic rays (GCR), solar energetic particles and particles in SAA (South Atlantic Anomaly). GCR are modulated by solar activity, at the period of solar minimum activity, GCR intensity is at maximum and the main contributor for space radiation is GCR. At present for space radiation measurements conducted by JSC (Johnson Space Center) SRAG (Space Radiation Analysis Group), the preferred active dosimeter sensitive to all LET (Linear Energy Transfer) is the tissue equivalent proportional counter (TEPC); the preferred passive dosimeters are thermoluminescence dosimeters (TLDs) and optically stimulated luminescence dosimeters (OSLDs) sensitive to low LET as well as CR-39 plastic nuclear track detectors (PNTDs) sensitive to high LET. For the method using passive dosimeters, radiation quantities for all LET can be obtained by combining radiation results measured with TLDs/OSLDs and CR-39 PNTDs. TEPC, TLDs/OSLDs and CR-39 detectors were used to measure the radiation field for the ISS (International Space Station) - Expedition 18-19/ULF2 space mission which was conducted from 15 November 2008 to 31 July 2009 - near the period of the recent solar minimum activity. LET spectra (differential and integral fluence, absorbed dose and dose equivalent) and radiation quantities were measured for positions TEPC, TESS (Temporary Sleeping Station, inside the polyethylene lined sleep station), SM-P 327 and 442 (Service Module - Panel 327 and 442). This paper presents radiation LET spectra measured with TEPC and CR-39 PNTDs and radiation dose measured with TLDs/OSLDs as well as the radiation quantities combined from results measured with passive dosimeters.

Zhou, Dazhuang↗

Implementation of a Smeared Crack Band Model in a Micromechanics Framework

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

Progressive Failure of a Unidirectional Fiber-Reinforced Composite Using the Method of Cells: Discretization Objective Computational Results

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

Enhanced detection limits for radiokrypton analysis

In this paper, we present a method for improving detection limits of Atom Trap Trace Analysis for the krypton radioisotopes 85 Kr and 81 Kr. For the case of 85 Kr this work demonstrates that systematic use of isotopically depleted gas for calibration and extended conditioning of the instrument results in a detection limit of 900 85 Kr atoms per 11 μl of Kr gas, equivalent to a 85 Kr/Kr isotopic abundance of 3 x 10 -15 . This improvement of roughly two orders of magnitude over previously reported limits will help to expand the reach of radiokrypton dating towards longer age ranges where most of the radioisotopes have decayed. Additionally, the method offers an opportunity to investigate radiokrypton production via spontaneous fission within naturally occurring minerals to understand potential underground production of these isotopes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Comparison and calibration of dose delivered by 137 Cs and x-ray irradiators in mice

Objective. The Office of Radiological Security, U.S. Department of Energy's National Nuclear Security Administration, is implementing a radiological risk reduction program which seeks to minimize or eliminate the use of high activity radiological sources, including 137 Cs, by replacing them with non-radioisotopic technologies, such as x-ray irradiators. The main goal of this paper is to evaluate the equivalence of the dose delivered by gamma- and x-ray irradiators in mice using experimental measurements and Monte Carlo simulations. We also propose a novel biophantom as an in situ dose calibration method. Approach. We irradiated mouse carcasses and 3D-printed mouse biophantoms in a 137 Cs irradiator (Mark I-68) and an x-ray irradiator (X-Rad320) at three voltages (160 kVp, 225 kVp and 320 kVp) and measured the delivered radiation dose. A Geant4-based Monte Carlo model was developed and validated to provide a comprehensive picture of gamma- and x-ray irradiation in mice. Main Results. Our Monte Carlo model predicts a uniform dose delivered in soft-tissue for all the explored irradiation programs and in agreement with the absolute dose measurements. Our Monte Carlo model shows an energy-dependent difference between dose in bone and in soft tissue that decreases as photon energy increases. Dose rate depends on irradiator and photon energy. We observed a deviation of the measured dose from the target value of up to –9% for the Mark I-68, and up to 35% for the X-Rad320. The dose measured in the 3D-printed phantoms are equivalent to that in the carcasses within 6% uncertainty. Significance. Our results suggest that 320 kVp irradiation is a good candidate to substitute 137 Cs irradiation barring a few caveats. There is a significant difference between measured and targeted doses for x-ray irradiation that suggests a strong need for in situ calibration, which can be achieved with 3D-printed mouse biophantoms. A dose correction is necessary for bone doses, which can be provided by a Monte Carlo calculation. Lastly, the biological implications of the differences in dose rates and dose per photon for the different irradiation methods should be carefully assessed for each small-animal irradiation experiment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Time-resolved electrical potential pump – X-ray photoelectron spectroscopy probe developments for investigating dynamic processes occurring at electrochemical interfaces

Electrode–electrolyte interfaces are of critical importance in several fields, including renewable energy, corrosion, and environmental chemistry. However, investigating these interfaces under operational conditions poses considerable challenges due to the limitations of the instrumentation employed. While recent advancements in in situ and operando techniques have enhanced our comprehension of the steady-state properties of solid-liquid interfaces, the dynamic behaviors of these systems remain inadequately explored. This study introduces a time-resolved X-ray photoelectron spectroscopy (XPS) technique designed to capture transient reaction intermediates and charging dynamics at electrified interfaces. The presented proof-of-principle study demonstrates that electrochemical processes, represented by an equivalent electrical circuit (EEC) model, can be probed and understood using square wave voltage pulses of a potentiostat synchronized to the modified data acquisition of an XPS setup. This method offers a valuable alternative to traditional pump–probe techniques, facilitating the investigation of a broader range of electrochemical systems. A dedicated software package for analyzing time- and energy-resolved XPS with a focus on extracting parameters of the EEC is geared towards benchmarking different EECs in future real-world electrochemical experiments.

Electrochemistry↗

Cradle-to-gate life cycle assessment of advanced composite panels incorporating CO 2 -derived multi-walled carbon nanotubes and hemp fiber for sustainable building applications

Advanced composite panels represent a promising pathway to reducing carbon emissions in the construction industry, yet comprehensive environmental impact assessments remain limited. Here, in this study, we conduct a life cycle assessment (LCA) to evaluate the environmental impacts of innovative composite panels produced from multi-walled carbon nanotubes (MWCNTs), hemp fiber (HF), recycled carbon fiber (rCF), and recycled polypropylene (PP), exploring their potential as baseline structural equivalents to conventional gypsum board. MWCNTs and HF play a critical role in sequestering carbon during raw material production, while the recycling processes for CF and PP generally require less energy compared to virgin material production. The LCA evaluates environmental performance using the TRACI 2.1 method, covering global warming potential (GWP), ozone depletion, smog formation, acidification, eutrophication, carcinogenic and non-carcinogenic effects, respiratory impacts, ecotoxicity, and fossil fuel depletion. Compositional variations—resin type (virgin vs. recycled), rCF content (9–29 wt%), and HF content (10–30 wt%)—are introduced for sensitivity and hotspot analyses. Results demonstrate that, when compared on the basis of preliminary structural equivalence, increasing recycled PP, rCF, and HF content can significantly reduce global warming potential compared to gypsum board. Beyond carbon reduction, the composite panels show trade-offs across other environmental categories. With the growing demand for composite materials in interior panels, ceiling systems, and exterior claddings, these findings highlight the environmental benefits and potential trade-offs of the proposed composites, establishing a foundational framework to support their continued development toward full building-system integration.

Advanced composite manufacturing↗

Wing Shape Sensing from Measured Strain

A new two-step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. In the first step, a measured strain is fitted using a piecewise least-squares curve fitting method together with the cubic spline technique. These fitted strains are integrated twice to obtain deflection data along the fibers. In the second step, computed deflection along the fibers are combined with a finite element model of the structure in order to interpolate and extrapolate the deflection and slope of the entire structure through the use of the System Equivalent Reduction and Expansion Process. The theory is first validated on a computational model, a cantilevered rectangular plate wing. The theory is then applied to test data from a cantilevered swept-plate wing model. Computed results are compared with finite element results, results using another strain-based method, and photogrammetry data. For the computational model under an aeroelastic load, maximum deflection errors in the fore and aft, lateral, and vertical directions are -3.2 percent, 0.28 percent, and 0.09 percent, respectively; and maximum slope errors in roll and pitch directions are 0.28 percent and -3.2 percent, respectively. For the experimental model, deflection results at the tip are shown to be accurate to within 3.8 percent of the photogrammetry data and are accurate to within 2.2 percent in most cases. In general, excellent matching between target and computed values are accomplished in this study. Future refinement of this theory will allow it to monitor the deflection and health of an entire aircraft in real time, allowing for aerodynamic load computation, active flexible motion control, and active induced drag reduction..

system equivalent reduction and expansion process↗

An implementation of the QMR method based on coupled two-term recurrences

The authors have proposed a new Krylov subspace iteration, the quasi-minimal residual algorithm (QMR), for solving non-Hermitian linear systems. In the original implementation of the QMR method, the Lanczos process with look-ahead is used to generate basis vectors for the underlying Krylov subspaces. In the Lanczos algorithm, these basis vectors are computed by means of three-term recurrences. It has been observed that, in finite precision arithmetic, vector iterations based on three-term recursions are usually less robust than mathematically equivalent coupled two-term vector recurrences. This paper presents a look-ahead algorithm that constructs the Lanczos basis vectors by means of coupled two-term recursions. Implementation details are given, and the look-ahead strategy is described. A new implementation of the QMR method, based on this coupled two-term algorithm, is described. A simplified version of the QMR algorithm without look-ahead is also presented, and the special case of QMR for complex symmetric linear systems is considered. Results of numerical experiments comparing the original and the new implementations of the QMR method are reported.

Freund, Roland W.↗