Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave model”

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 325 records · Page 18

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.

Geochemistry & Geophysics↗

Coordinated Satellite Observations of the Very Low Frequency Transmission Through the Ionospheric D Layer at Low Latitudes, Using Broadband Radio Emissions from Lightning

Both ray theory and full-wave models of very low frequency transmission through the ionospheric D layer predict that the transmission is greatly suppressed near the geomagnetic equator. We use data from the low-inclination Communication/Navigation Outage Forecast System satellite to test this semiquantitatively, for broadband very low frequency emissions from lightning. Approximate ground-truthing of the incident wavefields in the Earth-ionosphere waveguide is provided by the World Wide Lightning Location Network. Observations of the wavefields at the satellite are provided by the Vector Electric Field Instrument aboard the satellite. The data set comprises whistler observations with the satellite at magnetic latitudes<26deg. Thus, our conclusions, too, must be limited to the near-equatorial region and are not necessarily predictive of midlatitude whistler properties. We find that in most broadband recordings of radio waves at the satellite, very few of the lightning strokes result in a detectable radio pulse at the satellite. However, in a minority of the recordings, there is enhanced transmission of very low frequency lightning emissions through the D layer, at a level exceeding model predictions by at least an order of magnitude. We show that kilometric-scale D-layer irregularities may be implicated in the enhanced transmission. This observation of sporadic enhancements at low magnetic latitude, made with broadband lightning emissions, is consistent with an earlier review of D-layer transmission for transmission from powerful man-made radio beacons.

Jacobson, Abram R.↗

Comparison of techniques for approximating ocean bottom topography in a wave-refraction computer model

A study of the effects of using different methods for approximating bottom topography in a wave-refraction computer model was conducted. Approximation techniques involving quadratic least squares, cubic least squares, and constrained bicubic polynomial interpolation were compared for computed wave patterns and parameters in the region of Saco Bay, Maine. Although substantial local differences can be attributed to use of the different approximation techniques, results indicated that overall computed wave patterns and parameter distributions were quite similar.

Poole, L. R.↗

Modeling the Complete Planetary Subsurface Radio Remote Sensing Problem

Ground penetrating radar (GPR) methods are in the forefront of the search for subsurface water on Mars largely because of their platform versatility. In theory these instruments can operate from orbit, on the ground, and any altitude in between. Recent review papers and white papers make the strong case for using GPR methods in subsurface water exploration. But important questions have been raised recently about effect of unknown parameters on the ability to successfully obtain planetary GPR measurements. The almost completely unknown lower ionosphere on Mars may cause serious signal absorption problems for GPR instruments on orbiting platforms, and lossy upper layers of the ground may reduce the returned signal for any platform. What is clearly needed is a minimum approximations, full-wave model of the complete GPR problem, including ionospheric dispersion and absorption, surface transmission, and subsurface scattering. We are developing such a model that is as general purpose as possible, allowing arbitrary ionospheric parameters, surface roughness, and subsurface inhomogeneities. We present some of the details of this model, and highlight some of its capabilities with numerical examples. Additional information is contained in the original extended abstract.

Cummer, S. A.↗

Optical system storage design with diffractive optical elements

Optical data storage systems are gaining widespread acceptance due to their high areal density and the ability to remove the high capacity hard disk from the system. In magneto-optical read-write systems, a small rotation of the polarization state in the return signal from the MO media is the signal which must be sensed. A typical arrangement used for detecting these signals and correcting for errors in tracking and focusing on the disk is illustrated. The components required to achieve these functions are listed. The assembly and alignment of this complex system has a direct impact on cost, and also affects the size, weight, and corresponding data access rates. As a result, integrating these optical components and improving packaging techniques is an active area of research and development. Most designs of binary optic elements have been concerned with optimizing grating efficiency. However, rigorous coupled wave models for vector field diffraction from grating surfaces can be extended to determine the phase and polarization state of the diffracted field, and the design of polarization components. A typical grating geometry and the phase and polarization angles associated with the incident and diffracted fields are shown. In our current stage of work, we are examining system configurations which cascade several polarization functions on a single substrate. In this design, the beam returning from the MO disk illuminates a cascaded grating element which first couples light into the substrate, then introduces a quarter wave retardation, then a polarization rotation, and finally separates s- and p-polarized fields through a polarization beam splitter. The input coupler and polarization beam splitter are formed in volume gratings, and the two intermediate elements are zero-order elements.

