3D Wave Research Articles

Reproducibility of Cardiac Multifrequency MR Elastography in Assessing Left Ventricular Stiffness and Viscosity.

Cardiac magnetic resonance elastography (MRE) shows promise in assessing the mechanofunctional properties of the heart but faces clinical challenges, mainly synchronization with cardiac cycle, breathing, and external harmonic stimulation. To determine the reproducibility of in vivo cardiac multifrequency MRE (MMRE) for assessing diastolic left ventricular (LV) stiffness and viscosity. Prospective. This single-center study included a total of 28 participants (mean age, 56.6 ± 23.0 years; 16 male) consisting of randomly selected healthy participants (mean age, 44.6 ± 20.1 years; 9 male) and patients with aortic stenosis (mean age, 78.3 ± 3.8 years; 7 male). 3 T, 3D multifrequency MRE with a single-shot spin-echo planar imaging sequence. Each participant underwent two cardiac MMRE examinations on the same day. Full 3D wave fields were acquired in diastole at frequencies of 80, 90, and 100 Hz during a total of three breath-holds. Shear wave speed (SWS) and penetration rate (PR) were reconstructed as a surrogate for tissue stiffness and inverse viscous loss. Epicardial and endocardial ROIs were manually drawn by two independent readers to segment the LV myocardium. Shapiro-Wilk test, Bland-Altman analysis and intraclass correlation coefficient (ICC). P-value <0.05 were considered statistically significant. Bland-Altman analyses and intraclass correlation coefficients (ICC = 0.96 for myocardial stiffness and ICC = 0.93 for viscosity) indicated near-perfect test-retest repeatability among examinations on the same day. The mean SWS for scan and re-scan diastolic LV myocardium were 2.42 ± 0.24 m/s and 2.39 ± 0.23 m/s; the mean PR were 1.24 ± 0.17 m/s and 1.22 ± 0.14 m/s. Inter-reader variability showed good to excellent agreement for myocardial stiffness (ICC = 0.92) and viscosity (ICC = 0.85). Cardiac MMRE is a promising and reproducible method for noninvasive assessment of diastolic LV stiffness and viscosity. 2 TECHNICAL EFFICACY: 1.

A Multiscale Mathematical Model for the Fetal Blood Circulation of the Second Half of Pregnancy.

Doppler ultrasound is a commonly used method to assess hemodynamics of the fetal cardiovascular system and to monitor the well-being of the fetus. Indices based on the velocity profile are often used for diagnosis. However, precisely linking these indices to specific underlying physiology factors is challenging. Several influences, including wave reflections, fetal growth, vessel stiffness, and resistance distal to the vessel, contribute to these indices. Understanding these data is essential for making informed clinical decisions. Mathematical models can be used to investigate the relation between velocity profiles and physiological properties. This study presents a mathematical model designed to simulate velocity wave propagation throughout the fetal cardiovascular system, facilitating the assessment of factors influencing velocity-based indices. The model combines a one-fiber model of the heart with a 1D wave propagation model describing the larger vessels of the circulatory system and a lumped parameter model for the microcirculation. Fetal growth from 20 to 40 weeks of gestational age is incorporated by adjusting cardiac and circulatory parameter settings according to scaling laws. The model's results, including cardiac function, cardiac output distribution, and volume distribution, show a good agreement with literature studies for a growing healthy fetus from 20 to 40 weeks. In addition, Doppler indices are simulated in various vessels and agree with literature as well. In conclusion, this study introduces a novel closed-loop 0D-1D mathematical model that has been verified against literature studies. This model offers a valuable platform for analyzing factors influencing velocity-based indices in the fetal cardiovascular system.

Accuracy verification of the 2D spectral AIS method by the hull monitoring data

