Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Impulse Approximation”

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.

54 records · Page 3

Impulse response measurements of fast scintillator-based current mode detectors

In this study, experimental measurements were performed to characterize the impulse response of seven different gamma ray detector prototypes being considered to be fielded as a current mode diagnostic for neutron-diagnosed subcritical experiments at the Nevada National Security Site (NNSS). For such experiments, the impulse response of the fielded detector is required to have a full-width half-maximum (FWHM) of approximately 5 ns or less and a low amplitude tail. Each of the detector prototypes evaluated in this work were selected based on the potential to meet this performance requirement. An impulse of bremsstrahlung x-rays created from 2 MeV electrons from the NNSS’s Transformational Diagnostics and Imaging, Los Alamos, linear accelerator was measured with each detector. The measured impulse response was evaluated for each detector configuration. A comparison of these impulse responses revealed a relationship between fast timing performance and scintillator material, geometry, surface reflectivity conditions, and photomultiplier tube (PMT) selection. The detector configuration yielding an impulse response with the shortest FWHM of 3.82 ± 0.11 ns and possessing a low amplitude tail was a 5-inch (12.7-cm) diameter, 5-inch (12.7-cm) height cylindrical EJ-399-17-VI liquid scintillator in a black-painted housing coupled to an Adit D798B 5-inch (12.7-cm) PMT by means of a 0.079-inch (0.2-cm) EJ-560 optical coupling pad.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Entrainment, Detrainment, and Dilution of Dry and Moist Atmospheric Thermals

Here this study examines the entrainment, detrainment, and dilution of dry and moist (cloud) atmospheric thermals in large-eddy simulations. In a neutrally stable environment (with respect to dry dynamics), moist thermals have an increase in radius R with thermal height z t (α ≡ dR/dz t ) about 4 times smaller compared to dry thermals when density stratification is considered and ~2.4 times smaller without density stratification (i.e., applying the Boussinesq approximation). An analytic expression relating α to several dimensionless parameters is derived from the thermal impulse–circulation relation to clarify the factors impacting α. This expression shows that the difference in buoyancy structure between moist and dry thermals, with buoyancy concentrated in the central cores of moist thermals owing to latent heating, explains their smaller spreading rates. Individual contributions of entrainment and detrainment are analyzed using a direct parcel-based approach in the simulations. Moist thermals have similar fractional detrainment but much smaller fractional entrainment rates compared to dry thermals, consistent with the differences in α. Despite having smaller α, moist thermals are similarly dilute (quantified by a passive tracer) as dry thermals because of their greater mixing efficiency with the environment. Thus, moist thermals are substantially dilute but expand much less in size/volume as they rise compared to dry thermals. The α values for moist thermals in (dry) neutral and statically stable environments are similar, but fractional entrainment and especially detrainment rates are greater in the stable environment. Large detrainment rates are associated with a breakdown of the broader thermal vortex ring structure, especially with low environmental relative humidity, attributed in part to evaporation and buoyancy reversal.

54 ENVIRONMENTAL SCIENCES↗

Review of analytical studies on seismic fluid-structure interaction of base-supported cylindrical tanks

Analytical solutions for seismic fluid-structure interaction (FSI) of tanks are used by analysts and engineers for a number of purposes, including preliminary sizing, risk assessment, and verification of numerical models. Analytical solutions for seismic FSI of tanks were first developed in the 1930s and have been extended and/or modified in the following nine decades to accommodate different seismic inputs and boundary conditions. This report reviews published analytical solutions for rigid and flexible, base-supported, cylindrical tanks subjected to unidirectional horizontal motion of a small amplitude. The studies parsed the fluid-structure (FSI) response into an impulsive and a convective component, and addressed frequencies of lateral motions of the tank (i.e., impulsive frequencies), hydrodynamic pressures, frequencies of waves (i.e., convective frequencies), wave heights, and reactions at the base. The solutions were categorized into exact and approximate solutions based on the methodologies used. Herein, FSI responses are calculated and compared for a range of tank dimensions using the analytical solutions from the different studies. The responses presented here are normalized to be unitless, and can be used for tanks with different dimensions and mechanical properties, and subjected to different input motions. A sample steel water tank is analyzed using the normalized solutions for two earthquake ground motions with different frequency contents. Ignoring tank flexibility may significantly underestimate seismic responses of the tank and its contained fluid. The approximate solutions, which have been widely used for decades, may not be sufficiently accurate in the modern era. Calculation errors in the reviewed studies are identified and corrected, and differences are quantified for the sample tank.

