Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scale coupling”

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 613 records · Page 34

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Prediction of Long-Term Geochemical Change in Bentonite Based on the Interpretative THMC Model of the FEBEX In Situ Test

Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.

Zheng, Liange↗

The relationship between the galactic matter distribution, cosmic-ray dynamics, and gamma-ray production

The concept that cosmic-ray density is related to the matter to which cosmic rays are dynamically coupled through magnetic fields on the scale of galactic arm segments is considered with specific reference to gamma-ray astronomy. A model of the galactic-matter and cosmic-ray distributions is developed on the basis of 21-cm radio surveys and recent observations of the 2.6-mm CO emission line by assuming that cosmic rays are galactic in origin, their column density is proportional to the total interstellar-gas column density, their scale height is considerably larger than that of the matter, and the Galaxy is a spiral with an arm/interarm density ratio of about 3 to 1. It is found that there is a good correlation between the observed gamma-ray intensity and that predicted on the basis of essentially complete coupling of cosmic rays to the best estimate of atomic and molecular hydrogen in the Galaxy. Individual maxima observed in gamma-radiation from the central region of the Galaxy are shown to be well correlated with those predicted to result from certain spiral-arm tangents if all the matter is assumed to be modulated in a particular spiral-arm segment pattern.

Kniffen, D. A.↗

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗

Thermal confinement and transport in spherical tokamaks: a review

Here, we review the thermal plasma confinement and transport properties observed and predicted in low aspect ratio tokamaks, or spherical tokamaks (STs), which can depart significantly from those observed at higher aspect ratio. In particular, thermal energy confinement scalings show a strong, near linear dependence of energy confinement time on toroidal magnetic field, while the dependence on plasma current is more modest, the opposite of what is seen at higher aspect ratio. STs have revealed a very strong improvement in normalized confinement with decreasing collisionality, much stronger than at higher aspect ratio, which bodes well for an ST-based fusion pilot plant should this trend continue at an even lower collisionality than has already been accessed. These differences arise because of fundamental differences in transport in STs due to the more extreme toroidicity (i.e. reduced region of bad curvature), and to the relatively larger $E_r \times B$ shearing rates, both of which can suppress electrostatic drift wave instabilities at both ion and electron gyroradius scales. In addition, electromagnetic effects are much stronger in STs because they operate at high β T . Gyrokinetic (GK) studies, coupled with low- and high-k turbulence measurements, have shed light on the underlying physics controlling transport. At lower β T , both ion- and electron-scale electrostatic drift turbulence may be responsible for transport. At higher β T , microtearing, kinetic ballooning, and hybrid trapped electron/kinetic ballooning modes increasingly play a role, and they have a much stronger impact in the core of ST plasmas than at higher aspect ratio. Flow shear affects the balance between ion- and electron-scale modes. Non-linear GK simulations find regimes where the electron heat flux decreases with decreasing collisionality, consistent with the experimental global normalized confinement scaling. The ST is unique in that the relatively low toroidal magnetic field allows for localized measurements of electron-scale turbulence, and this coupled with turbulence measurements at ion-scales has facilitated detailed comparisons with GK simulations. These data have provided compelling evidence for the presence of ion temperature gradient and electron temperature gradient turbulence in some plasmas, and direct experimental support for the impact of experimental actuators like rotation shear, density gradient and magnetic shear on turbulence and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron Acceleration and Heating during Magnetic Reconnection in the Earth's Quasi-parallel Bow Shock

We perform a 2.5-dimensional particle-in-cell simulation of a quasi-parallel shock, using parameters for the Earth's bow shock, to examine electron acceleration and heating due to magnetic reconnection. The shock transition region evolves from the ion-coupled reconnection dominant stage to the electron-only reconnection dominant stage, as time elapses. The electron temperature enhances locally in each reconnection site, and ion-scale magnetic islands generated by ion-coupled reconnection show the most significant enhancement of the electron temperature. The electron energy spectrum shows a power law, with a power-law index around 6. We perform electron trajectory tracing to understand how they are energized. Some electrons interact with multiple electron-only reconnection sties, and Fermi acceleration occurs during multiple reflections. Electrons trapped in ion-scale magnetic islands can be accelerated in another mechanism. Islands move in the shock transition region, and electrons can obtain larger energy from the in-plane electric field than the electric potential in those islands. These newly found energization mechanisms in magnetic islands in the shock can accelerate electrons to energies larger than the achievable energies by the conventional energization due to the parallel electric field and shock drift acceleration. This study based on the selected particle analysis indicates that the maximum energy in the nonthermal electrons is achieved through acceleration in ion-scale islands, and electron-only reconnection accounts for no more than half of the maximum energy, as the lifetime of sub-ion-scale islands produced by electron-only reconnection is several times shorter than that of ion-scale islands.