Accuracy validation of the 2D spectral AIS method (AIS method using two-dimensional wave spectrum; 2D-AIS method) using hull monitoring data of a large ore carrier is discussed. The AIS method predicts the response of actual ships using AIS-based ship position data, hindcast wave data, and response characteristics data. The conventional AIS method uses representative parameters(H, T, θ) from wave spectrum modeling. The modeling involves many uncertainties. By using the 2D wave spectrum as it is, the effect of uncertainties associated with wave spectrum modeling is investigated. Compared to the analysis of the monitoring data conducted in this study, it was found that the 2D-AIS method improves the accuracy to about 15 %, while the conventional AIS method shows an error of up to about 35 % compared to the actual measurements in the long-term maximum expected values. This enhancement is attributed to the inability of the modeled wave spectrum to accurately represent the two-peak wave spectrum observed in extreme sea states, leading to a deficiency in reproducing wave frequency components corresponding to peak frequencies of the response. The findings underscore the effectiveness of the 2D-AIS method in refining long-term hull stress estimations, thus offering valuable insights for enhancing the safety and performance of hull structures.

Design of 3D vapor chamber for thermal management of permanent magnet synchronous motors

The vapor chamber is a high thermal conductivity device with significant potential for enhancing the power density of permanent magnet synchronous motors. Practical installation requirements and the heat curing process of thermal conductive adhesive impose strict demands on vapor chambers, such as precise 3D profiles and resistance to expansion. This paper presents a novel 3D wave vapor chamber with high thermal conductivity, a design not previously reported. Special internal support structures are engineered to ensure the vapor chamber with performances of thermal conductivity, anti-expansion, and bendability. Simulation results show an anti-expansion deformation of 6.28 % at 125 °C. Heat transfer experiments demonstrate that the 3D vapor chamber achieves an effective thermal conductivity of 4426.73 W/m·K and a thermal resistance of 0.136 °C/W under a 100 W heating load. Furthermore, the 3D vapor chamber-based cooling strategy manages a heat load of 143.93 W at a cutoff temperature of 90 °C, marking an 82.1 % improvement over traditional cooling methods.

A stable decoupled perfectly matched layer for the 3D wave equation using the nodal discontinuous Galerkin method

In outdoor acoustics, the calculations of sound propagating in air can be computationally heavy if the domain is chosen large enough for the sound waves to decay. The computational cost is lowered by strategically truncating the computational domain with an efficient boundary treatment. One commonly used boundary treatment is the perfectly matched layer (PML), which dampens outgoing waves without polluting the computed solution in the inner domain. The purpose of this study is to propose and assess a new perfectly matched layer formulation for the 3D acoustic wave equation, using the nodal discontinuous Galerkin finite element method. The formulation is based on an efficient PML formulation that can be decoupled to further increase the computational efficiency and guarantee stability without sacrificing accuracy. This decoupled PML formulation is demonstrated to be long-time stable, and an optimization procedure for the damping functions is proposed to enhance the performance of the formulation.

Modified three-dimensional nonlinear Schrödinger equation in finite water depth for gravity waves with influence of slowly varying currents

A new modified three-dimensional (3D) nonlinear Schrödinger equation (3DMNLSE) is derived for gravity waves in the presence of wind, dissipation, and two-dimensional slowly varying currents, which include transverse and longitudinal currents in finite water depth. The effect of currents on modulational instability (MI) is investigated. The divergence (convergence) effect of longitudinal favorable (adverse) currents decreases (increases) the MI growth rate and region of instability while suppressing (enhancing) the formation of freak waves. Meanwhile, the transverse currents hardly affect the occurrence of MI. Furthermore, some results from the simulations with the space evolution of 3DMNLSE are presented. The results show the ubiquitous occurrence of freak waves in 3D wave fields under certain sets of initial conditions. We demonstrate that larger waves can be triggered when a weakly modulated wave train enters a region of adverse currents. The maximum amplitude of a freak wave depends on the ratio of the current velocity to the wave phase velocity.

Hirota [formula omitted]-operator forms, multiple soliton waves, and other nonlinear patterns of a 2D generalized Kadomtsev–Petviashvili equation

