Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RESPONSE FUNCTIONS”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Predicting System Response at Unmeasured Locations

Traditional techniques to derive dynamic specification for components have a great deal of uncertainty. One of the major sources of uncertainty is that the number of response measurements in the operational system environment is insufficient to determine the component motion. This inadequacy is due to logistical limitations for data recording in field testing and space limitations for accelerometers, strain gages and associated wiring. Available measurements are often some distance from the component and therefore do not represent component motion. Typical straight-line envelopes of these unrepresentative measurements guarantee an increase in the uncertainty. In this paper multiple methods are attempted to expand a sparse set of field test measurements on a system to responses of interest that cannot be measured in the field due to the limitations. Proof of concept is demonstrated on the Modal Analysis Test Vehicle (MATV). The responses of interest, known as “truth responses”, are measured in a system vibration environment along with an optimized sparse set of 30 field responses. Methods to expand the field responses to the truth responses are demonstrated by comparing the acceleration spectral density of the expanded response to the measured response. Two methods utilize a validated finite element model of the MATV. One is developed from purely experiment based frequency response functions of a laboratory pre-test. Finally, these approaches are designed to drastically reduce the uncertainty of the component in-service motion as a basis for developing specifications that are guaranteed to be conservative with a known (instead of unknown) conservatism.

42 ENGINEERING↗
DOE OSTI2024

A Simplified Method for Predicting Shaker Voltage in IMMATs

Impedance Matched Multi-Axis Tests (IMMATs) can replicate in-service vibration induced stress more accurately than single axis shaker table tests as they can better match a part’s operational boundary conditions and excite it in multiple degrees of freedom simultaneously. Here, the shakers used in IMMATs are less powerful than shaker tables, so shaker force limits can be exceeded during tests if they are not placed adequately for the desired environment. The ability to predict shaker voltage and force before performing a test is, therefore, helpful in selecting shaker locations so that their limits are not exceeded. In this study, electrodynamic shakers were modeled as discrete electromechanical systems, and the shaker parameters were chosen to match experimentally obtained acceleration/voltage frequency response functions (FRFs). These models were coupled to a finite element model of the device under test (DUT) via dynamic substructuring, and the substructured model was demonstrated to accurately predict shaker voltage as well as the error in reproducing the environment at multiple accelerometer locations. A simple method called the FRF Multiplication method, in which the FRF of the substructured system is approximated as the product of two separate FRFs of the shaker and DUT respectively, was proposed and applied to the same system, yielding similar voltage and error predictions to those obtained using substructuring. Simple case studies were presented to explore the applicability of the proposed method, and it was demonstrated to have similar accuracy to the substructuring method in a range of cases. Additionally, we showed that while it was not possible to derive a unique model of the shakers from acceleration/voltage FRFs alone, the models that could be obtained were sufficient to predict test error almost perfectly and shaker voltage with less than 40 percent error.

42 ENGINEERING↗
DOE OSTI2026

Investigation of irradiation damage and heat deposition: a comparative analysis for HEU-to-LEU conversion in HFIR

The planned conversion of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel requires detailed evaluation of experiment-relevant parameters to ensure continued performance for materials testing and isotope production. Here, this study presents the first comprehensive assessment of displacements per atom (dpa) and heat deposition rates in target materials within the HFIR flux trap with both HEU and candidate LEU core configurations. Seven analyses were conducted to evaluate key performance metrics, including fast neutron flux distribution, cross section response functions, cross section data, and local dpa and heat deposition rates using mesh- and cell-based tallies. Simulations employed Shift, Monte Carlo N-Particle (MCNP), and the HIFR Controller (HFIRCON) tool suite for high-fidelity transport and depletion modeling. The LEU designs—using U 3 Si 2 -Al dispersion fuel and operating at 95 MW—were compared to the current 85 MW HEU configuration. Results show that while the candidate LEU cores exhibit higher dpa rates due to a harder spectrum and extended cycle lengths, they also demonstrate reduced heat deposition rates in irradiation experiments, primarily due to increased gamma self-shielding from higher 238 U content in the core. These findings confirm that LEU conversion can maintain HFIR’s materials irradiation capabilities but may require redesigning existing experimental hardware.

HEU↗
DOE OSTI2022

Receptance coupling substructure analysis and chatter frequency-informed machine learning for milling stability

This paper describes a milling stability identification approach that simultaneously considers: physics-based models for the tool tip frequency response functions and stability predictions; the binary result from a milling test (automatically labeled as stable or unstable based on frequency content); chatter frequency when an unstable result is obtained; and user risk tolerance. The algorithm applies probabilistic Bayesian machine learning with adaptive, parallelized Markov Chain Monte Carlo sampling to update the probability of stability with each milling test. Furthermore, the result is a robust solution for rapid convergence to optimized milling parameters for maximum metal removal rate using all available information.

