Validated simulations of dynamic crack propagation in single crystals using EFEM and XFEM
Abstract not provided
Zeng, Q.↗
SEARCH · Search NASA
Results for “PHASE VELOCITY”
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.
Liquid structure and temperature invariance of sound velocity in supercooled Bi melt
Structural rearrangement of liquid Bi in the vicinity of the melting point has been proposed due to the unique temperature invariant sound velocity observed above the melting temperature, the low symmetry of Bi in the solid phase and the necessity of overheating to achieve supercooling. The existence of this structural rearrangement is examined by measurements on supercooled Bi. The sound velocity of liquid Bi was measured into the supercooled region to high accuracy and it was found to be invariant over a temperature range of ∼60°, from 35° above the melting point to ∼25° into the supercooled region. The structural origin of this phenomenon was explored by neutron diffraction structural measurements in the supercooled temperature range. These measurements indicate a continuous modification of the short range order in the melt. The structure of the liquid is analyzed within a quasi-crystalline model and is found to evolve continuously, similar to other known liquid pnictide systems. The results are discussed in the context of two competing hypotheses proposed to explain properties of liquid Bi near the melting: (i) liquid bismuth undergoes a structural rearrangement slightly above melting and (ii) liquid Bi exhibits a broad maximum in the sound velocity located incidentally at the melting temperature.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Drag Force and Superfluidity in the Supersolid Stripe Phase of a Spin–Orbit-Coupled Bose–Einstein Condensate
The phase diagram of a spin–orbit-coupled two-component Bose gas includes a supersolid stripe phase, which is featuring density modulations along the direction of the spin–orbit coupling. This phase has been recently found experimentally [31]. In the present work, we characterize the superfluid behavior of the stripe phase by calculating the drag force acting on a moving impurity. Because of the gapless band structure of the excitation spectrum, the Landau critical velocity vanishes if the motion is not strictly parallel to the stripes, and energy dissipation takes place at any speed. Moreover, due to the spin–orbit coupling, the drag force can develop a component perpendicular to the velocity of the impurity. Finally, by estimating the time over which the energy dissipation occurs, we find that for slow impurities, the effects of friction are negligible on a time scale up to several seconds, which is comparable with the duration of a typical experiment.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Gravitational waves from a very strong electroweak phase transition
We investigate the production of a stochastic background of gravitational waves in the electroweak phase transition. We consider extensions of the Standard Model which can give very strongly first-order phase transitions, such that the transition fronts either propagate as detonations or run away. To compute the bubble wall velocity, we estimate the friction with the plasma and take into account the hydrodynamics. We track the development of the phase transition up to the percolation time, and we calculate the gravitational wave spectrum generated by bubble collisions, magnetohydrodynamic turbulence, and sound waves. For the kinds of models we consider, we find parameter regions for which the gravitational waves are potentially observable at the planned space-based interferometer eLISA. In such cases, the signal from sound waves is generally dominant, while that from bubble collisions is the least significant of them. Since the sound waves and turbulence mechanisms are diminished for runaway walls, the models with the best prospects of detection at eLISA are those which do not have such solutions. In particular, we find that heavy extra bosons provide stronger gravitational wave signals than tree-level terms.
79 ASTRONOMY AND ASTROPHYSICS↗
SCATTERING OF THE f-MODE BY SMALL MAGNETIC FLUX ELEMENTS FROM OBSERVATIONS AND NUMERICAL SIMULATIONS
The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models that differ in radius and total magnetic flux. A quiet-Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m = 0 and m = {+-}1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p{sub 1}-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the ensemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase shifts with the azimuthal order m appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.
79 ASTRONOMY AND ASTROPHYSICS↗
Quasar Rain: The Broad Emission Line Region as Condensations in the Warm Accretion Disk Wind
The origin of the broad emission line region (BELR) in quasars and active galactic nuclei is still unclear. I propose that condensations form in the warm, radiation-pressure-driven, accretion disk wind of quasars creating the BEL clouds and uniting them with the other two manifestations of cool (∼10{sup 4} K) gas in quasars, the low ionization phase of the warm absorbers (WAs) and the clouds causing X-ray eclipses. The cool clouds will condense quickly (days to years), before the WA outflows reach escape velocity (which takes months to centuries). Cool clouds form in equilibrium with the warm phase of the wind because the rapidly varying X-ray quasar continuum changes the force multiplier, causing pressure waves to move gas into stable locations in pressure–temperature space. The narrow range of two-phase equilibrium densities may explain the (luminosity){sup 1/2} scaling of the BELR size, while the scaling of cloud formation timescales could produce the Baldwin effect. These dense clouds have force multipliers of order unity and so cannot be accelerated to escape velocity. They fall back on a dynamical timescale (months to centuries), producing an inflow that rains down toward the central black hole. As they soon move at Mach ∼10–100 with respect to the WA outflow, these “raindrops” will be rapidly destroyed within months. This rain of clouds may produce the elliptical BELR orbits implied by velocity-resolved reverberation mapping in some objects and can explain the opening angle and destruction timescale of the narrow “cometary” tails of the clouds seen in X-ray eclipse observations. Some consequences and challenges of this “quasar rain” model are presented, along with several avenues for theoretical investigation.
79 ASTRONOMY AND ASTROPHYSICS↗
Identifying and discriminating phase transitions along decaying shocks with line imaging Doppler interferometric velocimetry and streaked optical pyrometry [Identify phase transition along decaying shocks with line imaging VISAR and streaked optical pyrometry]
Ultrafast line-imaging velocity interferometer system for any reflector and streaked optical pyrometry are now commonly used to obtain high precision equation of state and electronic transport data under dynamic compression at major high energy density science facilities. In this work, we describe a simple way to improve distinguishing phase transformation signatures from other signals when monitoring decaying shock waves. The line-imaging capability of these optical diagnostics offers additional supporting evidence to the assignment of particular anomalies—such as plateaus or reversals—to the occurrence of a phase transition along the Hugoniot. Finally, we illustrate the discussion with two example datasets collected during laser driven shock compression of quartz and stishovite.
47 OTHER INSTRUMENTATION↗
Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon
The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.
Science & Technology - Other Topics↗
Hydrodynamics of a multicomponent vesicle under strong confinement
Herein we numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.
97 MATHEMATICS AND COMPUTING↗
Di-Higgs production in the 4b channel and gravitational wave complementarity
We present a complementarity study of gravitational waves and double Higgs production in the 4 b channel, exploring the gauge singlet scalar extension of the SM. This new physics extension serves as a simplified benchmark model that realizes a strongly first-order electroweak phase transition necessary to generate the observed baryon asymmetry in the universe. In calculating the signal-to-noise ratio of the gravitational waves, we incorporate the effect of the recently discovered significant suppression of the gravitational wave signals from sound waves for strong phase transitions, make sure that supercooled phase transitions do complete and adopt a bubble wall velocity that is consistent with a successful electroweak baryogenesis by solving the velocity profiles of the plasma. The high-luminosity LHC sensitivity to the singlet scalar extension of the SM is estimated using a shape-based analysis of the invariant 4 b mass distribution. We find that while the region of parameter space giving detectable gravitational waves is shrunk due to the new gravitational wave simulations, the qualitative complementary role of gravitational waves and collider searches remain unchanged.
79 ASTRONOMY AND ASTROPHYSICS↗
Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements
Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Determination of the Flow Rate of Coolant in Pipelines with Unstabilized Flow of Fluid in the Area Following a Bend
Questions related to the determination of the flow rate of coolant in the pipelines of electric plants are considered. Amethod of determining the flow of single-phase coolant for the case of an asymmetric unstabilized velocity profile that is present following a bend in a pipeline is proposed. Results from the use of the proposed method in novel experiments on the testing unit of the Élektrogorsk Scientific Research Center and in the steam pipelines of power unit No. 3 of the Balakovo nuclear power plant are presented.
42 ENGINEERING↗
Visualizing pectin polymer-polymer entanglement produced by interfacial water movement
In this report, we investigated the physical conditions for creating pectin polymer-polymer (homopolymer) entanglement. The potential role of water movement in creating pectin entanglement was investigated by placing water droplets-equivalent to the water content of two gel phase films-between two glass phase films and compressing the films at variable probe velocities. Slow probe velocity (0.5 mm/sec) demonstrated no significant debonding. Furthermore, corresponding videomicroscopy demonstrated an occasional water bridge, but no evidence of stranding or polymer entanglement. In contrast, fast probe velocity (5 mm/sec) resulted in 1) an increase in peak adhesion strength, 2) a progressive debonding curve, and 3) increased work of cohesion (p < .001). Corresponding videomicroscopy demonstrated pectin stranding and delamination between pectin films. Scanning electron microscopy images obtained during pectin debonding provided additional evidence of both stranding and delamination. We conclude that water movement can supply the motive force for the rapid chain entanglement between pectin films.
59 BASIC BIOLOGICAL SCIENCES↗
Impact damage of narrow silicon carbide (SiC) ceramics with and without environmental barrier coatings (EBCs) by various foreign object debris (FOD) simulants
The leading and trailing edges of turbine airfoils are highly susceptible to impact by foreign object debris (FOD) entrained in the gas stream of turbine engines. In this study, a narrow specimen geometry is implemented to study the FOD impact behavior of low curvature airfoil edges. Silicon carbide (SiC) ceramic specimens with and without an environmental barrier coating (EBC) are used as target materials. The EBC consists of an air plasma sprayed mullite topcoat (~ 166 µm) and silicon bond coat (~ 58 µm). Spherical (∅ 1.5 mm) FOD simulants of silicon nitride (Si3N4), partially stabilized zirconia (PSZ), and steel are used as impactors. Additionally, a modified gas gun is used to generate impact damage at normal incidence and velocities ranging between 300 to 400 m/s. In situ phase contrast X-ray radiographs, which render 2D representations of the internal transient damage, are captured during the impact experiments. Impactor hardness and toughness are found to primarily influence damage in the coating. Relative to wide specimens, the narrow specimens experience a greater level of cracking and this is postulated to result from reduced self-confinement. Irrespective of specimen geometry and impactor material type, limited cracking (i.e. enhanced shielding) is observed for the coated condition.
36 MATERIALS SCIENCE↗
GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design
GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Saturn's Weather‐Driven Aurorae Modulate Oscillations in the Magnetic Field and Radio Emissions
Abstract The Cassini spacecraft revealed that Saturn's magnetic field displayed oscillations at a period originally thought to match the planetary rotation rate but later found not to. One of many proposed theories predicts that a polar twin‐cell neutral weather system drives this variation, producing observable differences in flows within Saturn's ionosphere. Here, using spectral observations of auroral emission lines taken by the Keck Observatory's Near Infrared Echelle Spectrograph (Keck‐NIRSPEC) in 2017, we derive ion line‐of‐sight velocity maps after grouping spectra into rotational quadrants matching phases of the planetary magnetic field. We measure 0.5 km s −1 wind systems in the ionosphere consistent with predicted neutral twin‐vortex flow patterns. These findings demonstrate that neutral winds in Saturn's polar regions cause the rotational period, as determined via the magnetic field, to exhibit differences and time variabilities relative to the planet's true period of rotation in a process never before seen within planetary atmospheres.
Chowdhury, M. N.↗
Closed‐loop hierarchical control strategies for connected and autonomous hybrid electric vehicles with random errors
This study presents a closed‐loop hierarchical control strategy for a group of connected and autonomous hybrid electric vehicles (HEVs) with the purpose of optimising fuel economy while guaranteeing traffic mobility and vehicle safety. In the hierarchical control architecture, a decentralised stochastic model predictive control (SMPC)‐based higher level controller incorporating signal phase and timing information is formulated to optimise the velocity profile of each vehicle, and an adaptive equivalent consumption minimisation strategy based lower level controller is employed for the energy management control of the HEVs. Creatively, random errors of the control variable are considered in the SMPC framework. The errors are discretised and modelled as a Markov process and the state transition matrix of the errors is generated randomly to capture the error transition dynamics. To solve the SMPC problem more efficiently, a scenario‐based SMPC is employed. Moreover, the propulsion and recuperation efficiencies of the lower level controller are calculated at each time step with measurable variables, and fed back to the higher level controller for velocity optimisation in next time step. Simulation results validate the control effectiveness and advantages of the proposed control architecture.
Qiu, Shaolin↗
Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state
Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y , as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.
36 MATERIALS SCIENCE↗