Contact Wave Research Articles

New adaptive low-dissipation central-upwind schemes

We introduce new second-order adaptive low-dissipation central-upwind (LDCU) schemes for the one- and two-dimensional hyperbolic systems of conservation laws. The new adaptive LDCU schemes employ the recently proposed LDCU numerical fluxes computed using the point values reconstructed with the help of adaptively selected nonlinear limiters. To this end, we use a smoothness indicator to detect “rough” parts of the computed solution, where the piecewise linear reconstruction is performed using an overcompressive limiter, which leads to extremely sharp resolution of shock and contact waves. In the “smooth” areas, we use a more dissipative limiter to prevent appearance of artificial kinks and staircase-like structures there. In order to avoid oscillations, we perform the reconstruction in the local characteristic variables obtained using the local characteristic decomposition. We use a smoothness indicator from Löhner (1987) [34] and apply the developed schemes to the one- and two-dimensional Euler equations of gas dynamics. The obtained numerical results clearly demonstrate that the new adaptive LDCU schemes outperform the original ones.

Advances in Cardiovascular Wearable Devices.

Cardiovascular diseases are a leading cause of death worldwide. They mainly include coronary artery disease, rheumatic heart disease, andcerebrovascular disease, and. Cardiovascular diseases can be better managed and diagnosed using wearable devices. Wearable devices, in comparison to traditional cardiovascular diagnostic tools, are not only inexpensive but also have the potential to provide continuous real-time monitoring. This paper reviews some of the recent advances in cardiovascular wearable devices. It discusses traditional implantable devices for cardiovascular diseases as well as wearable devices. The different types of wearable devices are categorized based on different technologies, namely using galvanic contact, photoplethysmography (PPG), and radio frequency (RF) waves. It also highlights the use of artificial intelligence (AI) in cardiovascular disease diagnostics as well as future perspectives on cardiovascular devices.

Frictional Properties of Simulated Fault Gouges Subject to Normal Stress Oscillation and Implications for Induced Seismicity

AbstractUnder critical conditions where experimental fault slip exhibits self‐sustained oscillation, effects of normal stress oscillation (NSO) on fault strength and stability remain poorly understood, as do potential effects of NSO on natural and induced seismicity. In this study, we employed double direct‐shear testing to investigate the frictional behavior of a synthetic, slightly velocity‐weakening (SVW) fault gouge (characterized by self‐sustained oscillation under quasi‐static shear loading), when subjected to NSO at different amplitudes (5%–20% of 5 MPa) and frequencies (0.001–1 Hz). During the experiment, fault displacement and gouge layer thickness were measured. Transmitted ultrasonic waves were also employed to probe grain contact states within the gouge layer. Our results show that fault weakening and unstable slip can be triggered at NSO frequencies ranging from 0.03 to 0.1 Hz and amplitudes exceeding 5%. Interestingly, an amplified shear stress drop and weakening effect were observed when the NSO frequency fell in 0.05–0.1 Hz. Analysis of transmitted ultrasonic waves in tests on the SVW gouge revealed fault dilation, accompanied by unstable slip and weakening. By extending an existing microphysical model (the “Chen‐Niemeijer‐Spiers [CNS]” model), to account for elastic effects of NSO on gouge microstructure and grain contact state, the mechanical and wave data obtained in our experiments on the SVW gouge was reproduced, suggesting an approach for modeling fault instability under upper crustal (SVW) conditions where normal stress is perturbed by subsurface operations, such as periodic gas storage stimulation of reservoir formations.

Stability of viscous contact wave for the full compressible Navier–Stokes–Korteweg equations with large perturbation