Kostuk, Raymond K.↗

Spin Waves and Magnetic Exchange Hamiltonian in CrSBr

Abstract CrSBr is an air‐stable two‐dimensional (2D) van der Waals semiconducting magnet with great technological promise, but its atomic‐scale magnetic interactions—crucial information for high‐frequency switching—are poorly understood. An experimental study is presented to determine the CrSBr magnetic exchange Hamiltonian and bulk magnon spectrum. The A ‐type antiferromagnetic order using single crystal neutron diffraction is confirmed here. The magnon dispersions are also measured using inelastic neutron scattering and rigorously fit the excitation modes to a spin wave model. The magnon spectrum is well described by an intra‐plane ferromagnetic Heisenberg exchange model with seven nearest in‐plane exchanges. This fitted exchange Hamiltonian enables theoretical predictions of CrSBr behavior: as one example, the fitted Hamiltonian is used to predict the presence of chiral magnon edge modes with a spin‐orbit enhanced CrSBr heterostructure.

2D materials↗

Nonlinear Dynamics of Cantilever-Sample Interactions in Atomic Force Microscopy

The interaction of the cantilever tip of an atomic force microscope (AFM) with the sample surface is obtained by treating the cantilever and sample as independent systems coupled by a nonlinear force acting between the cantilever tip and a volume element of the sample surface. The volume element is subjected to a restoring force from the remainder of the sample that provides dynamical equilibrium for the combined systems. The model accounts for the positions on the cantilever of the cantilever tip, laser probe, and excitation force (if any) via a basis set of set of orthogonal functions that may be generalized to account for arbitrary cantilever shapes. The basis set is extended to include nonlinear cantilever modes. The model leads to a pair of coupled nonlinear differential equations that are solved analytically using a matrix iteration procedure. The effects of oscillatory excitation forces applied either to the cantilever or to the sample surface (or to both) are obtained from the solution set and applied to the to the assessment of phase and amplitude signals generated by various acoustic-atomic force microscope (A-AFM) modalities. The influence of bistable cantilever modes of on AFM signal generation is discussed. The effects on the cantilever-sample surface dynamics of subsurface features embedded in the sample that are perturbed by surface-generated oscillatory excitation forces and carried to the cantilever via wave propagation are accounted by the Bolef-Miller propagating wave model. Expressions pertaining to signal generation and image contrast in A-AFM are obtained and applied to amplitude modulation (intermittent contact) atomic force microscopy and resonant difference-frequency atomic force ultrasonic microscopy (RDF-AFUM). The influence of phase accumulation in A-AFM on image contrast is discussed, as is the effect of hard contact and maximum nonlinearity regimes of A-AFM operation.

Cantrell, John H.↗

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling Greenland Outlet Glaciers