Solar magnetic reconnection↗

Modeling of crystallization fouling on a horizontal-tube falling-film evaporator for thermal desalination

Scale formation within horizontal-tube falling-film evaporators is an important issue for thermal desalination, because of the deterioration of heat transfer performance and added maintenance costs. Predicting the fouling process is crucial for desalination plant design, operation, and maintenance. In this paper, a model is developed to predict the crystallization fouling process during seawater falling-film flow over a horizontal tube bundle. The scaling is estimated based on deposition theory, with couple heat transfer including the in-tube steam condensation, conduction through the tube wall and the scale layer, and falling-film evaporation. The spatial and temporal variations of temperature, heat transfer coefficient, and scale thickness are predicted, and the effects of variations in process parameters of steam, feed seawater, and tube properties on the scale thickness, evaporation rate, and heat transfer coefficient for falling-film evaporation are analyzed. Comparisons to existing experimental data show the scaling layer thickness falls in the general range of reported measurements, with a slight overestimation due to neglecting the delay of scaling onset under real conditions. The scale layer thickness increases dramatically as the in-tube steam pressure is increased or the seawater flow rate is decreased. In general, the scaling layer becomes thicker on the lower tubes in the tube bundle, due to the heat and mass transfer of the falling film. The effect of tube material on scaling appears mainly dependent on thermal conductivity; thus, using polymer rather than a stainless steel tube decreases scaling and evaporation rate 96% and 88%, respectively. Finally, this work has potential to guide thermal desalination plant design, material selection, and operating parameter optimization.

42 ENGINEERING↗

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗

SCM overview & the EBR-II shutdown heat removal tests validation.

Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Modulation of Cloud Optical Properties by Vertical Circulations Associated with a Jet Streak Exit Region: The November 26 FIRE Cirrus Case Study

The period from 18 UTC 26 November 1991 to roughly 23 UTC 26 November 1991 has become a focal case study of the FIRE (First International Satellite Cloud Climatology Regional Experiment) Cirrus-11 field campaign. The middle and upper tropospheric cloud data that were collected have allowed FIRE scientists to learn a great deal about the morphological structure and microphysical and radiative characteristics of the mid-latitude cirrus that occurred during that time. An important component of this effort in determining the synoptic scale forcing that existed during this time. By forcing, we mean the coupling between the background vertical air motions and the large scale moisture budget that initiated and maintained cirrus cloud in the study region. Defining the synoptic scale forcing is one of the stated scientific objectives of the FIRE program.

Mace, Gerald G.↗

Gravitational wave probes of massive gauge bosons at the cosmological collider

We extend the reach of the "cosmological collider" for massive gauge boson production during inflation from the CMB scales to the interferometer scales. Considering a Chern-Simons coupling between the gauge bosons and the pseudoscalar inflaton, one of the transverse gauge modes is efficiently produced and its inverse decay leaves an imprint in the primordial scalar and tensor perturbations. We study the correlation functions of these perturbations and derive the updated constraints on the parameter space from CMB observables. We then extrapolate the tensor power spectrum to smaller scales consistently taking into account the impact of the gauge field on inflationary dynamics. Our results show that the presence of massive gauge fields during inflation can be detected from characteristic gravitational wave signals encompassing the whole range of current and planned interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗

Ground test of a large scale 'D' vented thrust deflecting nozzle

Future V/STOL aircraft will require efficient techniques for changing the thrust vector from the vertical direction for VTOL operation to the horizontal direction for conventional flight. Most V/STOL concepts utilize thrust vectoring nozzles to provide this variation in the thrust vector direction. An experimental test program was initiated to demonstrate the capabilities of a large scale 'D' vented thrust deflecting system coupled with a high bypass ratio turbofan engine. Data were obtained for a 'D' vented nozzle mounted behind a YTF-34-F5 turbofan engine. Preliminary data are presented for a variety of test conditions. Attention is given to aspects of 'D' vented nozzle design, the test apparatus, engine-nozzle compatibility, exit area variation, longitudinal vectoring performance, nozzle temperature distribution, and large scale - small scale comparisons.