This study focuses on the stability of the viscous contact wave for the one-dimensional full Navier–Stokes–Korteweg equations with density-temperature dependent transport coefficients and large perturbation. Our findings demonstrate a Nishida–Smoller type result, indicating that the solution remains stable under large perturbation as long as γ−1 is sufficiently small. Notably, the smallness of the capillary coefficient is unnecessary. We then employ the initial layer analysis technique to investigate the asymptotic behaviour in the Wk,p norm. We show that the capillary term has a smoothing effect, which implies that the strong solution is indeed a smooth one. Our results represent an improvement over those previously reported in Chen and Sheng (2019 Nonlinearity 32 395–444). Furthermore, by applying the method in this study to the isothermal case, we can achieve a better outcome than Germain and LeFloch (2016 Commun. Pure Appl. Math. 69 3–61).

A modified virtual time reversal method for enhancing monitoring sensitivity of bolt preloads based on ultrasonic guided waves

Bolted joints are prone to loosening due to inaccurate initial preloads and time-varying service conditions. Real-time monitoring of bolt preloads is therefore an important means to warrant the reliability and safety of a bolted structure, though it remains a daunting task to implement precise monitoring of early bolt loosening. A modified virtual time reversal (MVTR) method is developed in this study, aimed at enhanced sensitivity of bolt preload monitoring. To begin with, the frequency response function (FRF) of an intact bolted structure is ascertained using the excitation wave and the received wave signal through the bolt. Then a received wave signal captured from the bolted structure under monitoring is time-reversed, and a refocused signal is calculated in virtue of the time-reversed signal and the FRF in the intact state. Bolt tightness indices are established by the peak amplitude or energy of the refocused wave packet. To facilitate interpretation of the principle of the proposed MVTR, a semi-analytic model is developed, which combines the waveform superposition method and finite element method. With the model, the effects of interface contact length, wave mode and signal frequency on the tightness indices are investigated. The effectiveness of the proposed MVTR is experimentally verified using three representative bolted lap structures featuring single bolt or multiple bolts. The results confirm the findings from the semi-analytic model, demonstrating that the MVTR enhances monitoring sensitivity compared to prevailing approaches using linear or nonlinear guided wave features.

An adaptive mesh refinement method based on a characteristic‐compression embedded shock wave indicator for high‐speed flows

AbstractNumerical simulation of high‐speed flows often needs a fine grid for capturing detailed structures of shock or contact wave, which makes high‐order discontinuous Galerkin methods (DGMs) extremely costly. In this work, a characteristic‐compression based adaptive mesh refinement (AMR, h‐adaptive) method is proposed for efficiently improving resolution of the high‐speed flows. In order to allocate computational resources to needed regions, a characteristic‐compression embedded shock wave indicator is developed on incompatible grids and employed as the criterion for AMR. This indicator applies the admissible jumps of eigenvalues to measure the local compression of homogeneous characteristic curves, and theoretically can capture regions of characteristic‐compression which contain structures of shock, contact waves and vortices. Numerical results show that the proposed h‐adaptive DGM is robust, efficient and high‐resolution, it can capture dissipative shock, contact waves of different strengths and vortices with low noise on a rather coarse grid, and can significantly improve resolution of these structures through mild increase of computational resources as compared with the residual‐based h‐adaptive method.

Recognition of Tunnel Fracture Zones in Seismic Waves and Ground-Penetrating Radar Data

Fracture zones in front of tunnel faces can easily cause falling blocks and landslides during the construction process. Using seismic waves and ground-penetrating radar (GPR) data, we extracted the features of fracture zones and achieved the advanced prediction of tunnel fracture zones. The energy variation in the reflected waves propagated by seismic waves at interfaces with different impedances of contact waves was found to manifest as positive and negative reflections, and the amplitude of reflected signals within the fracture zone areas thus increased. We designed a superimposed velocity spectrum, divided the areas of variation in wave velocity, and constructed the three-dimensional spatial distribution of the tunnel fracture zones. Based on the phase change, increase in amplitude, and increase in the center-frequency characteristics of the one-dimensional time waveform of the electromagnetic waves in the fault zone area (A-scan), we located the characteristic points of the fracture zones and observed the occurrence of in-phase axis misalignment in two-dimensional scanning (B-scan). We then implemented the identification of fracture zones. This method predicted the fractured area in the rock surrounding the Liangwangshan Tunnel, and during the tunnel excavation, the fracture zones appeared in the recognition area.