42 ENGINEERING↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ExaCA v2.0: A versatile, scalable, and performance portable cellular automata application for additive manufacturing solidification

The previously established ExaCA software for performance portable alloy grain structure simulation has been updated to better represent the solidification behavior during complex alloy processing conditions, such as those encountered during metal additive manufacturing (AM), and for improved performance and scalability. Here, an extension to the time–temperature history input data format and the core ExaCA algorithm to include an arbitrary number of melting and solidification events yielded improved prediction of texture for various melt pool geometries, expanding the range of AM-relevant conditions that can be accurately simulated. Improved heat transport process simulation coupling, including the creation of large raster datasets from single track time–temperature history data and in-memory coupling with the new, performance portable finite difference code Finch, were also demonstrated in example studies on the effect of multilayer AM microstructure predictions on hatch spacing and cell size, respectively. Additional new features are detailed and demonstrated, including the ability to perform simulations using various interfacial response function forms, execute simulations on state-of-the-art hardware, improved usability through post-processing versatility, and improved strong and weak scaling performance. The performance, physics, and versatility improvements demonstrated here will further enable large-scale studies on AM process–microstructure relationships that were not previously possible. Furthermore, the usability improvements and ability to run coupled AM process–microstructure simulations using the Finch-ExaCA workflow will facilitate broader use of this open-source software by the computational materials community.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Serpent neutronics model of Wendelstein 7-X for 14.1 MeV neutrons

In this work, a Serpent 2 neutronics model of the Wendelstein 7-X (W7-X) stellarator is prepared, and an response function for the Scintillating-Fibre neutron detector (SciFi) is calculated using the model. The neutronics model includes the simplified geometry for the key components of the stellarator itself as well as the torus hall. The objective of the model is to assess the 14.1 MeV neutron flux from deuteron-triton fusions in W7-X, where the neutrons are modelled only until they have slowed down to 1 MeV energy. The key messages of this article are: demonstration of unstructured mesh geometry usage for stellarators, W7-X in particular; technical documentation of the model and first insights in fast neutron behaviour in W7-X, especially related to the SciFi: the model indicates that the superconducting coils are the strongest scatterers and block neutrons from large parts of the plasma. The back-scattering from e.g. massive steel support structures is found to be small. The SciFi will detect neutrons from an extended plasma volume in contrast to having an effective line-of-sight.

14 MeV↗
DOE OSTI2024

Application of a Geochemically Informed Leak Detection (GILD) Model to CO 2 Injection Sites on the United States Gulf Coast

The Gulf Coast region possesses great potential for CO 2 enhanced oil recovery (EOR) and CO 2 storage. A geochemically informed leak detection (GILD) model has been applied to CO 2 injection sites on the Gulf Coast with considerations of measurement variability. The Jasper aquifer in Montgomery County, Texas, was chosen to demonstrate the method. Based on background data from wells in the area, combinations of mineral and fluid compositions were used to create 23 scenarios for the geochemical model. The output from the geochemical model was used to identify sensitive monitoring species, and response functions were generated for these as a function of the CO 2 leakage concentration. The sources of measurement variability for background conditions were characterized from the Jasper aquifer background data, and then normalized using the coefficient of variation of each species across the monitoring wells. Bayesian belief network (BBN) models were constructed, and measurement variability of different levels were added to compare leak detection probabilities. Increasing measurement variability decreased the power to detect a leak of a given size. For a moderately high CO 2 concentration of 0.2 mol/kg, the probability of detecting this leakage effect using pH as the monitoring variable in an aquifer with calcite decreases from 98% (no measurement variability) to 61% (medium variability) to 33% (high variability). The loss in power of the sampling protocol with increasing measurement variability is similar in magnitude when Ca 2+ or HCO 3 - is used as the monitoring parameter, but only for aquifers with calcite.

58 GEOSCIENCES↗
DOE OSTI2021

Borderline first-order magnetic phase transition in AlFe 2 B 2

The thermal evolution of lattice parameters coupled with heat capacity data provide insight into tailorable magnetism-structure attributes in the orthorhombic compound AlFe 2 B 2 that was synthesized with and without small additions of gallium. Temperature-dependent X-ray powder diffraction experiments conducted through the magnetic phase transition reveal that the a- and b-parameters of both samples increase with increasing temperature while the c-parameter decreases. While a weak volumetric thermal expansion is noted over a range of temperatures well below and above the magnetic phase transition, anomalous behavior was observed within the phase transition region itself to reveal a magnetostructural phase transition with borderline first-order character in the Ga-modified sample but of more second-order character in the Ga-free sample. It is established that the nearest-neighbor Fe-Fe interatomic distance within the (ab)-plane plays a dominant role in influencing the magneto-functional response of these compounds. The magnetocaloric properties are discussed in the context of temperature-induced changes of the interatomic bonding that are influenced by the hypothesized presence of iron antisite defects in the AlFe 2 B 2 lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modeling multiple scattering transient of an ultrashort laser pulse by spherical particles

