Generalized continuum mixture theory for multi-material shock physics
Not Available
Schumacher, Shane C.↗
SEARCH · Search NASA
Results for “generalized continuum”
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.
Generalized continuum theory for nematic elastomers: Non-affine motion and characteristic behavior
Not provided.
Materials Science↗
Mechanical analysis of heterogeneous materials with higher-order parameters
Abstract Even though heterogeneous porous materials are widely used in a variety of engineering and scientific fields, such as aerospace, energy-storage technology, and bio-engineering, the relationship between effective material properties of porous materials and their underlying morphology is still not fully understood. To contribute to this knowledge gap, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we use the second-order scheme that is an extension of the first-order computational homogenization framework, where a generalized continuum enables us to introduce length scale into the material constitutive law and capture both pore size and pore distribution. Through several numerical case studies with different combinations of porosity, pore shapes, and distributions, we systematically studied the relationship between the underlying morphology and effective mechanical properties. The results highlight the necessity of higher-order homogenization in understanding the mechanical properties and reveal that higher-order parameters are required to capture the role of realistic pore morphologies on effective mechanical properties. Furthermore, for specific pore shapes, higher-order parameters exhibit dominant influence over the first-order continuum.
42 ENGINEERING↗
Discovery of a z ∼ 0.8 ultra steep spectrum radio halo in the MeerKAT-South Pole Telescope Survey
Radio haloes are diffuse synchrotron sources that trace the turbulent intracluster medium (ICM) of galaxy clusters. However, their origin remains unknown. Two main formation models have been proposed: the hadronic model, in which relativistic electrons are continuously injected by cosmic-ray protons; and the leptonic turbulent re-acceleration model, where cluster mergers re-energize electrons in situ. A key discriminant between the two models would be the existence of ultra-steep spectrum radio haloes (USSRHs), which can only be produced through turbulent re-acceleration. Here, we report the discovery of an USSRH in the galaxy cluster SPT-CLJ2337–5942 at redshift $z = 0.78$ in the MeerKAT-South Pole Telescope 100 deg$^2$ UHF (0.58–1.09 GHz) survey. This discovery is noteworthy for two primary reasons: it is the highest redshift USSRH system to date; and the close correspondence of the radio emission with the thermal ICM as traced by Chandra X-ray observations, further supporting the leptonic re-acceleration model. The halo is underluminous for its mass, consistent with a minor merger origin, which produces steep-spectrum, lower luminosity haloes. This result demonstrates the power of wide-field, high-fidelity, low-frequency ($\lesssim 1$ GHz) surveys like the MeerKAT-SPT 100 deg$^2$ programme to probe the origin and evolution of radio haloes over cosmic time, ahead of the Square Kilometre Array.
X-rays: galaxies↗
Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations
In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.
97 MATHEMATICS AND COMPUTING↗
Diagnosis of 3D magnetic field and mode composition in MHD turbulence with Y -parameter
ABSTRACT Magnetic fields are crucial in numerous astrophysical processes within the interstellar medium (ISM). However, the detailed determination of magnetic field geometry is notoriously challenging. Based on the modern magnetohydrodynamic (MHD) turbulence theory, we introduce a novel statistical technique, the ‘Y-parameter’, to decipher the magnetic field inclination in the ISM and identify dominant turbulence modes. The Y-parameter, calculated as the ratio of anisotropies of different Stokes parameter combinations, displays contrasting trends with the mean-field inclination angle in Alfvénic and compressible turbulence modes. A Y-parameter value around 1.5 ± 0.5 provide a statistical boundary to determine the dominant MHD turbulence modes. We have discovered specific correlations between the Y-parameter value and the inclination angle that unveil the dominant turbulence mode. This methodology, when applied to future radio polarization surveys such as LOFAR and SKA, promises to significantly enhance our knowledge of 3D magnetic field in the ISM and improve our understanding of interstellar turbulence.
79 ASTRONOMY AND ASTROPHYSICS↗
Multi-scale VLBI observations of the candidate host galaxy of GRB 200716C
We present the discovery and the subsequent follow up of radio emission from SDSS J130402.36+293840.6 (J1304+2938), the candidate host galaxy of the gamma-ray burst (GRB) GRB 200716C. The galaxy is detected in the RACS (0.89 GHz), the NVSS, the Apertif imaging survey, and the FIRST (1.4 GHz), the VLASS (3 GHz), and in public LOFAR (130–170 MHz), WISE (3.4–22 μm), and SDSS (z, i, r, g, u filters) data. The luminosity inferred at 1.4 GHz is (5.1 ± 0.2) × 10 30 erg s -1 Hz -1 . To characterise the emission and distinguish between different components within the galaxy, we performed dedicated, high-sensitivity and high-resolution observations with the European VLBI Network (EVN) + e-MERLIN at 1.6 and 5 GHz. We did not detect any emission from a compact core, suggesting that the presence of a radio-loud active galactic nucleus (AGN) is unlikely, and therefore we ascribe the emission observed in the public surveys to star-forming regions within the galaxy. We confirm and refine the redshift estimate, z = 0.341 ± 0.004, with a dedicated Telescopio Nazionale Galileo (TNG) spectroscopic observation. Finally, we compiled a list of all the known hosts of GRB afterglows detected in radio and computed the corresponding radio luminosity: if GRB 200716C belongs to J1304+2938, this is the third most radio-luminous host of a GRB, implying one of the highest star-formation rates (SFRs) currently known, namely SFR ~ 324±61 M ⊙ yr -1 . On the other hand, through the analysis of the prompt emission light curve, recent works suggest that GRB 200716C might be a short-duration GRB located beyond J1304+2938 and gravitationally lensed by an intermediate-mass black hole (IMBH) hosted by the galaxy. Neither the public data nor our Very Long Baseline Interferometry (VLBI) observations can confirm or rule out the presence of an IMBH acting as a (milli-)lens hosted by the galaxy, a scenario still compatible with the set of radio observations presented in this work.
79 ASTRONOMY AND ASTROPHYSICS↗
Towards the design of nature-inspired materials: Impact of complex pore morphologies via higher-order homogenization
Even though the development of novel materials that mimic nature is widely used in a variety of engineering and scientific fields, the relationship between effective material properties and underlying, often complex pore morphology is still not fully understood. To address this knowledge gap and accelerate the development of novel nature-inspired materials, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we create unique pore morphologies with varying levels of complexity that serve as a more realistic representation of natural materials. Here, we then use the second-order homogenization method to capture the role of pore size, shape, orientation, and distribution on effective properties. By creating different pore morphologies, we systematically studied the relationship between morphology and effective mechanical properties. The results highlight the necessity of higher-order parameters to fully capture the role of realistic pore morphologies on effective mechanical properties and provide a path forward in the design of nature-inspired materials.
36 MATERIALS SCIENCE↗
Measurement of the anisotropy power spectrum of the radio synchrotron background
We present the first targeted measurement of the power spectrum of anisotropies of the radio synchrotron background, at 140 MHz, where it is the overwhelmingly dominant photon background. This measurement is important for understanding the background level of radio sky brightness, which is dominated by steep-spectrum synchrotron radiation at frequencies below ν ~ 0.5 GHz and has been measured to be significantly higher than that produced by known classes of extragalactic sources and most models of Galactic halo emission. We determine the anisotropy power spectrum on scales ranging from 2° to 0.2 arcmin with Low-Frequency Array observations of two 18-deg 2 fields – one centred on the Northern hemisphere’s coldest patch of radio sky where the Galactic contribution is smallest and the other offset from that location by 15°. We find that the anisotropy power is higher than that attributable to the distribution of point sources above 100 μJy in flux. Furthermore, this level of radio anisotropy power indicates that if it results from point sources, those sources are likely at low fluxes and incredibly numerous, and likely clustered in a specific manner.
79 ASTRONOMY AND ASTROPHYSICS↗
Young but fading radio sources: searching for remnants among compact steep-spectrum radio sources
The incidence of young but fading radio sources provides important information on the life cycle of radio emission in radio-loud active galactic nuclei. Despite its importance for constraining the models of radio source evolution, there are no systematic studies of remnants in complete samples of young radio sources. We report results of the study of 18 compact steep-spectrum (CSS) radio sources, selected from the statistically complete B3-VLA (Very Large Array) CSS sample, characterized by a steep optically thin spectrum (α ≥ 1.0) and no core detection in earlier studies. Our deep multifrequency VLA, pc-scale Very Long Baseline Array, and enhanced Multi Element Remotely Linked Interferometer Network (e-MERLIN) observations allowed us to locate the core component in 10 objects. In three CSS sources, there is no clear evidence of present-time active regions, suggesting that they are likely in a remnant phase. Among sources with core detection, we find three objects that have no clear active regions (hotspots) at the edges of the radio structure, suggesting that the radio emission may have just restarted. Our results support a power-law distribution of the source ages, although the poor statistics prevents us from setting solid constraints on the percentage of remnants and restarted sources in subpopulations of radio sources.
79 ASTRONOMY AND ASTROPHYSICS↗
Diffuse sources, clustering, and the excess anisotropy of the radio synchrotron background
ABSTRACT We present the largest low frequency (120 MHz) arcminute resolution image of the radio synchrotron background (RSB) to date, and its corresponding angular power spectrum of anisotropies (APS) with angular scales ranging from 3° to 0.3 arcmin. We show that the RSB around the north celestial pole has a significant excess anisotropy power at all scales over a model of unclustered point sources based on source counts of known source classes. This anisotropy excess, which does not seem attributable to the diffuse Galactic emission, could be linked to the surface brightness excess of the RSB. To better understand the information contained within the measured APS, we model the RSB varying the brightness distribution, size, and angular clustering of potential sources. We show that the observed APS could be produced by a population of faint clustered point sources only if the clustering is extreme and the size of the Gaussian clusters is ≲1 arcmin. We also show that the observed APS could be produced by a population of faint diffuse sources with sizes ≲1 arcmin, and this is supported by features present in our image. Both of these cases would also cause an associated surface brightness excess. These classes of sources are in a parameter space not well probed by even the deepest radio surveys to date.
79 ASTRONOMY AND ASTROPHYSICS↗
Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)
Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A More Transparent Infrared Window
The infrared window region (780–1,250 cm –1 , 12.8 to 8.0 μm) is of great importance to Earth's climate due to its high transparency and thermal energy. We present here a new investigation of the transparency of this spectral region based on observations by interferometers of downwelling surface radiance at two DOE Atmospheric Radiation Measurement program sites. We focus on the dominant source of absorption in this region, the water vapor continuum, and derive updated values of spectral absorption coefficients for both the self and foreign continua. Our results show that the self continuum is too strong in the previous version of Mlawer-Tobin_Clough-Kneizys-Davies (MT_CKD) water vapor continuum model, a result that is consistent with other recent analyses, while the foreign continuum is too weak in MT_CKD. In general, the weaker self continuum derived in this study results in an overall increase in atmospheric transparency in the window, although in atmospheres with low amounts of water vapor the transparency may slightly decrease due to the increase in foreign continuum absorption. These continuum changes lead to a significant decrease in downwelling longwave flux at the surface for moist atmospheres and a modest increase in outgoing longwave radiation. The increased fraction of surface-leaving radiation that escapes to space leads to a notable increase (~5–10%) in climate feedback, implying that climate simulations that use the new infrared window continuum will show somewhat less warming than before. This study also points out the possibly important role that aerosol absorption may play in the longwave radiative budget.
58 GEOSCIENCES↗
Data/code for Mlawer et al. 2024 (A more transparent infrared window)
The infrared window region (780–1,250 cm –1 , 12.8 to 8.0 μm) is of great importance to Earth's climate due to its high transparency and thermal energy. We present here a new investigation of the transparency of this spectral region based on observations by interferometers of downwelling surface radiance at two DOE Atmospheric Radiation Measurement program sites. We focus on the dominant source of absorption in this region, the water vapor continuum, and derive updated values of spectral absorption coefficients for both the self and foreign continua. Our results show that the self continuum is too strong in the previous version of Mlawer-Tobin_Clough-Kneizys-Davies (MT_CKD) water vapor continuum model, a result that is consistent with other recent analyses, while the foreign continuum is too weak in MT_CKD. In general, the weaker self continuum derived in this study results in an overall increase in atmospheric transparency in the window, although in atmospheres with low amounts of water vapor the transparency may slightly decrease due to the increase in foreign continuum absorption. These continuum changes lead to a significant decrease in downwelling longwave flux at the surface for moist atmospheres and a modest increase in outgoing longwave radiation. The increased fraction of surface-leaving radiation that escapes to space leads to a notable increase (~5–10%) in climate feedback, implying that climate simulations that use the new infrared window continuum will show somewhat less warming than before. This study also points out the possibly important role that aerosol absorption may play in the longwave radiative budget.
54 ENVIRONMENTAL SCIENCES↗
Noncollinear phases in moiré magnets
Significance Moiré patterns, which result when two or more two-dimensional materials with incommensurate or rotated lattices are layered together, create controllable electronic bands that have recently been shown to induce a tremendous wealth of physical phenomena. Here we initiate the theoretical study of moiré patterns’ influence on magnetic states of localized spins. We construct a general formalism using continuum field theory and present thorough analyses in twisted bilayers of antiferromagnets and also ferromagnets, which are within the experimental reach of Van der Waals materials. Two-dimensional magnets, well known for their strong intrinsic spin fluctuations, may serve as a new platform for moiré effects and open the door to a large class of novel phenomena that were once unimaginable.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al
Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.
36 MATERIALS SCIENCE↗
On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity
The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.
36 MATERIALS SCIENCE↗
A Generalizable Approach to Direct the Self-Assembly of Functional Blue-Phase Liquid Crystals
Blue phases (BPs) are soft and stimuli-responsive photonic crystals that are interesting for sensing, display, and lasing applications. Polycrystallinity has, however, limited both the study and applicability of BPs. Continuum simulations, which lack molecule-specific details, predict that striped chemical patterns, consisting of alternating homeotropic and planar regions, can be used to nucleate single crystals of BPs. Here it is experimentally demonstrated that, independent of the chemistry and complexity of the BP forming material, chemical patterns direct the self-assembly of BPs; these results indicate that a general, thermodynamics-based continuum description of BP self-assembly is adequate for a broad range of materials. When the pattern periodicity equals the BPII unit cell size, single crystals with (100) orientation form for all studied materials. Chemical patterns also promote the growth of BPI crystals with (110) orientation. Importantly, the self-assembly of a photopolymerizable BP is directed and thermally stable, UV-polymerized single crystals are prepared. As a result, the ability to create arbitrarily large, polymeric photonic single crystals with uniform lattice orientation, may have broad implications for optical applications.
36 MATERIALS SCIENCE↗