In the present paper, extensive research is conducted on a 2D generalized Kadomtsev–Petviashvili (2D-gKP) equation which models water waves with long wavelengths. The study starts by establishing Hirota D-operator forms of the governing equation using the Bell polynomial method (BPM). Through checking the three-soliton condition for the 2D-gKP equation, its integrability is then examined, and as a consequence, single-, double-, and (special) triple-soliton waves are derived from the modified Hirota method. In the end, after deriving lump and breather waves of the 2D-gKP equation through ansatz methods, a theoretical discussion along with its proof is presented. Attempts have been made to elucidate the physical features of nonlinear waves by displaying such waves in 2D and 3D postures. The paper’s findings contribute certainly to research on nonlinear waves of KP-type equations.

Detection of natural pulse waves (PWs) in 3D using high frame rate imaging for anisotropy characterization

IntroductionNumerous studies have shown that natural mechanical waves have the potential to assess the elastic properties of the myocardium. When the Aortic and Mitral valves close, a shear wave is produced, which provides insights into tissue stiffness. In addition, the Atrial Kick (AK) generates a wave similar to Pulse Waves (PWs) in arteries, providing another way to assess tissue stiffness. However, tissue anisotropy can also impact PW propagation, which is currently underexplored. This study aims to address this gap by investigating the impact of anisotropy on PW propagation in a phantom.MethodsTube phantoms were created using Polyvinyl Alcohol (PVA). Anisotropy was induced between two sets of two freeze-thaw cycles by stretching and twisting the material. The study first tests and validates the procedure of making helical anisotropic vessel phantoms using the shear wave imaging technique (by estimating the shear wave speed at different probe angles). Using plane wave ultrasound tomography synchronized with a peristaltic pump, 3D high frame rate imaging is performed and used to detect the 3D propagation pattern of PW for each manufactured vessel phantom. Finally, the study attempts to extract the anisotropic coefficient of the vessel using pulse wave propagation angle.ResultsThe Shear wave imaging results obtained for the isotropic vessel show very similar values for each probe angle. On the contrary, the results obtained for the axial anisotropy vessel show a region with a higher shear wave speed at about 0°, corresponding to the long axis of the vessel. Finally, the results obtained for the helical anisotropy depicted increasing shear wave velocity value from −20° to 20°. For the axial phantom, the wavefront of the pulse wave is perpendicular to the long axis of the vessel, while oriented for the helical anisotropic vessels phantom. The pulse wave propagation angle increased with the number of twists made during the vessel manufacturing.DiscussionThe results show that anisotropy can be induced in PVA vessel phantoms by stretching and twisting the material in freeze-thaw cycles. The findings also suggest that vessel anisotropy affects pulse wave propagation angles. Estimating the pulse wave propagation angle may be interesting in characterizing tissue anisotropy in organs where such waves are naturally present.

Deep neural helmholtz operators for 3D elastic wave propagation and inversion

Summary Numerical simulations of seismic wave propagation in heterogeneous 3D media are central to investigating subsurface structures and understanding earthquake processes, yet are computationally expensive for large problems. This is particularly problematic for full waveform inversion, which typically involves numerous runs of the forward process. In machine learning there has been considerable recent work in the area of operator learning, with a new class of models called neural operators allowing for data-driven solutions to partial differential equations. Recent works in seismology have shown that when neural operators are adequately trained, they can significantly shorten the compute time for wave propagation. However, the memory required for the 3D time domain equations may be prohibitive. In this study, we show that these limitations can be overcome by solving the wave equations in the frequency domain, also known as the Helmholtz equations, since the solutions for a set of frequencies can be determined in parallel. The 3D Helmholtz neural operator is 40 times more memory-efficient than an equivalent time-domain version. We employ a Helmholtz neural operator for 2D and 3D elastic wave modeling, achieving two orders of magnitude acceleration compared to a baseline spectral element method. The neural operator accurately generalizes to variable velocity structures and can be evaluated on denser input meshes than used in the training simulations. We also show that when solving for wavefields strictly on the surface, the accuracy can be significantly improved via a graph neural operator layer. In leveraging automatic differentiation, the proposed method can serve as an alternative to the adjoint-state approach for 3D full-waveform inversion, reducing the computation time by a factor of 350.

