Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave setup”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗

Experimental demonstration of externally driven millimeter-wave particle accelerator structure

We report the experimental demonstration of a mm-wave electron accelerating structure powered by a high-power rf source. We demonstrate reliable coupling of an unprecedented rf power—up to 575kW into the mm-wave accelerator structure using a quasi-optical setup. This standing wave accelerating structure consists of a single-cell copper cavity and a Gaussian to TM 01 mode converter. The accelerator structure is powered by 110 GHz, 10-ns long rf pulses. These pulses are chopped from 3 microsecond pulses from a gyrotron oscillator using a laser-driven silicon switch. We show an unprecedented high gradient up to 230MV/m that corresponds to a peak surface electric field of more than 520 MV/m. We have achieved these results after conditioning the cavity with more than 100,000 pulses. We also report preliminary measurements of rf breakdown rates, which are important for understanding rf breakdown physics in the millimeter-wave regime. Furthermore, these results open up many frontiers for applications not only limited to the next generation particle accelerators but also x-ray generation, probing material dynamics, and nonlinear light-matter interactions at mm-wave frequency.

43 PARTICLE ACCELERATORS↗

Geometry Optimization of Cable-Based Actuation for Small-Scale Model Testing of a Floating Marine Turbine: Preprint

This research aims to apply combined wave and tidal current loads to a small-scale floating marine current turbine in a wave tank, where an actuation system applies hydrodynamic and mooring forces on the hardware based on results from a simulation. We use a real-time hybrid test setup with physical wave forcing from the wave tank and simulated current and mooring forces implemented through a tensioned cable array. For this paper, we developed an optimization algorithm that adjusts the hardware geometry of the cable array to achieve more efficient tension allocation across the cables for the loads that need to be actuated. By adjusting the points where the cables are attached to the floating platform and the angles between the platform and the winches, an optimal cable geometry can be found to minimize tension variations in the lines, maintain the desired pretension, and prevent excessive tensions. We present the optimization problem formulation, the actuation system evaluation approach, and optimization results that show effective cable actuation setups that are being considered for implementation in the wave tank tests.

cable-based actuation↗

Highly Accurate Experimental Heave Decay Tests with a Floating Sphere: A Public Benchmark Dataset for Model Validation of Fluid–Structure Interaction

Highly accurate and precise heave decay tests on a sphere with a diameter of 300 mm were completed in a meticulously designed test setup in the wave basin in the Ocean and Coastal Engineering Laboratory at Aalborg University, Denmark. The tests were dedicated to providing a rigorous benchmark dataset for numerical model validation. The sphere was ballasted to half submergence, thereby floating with the waterline at the equator when at rest in calm water. Heave decay tests were conducted, wherein the sphere was held stationary and dropped from three drop heights: a small drop height, which can be considered a linear case, a moderately nonlinear case, and a highly nonlinear case with a drop height from a position where the whole sphere was initially above the water. The precision of the heave decay time series was calculated from random and systematic standard uncertainties. At a 95% confidence level, uncertainties were found to be very low—on average only about 0.3% of the respective drop heights. Physical parameters of the test setup and associated uncertainties were quantified. A test case was formulated that closely represents the physical tests, enabling the reader to do his/her own numerical tests. The paper includes a comparison of the physical test results to the results from several independent numerical models based on linear potential flow, fully nonlinear potential flow, and the Reynolds-averaged Navier–Stokes (RANS) equations. A high correlation between physical and numerical test results is shown. The physical test results are very suitable for numerical model validation and are public as a benchmark dataset.

42 ENGINEERING↗

Exploration of Water Jets in Supersonic Crossflow Using X-Ray Diagnostics

Here, near-field structures of pure- and aerated-liquid jets injected into Mach 2 crossflow environments were experimentally investigated at the 7-BM beamline of the Advanced Photon Source at Argonne National Laboratory using high-speed shadowgraph imaging, high-speed x-ray imaging, and pathlength-integrated x-ray fluorescence. Liquid was flush injected into a blowdown supersonic wind tunnel. Test section top and side windows were fitted with a polycarbonate plate for high x-ray transmittance. An axisymmetric aerated-liquid injector was fitted with an exchangeable nozzle adaptor to generate a liquid jet. Water and nitrogen, the gas used for aerated-liquid jets, were doped with specially selected dopants to facilitate x-ray measurements. Liquid column structures for a pure-liquid jet, typically masked by dense droplet clouds in a supersonic crossflow environment, were visualized with the present x-ray imaging setup. Surface-wave formation, movement, and transformation on the column windward surface of the pure-liquid jet, along with column deformation and breakup processes, were qualitatively characterized. We measured time-averaged line-of-sight (LOS) liquid mass distributions within the liquid jet near fields. Cross-sectional liquid column contours of pure-liquid jets at various deformation stages were reconstructed from LOS liquid mass distributions. The pure-liquid jet deformed liquid column reached a maximum width of 1.47 times the orifice diameter before exhibiting rapid mass removal. We also observed the plume-crossing phenomenon in the aerated-liquid jet.