The objective of this project was to develop simple yet realistic models of Greenland outlet glaciers to better understand ongoing changes and to identify possible causes for these changes. Several approaches can be taken to evaluate the interaction between climate forcing and ice dynamics, and the consequent ice-sheet response, which may involve changes in flow style. To evaluate the icesheet response to mass-balance forcing, Van der Veen (Journal of Geophysical Research, in press) makes the assumption that this response can be considered a perturbation on the reference state and may be evaluated separately from how this reference state evolves over time. Mass-balance forcing has an immediate effect on the ice sheet. Initially, the rate of thickness change as compared to the reference state equals the perturbation in snowfall or ablation. If the forcing persists, the ice sheet responds dynamically, adjusting the rate at which ice is evacuated from the interior to the margins, to achieve a new equilibrium. For large ice sheets, this dynamic adjustment may last for thousands of years, with the magnitude of change decreasing steadily over time as a new equilibrium is approached. This response can be described using kinematic wave theory. This theory, modified to pertain to Greenland drainage basins, was used to evaluate possible ice-sheet responses to perturbations in surface mass balance. The reference state is defined based on measurements along the central flowline of Petermann Glacier in north-west Greenland, and perturbations on this state considered. The advantage of this approach is that the particulars of the dynamical flow regime need not be explicitly known but are incorporated through the parameterization of the reference ice flux or longitudinal velocity profile. The results of the kinematic wave model indicate that significant rates of thickness change can occur immediately after the prescribed change in surface mass balance but adjustments in flow rapidly diminish these rates to a few cm/yr at most. The time scale for adjustment is of the order of a thousand years or so.

vanderVeen, Cornelis↗

Triaxiality and the nature of low-energy excitations in Ge 76

The deformation properties of the low-lying states in 76 Ge have been investigated following a safe energy Coulomb excitation measurement with the GRETINA tracking array and CHICO2 heavy-ion counter at the ATLAS accelerator facility at Argonne National Laboratory. A comprehensive set of transition and static E2 matrix elements were extracted from the measured differential Coulomb cross-sections, and compared with results of configuration interaction shell-model calculations and computations carried out within the framework of the generalized triaxial rotor model. The remarkable agreement between the calculated and experimental data supports a near-maximum triaxial deformation for the ground state of 76 Ge. Additionally, the degree of softness of the asymmetry in 76 Ge and 76 Se was investigated using rotational invariants generated from configuration interaction shell-model wave functions computed with the jj44b and JUN45 effective interactions. The resulting invariants are shown to be consistent with a stiff triaxial deformation in 76 Ge and a predominantly soft triaxial potential for 76 Se, in agreement with the conclusions of recent works by this collaboration.

59 ≤ A ≤ 89↗

Quantum magic and multipartite entanglement in the structure of nuclei

Motivated by the Gottesman-Knill theorem, we present a detailed study of the quantum complexity of p -shell and s d -shell nuclei. Valence-space nuclear shell-model wave functions generated by the code are mapped to qubit registers using the Jordan-Wigner mapping (12 qubits for the p shell and 24 qubits for the s d shell), from which measures of the many-body entanglement ( n -tangles) and magic (nonstabilizerness) are determined. While exact evaluations of these measures are possible for nuclei with a modest number of active nucleons, Monte Carlo simulations are required for the more complex nuclei. The broadly applicable Pauli-string I ̂ Z ̂ exact (PSIZe) Markov chain Monte Carlo (MCMC) technique is introduced to accelerate the evaluation of measures of magic in deformed nuclei (with hierarchical wave functions), by factors of ≈ 8 for some nuclei. Significant multinucleon entanglement is found in the s d shell, dominated by proton-neutron configurations, along with significant measures of magic. This is evident not only for the deformed states, but also for nuclei on the path to instability via regions of shape coexistence and level inversion. These results indicate that quantum-computing resources will accelerate precision simulations of such nuclei and beyond. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spin dynamics in natural multiferroic pyroxene NaFeSi 2 ⁢O 6