Fourier-based beamforming for 3D plane wave imaging and application in vector flow imaging using selective compounding

Objective. Ultrafast ultrasound imaging using planar or diverging waves for transmission is a promising approach for efficient 3D imaging with matrix arrays. This technique has advantages for B-mode imaging and advanced techniques, such as 3D vector flow imaging (VFI). The computation load of the cross-beam technique is associated with the number of transmit anglesmand receive anglesn. The full velocity vector is obtained using the least square fashion. However, the beamforming is repeatedm × ntimes using a conventional time-domain delay-and-sum (DAS) beamformer. In the 3D case, the collection and processing of data from different beams increase the amount of data that must be processed, requiring more storage capacity and processing power. Furthermore, the large computation complexity of DAS is another major concern. These challenges translate into longer computational times, increased complexity in data processing, and difficulty in real-time applications.Approach. In response to this issue, this study proposes a novel Fourier domain beamformer for 3D plane wave imaging, which significantly increases the computational speed. Additionally, a selective compounding strategy is proposed for VFI, which reduces the beamforming process fromm × ntom(wheremandnrepresent the number of transmission and reception, respectively), effectively shortening the processing time. The underlying principle is to decompose the receive wavefront into a series of plane waves with different slant angles. Each slant angle can produce a sub-volume for coherent or selective compounding. This method does not rely on the assumption that the plane wave is perfect and the results show that our proposed beamformer is better than DAS in terms of resolution and image contrast. In the case of velocity estimation, for the Fourier-based method, only Tx angles are assigned in the beamformer and the selective compounding method produces the final image with a specialized Rx angle.Main results. Simulation studies andin vitroexperiments confirm the efficacy of this new method. The proposed beamformer shows improved resolution and contrast performance compared to the DAS beamformer for B-mode imaging, with a suppressed sidelobe level. Furthermore, the proposed technique outperforms the conventional DAS method, as evidenced by lower mean bias and standard deviation in velocity estimation for VFI. Notably, the computation time has been shortened by 40 times, thus promoting the real-time application of this technique. The efficacy of this new method is verified through simulation studies andin vitroexperiments and evaluated by mean bias and standard deviation. Thein vitroresults reveal a better velocity estimation: the mean bias is 2.3%, 3.4%, and 5.0% forvx,vy, andvz, respectively. The mean standard deviation is 1.8%, 1.7%, and 3.4%. With DAS, the evaluated mean bias is 9.8%, 4.6%, and 6.7% and the measured mean standard deviation is 7.5%, 2.5%, and 3.9%.Significance. In this work, we propose a novel Fourier-based method for both B-mode imaging and functional VFI. The new beamformer is shown to produce better image quality and improved velocity estimation. Moreover, the new VFI computation time is reduced by 40 times compared to conventional methods. This new method may pave a new way for real-time 3D VFI applications.

Structural Design and Horizontal Wave Force Estimation of a Wall-Climbing Robot for the Underwater Cleaning of Jackets