Theory and simulation of shock waves freely propagating through monoatomic non-Boltzmann gas

The effect of non-Boltzmann energy distributions on the free propagation of shock waves through a monoatomic gas is investigated via theory and simulation. First, the non-Boltzmann heat capacity ratio γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma $$\\end{document}, as a key property for describing shock waves, is derived from first principles via microcanonical integration. Second, atomistic molecular dynamics simulations resembling a shock tube setup are used to test the theory. The presented theory provides heat capacity ratios ranging from the well-known γ=5/3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma = 5/3$$\\end{document} for Boltzmann energy-distributed gas to γ→1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma \\rightarrow 1$$\\end{document} for delta energy-distributed gas. The molecular dynamics simulations of Boltzmann and non-Boltzmann driven gases suggest that the shock wave propagates about 9% slower through the non-Boltzmann driven gas, while the contact wave appears to be about 4% faster if it trails non-Boltzmann driven gas. The observed slowdown of the shock wave through applying a non-Boltzmann energy distribution was found to be consistent with the classical shock wave equations when applying the non-Boltzmann heat capacity ratio. These fundamental findings provide insights into the behavior of non-Boltzmann gases and might help to improve the understanding of gas dynamical phenomena.Graphical abstract

Global stability of viscous contact wave for compressible Navier–Stokes equations with temperature‐dependent transport coefficients and large data

We are concerned with the large‐time behavior of the viscous contact wave for a one‐dimensional compressible Navier–Stokes equations whose transport coefficients depend on the temperature. It is shown that if the adiabatic exponent satisfies a suitable inequality and strength is small enough, the unique solution global in time to ideal polytropic gas exists and asymptotically tends toward the viscous contact wave under large initial perturbation for any . Subsequently, we also prove that the large‐time asymptotic stability of the contact wave has a uniform convergence rate under initial data. Our results improve ‐rate in [F. Huang, Z. Xin, and T. Yang, Adv. Math. 219 (2008), 1246–1297] which studied the constant coefficients Navier–Stokes system. The proofs are given by the elementary energy estimate and anti‐derivative method.

Finite-volume two-step scheme for solving the shear shallow water model

<abstract><p>The shear shallow water (SSW) model introduces an approximation for shallow water flows by including the effect of vertical shear in the system. Six non-linear hyperbolic partial differential equations with non-conservative laws make up this system. Shear, contact, rarefaction, and shock waves are all admissible in this model. We developed the finite-volume two-step scheme, the so-called generalized Rusanov (G. Rusanov) scheme, for solving the SSW model. This method is split into two stages. The first one relies on a local parameter that permits control over the diffusion. In stage two, the conservation equation is recovered. Numerous numerical instances were taken into consideration. We clarified that the G. Rusanov scheme satisfied the C-property. We also compared the numerical solutions with those obtained from the Rusanov, Lax-Friedrichs, and reference solutions. Finally, the G. Rusanov technique may be applied for solving a wide range of additional models in developed physics and applied science.</p></abstract>

Stability and decay rate of viscous contact wave to one-dimensional compressible Navier-Stokes equations

Stability and decay rate of viscous contact wave to one-dimensional compressible Navier-Stokes equations

A novel flux splitting based on wave-particle splitting for ideal magnetohydrodynamics