The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gréhan (2000) is used to derive an energetic model of the scattering transient. It is shown that the spreading and extinction of the pulse can be decoupled from the transient of scattering, which allows to describe each phenomenon individually. The transient of scattering is modeled with the Lorenz-Mie Theory and thus is also valid for a relative refractive index lower than one, contrary to the Debye series expansion which does not converge close to the critical angle. This is made possible after the introduction of a new physical object, the Scattering Impulse Response Function (SIRF) which allows to detect the different modes of scattering transient, in time and direction. Here the present approach is more generic, as it enables to simulate clouds of air bubbles in water, which was not possible previously. Two different approaches are proposed within the Monte Carlo framework. The first is a pure Monte Carlo approach where the delay due to the scattering is randomly drawn at each event, while the second is based on the transport of the whole scattering signal. They are both embedded in the Monte Carlo code Scatter3D. Both models produce equivalent trends and are validated against published numerical results. They are then applied to the multiple scattering of ultra short pulse by a cloud of bubble in water in the forward direction. The pulse spread due to the propagation in water is computed for a wide range of traveled distances and pulse durations, and the optimal pulse duration is given to minimize the pulse spread at a given distance. The main result is that the scattered photons exit the turbid medium earlier than the ballistic photons and produce a double peak related to the refraction in the bubble. This demonstrates the possibility to develop new diagnostics to characterize dynamic bubbly flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Modeling the processes of soil moisture in regulating microbial and carbon-nitrogen cycling

Soil carbon (C) and nitrogen (N) cycles and their complex responses to hydro-climatic forcing have gained increasing attention. While the temperature effects have been intensively studied, soil moisture response functions (SMRFs) are not well documented for various microbial and enzymatic processes due to the difficulties in directly measuring and differentiating the moisture effects on various processes. In this work we extended our C-only Microbial-ENzyme Decomposition (MEND) model to the C-N coupled MEND model with flexible element stoichiometry. Our model calibration showed good agreement between simulated and observed C:N ratios in soil organic matter and microbial biomass, as well as the ammonium and nitrate concentrations. We show that the selection of SMRFs for specific biogeochemical processes could result in significant differences in model simulated microbial and C-N processes. In particular, it is essential to account for the soil moisture effects on microbial dormancy and resuscitation, as the changes in microbial physiology under favorable or stressful conditions will exert strong controls on soil C and N dynamics. We also advocate the utilization of dynamic (time-variant) data (e.g., litter input, N deposition, soil temperature and moisture), instead of time-invariant data, to drive model simulations and analyses. Dynamic forcing data (particularly dynamic soil moisture) better represent the real-world climate and environmental conditions, which could facilitate more realistic modeling and understanding of soil C and nutrient cycling in a changing world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Physics-informed Bayesian machine learning case study: Integral blade rotors

This paper provides a physics-informed Bayesian machine learning (PIBML) description and case study. The PIBML approach applies three physics-based models to establish the initial beliefs before testing to determine the probability of milling stability (or prior). These include: receptance coupling substructure analysis (RCSA) prediction for the tool tip frequency response functions; finite element software prediction of the mechanistic force model coefficients; and a spindle speed-dependent power law model for process damping. Testing was then performed to identify optimal stable machining conditions using an expected improvement in material removal rate criterion. The prior probability of stability was updated using the test results to determine the posterior probability of stability. The test results were compared to the parameter recommendations provided by the endmill manufacturer. A demonstration integral blade rotor was machined at the optimal stable machining conditions for 304 stainless steel and 6061-T6 aluminum. Finally, the disagreement between manufacturer recommendations and milling performance in both materials tested emphasizes the need for broad implementation of PIBML approaches to increase machining productivity and efficiency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Local modeling for FRF estimation with noisy input measurements

The frequency response function (FRF) is an essential means by which dynamic systems are qualified. In recent years, local modeling approaches have been extensively researched and shown to significantly outperform traditional FRF estimators. However, the standard local modeling approach assumes a perfectly-known system input, which results in biased FRF estimates in the presence of input noise. This paper derives a simple adjustment that can be used to improve FRF estimation for systems subjected to random excitation with noisy input data. This improvement can be implemented with little modification to standard local modeling algorithms and with little additional computational burden. The adjustment is coupled with a model selection procedure to avoid underfitting and overfitting. In conclusion, the methods presented in this paper are validated on a simulation, and they are shown to reduce bias due to input noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Increased damping through captured powder in additive manufacturing