42 ENGINEERING↗

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗

Angle-dependent strong-field ionization and fragmentation of carbon dioxide measured using rotational wave packets

In this work, we experimentally study the angle-dependent single ionization of carbon dioxide (CO 2 ) by linearly and circularly polarized pulses. The angle dependence of the ionization probability by linearly polarized pulses extracted from time-domain measurements on an impulsively excited rotational wave packet is compared with data obtained from a direct angle-scan measurement. The results from the measurement with linear and circular polarization are consistent with the adiabatic ionization approximation. Further, we extend the time-domain method to extract the dependence of the asymptotic momentum distribution of fragment ions on the orientation of the molecular axis, and apply it to investigate dissociative double ionization of CO 2 . We show that such measurements can directly test the validity of the axial recoil approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Monitoring the propagation of mechanical discontinuity using data-driven causal discovery and supervised learning

Mechanical wave transmission through a material is influenced by the mechanical discontinuity in the material. The propagation of embedded discontinuities can be monitored by analyzing the wave-transmission measurements recorded by a multipoint sensor system placed on the surface of the material. The proposed workflow monitors the propagation of mechanical discontinuity through three stages, namely initial, intermediate, and final stages, by using supervised learning followed by data-driven causal discovery. To the end, the workflow processes the multipoint waveform measurements resulting from a single impulse source, while considering the effects of wave attenuation, dispersion and multiple wave-propagation modes due to the discontinuity and material boundaries. Among various feature reduction techniques ranging from decomposition methods to manifold approximation methods, the features derived based on statistical parameterizations of the measured waveforms lead to reliable monitoring that is robust to changes in precision, resolution, and signal-to-noise ratio of the multipoint sensor measurements. The numbers of zero-crossing, negative-turning, and positive turning in the waveforms are the strongest causal signatures of the propagation of mechanical discontinuity. Higher order moments of the waveforms, such as variance, skewness and kurtosis, are also strong causal signatures of the propagation. Finally, the newly discovered causal signatures confirm that the statistical correlations and conventional feature rankings are not always statistically significant indicators of causality.

42 ENGINEERING↗

Impulse Noise Sampling at 9965 - 300-pound energetic experimant

On October 1, 2022, sound level measurements were taken at various locations throughout Kirtland Air Force Base (KAFB) and Southeastern Albuquerque. The purpose was to support sound propagation modeling predictions and sound regulations for public exposure during the detonation of an approximately 300-pound energetic experiment. Ground Zero was located on Range 7 of Sandia Thunder Range (06647). A total of 8 measurement locations were identified (e.g., 5 on KAFB and 3 in the Southeastern Albuquerque neighborhoods).

42 ENGINEERING↗

Attosecond response of molecules to impulsive ionization

When matter interacts with energetic radiation it can undergo sudden, or impulsive, ionization. This process can drive chemical change and occurs widely in space and planetary atmospheres, yet its comprehensive description challenges our current theoretical and computational capabilities as it requires advanced treatment of electron correlation and nonadiabatic dynamics beyond the Born–Oppenheimer approximation. Here, in this study, we measure the response of the para-aminophenol molecule to sudden ionization. Using attosecond X-ray absorption spectroscopy, we resolve the ultrafast dynamics of the ionized molecule with atomic precision. A subfemtosecond decay corresponds to states undergoing non-radiative decay, whereas few-femtosecond oscillatory signatures are associated with electronic wavepacket motion in stable cation states that later couple to nuclear motion. We compare our measurement with state-of-the-art computational modelling, qualitatively reproducing the observed response across multiple timescales. These results provide a benchmark for computational models of sudden ionization and ultrafast charge motion in matter.