A novel flux splitting method based on wave-particle splitting is developed for one-dimensional ideal magnetohydrodynamics. While ideal magnetohydrodynamics (MHD) equations are non-convex with non-homogeneous flux as opposed to their hydrodynamic counterparts, the present flux splitting methods cannot develop Riemann solver. The proposed approach based on wave-particle method referred as Advection Magnetic-direction Wave Particle Splitting (AMWPS) scheme splits the flux into four parts: the wave-like total pressure part, the wave-like normal magnetic part, the particle-like advection part and the particle-like tangential magnetic part. The particle-like parts can be solved by the scheme developed by Toro and Vazquez-Cendon (TV scheme) and we propose a novel approach to the solution of Riemann problem formed by the wave-like parts including total pressure sub-flux and normal magnetic sub-flux. As for the ordered wave foliation in MHD and other limitations in wave structure, we also analyze the relation between plasma β and specific heat ratio γ. The advantage of AMWPS is that it can sharply capture isolated, stationary Alfven wave discontinuity. Several one-dimensional and two-dimensional MHD problems have been tested to highlight the accuracy, positivity preservation and robustness of AMWPS scheme and comparative studies show that AMWPS significantly outperforms the Riemann solver named Harten, Lax and Leer for contact wave (HLLC) in most cases.

A low diffusion flux-split scheme for all Mach number flows

A low diffusion version of the Harten–Lax–Leer convective pressure split (HLL-CPS) scheme for resolving the shear layers and the flow features at low Mach numbers is presented here. The low diffusion HLL-CPS scheme is obtained by reconstructing the velocities at the cell interface with the face normal Mach number and a pressure function. Asymptotic analysis of the modified scheme shows a correct scaling of the pressure at low Mach numbers and a significant reduction in numerical dissipation. The robustness of the HLL-CPS scheme for strong shock is improved by reducing the contribution of the contact wave in the vicinity of the shock. The improvement in robustness for strong shock is demonstrated analytically through linear perturbation and matrix stability analyses. A set of numerical test cases are solved to demonstrate the efficacy of the proposed scheme over a wide range of Mach numbers.

A simple HLLE-type scheme for all Mach number flows

A simple HLLE-type scheme is proposed for all Mach number flows. In the proposed scheme, no extra wave structure is added in the HLLE scheme to resolve the shear wave while the contact wave is resolved by adding a wave structure similar to the HLLEM scheme. The resolution of the shear layers and the flow features at low Mach numbers is achieved by a velocity reconstruction method based on the face normal Mach number. Robustness against the numerical instabilities is achieved by scaling the velocity reconstruction function in the vicinity of shock with a multi-dimensional pressure sensor. The ability of the proposed scheme to resolve low Mach flow features is demonstrated through asymptotic analysis while the stability of the proposed scheme for strong shock is demonstrated through linear perturbation and matrix stability analyses. A set of numerical test cases are solved to show that the proposed scheme is free from numerical shock instability problems at high speeds and is capable of resolving the flow features at very low Mach numbers.

Damage localization using contact and non-contact narrow frequency band elastic wave generation

In this paper analyses of non-contact elastic wave generation based on an air-coupled transducer (ACT) are presented. Authors compared contact piezoelectric based (PZT) wave generation, with non-contact one using ACT for a narrow frequency band excitation. Elastic wave sensing was based in both cases on a scanning laser Doppler vibrometer. The problem of ACT slope angle for the effectiveness of elastic wave modes generation is analysed. It was shown that for an excitation frequency of 40 kHz it is possible to generate A0 wave mode both in the case of PZT and ACT. Both methods of wave generation are compared based on an analysis of damage localisation results. For this purpose wave irregularity mapping (WIM) algorithm was utilised together with two attenuation compensation methods (based on energy normalisation and weighting function). Specimens in the form of simple square GFRP panel and bonded GFRP aerospace structure with omega stiffeners were investigated. The authors investigated simulated damage (Teflon inserts) and real delamination. In this research very cheap and popular, low-frequency ACT was utilised giving the possibility of realisation of low-cost excitation method. ACTs connected with a laser vibrometry sensing allows of realisation of full non-contact damage localisation method. Authors presented also non typical but very useful application of SLDV allowing to visualise the effect of acoustic wave generation in the air and determination of the pattern of generated acoustic waves. This was based on the photo-acoustic effect.

Asymptotic stability of viscous contact wave to a radiation hydrodynamic limit model