47 OTHER INSTRUMENTATION↗

Performance Analysis of PIConGPU: Particle-in-Cell on GPUs using NVIDIA’s NSight Systems and NSight Compute

PIConGPU, Particle In Cell on GPUs, is an open source simulations framework for plasma and laser-plasma physics used to develop advanced particle accelerators for radiation therapy of cancer, high energy physics and photon science. While PIConGPU has been optimized for at least 5 years to run well on NVIDIA GPU-based clusters, there has been limited exploration by the development team of potential scalability bottlenecks using recently updated and new tools including NVIDIA’s NVProf tool and the brand-new NVIDIA NSight Suite (Systems and Compute) tools. PIConGPU is a highly optimized application that runs production jobs at scale on a system Oak Ridge Leadership Facility’s (OLCF) Summit supercomputer (using the full machine at 4600 nodes; at 98% of GPU utilization on all ~28000 NVIDIA Volta GPUs). PIConGPU has been selected as one of the the eight applications for OLCF’s coveted Center for Accelerated Application Readiness (CAAR) program aimed at the facility’s Frontier supercomputer (OLCF’s first exascale system to launch in 2021), to partner with our vendors (primary vendors: AMD and Cray/HPE) ensuring that Frontier will be able to perform large-scale science when it opens to users in 2022. To this effect, performance engineers on the PIConGPU team wanted to dive deep into the application to understand at the finest granularity, which portions of the code could be further optimized to exploit the hardware on Summit at it’s maximum potential and also to elucidate which key kernels should be tracked and optimized for the CAAR effort to port this code to Frontier. Any bottlenecks that are observed via performance profiling on Summit are likely to also impact scalability on the Frontier-dev system and the Frontier Early Access (EA) system. Additionally, the engineers wanted to take a closer look at the newest NVIDIA profiling tools which allows us to identify the most useful features on these tools and will provide an opportunity to compare it to new AMD and Cray’s performance analysis tool releases and provide feedback to our vendor partners on what features are most important and mission critical for CAAR efforts. The primary goal of this report is to focus on the evaluation of PIConGPU’s most time-intensive kernels using NVProf and NSight Suite. Three kernels, Current Deposition (also known as Compute Current), Particle Push (Move and Mark), and Shift Particles are known to be some of the most time-consuming kernels in PIConGPU. The Current xi Deposition kernel and Particle Push kernel both set up the particle attributes for running any physics simulation with PIConGPU, so it is crucial to improve the performance of these two kernels. In this report, we measure single GPU metrics for the three kernels, offer high level takeaways from the conducted analysis, and compare the profiling data from NSight Compute to that of NVProf. This analysis was performed using a grid size of 240 x 272 x 224, and 10 time steps with the Mid-November Figure of Merit (FOM) run setup. The Traveling Wave Electron Acceleration (TWEAC) science case used in this run is a representative science case for PIConGPU. This execution can also be used for baseline analysis on AMD MI50/ MI60 systems. As of the time of writing, the PIConGPU application has limited use for features of NSight Systems, so this report will mainly focus on insights garnered from NSight Compute. For this analysis, we run the “full” metric set available in NSight Compute version 2020.1.2 and use NSight Systems version 2020.3.1 to generate the application timeline.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Features of correlation measurements of the parameters of pulsed hyperspectral optical fields using an asymmetric interferometer