Spin dynamics in the natural mineral aegirine, NaFeSi 2⁢ O 6 , a member of the pyroxene family, was studied by elastic and inelastic neutron scattering. Magnetization and specific-heat measurements as well as single-crystal neutron diffraction maps, taken in the temperature range 2–20 K, confirm two successive magnetic transitions at 8.8 and 5.8 K, consistent with previous studies. The observed spin-wave excitations emerge from the incommensurate magnetic Bragg peaks corresponding to the propagation vector 𝑘 ICM = (0,0.77,0), and extend up to energies of about 1.5 meV. In the low-temperature helical phase, the spin dynamics of the Fe3+ ions is well described by a simple linear spin-wave model. The observed excitations can be modeled using a spin Hamiltonian that includes three primary exchange interactions—intrachain coupling 𝐽 = 0.142⁢(2) ⁢meV, and interchain couplings 𝐽 1 = 0.083⁢(1) ⁢meV and 𝐽 2 = 0.186⁢(1) meV—and an easy-plane anisotropy 𝐷 = 0.020⁢(6)⁢ meV. Furthermore, our results show that no single exchange interaction dominates the spin dynamics. The similar strengths of the intrachain and interchain couplings point to the fact that the magnetic interactions in aegirine are three dimensional rather than confined along one direction. As a result, the system cannot be considered quasi one dimensional, as previously suggested, and calls for further investigations.

Magnetic coupling↗

Reconstructing the reflectivity of liquid surfaces from grazing incidence X-ray off-specular scattering data

The capillary wave model of a liquid surface predicts both the X-ray specular reflection and the diffuse scattering around it. A quantitative method is presented to obtain the X-ray reflectivity (XRR) from a liquid surface through the diffuse scattering data around the specular reflection measured using a grazing incidence X-ray off-specular scattering (GIXOS) geometry at a fixed horizontal offset angle with respect to the plane of incidence. With this approach the entire Q z -dependent reflectivity profile can be obtained at a single, fixed incident angle. This permits a much faster acquisition of the profile than with conventional reflectometry, where the incident angle must be scanned point by point to obtain a Q z -dependent profile. The XRR derived from the GIXOS-measured diffuse scattering, referred to in this paper as pseudo-reflectivity, provides a larger Q z range compared with the reflectivity measured by conventional reflectometry. Transforming the GIXOS-measured diffuse scattering profile to pseudo-XRR opens up the GIXOS method to widely available specular XRR analysis software tools. Here the GIXOS-derived pseudo-XRR is compared with the XRR measured by specular reflectometry from two simple vapor–liquid interfaces at different surface tension, and from a hexadecyltrimethylammonium bromide monolayer on a water surface. For the simple liquids, excellent agreement (beyond 11 orders of magnitude in signal) is found between the two methods, supporting the approach of using GIXOS-measured diffuse scattering to derive reflectivities. Pseudo-XRR obtained at different horizontal offset angles with respect to the plane of incidence yields indistinguishable results, and this supports the robustness of the GIXOS-XRR approach. The pseudo-XRR method can be extended to soft thin films on a liquid surface, and criteria are established for the applicability of the approach.

36 MATERIALS SCIENCE↗

Probing the chiral magnetic wave with charge-dependent flow measurements in Pb-Pb collisions at the LHC

The Chiral Magnetic Wave (CMW) phenomenon is essential to provide insights into the strong interaction in QCD, the properties of the quark-gluon plasma, and the topological characteristics of the early universe, offering a deeper understanding of fundamental physics in high-energy collisions. Measurements of the charge-dependent anisotropic flow coefficients are studied in Pb-Pb collisions at center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV to probe the CMW. In particular, the slope of the normalized difference in elliptic (v 2 ) and triangular (v 3 ) flow coefficients of positively and negatively charged particles as a function of their event-wise normalized number difference, is reported for inclusive and identified particles. The slope ${r}_3^{\textrm{Norm}}$ is found to be larger than zero and to have a magnitude similar to ${r}_2^{\textrm{Norm}}$, thus pointing to a large background contribution for these measurements. Furthermore, ${r}_2^{\textrm{Norm}}$ can be described by a blast wave model calculation that incorporates local charge conservation. In addition, using the event shape engineering technique yields a fraction of CMW (f CMW ) contribution to this measurement which is compatible with zero. This measurement provides the very first upper limit for f CMW , and in the 10–60% centrality interval it is found to be 26% (38%) at 95% (99.7%) confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for the Chiral Magnetic Effect with charge-dependent azimuthal correlations in Xe–Xe collisions at s NN = 5.44 TeV