This paper is concerned with the large time behavior of the solutions for 1D radiation hydrodynamic limit model without viscosity and its asymptotic stability of the viscous contact discontinuity wave under the smallness assumption of the strength of the contact wave and initial perturbations. The present pressure includes a fourth-order term about the absolute temperature from radiation effect which brings the main difficulty. Furthermore, the dissipative of the system is weaker for the lack of viscosity. All these make the problem more challenging. The prove is mainly based on the energy method, including normal and radial directions energy estimates.

Small Knudsen rate of convergence to contact wave for the Landau equation

In this paper, we study the hydrodynamic limit to contact discontinuity of the compressible Euler system for the Landau equation with the physical Coulomb interaction as the Knudsen number ε>0 is vanishing. Our results are twofold. First, we construct the unique global-in-time solution to the Landau equation with suitably small initial data as ε>0 is sufficiently small. Second, we prove that the solution of the Landau equation converges to a local Maxwellian whose fluid quantities are the contact discontinuity of the corresponding Euler system uniformly away from the discontinuity as ε→0. Moreover, the uniform convergence rate is also obtained. The main idea is to introduce a new time-velocity weigh function to induce an extra quartic dissipation term for treating the large-velocity growth in the nonlinear estimates due to degeneration of the linearized Landau operator in the Coulomb case.

Decay rate to contact discontinuities for the one-dimensional compressible Navier–Stokes equations with a reacting mixture

In this paper, we investigate the nonlinear stability of contact waves for the Cauchy problem to the compressible Navier–Stokes equations for a reacting mixture in one dimension. If the corresponding Riemann problem for the compressible Euler system admits a contact discontinuity solution, it is shown that the contact wave is nonlinearly stable, while the strength of the contact discontinuity and the initial perturbation are suitably small. Especially, we obtain the convergence rate by using anti-derivative methods and elaborated energy estimates.

Preoperative strain ultrasound elastography can predict occult central cervical lymph node metastasis in papillary thyroid cancer: a single-center retrospective study.

To determine whether preoperative ultrasound elastography can predict occult central cervical lymph node metastasis (CCLNM) in patients with papillary thyroid cancer. This retrospective study included 541 papillary thyroid cancer patients with clinically negative lymph nodes prior to surgery between July 2019 and December 2021. Based on whether CCLNM was present on postoperative pathology, patients were categorized as CCLNM (+) or CCLNM (-). Preoperative clinical data, conventional ultrasound features, and ultrasound elastography indices were compared between the groups. Univariate and multivariate logistic regression analysis were performed to identify the independent predictors of occult CCLNM. A total of 36.60% (198/541) patients had confirmed CCLNM, while 63.40% (343/541) did not. Tumor location, bilaterality, multifocality, echogenicity, margin, shape, vascularity, capsule contact, extrathyroidal extension, aspect ratio, and shear wave elasticity parameters were comparable between the groups (all P > 0.05). Univariate analysis showed statistically significant differences between the two groups in age, sex, tumor size, calcification, capsule invasion, and strain rates ratio in strain ultrasound elastography (all P < 0.05). In multivariate logistic regression analysis, the independent predictors of occult CCLNM were age (OR = 0.975, 95% CI = 0.959-0.991, P = 0.002), sex (OR = 1.886, 95% CI = 1.220-2.915, P = 0.004), tumor size (OR = 1.054, 95% CI = 1.014-1.097, P = 0.008), and strain rates ratio (OR = 1.178, 95% CI = 1.065-1.304, P = 0.002). Preoperative strain ultrasound elastography can predict presence of occult CCLNM in papillary thyroid cancer patients and help clinicians select the appropriate treatment strategy.

Stability of viscous contact wave for the radiative and reactive gas with free boundary

The free boundary problem about the stability of viscous contact wave for the radiative and reactive gas is established by a basic energy method under the small perturbation. The present pressure includes a fourth order term about the absolute temperature from radiation effect as well as the ideal polytropic part, which brings the main difficulty to prove the asymptotic stability of the viscous contact wave.