This paper describes the retention of metal powder within additively manufactured structures to achieve increased damping. In this work, solid and hollow (powder filled) aluminum cylinders were fabricated in using laser powder bed fusion where the internal, captured powder in the hollow cylinders served as an energy dissipation mechanism. Impact testing was completed to compare the frequency response functions of three cylinder geometries: one solid cross-section and two hollow cross-sections with different inner diameters and, therefore, different powder volumes. Increased damping with larger inner diameter was observed. The damping was quantified using a structural damping model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic force and stability prediction for milling using feed rate scheduling software and time-domain simulation

This paper describes: 1) the use of feed rate scheduling software to predict the radial depth of cut variation for three-axis milling toolpaths and; 2) the use of the radial depth profile in a time-domain simulation to predict dynamic cutting forces. The time-domain simulation, which also includes the tool tip frequency response functions and force model (which relates the cutting force components to the chip geometry) as inputs, enables dynamic force profiles to be predicted and parameter combinations that cause chatter to be identified. A ramp geometry is selected that provides constantly varying radial depth and force predictions are completed at multiple axial depths for comparison to measured forces. Both stable and unstable (chatter) milling conditions were observed with good agreement between time-domain simulation and measurement results. The value of combining the feed rate scheduling software and time-domain simulation is demonstrated.

42 ENGINEERING↗
DOE OSTI2022

A machining digital twin for hybrid manufacturing

Hybrid manufacturing consisting of metal additively manufactured preforms and computer numerical control (CNC) machining has been established to be an effective method for high material use rates. However, hybrid manufacturing introduces unique challenges. Near-net shape designs are typically selected, which result in a smaller margin for part placement within the stock and stringent requirements for work coordinate system identification. Additionally, less stock material reduces the preform stiffness, which limits the material removal rates during machining. This paper demonstrates a digital twin for CNC machining of a wire arc additively manufactured preform that implements: 1) structured light scanning for stock model identification and tool path generation; 2) a fused filament fabrication apparatus to attach temporary fiducials and scan targets to the preform that enable coordinate system definition for both the CAM and CNC machine; 3) preform and tool tip frequency response function measurements to enable stable milling parameter selection; and 4) post-manufacturing measurements of geometry, surface finish, and structural dynamics to confirm designer intent. These efforts define key components of the machining digital twin for hybrid manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Comparison of gamma-ray spectral analysis methods for thick-walled UF 6 cylinders

Three data analysis methods were compared for gamma ray spectra collected by planar and coaxial high-purity germanium detectors (HPGe) on thick-walled uranium hexafluoride cylinders. The cylinder contents range from natural to low-enriched uranium (4.95% 235 U). Each spectrum was analyzed with the traditional 185.7 keV enrichment meter method, a 185.7/1001 keV ratio method, and the Fixed Energy Response Function Analysis with Multiple Efficiency (FRAM) software that utilizes low- and high-energy gamma rays and X-rays. When results for all cylinders were compared, the enrichment meter method and the 185.7/1001 keV peak ratio method consistently provided the most accurate data for all detector types. The error on the enrichment calculations analyzed with FRAM did not meet the international target value whereas those calculated with the enrichment meter and ratio method did. Enrichments calculated by FRAM were found to skew to higher values on natural uranium and to lower values on low-enriched uranium. FRAM calculated the least-precise enrichments; replicate analyses of single cylinders resulted in 235 U values that varied by as much as 19%. FRAM does have the significant advantages in that it does not require the use of calibration cylinders, infinite thickness criteria are not necessary, and it can be applied to any HPGe detector. Finally, if calibration cylinders are available, the enrichment meter and 185.7/1001 peak ratio are the recommended methods as they are the most accurate and precise.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Characterization of Multianode Photomultiplier Tubes for use in the CLAS12 RICH detector

In this work, we present results of the detailed study of several hundred Hamamatsu H12700 Multianode Photomultiplier Tubes (MaPMTs), characterizing their response to the Cherenkov light photons in the second Ring Imaging Cherenkov detector, a part of the CLAS12 upgrade at Jefferson Lab. The total number of pixels studied was 25536. The single photoelectron spectra were measured for each pixel at different high voltages and light intensities of the laser test setup. Using the same dedicated front-end electronics as in the first RICH detector, the setup allowed us to characterize each pixel’s properties such as gain, quantum efficiency, signal crosstalk between neighboring pixels, and determine the signal threshold values to optimize their efficiency to detect Cherenkov photons. A recently published state-of-the-art mathematical model, describing photon detector response functions measured in low light conditions, was extended to include the description of the crosstalk contributions to the spectra. The database of extracted parameters will be used for the final selection of the MaPMTs, their arrangement in the new RICH detector, and the optimization of the operational settings of the front-end electronics. The results show that the characteristics of the H12700 MaPMTs satisfy our requirements for the position-sensitive single photoelectron detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