Currently, divers face significant safety risks when cleaning marine organisms from the steel structures of offshore underwater platform jackets. Consequently, utilizing robots instead of divers to carry out underwater biofouling removal operations will be an important development direction for the underwater maintenance of offshore platforms in the future. In this study, a wall-climbing robot was designed to clean marine organisms from the underwater surface of a platform jacket leg. The overall structure of the underwater cleaning wall-climbing robot is introduced, including the cleaning actuator and the variable curvature-adapted connecting rod mechanism. The corresponding relationship between the variable curvature-adapted connecting rod mechanism and the jacket leg is analyzed in detail. The variable curvature-adapted connecting rod mechanism was optimized using a genetic algorithm to ensure that the underwater cleaning wall-climbing robot can adapt to a minimum diameter of 1 m for the jacket leg. By drawing on Airy wave theory and random wave theory, the Airy wave parameters for waves were analyzed under different sea conditions, considering practical application scenarios. By using Fluent software 2022, a 2D numerical wave tank was constructed to simulate waves under various sea conditions, and the wave surface shapes for different sea states were determined. By building on the Morison equation, a method for calculating the horizontal wave forces on the underwater cleaning wall-climbing robot using the equivalent area and equivalent volume is proposed. By using the two aforementioned methods, the horizontal wave forces on the underwater cleaning wall-climbing robot under specific sea states were determined. The horizontal wave forces of the underwater cleaning wall-climbing robot under different sea conditions were analyzed and simulated in a 3D numerical wave tank. By comparing the theoretical analysis results with the numerical simulation results, where the maximum difference at the extreme points is approximately 11%, the feasibility of the proposed horizontal wave force estimation method was verified.

Automatic recognition of different 3D soliton wave types using deep learning methods

Automatic recognition of different 3D soliton wave types using deep learning methods

An innovative 3D wave channel design for boosting the dynamic response of proton exchange membrane fuel cell stack

The dynamic response of proton exchange membrane fuel cell (PEMFC) is a crucial metric for evaluating these systems in vehicular applications. The forced convection effect can significantly enhance the dynamic response capabilities of proton exchange membrane fuel cells. In this study, we propose the application of the forced convection effect within a parallel flow field. This approach retains the advantages of low pressure drop at the inlets and outlets of the parallel flow field and facilitates rapid distribution of reaction gases, while simultaneously enhancing the dynamic response of actual PEMFC stacks without compromising drainage efficiency. we propose an innovative 3D wave channel design for PEMFC. Specifically, we design and assemble a stack optimized for performance. The proposed design is tested against the traditional 2D straight channel stack, with the 3D wave channel stack demonstrating an increase in peak power density by 8.05% to 33.92% across a relative humidity range of 20% to 100%. This novel design outperforms the 2D straight channel design in terms of gas distribution uniformity, and enhancedmass transfer to the diffusion layer, as well as reduction of concentration polarization in the PEMFC. It also enhances self-humidification under low humidity conditions. Most importantly, the 3D wave channel stack exhibits a superior dynamic response, with voltage fluctuations after current loading being 74.36% to 91.15% lower than those of the 2D straight channel stack. Furthermore, the 3D wave channel stack maintains superior voltage stability under starvation conditions. The findings from this study contribute significantly to optimizing the dynamic response capability of PEMFC systems, highlighting the potential of the 3D wave channel design in enhancing PEMFC performance.

Identification of bored pile defects utilizing torsional low strain integrity test: Theoretical basis and numerical analysis

The torsional low strain integrity test (TLSIT), known for its advantages such as a smaller detection blind zone, improved identification of shallowly buried defects, stable phase velocity for signal interpretation, and better adaptability for existing pile testing. However, it lacks a comprehensive understanding of the authentic three-dimensional (3D) strain wave propagation mechanism, particularly wave reflection and transmission at defects. To address this gap, a novel 3D theoretical framework is introduced in this context to model the authentic 3D wave propagation during the TLSIT. The proposed approach is validated by comparing its results with those obtained from 3D finite element method (FEM) simulations and simplified 1D (one-dimensional) and 3D analytical solutions. Additionally, a parametric study is conducted to enhance insights into the formation mechanism of high-frequency interference observed during the TLSIT. Finally, a defect identification study is performed to provide guidance for interpreting the wave spectrum in terms of defect characteristics.

Experimental investigation of different sub-regimes in horizontal stratified and intermittent gas-liquid two-phase flow: Flow map and analysis of pressure drop fluctuations

