Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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 703 records · Page 39

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗

Theoretical investigations on Kerr and Faraday rotations in topological multi-Weyl Semimetals

Motivated by the recent proposal of giant Kerr rotation in WSMs, we investigate the Kerr and Faraday rotations in time-reversal broken multi-Weyl semimetals (mWSMs) in the absence of an external magnetic field. Using the framework of Kubo response theory, we find that both the longitudinal and transverse components of the optical conductivity in mWSMs are modified by the topological charge ( n n ). Engendered by the optical Hall conductivity, we show in the thin film limit that, while the giant Kerr rotation and corresponding ellipticity are independent of n n , the Faraday rotation and its ellipticity angle scale as n n and n^2 n 2 , respectively. In contrast, the polarization rotation in semi-infinite mWSMs is dominated by the axion field showing n n dependence. In particular, the magnitude of Kerr (Faraday) angle decreases (increases) with increasing n n in Faraday geometry, whereas in Voigt geometry, it depicts different n n -dependencies in different frequency regimes. The obtained results on the behavior of polarization rotations in mWSMs could be used in experiments as a probe to distinguish single, double, and triple WSMs, as well as discriminate the surfaces of mWSMs with and without hosting Fermi arcs.

36 MATERIALS SCIENCE↗

Formulation and Performance Evaluation of Epoxy Sealant Systems for Double-Shell Tank Bottom Refurbishment

The performance of epoxy sealants used in the refurbishment of double-shell tank (DST) systems requires balancing processability, thermomechanical stability, and adhesion to cementitious substrates. This study incorporates Heloxy 8 as a reactive diluent into Westlake 862 epoxy to tailor workability and cured-state properties. Rheological time-sweep analysis demonstrates that increasing the diluent content significantly reduces complex viscosity and extends workability, thereby improving pumpability and flow for large-area applications. However, the targeted 2-hour processing window is not fully achieved. Differential scanning calorimetry (DSC) confirms that all formulations cure at room temperature to glass transition temperatures ( T g ) at least 20 °C above the maximum DST operating temperature (27 °C), thereby ensuring service in the glassy regime. Dynamic mechanical analysis (DMA) reveals formulation-dependent reductions in tan delta and increases in storage modulus, indicating increasingly elastic and mechanically stable networks with diluent incorporation. Pull-off adhesion testing shows that modified formulations (70–90% Westlake epoxy) exhibit significantly higher adhesion strengths than the unmodified system. Grout cohesive failure indicates that interfacial bonding exceeds substrate strength. Collectively, these results demonstrate that controlled reactive diluent incorporation enables optimization of processing behavior, interfacial adhesion, and thermomechanical performance, supporting the suitability of the modified epoxy systems as durable sealant layers for cementitious barrier applications in hazardous waste containment infrastructure.

Differential scanning calorimetry↗

A Unified Interpretation of Variability in Precipitation Isotope Ratios

Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.

54 ENVIRONMENTAL SCIENCES↗

Dependence of alpha-particle-driven Alfven eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conven- tional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfven eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently- developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvenic or nearly sub-Alfvenic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Visualizing Anisotropic Oxygen Diffusion in Ceria under Activated Conditions

Oxygen reactivity plays a key role in the performance of ceria-based catalysts. Aberration-corrected transmission electron microscopy and molecular dynamics simulations were used to study the oxygen atom diffusion in ceria under activated conditions. Reactive oxygen atom and its real-time diffusion were visualized. Additionally, the interplay between cerium and oxygen atoms originating from a Coulomb interaction was revealed by the out-of-plane buckling of cerium atoms associated with oxygen transport. Anisotropic oxygen atom diffusion that depends on crystal orientations was discovered, demonstrating a preferential [001] crystallographic diffusion pathway. These findings reveal prospects for applications of anisotropic orientation-relevant fluorite-structured oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Singularities in nearly uniform one-dimensional condensates due to quantum diffusion