Differences in correlation measurements of the parameters of pulsed hyperspectral optical fields using symmetric and asymmetric interferometers are considered. It is shown analytically that the resulting cross-correlation function is sensitive to phase perturbations in the original wave field. The considered setup, which contains a telescopic reflective 4f system of parabolic mirrors in one arm, demonstrates that in the case of an asymmetric interferometer, the presence of aberrations leads to degradation of the reconstructed image, whereas in the case of symmetric interferometers these aberrations do not affect the result. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fast Dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object's luminosity. We describe a dissipation mechanism that works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields B 1 and B 2 that propagate in opposite directions along the background magnetic field B 0 and collide. The collision forms a thin current sheet perpendicular to B 0 , which absorbs the incoming waves. The current sheet is sustained by an electric field E breaking the magnetohydrodynamic condition E < B and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude A = B 1 / B 0 = B 2 /B 0 propagating in a uniform B 0 . The mechanism is activated when A > 1/2. It dissipates a large fraction of the wave energy, † = (2A - 1)/A 2 , reaching 100% when A = 1. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate A of main interest, the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at A $\gg$ 1. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗

Experimental study of high-energy single-pulse nanosecond discharges in pin-to-pin configuration

Single-pulse pin-to-pin nanosecond pulsed discharges in ambient air (gap distances of 3–7 mm) were investigated at high pulse energies (∼20–30 mJ deposited energy per 11 ns full width at half maximum pulse). A Michelson interferometry setup (532 nm continuous-wave laser) was employed to record time-resolved interferograms of the discharge, enabling spatially resolved calculation of the electron number density. An intensified charge-coupled device camera was used to capture the spatiotemporal evolution of the discharge (streamer formation and spark channel development), and Coherent Anti-Stokes Raman Spectroscopy (CARS) was used to measure post-discharge N 2 vibrational temperatures. Discharge current and voltage were monitored with a back-current shunt. The discharge was initiated with simultaneous cathode-directed and anode-directed streamers that bridge the gap within ∼1 ns and form a luminous, filamentary plasma channel immediately after breakdown. Laser interferometry measurements showed peak electron number densities of the order of 10 1 7–10 18 cm −3 , occurring about 15–20 ns after pulse arrival at the discharge gap. Shorter gap discharges yielded higher peak electron densities, consistent with the higher energy density in the smaller gaps. Spatially, the electron density was highest near the electrodes and decreased toward the midgap region, with 5 and 7 mm gaps exhibiting a pronounced drop in the central channel. CARS indicated initial vibrational temperatures of approximately 4000–6000 K in the spark core of ∼50 ns after the discharge onset, decaying on a ∼500 ns timescale as the plasma cooled and recombined.

Coherent anti-Stokes Raman spectroscopy↗

High resolution time- and angle-resolved photoemission spectroscopy with 11 eV laser pulses

Performing time- and angle-resolved photoemission (tr-ARPES) spectroscopy at high momenta necessitates extreme ultraviolet laser pulses, which are typically produced via high harmonic generation (HHG). Despite current advances, HHG-based setups still require large pulse energies (from hundreds of μ J to mJ) and their energy resolution is limited to tens of meV. Here, we present a novel 11 eV tr-ARPES setup that generates a flux of 5 × 10 10 photons/s and achieves an unprecedented energy resolution of 16 meV. It can be operated at high repetition rates (up to 250 kHz) while using input pulse energies down to 3 µ J. We demonstrate these unique capabilities by simultaneously capturing the energy and momentum resolved dynamics in two well-separated momentum space regions of a charge density wave material ErTe 3 . This novel setup provides the opportunity to study the non-equilibrium band structure of solids with exceptional energy and time resolutions at high repetition rates.

47 OTHER INSTRUMENTATION↗

Electronic and Vibrational Coherence in Heterogeneous Electron Transfer (Final Report)

Interfacial electron transfer reactions constitute key physical phenomena central to a variety of energy related transport and conversion processes such as catalysis, photocatalysis, photovoltaics, energy storage, molecular electronics, etc. Heterogeneous material systems like organic/inorganic interfaces are of particular interest because they provide great potential to tailor properties according to application specific requirements. While electron transfer dynamics has been studied intensively in heterogeneous systems many questions remain unanswered. The significance of electronic-vibrational coupling and coherence is one aspect that has yet to be fully explored. The objective of this project is to develop and implement new ultrafast spectroscopic techniques based on pump four-wave mixing spectroscopy (pump-FWM) for studying the dynamics of vibrational coupling and coherence in heterogeneous electron transfer (HET) at molecule/semiconductor interfaces. The implementation of the spectroscopic setup allows to explore different four-wave mixing techniques and thus to identify a suitable approach to study the significance of vibrational as well as electronic coherence in HET.

14 SOLAR ENERGY↗

Multidimensional four-wave-mixing spectroscopy with squeezed light