Horizontal stratified and intermittent gas-liquid two-phase flows exhibit several sub-regimes. Identifying these correctly would enable the development of more robust predictive models. This study reports on experimental investigation, based on observations and the collection of time series of pressure drop in a 40 mm ID pipe. Nine different sub-regimes (SS, 2D wave, 3D wave, RW, ED+RW, PS+RW, plug, LAS, HAS) were revealed. An original flow pattern map for this diameter is proposed including transitory regions between the sub-regime. Comparison of experimental observations with existing flow maps has highlighted the effect of pipe diameter on transition of sub-regimes. As reported in literature the sub-regimes, could be identified by direct visualization of the pressure drop time series obtained from differential pressure sensor a swell as using Probability Density Function. Furthermore, statistical analysis of standard deviation as function of gas superficial velocity enabled detection of the transition of various sub-regimes. Finally, a space feature based on standard deviation and mixture Froude number is proposed as an identification tool between the different sub-regimes investigated.

Temporal wave dynamics, phase portrait and qualitative analysis of the time-dependent (2+1)-dimensional Zakharov-Kuznetsov equation

In this article, we analyze the effect of time-dependent coefficients and the complex wave dynamics of the (2+1)-dimensional Zakharov-Kuznetsov (ZK) equation. This equation provides a detailed, insightful, and realistic description of space physics, plasma physics, controlled fusion, and nonlinear sciences. The wave solutions are established using the generalized Kudryashov, modified simple equation, and modified sine-Gordon expansion techniques and are illustrated by graphical depictions, which provide valuable insight into understanding the complex dynamics of waves across different physical systems. Exact solitary wave solutions offer a dependable approach to investigating the behavior of a system under particular conditions and facilitating a comprehensive understanding of its dynamics. We also conduct a stability analysis and present the phase portrait of the solutions, which are useful in various fields, including physics, plasma physics, chemistry, biology, economics, and sociology. We ascertain that the profiles of 3D and 2D soliton-shaped waves are significantly affected by dynamic changes in coefficients, wave velocity, and associated model parameters. This research could help clarify the dynamics of intricate systems, paving to a better understanding and analysis of the temporal aspects of various phenomena.

What Can Hurricane Sam (2021) Tell Us About Extreme Ocean Waves Under Tropical Cyclones?