Charge-dependent two- and three-particle correlations measured in Xe–Xe collisions at $\sqrt{s_{NN}}$ = TeV are presented. Results are obtained for charged particles in the pseudorapidity range |η| and transverse momentum interval 02. ≤ p T < GeV/c for different collision centralities. The three-particle correlator γ αβ ≡ $\langle$cos(φ α + φ β - 2Ψ a )$\rangle$, calculated for different combinations of charge sign α and β, is expected to be sensitive to the presence of the Chiral Magnetic Effect (CME). Its magnitude is similar to the one observed in Pb–Pb collisions in contrast to a smaller CME signal in Xe–Xe collisions than in Pb–Pb collisions predicted by Monte Carlo (MC) calculations including a magnetic field induced by the spectator protons. These observations point to a large non-CME contribution to the correlator. Furthermore, the charge dependence of γ αβ can be described by a blast wave model calculation that incorporates background effects and by the Anomalous Viscous Fluid Dynamics model with values of the CME signal consistent with zero. The Xe–Xe and Pb–Pb results are combined with the expected CME signal dependence on the system size from the MC calculations including a magnetic field to obtain the fraction of CME contribution in γ αβ , $f$ CME . The CME fraction is compatible with zero for the 30% most central events in both systems and then becomes positive. This yields an upper limit of 2% (3%) and 25% (32%) at 95% (99.7%) confidence level for the CME signal contribution to γ αβ in the 0–70% Xe–Xe and Pb–Pb collisions, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of 3 Λ H production in Pb–Pb collisions at √ s NN = 5.02 TeV

The first measurement of 3 Λ H and $^{3}_{\overline{Λ}}$$\overline{H}$ differential production with respect to transverse momentum and centrality in Pb–Pb collisions at √ s NN = 5.02 TeV is presented. The 3 Λ H has been reconstructed via its two-charged-body decay channel, i.e., 3 Λ H → 3 He → π - . A Blast-Wave model fit of the p T -differential spectra of all nuclear species measured by the ALICE collaboration suggests that the kinetic freeze-out surface is consistent with that of other nuclei. The ratio between the integrated yields of 3 Λ H and 3 He is compared to predictions from the statistical hadronisation model and the coalescence model, with the latter being favoured by the presented measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ponderomotive force driven density modifications parallel to B on the LAPD

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Local measurement of bulk thermal diffusivity using photothermal radiometry

Here, an experimental methodology using photothermal radiometry is developed for the accurate measurement of bulk thermal diffusivity of nuclear fuels and materials irradiated to high doses. Under these conditions, nuclear fuels, such as uranium oxide, and moderator materials, such as graphite, become friable, which requires characterization techniques that can accommodate irregularly shaped fragments. Photothermal radiometry, a good candidate for this application, involves locally heating a sample by using a laser and measuring the temperature field by monitoring blackbody radiation. The interaction volume for this study, less than a millimeter, is carefully chosen to sample a statistically significant number of large-scale structural features, such as pores and gas filled bubbles, and is small enough that the sample fragments can be treated as a thermal half-space. The thermal diffusivity standards considered in this study cover a range of thermal diffusivities representative of both fresh and spent nuclear fuels. We also consider a sample having a porous microstructure representative of large-scale structures found in materials irradiated to high doses. Our measurement methodology circumvents complex thermal wave models that address optical diffraction, nonlinear transfer function associated with blackbody radiation, and finite sample size effects. Consequently, the large measurement uncertainty associated with modeling these effects can be avoided. While the emphasis here is on nuclear fuels and materials, this measurement approach is well suited to measure thermal transport in a variety of technologically important materials associated with advanced synthesis techniques. Examples range from small, exotic single crystals grown using hydrothermal growth techniques to additively manufactured components having complex geometries.

36 MATERIALS SCIENCE↗