Driver, Taran [SLAC National Accelerator Laborator↗

Sound propagation in realistic interactive 3D scenes with parameterized sources using deep neural operators

We address the challenge of acoustic simulations in three-dimensional (3D) virtual rooms with parametric source positions, which have applications in virtual/augmented reality, game audio, and spatial computing. The wave equation can fully describe wave phenomena such as diffraction and interference. However, conventional numerical discretization methods are computationally expensive when simulating hundreds of source and receiver positions, making simulations with parametric source positions impractical. To overcome this limitation, we propose using deep operator networks to approximate linear wave-equation operators. This enables the rapid prediction of sound propagation in realistic 3D acoustic scenes with parametric source positions, achieving millisecond-scale computations. By learning a compact surrogate model, we avoid the offline calculation and storage of impulse responses for all relevant source/listener pairs. Our experiments, including various complex scene geometries, show good agreement with reference solutions, with root mean squared errors ranging from 0.02 to 0.10 Pa. Notably, our method signifies a paradigm shift as—to our knowledge—no prior machine learning approach has achieved precise predictions of complete wave fields within realistic domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Combined Optical and Radio-Frequency Measurements of a Lightning Megaflash by the FORTE Satellite

The optical and VHF instrumentation on the Fast On-Orbit Recording of Transient Events (FORTE) satellite is used to document the combined phenomenology evolution of a lightning “megaflash” – mesoscale lightning that propagates laterally over exceptional distances. We identify a FORTE flash whose maximum extent was 82 km and inferred length over multiple distinct branches exceeded 100 km. This flash lasted 1.2 s and produced 250 optical and 591 RF events. We find that the channel development mapped by FORTE’s pixelated lightning imager (LLS) occurred at a typical speed of 2.6x10 5 m s -1 and was accompanied by sustained periods of VHF emission that could individually exceed 100 ms in duration. The impulsive IC events generated by the flash indicate that this development occurred at altitudes between 3 and 8 km. Four +CG strokes were identified in the VHF waveform data that are responsible for two of the three highly-radiant LLS groups (two of the +CGs were not as optically bright as the others). These strokes occurred at different locations throughout the flash footprint with the most distant strokes separated by approximately 50 km. These space-based observations match previous observations of megaflashes as well as ground-based measurements of negative leader development during “spider” lightning, suggesting that FORTE is sensing the same phenomena.

Peterson, Michael Jay↗

Discretization Writeup for Grey Flux-Limited Radiation Diffusion

This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetohydrodynamic Wave Mode Conversion in a Pseudostreamer Structure

The interaction of a fast magnetoacoustic wave with a magnetic null point in a pseudostreamer magnetic configuration is modeled. The incident fast wave is excited by an impulsive velocity driver located in the corona outside the Alfvén-acoustic equipartition layer. The fast wave front experiences refraction and wraps around the null point, while the front segment approaching the null point along the magnetic bisector passes through it. In the equipartition layer surrounding the null point, part of the fast wave energy is converted into slow magnetoacoustic waves. Slow waves propagate outward from the null point along the magnetic separatrices. The fronts of slow waves exhibit a characteristic chevron shape. The fast-to-slow mode conversion efficiency increases with the width of the initial perturbation, following a power-law dependence with an index of approximately 1.06. For initial wave amplitudes below about 10% of the local Alfvén speed, the dependence of the conversion efficiency on the initial amplitude is weak. After interacting with the null point, the amplitudes of the fast waves decrease significantly owing to energy conversion and the spreading of wave energy over a broad spatial domain. In contrast, the energy of slow waves remains concentrated in the vicinity of the magnetic separatrices, resulting in less pronounced geometric attenuation of these waves. These findings are consistent with observational reports of the interaction of a coronal fast wave with a pseudostreamer.

Zhong, Yu [Univ. of Warwick, Coventry (United King↗

Structural modeling and dynamic analysis of condensate storage tanks in nuclear power plants

Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Real-Time Seismic System for Monitoring, Imaging, and Characterization (RT-SEISMIC)

The goal of this Laboratory Directed Research and Development (LDRD) project was to develop a borehole seismic source and sensor array to enable real-time seismic imaging at scales and conditions relevant to the energy industry including both fossil-energy and geothermal. In FY21 and FY22, we designed, built, and tested both a prototype impulse source module for generating seismic energy and a sensing module for recording ground motions generated by the source module array. A pneumatically driven vibratory source was also designed. The source modules were fabricated with all high temperature components and the team has worked to incorporate the current RT-SEISMIC electronics design into a commercially available, high temperature silicon-on-insulator chip integrated circuit. Several issues were identified during fabrication and lab testing that led to redesign of several system components and subsequent retesting. The final round of testing showed that while metal/graphite-based seals worked quite well for static seals, they were unable to provide an adequate gas seal for dynamic, reciprocating part movements which necessitated a final redesign using Kalrez. This change will result in a continuous temperature rating of approximately 275 degrees C for the system. While a field test of the RT-SEISMIC system was targeted in FY22, due to the extended lab testing and redesign efforts, field testing was not achieved. As a result of this LDRD investment, several sponsors have expressed interest in RT-SEISMIC and we expect to continue towards a field demonstration of the full system in the future.

58 GEOSCIENCES↗

Exclusive $J/ψ$ photoproduction on nuclei

Motivated by the recent experimental developments, the Pom-CQM model of the γ+N → J/ψ+N reaction of Lee et al. [Eur. Phys. J. A. 58, 252 (2022)] and Sakinah et al. [Phys. Rev. C. 109, 065204 (2024)] has been applied to predict the exclusive J/ψ photo-production on nuclei (A). Within the multiple scattering theory, the calculations have been performed by including the impulse amplitude $T^{\textrm{IMP}}_{J/ψA,γA}$ and the J/ψ-nucleus final state interaction (FSI) amplitude $T^{\textrm{FSI}}_{J/ψA,γA}$. For the deuteron target, $T^{\textrm{IMP}}_{J/ψd,γd}$ is calculated exactly using the wave function generated from the realistic nucleon-nucleon potentials. Here, it is found that, near the threshold region, the J/ψ photoproduction cross sections depend sensitively on the d-state of the deuteron wave function. The FSI amplitude $T^{\textrm{FSI}}_{J/ψA,γA}$ is calculated using the first-order optical potential constructed from the J/ψ-N scattering amplitude generated from the employed Pom-CQM model. It turns out that the FSI has significant effects in the large momentum-transfer region. By using the conventional fixed scatter approximation (FSA) and the nuclear form factors from the variational Monte-Carlo (VMC) calculations of Lonardoni et al. [Phys. Rev.C. 96, 024326 (2017)], the cross sections of the J/ψ photo-production on 4 He, 16 O, and 40 Ca are also predicted for future experimental investigations at JLab and EIC.

Kim, Sang-Ho [Soongsil University, Seoul (Korea, R↗

The response of accident tolerant fuel cladding to LOCA burst testing: A comparative study of leading concepts

Accident tolerant fuel claddings seek to improve safety margins during loss-of-coolant accident (LOCA) scenarios by reducing cladding steam oxidation rates, hydrogen production, delaying or mitigating cladding burst, and reducing cladding deformation. After rapid and extensive progress in recent years, Cr coated Zr-based alloys, ferritic FeCrAl, and SiC/SiC ceramic matrix composites (CMCs) have emerged as the most promising cladding concepts, but there are few studies comparing the accident tolerance of all three concepts simultaneously. Here, in this study, LOCA burst testing has been conducted on Cr high impulse power magnetron sputtered Zry-4 (Cr/Zry-4), C26M (FeCrAl), and SiC/SiC SiGA (CMC). Experimental observations indicate Cr/Zry-4 burst at higher temperatures than Zry-4 however the effectiveness of the coating on mitigating ZrO 2 formation via Cr 2 O 3 formation depended on physical proximity to the burst opening. No dramatic improvement in diametrical strain nor opening geometry was observed. C26M burst at higher temperatures than Zry-4 at lower pressures, exhibiting small amounts of strain and smaller openings, and formed slow-growing Al 2 O 3 regardless of location. The CMCs did not burst or show any macroscopic signs of deformation, all the way up to an internal cladding pressure of 15.2 MPa. High temperature expanding Nb plug testing indicated that CMC room temperature strength may persist up to 1650°C, with an approximately 40% loss of strength at 1900°C. CMC steam oxidation was also compared to reference chemically vapor deposited SiC up to 1700°C, and CMCs were found to oxidize at slightly accelerated rates that were still orders of magnitude lower than bare Zry-4. Projected H 2 generation during LOCA burst testing showed that when accounting for cladding burst exposing inner diameters to steam, Cr coated Zry-4 reduced H 2 generation by a factor of two compared with bare Zry-4, while both C26M and SiC/SiC reduced H 2 generation by orders of magnitude relative to bare Zry-4.

36 MATERIALS SCIENCE↗

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