AbstractWe investigate the ocean wave field under Hurricane Sam (2021). Whilst measurements of waves under Tropical Cyclones are rare, an unusually large number of quality in situ and remote measurements are available in that case. First, we highlight the good consistency between the wave spectra provided by the Surface Waves Investigation and Monitoring (SWIM) instrument onboard the China‐France Oceanography SATellite, the in situ spectra measured by National Data Buoy Center buoys, and a saildrone. The impact of strong rains on SWIM spectra is then further investigated. We show that whereas the rain definitely affects the normalized radar cross section, both the innovative technology (beam rotating scanning geometry) and the post‐processing processes applied to retrieve the 2D wave spectra ensure a good quality of the resulting wave spectra, even in heavy rain conditions. On this basis, the satellite, airborne and in situ observations are confronted to the analytical model proposed by Kudryavtsev et al. (2015, https://doi.org/10.1002/2015JC011284). We show that an extended fetch mechanism may be invoked to explain the large significant wave height observed in the right front quadrant of Hurricane Sam.

Patchwork structure of continental lithosphere captured in 3D body wave images of its anisotropic fabrics

This paper presents an overview of research conducted for more than five decades around Vladislav Babuška and collaborators on large-scale seismic anisotropy in tectonically different regions of continental lithosphere in Europe. A wide range of independent data sets and methods are covered. It also briefly touches laboratory measurements of velocity anisotropy on rock samples from the crust and the upper mantle, and emphasizes the importance of considering anisotropy in studies of the Earth structure. The anisotropy is responsible for even larger velocity variations than those due to composition of the most abundant upper mantle rocks (peridotites). The large-scale in-situ measurements of the upper mantle anisotropy capture fabrics of the mantle lithosphere, and enables mapping lateral changes in its structure. The joint inversion/interpretation of the teleseismic body-wave anisotropic parameters, such as variations of directional terms of relative travel time residuals of P waves, shear-wave splitting or the coupled anisotropic-isotropic teleseismic P-wave tomography, image the continental lithosphere as a mosaic of anisotropic domains. Each of the domains has its own thickness and fossil fabric characterized by tilted symmetry axes. We map boundaries of the domains in dependence on the fabric changes. The boundaries can be either narrow and steep or broader and inclined, with an offset relative to boundaries of the related crustal bocks, which can reach several tens of kilometres. This overview presents the European lithosphere-asthenosphere boundary (LAB) and shows examples of anisotropic fabrics of the mantle lithosphere domains and their boundaries in different parts of the European plate.

Topology optimization of blazed gratings under conical incidence

A topology optimization method is presented and applied to a blazed diffraction grating in reflection under conical incidence. This type of grating is meant to disperse the incident light on one particular diffraction order, and this property is fundamental in spectroscopy. Conventionally, a blazed metallic grating is made of a sawtooth profile designed to work with the ±1st diffraction order in reflection. In this paper, we question this intuitive triangular pattern and look for optimal opto-geometric characteristics using topology optimization based on finite element modelling of Maxwell’s equations. In practical contexts, the grating geometry is mono-periodic, but it is enlightened by a 3D plane wave with a wave vector outside of the plane of invariance. Consequently, this study deals with the resolution of direct and inverse problems using the finite element method in this intermediate state between 2D and 3D: the so-called conical incidence. A multi-wavelength objective is used in order to obtain a broadband blazed effect. Finally, several numerical experiments are detailed. Our numerical results show that it is possible to reach a 98% diffraction efficiency on the −1st diffraction order if the optimization is performed on a single wavelength, and that the reflection integrated over the [400,1500] nm wavelength range can be 29% higher in absolute terms, 56% in relative terms, than that of the sawtooth blazed grating when using a multi-wavelength optimization criterion (from 52% to 81%).

3D wave dispersion analysis of graphene platelet-reinforced ultra-stiff double functionally graded nanocomposite sandwich plates with metamaterial honeycomb core layer

This research addresses the three-dimensional thermomechanical wave propagation behavior in sandwich composite nanoplates with a metamaterial honeycomb core layer and double functionally graded (FG) ultra-stiff surface layers. Due to its potential for high-temperature applications, pure nickel (Ni) is preferred for the honeycomb core layer, and an Al2O3/Ni ceramic-metal matrix is preferred for the surface layers. The functional distribution of graphene platelets (GPLs) in three different patterns, Type-U, Type-X, and Type-O, in the metal-ceramic matrix with a power law distribution provides double-FG properties to the surface layers. The mechanical and thermal material characteristics of the core and surface layers, as well as the reinforcing GPLs, are temperature-dependent. The pattern of temperature variation over the plate thickness is considered to be nonlinear. The sandwich nanoplate’s motion equations are obtained by combining the sinusoidal higher-order shear deformation theory (SHSDT) with nonlocal integral elasticity and strain gradient elasticity theories. The wave equations are established by using Hamilton’s principle. Parametric simulations and graphical representations are performed to analyze the effects of honeycomb size variables, wave number, the power law index, the GPL distribution pattern, the GPL weight ratio, and the temperature rise on three-dimensional wave propagation in an ultra-stiff sandwich plate. The results of the analysis reveal that the 3D wave propagation of the sandwich nanoplate can be significantly modified or tuned depending on the desired parameters and conditions. Thus, the proposed sandwich structure is expected to provide essential contributions to radar/sonar stealth applications in air, space, and submarine vehicles in high or low-temperature environments, protection of microelectromechanical devices from high noise and vibration, soft robotics applications, and wearable health and protective equipment applications.