Rosenberg, E. W.↗

Microwave Remote Sensing and the Cold Land Processes Field Experiment

The Cold Land Processes Field Experiment (CLPX) has been designed to advance our understanding of the terrestrial cryosphere. Developing a more complete understanding of fluxes, storage, and transformations of water and energy in cold land areas is a critical focus of the NASA Earth Science Enterprise Research Strategy, the NASA Global Water and Energy Cycle (GWEC) Initiative, the Global Energy and Water Cycle Experiment (GEWEX), and the GEWEX Americas Prediction Project (GAPP). The movement of water and energy through cold regions in turn plays a large role in ecological activity and biogeochemical cycles. Quantitative understanding of cold land processes over large areas will require synergistic advancements in 1) understanding how cold land processes, most comprehensively understood at local or hillslope scales, extend to larger scales, 2) improved representation of cold land processes in coupled and uncoupled land-surface models, and 3) a breakthrough in large-scale observation of hydrologic properties, including snow characteristics, soil moisture, the extent of frozen soils, and the transition between frozen and thawed soil conditions. The CLPX Plan has been developed through the efforts of over 60 interested scientists that have participated in the NASA Cold Land Processes Working Group (CLPWG). This group is charged with the task of assessing, planning and implementing the required background science, technology, and application infrastructure to support successful land surface hydrology remote sensing space missions. A major product of the experiment will be a comprehensive, legacy data set that will energize many aspects of cold land processes research. The CLPX will focus on developing the quantitative understanding, models, and measurements necessary to extend our local-scale understanding of water fluxes, storage, and transformations to regional and global scales. The experiment will particularly emphasize developing a strong synergism between process-oriented understanding, land surface models and microwave remote sensing. The experimental design is a multi-sensor, multi-scale (1-ha to 160,000 km ^ {2}) approach to providing the comprehensive data set necessary to address several experiment objectives. A description focusing on the microwave remote sensing components (ground, airborne, and spaceborne) of the experiment will be presented.

Kim, Edward J.↗

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of a Reactor Model for Chemical Conversion of Lunar Regolith

Lunar regolith will be used for a variety of purposes such as oxygen and propellant production and manufacture of various materials. The design and development of chemical conversion reactors for processing lunar regolith will require an understanding of the coupling among the chemical, mass and energy transport processes occurring at the length and time scales of the overall reactor with those occurring at the corresponding scales of the regolith particles. To this end, a coupled transport model is developed using, as an example, the reduction of ilmenite-containing regolith by a continuous flow of hydrogen in a flow-through reactor. The ilmenite conversion occurs on the surface and within the regolith particles. As the ilmenite reduction proceeds, the hydrogen in the reactor is consumed, and this, in turn, affects the conversion rate of the ilmenite in the particles. Several important quantities are identified as a result of the analysis. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the time for hydrogen to diffuse into the pores of the regolith particles and the chemical reaction time. The paper investigates the relationships between these quantities and their impact on the regolith conversion. Application of the model to various chemical reactor types, such as fluidized-bed, packed-bed, and rotary-bed configurations, are discussed.

Hedge, uday↗

Development of a Reactor Model for Chemical Conversion of Lunar Regolith

Lunar regolith will be used for a variety of purposes such as oxygen and propellant production and manufacture of various materials. The design and development of chemical conversion reactors for processing lunar regolith will require an understanding of the coupling among the chemical, mass and energy transport processes occurring at the length and time scales of the overall reactor with those occurring at the corresponding scales of the regolith particles. To this end, a coupled transport model is developed using, as an example, the reduction of ilmenite-containing regolith by a continuous flow of hydrogen in a flow-through reactor. The ilmenite conversion occurs on the surface and within the regolith particles. As the ilmenite reduction proceeds, the hydrogen in the reactor is consumed, and this, in turn, affects the conversion rate of the ilmenite in the particles. Several important quantities are identified as a result of the analysis. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the time for hydrogen to diffuse into the pores of the regolith particles and the chemical reaction time. The paper investigates the relationships between these quantities and their impact on the regolith conversion. Application of the model to various chemical reactor types, such as fluidized-bed, packed-bed, and rotary-bed configurations, are discussed.

Hegde, U.↗

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