Dissipative systems often exhibit wavelength-dependent loss rates. One prominent example is Rydberg polaritons formed by electromagnetically induced transparency, which have long been a leading candidate for studying the physics of interacting photons and also hold promise as a platform for quantum information. In this system, dissipation is in the form of quantum diffusion, i.e., proportional to k 2 (k being the wavevector) and vanishing at long wavelengths as k→0. Here, we show that one-dimensional condensates subject to this type of loss are unstable to long-wavelength density fluctuations in an unusual manner: after a prolonged period in which the condensate appears to relax to a uniform state, local depleted regions quickly form and spread ballistically throughout the system. We connect this behavior to the leading-order equation for the nearly uniform condensate—a dispersive analog to the Kardar-Parisi-Zhang equation—which develops singularities in finite time. Furthermore, we show that the wavefronts of the depleted regions are described by purely dissipative solitons within a pair of hydrodynamic equations, with no counterpart in lossless condensates. We close by discussing conditions under which such singularities and the resulting solitons can be physically realized.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Mountaintop View Requires Minimal Sorting: A Faster Contour Tree Algorithm

Consider a scalar field f : M → R, where M is a triangulated simplicial mesh in R d . A level set, or contour, at value v is a connected component of f –1 (v). As v is changed, these contours change topology, merge into each other, or split. Contour trees are concise representations of f that track this contour behavior. The vertices of these trees are the critical points of f, where the gradient is zero. The edges represent changes in the topology of contours. It is a fundamental data structure in data analysis and visualization, and there is significant previous work (both theoretical and practical) on algorithms for constructing contour trees. Suppose M has n vertices, N facets, and t critical points. A classic result of Carr, Snoeyink, and Axen (2000) gives an algorithm that takes O(n log n+Nα(N)) time (where α(·) is the inverse Ackermann function). A further improvement to O(t log t + N) time was given by Chiang et al. All these algorithms involve a global sort of the critical points, a significant computational bottleneck. Unfortunately, lower bounds of Ω(t log t) also exist. We present the first algorithm that can avoid the global sort and has a refined time complexity that depends on the contour tree structure. Intuitively, if the tree is short and fat, we get significant improvements in running time. For a partition of the contour tree into a set of descending paths, P, our algorithm runs in O($\Sigma$ pϵP |p| log |p| + tα(t) + N). This is at most O(t log D + N), where D is the diameter of the contour tree. Moreover, it is O(tα(t) + N) for balanced trees, a significant improvement over the previous complexity. Our algorithm requires numerous ideas: partitioning the contour tree into join and split trees, a local growing procedure to iteratively build contour trees, and the use of heavy path decompositions for the time complexity analysis. There is a crucial use of a family of binomial heaps to maintain priorities, ensuring that any comparison made is between comparable nodes of the contour tree. We also prove lower bounds showing that the $\Sigma$ pϵP |p| log |p| complexity is inherent to computing contour trees.

97 MATHEMATICS AND COMPUTING↗

Numerically Relevant Timescales in the MG2 Microphysics Model

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. Here, we hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate.

58 GEOSCIENCES↗

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗

Linear Viscoelastic Properties of Putative Cyclic Polymers Synthesized by Reversible Radical Recombination Polymerization (R3P)

Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hygrothermal Performance of Bio-Based Materials in Residential Building Wall Envelope Systems

In the United States, building envelope systems contribute approximately 28% of building energy consumption, with walls being a primary factor. As the construction industry explores new approaches to lower the energy intensity needed to manufacture materials and, in turn, reduce energy consumption, bio-based materials are gaining attention as a promising alternative. While prior studies have demonstrated that bio-based materials can decrease material energy intensity, their long-term hygrothermal behavior within a whole residential wall system is not yet fully understood. This study aims to evaluate the potential of replacing Oriented Strand Board (OSB) and drywall in residential wall systems with bio-based alternatives. The research focuses on two bio-based materials, assessing their hygrothermal performance under Chicago climate condition (Climate Zone 5). To achieve that, a series of laboratory experiments were conducted to measure key properties of the bio-based materials such as thickness, dimensions, density, thermal conductivity (as a function of temperature and humidity), moisture-dependent permeance and water vapor permeability (WVP), and sorption isotherms. These properties were then used in WUFI® simulations to predict the moisture durability of a standard residential wall system in Climate Zone 5. Results show that substituting OSB and drywall with bio-based materials can achieve acceptable moisture durability, effectively mitigating risks of mold growth and structural damage over time.