A set of nonlinear spectroscopic measurements with quantum light is presented here. It involves a strong pump with frequency ω pu and a weak probe atω pr interacting with a solid state target to generate a conjugated beam ω c =2ω pu -ω pr via nondegenerate four-wave mixing. The $\chi^{(3)}$ susceptibility can be measured by the noise spectra of the intensity difference of the squeezed beams. We discuss three spectroscopic setups based on squeezed light: first, single four-wave mixing in a single crystal; second, cascading scheme involving two crystals; third, an $\textit{SU}$(1, 1) interferometer based on two separate four-wave mixing processes. We further investigate the microscopic noise and optical losses in all three setups. Simulations are presented for silicon-vacancy color centers in diamond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗

Controlled Injection of Electrons for Improved Performance of Laser-Wakefield Acceleration

Our lab performed a series of experiments on controlled injection of LWFA with two extreme laser pulses: one laser pulse was used to drive LWFA. In contrast, the other precisely controls the injection of electrons into the wakefield plasma wave. We have demonstrated an efficient reduction of electron beam emittance by optical injection and novel spatio-temporal manipulation of the electron beam by the extreme laser-plasma waves. Moreover, the optical injection process involves unprecedented new physics, which will benefit studies in astrophysics, inertial confinement fusion, laser-driven x-ray light sources, and ultrafast science. We further improved the performance of LWFA by a better understanding and control of the extreme laser-plasma wave dynamics for optical injection. Three novel mechanisms were identified in the optical injection process: (1) collision of two laser pulses, (2) collision of a laser pulse with a wakefield plasma wave, and (3) collision of two laser wakefields. All three methods improved electron-beam energy spread, emittance, stability, and reproducibility. First, we identified, compared, and independently controlled these different injection mechanisms by controlling the polarization of the injector laser pulse. Second, we manipulated the spatiotemporal splitting of LWFA electron beams to reduce their emittance. Third, we applied the same laser-plasma-wave dynamics and the same experimental setup for improved performance of ion acceleration, which has also been predicted by the recent simulations.

43 PARTICLE ACCELERATORS↗

TRUST Nonlinear Dynamics Testbed Assessment

The following assessment evaluates the efficacy of the control script for carrying out a linear signal to generate a linear mechanical response of the system: the TRUST nonlinear dynamics (TRUSTND) testbed. The hardware has three main components: the controller (NI PXIe 8861) embedded in the National Instruments chassis (NI PXIe 1092), the signal amplifier (The Modal Shop Linear Power Amplifier 2050E09) to amplify the output of a custom LabVIEW script, and a shaker (The Modal Shop Electrodynamic Exciter 2075E) where wave spring specimens live and are tested within the attached aluminum testbed (center mass). Different kinds of signals (single-tone, swept frequency and white noise) were sent through this hardware in the E-1 lab space at TA-53. First, the linearity of the setup was tested by specifying a sine wave in the control script and checking the quality of the oscillations in a wave spring (McMaster-Carr 9714K19 [1]). When passing a sine wave through the amplifier was succeeded, the frequency response of the wave spring was tested with white noise for comparison to finite-element model predictions. The testbed outputs a repeatable, linear response to the input excitation when used with the wave spring. These measurements reduce the uncertainty associated with the testbed and test procedure to better characterize uncertainty associated with material behavior (SX358 foam samples) in future tests

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermo-acoustoelastic effect of Rayleigh wave: Theory and experimental verification

Previous studies showed that the thermally-induced ultrasonic bulk wave velocity change could be used to measure acoustoelastic coefficients and third-order elastic constants of elastic materials. This method is naturally immune from the ambient temperature effect and has improved sensitivity and a simpler test setup than the conventional acoustoelastic test. However, Rayleigh wave is preferred for thick components or structures with only one accessible surface. In this work, the thermo-hyperelastic constitutive equation, along with acoustoelastic theory, is used to derive the expression of the thermo-acoustoelastic coefficient (TAEC) of Rayleigh wave. The numerical relationship between the TAEC of Rayleigh wave and Murnaghan constants ( l, m and n ) are given for common metals. Further, the TAEC expressions for Rayleigh wave and shear wave are similar, and both are dominated by the constant m . The TAEC of Rayleigh wave was measured on an aluminum 6061 specimen using the thermal modulation experiment in a temperature range of 22 ~ 35°C. The measured TAEC shows good agreement with the theoretical calculation. Then the third-order elastic constants were calculated based on TAECs of bulk waves and Rayleigh wave.

42 ENGINEERING↗