Palani, Hevar [ORNL] (ORCID:0000000220211994)↗

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗

Surrogate approximation of the Grad–Shafranov free boundary problem via stochastic collocation on sparse grids

In magnetic confinement fusion devices, the equilibrium configuration of a plasma is determined by the balance between the hydrostatic pressure in the fluid and the magnetic forces generated by an array of external coils and the plasma itself. The location of the plasma is not known a priori and must be obtained as the solution to a free boundary problem. The partial differential equation that determines the behavior of the combined magnetic field depends on a set of physical parameters (location of the coils, intensity of the electric currents going through them, magnetic permeability, etc.) that are subject to uncertainty and variability. The confinement region is in turn a function of these stochastic parameters as well. In this work, we consider variations on the current intensities running through the external coils as the dominant source of uncertainty. This leads to a parameter space of dimension equal to the number of coils in the reactor. With the aid of a surrogate function built on a sparse grid in parameter space, a Monte Carlo strategy is used to explore the effect that stochasticity in the parameters has on important features of the plasma boundary such as the location of the x-point, the strike points, and shaping attributes such as triangularity and elongation. Furthermore, the use of the surrogate function reduces the time required for the Monte Carlo simulations by factors that range between 7 and over 30.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Corrosion Response of Ni-19Cr-5Fe in Molten NaCl-MgCl 2 Salt

Ni-based alloys are strong candidates for use in high-temperature molten salt reactors due to their superior corrosion resistance and mechanical stability compared to stainless steels. In this study, the static corrosion behavior of a Ni-19Cr-5Fe model alloy was systematically evaluated in a purified molten NaCl-MgCl 2 salt at 700°C for 30, 240, and 500 h. Post-exposure analyses were conducted to assess microstructural evolution, corrosion depth, and elemental depletion profiles. Corrosion rates, quantified by chromium depletion depth, followed an inverse power-law trend with increasing exposure time, indicating diffusion-limited kinetics. This trend is attributed to limitations in mass transport in the static salt and to the progressive local depletion of reactive chromium species. Coupled electron backscatter diffraction and energy dispersive X-ray spectroscopy analysis further revealed that the grain boundary character significantly influences corrosion susceptibility: high-angle grain boundaries exhibited pronounced Cr depletion and pitting, while low-angle and Σ3 boundaries remained comparatively resistant. These results offer valuable insight into the role of microstructure in corrosion processes and reinforce the importance of time-dependent material evaluation in molten salt environments relevant to advanced reactor designs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Standoff Infrared Measurements of Chemical Plume Dynamics in Complex Terrain Using a Combination of Active Swept-ECQCL Laser Spectroscopy with Passive Hyperspectral Imaging

Chemical plume detection and modeling in complex terrain present numerous challenges. We present experimental results from outdoor releases of two chemical tracers (sulfur hexafluoride and Freon-152a) from different locations in mountainous terrain. Chemical plumes were detected using two standoff instruments collocated at a distance of 1.5 km from the plume releases. A passive long-wave infrared hyperspectral imaging system was used to show time- and space-resolved plume transport in regions near the source. An active infrared swept-wavelength external cavity quantum cascade laser system was used in a standoff configuration to measure quantitative chemical column densities with high time resolution and high sensitivity along a single measurement path. Both instruments provided chemical-specific detection of the plumes and provided complementary information over different temporal and spatial scales. The results show highly variable plume propagation dynamics near the release points, strongly dependent on the local topography and winds. Effects of plume stagnation, plume splitting, and plume mixing were all observed and are explained based on local topographic and wind conditions. Measured plume column densities at distances ~100 m from the release point show temporal fluctuations over ~1 s time scales and spatial variations over ~1 m length scales. The results highlight the need for high-speed and spatially resolved measurement techniques to provide validation data at the relevant spatial and temporal scales required for high-fidelity terrain-aware microscale plume propagation models.

47 OTHER INSTRUMENTATION↗