Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave modeling”

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 1,369 records · Page 76

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

Parameter-Varying Hydrodynamic Model of a Single Vane in a Variable-Geometry Oscillating Surge Wave Energy Converter: Preprint

This paper presents a study on the time-varying hydrodynamic modeling of an individual flap of a variable geometry oscillating surge wave energy converter (VGOSWEC). The WEC design incorporates controlled surfaces that can modify their orientation relative to the wave motion, reducing hydrodynamic pressure and related loads. This research focuses on characterizing the behavior of the oscillating WEC using a simplified model and three methods for achieving a continuous time-variant model: discrete hydrodynamic parameters, interpolation of hydrodynamic parameters, and a fitting function. The results of this study contribute to the understanding of time-varying hydrodynamic effects in variable geometry oscillating WECs. The findings provide insights into the potential for reducing structural loads and improving the overall performance of such devices. Further research and development in this area could lead to advancements in WEC technologies, enabling their integration into the competitive energy market.

HYDRO ENERGY,TIDAL AND WAVE POWER↗

Crustal Structure of the Mesopotamian Plain, East of Iraq

The crustal structure of Iraq was investigated through analyzing teleseismic data from 12 new seismic stations. Three seismic stations are located within the Zagros Fold-Thrust Belt, and nine are within the Mesopotamian Plain. Joint inversion of P wave receiver function and Rayleigh wave dispersion data were employed to resolve S wave velocity structure models beneath each station. Combining these models with available Moho depths within and around the study area reveals that the Moho depth is smoothly increasing from the Arabian platform toward the Zagros Mountains. Further, an exceptional deeper root is observed near the eastern edge of the Mesopotamian Plain where sedimentary pile is the thickest. This root was interpreted as a structure inherited from the successive Mesozoic rifting of the NE Arabian platform enhanced by sedimentary loading and progressive subsidence. A low-velocity uppermost mantle was observed beneath the Arabian Foredeep attesting the existence of a low-strength lithospheric mantle beneath the region southwest from the Zagros deformation front. The weak rheology of the uppermost mantle may have allowed local sinking of the crust and deepening of the Moho boundary due to vertical loading and the Late Miocene lateral contraction.

58 GEOSCIENCES↗

Numerical simulation of magnetohydrodynamic shock propagation in the corona

Recent developments in the field of numerical simulation models for the study of shock wave propagation in the corona are presented. These models are based on gasdynamic (GD) and ideal (that is, dissipationless, except at shocks) magnetohydrodynamic (MHD) theories. The characteristics and physical interpretations of the results derived from these models are discussed in some detail. The most significant physical results obtained to date are provided by the two-dimensional non-planar, time-dependent, MHD numerical simulation model. In this model, the non-linear interaction among the three essential MHD waves, i.e., fast-, slow-, and Alfven waves are demonstrated. Finally, the physical relevance of these numerical simulation models in relation to observed solar activity is presented.

Wu, S. T.↗

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The effects of profiles on supersonic jet noise

The effect of velocity profiles on supersonic jet noise are studied by using stability calculations made for a shock-free coannular jet, with both the inner and outer flows supersonic. The Mach wave emission process is modeled as the noise generated by the large scale turbulent structures or the instability waves in the mixing region. Both the vortex-sheet and the realistic finite thickness shear layer models are considered. The stability calculations were performed for both inverted and normal velocity profiles. Comparisons are made with the results for an equivalent single jet, based on equal thrust, mass flow rate and exit area to that of the coannular jet. The advantages and disadvantages of these velocity profiles as far as noise radiation is concerned are discussed. It is shown that the Rayleigh's model prediction of the merits and demerits of different velocity profiles are in good agreement with the experimental data.

Tiwari, S. N.↗

Applicability of the single equivalent point dipole model to represent a spatially distributed bio-electrical source

Although the single equivalent point dipole model has been used to represent well-localised bio-electrical sources, in realistic situations the source is distributed. Consequently, position estimates of point dipoles determined by inverse algorithms suffer from systematic error due to the non-exact applicability of the inverse model. In realistic situations, this systematic error cannot be avoided, a limitation that is independent of the complexity of the torso model used. This study quantitatively investigates the intrinsic limitations in the assignment of a location to the equivalent dipole due to distributed electrical source. To simulate arrhythmic activity in the heart, a model of a wave of depolarisation spreading from a focal source over the surface of a spherical shell is used. The activity is represented by a sequence of concentric belt sources (obtained by slicing the shell with a sequence of parallel plane pairs), with constant dipole moment per unit length (circumferentially) directed parallel to the propagation direction. The distributed source is represented by N dipoles at equal arc lengths along the belt. The sum of the dipole potentials is calculated at predefined electrode locations. The inverse problem involves finding a single equivalent point dipole that best reproduces the electrode potentials due to the distributed source. The inverse problem is implemented by minimising the chi2 per degree of freedom. It is found that the trajectory traced by the equivalent dipole is sensitive to the location of the spherical shell relative to the fixed electrodes. It is shown that this trajectory does not coincide with the sequence of geometrical centres of the consecutive belt sources. For distributed sources within a bounded spherical medium, displaced from the sphere's centre by 40% of the sphere's radius, it is found that the error in the equivalent dipole location varies from 3 to 20% for sources with size between 5 and 50% of the sphere's radius. Finally, a method is devised to obtain the size of the distributed source during the cardiac cycle.

NASA Discipline Cardiopulmonary↗

Revisiting Theoretical Limits for One-Degree-of-Freedom Wave Energy Converters: Preprint

This work revisits the theoretical limits of one-degree of freedom wave energy converters (WECs). This paper considers the floating sphere used in the Ocean Energy Systems (OES) Task 10 Wave Energy Converter modeling and verification effort for analysis. Analytical equations are derived to determine bounds on the motion amplitude, time-averaged power, and power take-off (PTO) force. A unique result was found that shows the time-averaged power absorbed by a WEC can be defined solely by the inertial properties and the radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force amplitude was derived which has provided upper and lower bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression can only provide a lower bound as there is theoretically no upper bound. These bounds are used to compare the performance of a floating sphere if it were to extract energy using surge or heave motion. The analysis shows that because of the differences in hydrodynamic coefficients of each oscillating mode there will be different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on power absorption while also utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate the time-averaged power with a nonideal PTO is modified by the mechanical-to-electrical efficiency and the ratio of the PTO spring and damping coefficients. These were separated in the expression which allows for limits to be set on the possible values of PTO coefficients to ensure a net flow of power to the grid.

50 EE - Wind and Water Power Program - Water (EE-4↗

Galactic civilizations: Population dynamics and interstellar diffusion

The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the environment. Dispersal is fundamentally a diffusion process; a density-dependent diffusivity describes interstellar emigration. Two models are considered: the first describing zero population growth (ZPG), and the second which also includes local growth and saturation of a planetary population, and for which an asymptotic traveling wave solution is found.

Newman, W. I.↗

A semikinetic model for early stage plasmasphere refilling. II - Effects of wave-particle interactions

The paper treats the early stages of plasmasphere refilling along an initially depleted L = 4 magnetic flux tube through a semikinetic model, with special attention given to the effects of wave-particle interactions in which stochastic diffusion of ions in perpendicular velocity due to equatorially concentrated electromagnetic ion cyclotron waves plays a central role. Characteristic individual ion trajectories are examined, as well as the devolution of bulk parameters and ion distribution functions when equal 'polar wind' streams are injected at the northern and southern ionospheres. In the ion trajectories, relatively modest and realistic perpendicular electric field power levels are found to lead to decreased mirror latitudes, substantial acceleration, and equatorial entrainment of these ions. A substantial equatorial density depletion is also found, in agreement with the results of Olsen (1992).

Lin, J.↗

Mach wave emission from a high temperature supersonic jet

The paper considers the compressible Rayleigh equation as a model for the Mach wave emission mechanism associated with high temperature supersonic jets. Solutions to the compressible Rayleigh equation reveal the existence of several families of supersonically convecting instability waves. These waves directly radiate noise to the jet far field. The predicted noise characteristics are compared to previously acquired experimental data for an axisymmetric Mach 2 fully pressure balanced jet operating over a range of jet operating total temperatures from ambient to 1370 K. The results of this comparison show that the first order supersonic instability wave and the Kelvin-Helmholtz first, second, and third order modes have directional radiation characteristics that are in agreement with observed data. The assumption of equal initial amplitudes for all of the waves leads to the conclusion that the flapping mode of instability dominates the noise radiation process of supersonic jets. At a jet temperature of 1370 K, supersonic instability waves are predicted to dominate the noise radiated at high frequency at narrow angles to the jet axis.

Seiner, J. M.↗

Lightweight, Flexible Electromagnetic Shielding Composite Films Reinforced with Recycled Carbon Fibers and Carbon Nanofillers

Lightweight, flexible polymer composites are widely adopted in modern electronic devices and systems for high-efficiency electromagnetic interference (EMI) shielding. Reinforcing polymer composites with conductive fillers offers a promising alternative to conventional metal-based shielding material thanks to their low density, tunable properties, and flexibility. Herein, lightweight, flexible ultra-high molecular weight polyethylene composite films reinforced with recycled carbon fibers (rCFs) and carbon-based nanofillers, including graphene nanoplatelets (GNPs) and carbon nanotubes, are reported for EMI shielding applications. The incorporation of these nanofillers significantly reduces the required content of rCFs and improves processability while maintaining comparable shielding effectiveness. The addition of these nanofillers with rCFs enhances the EMI shielding effectiveness by up to 15 dB. Incorporating 1 wt% GNPs can replace 5 wt% rCFs and achieve comparable EMI shielding performance, while 5 wt% of either nanofillers can substitute for 10 wt% rCFs. Numerical modeling of electromagnetic wave transmission reveals that increasing nanofiller concentrations enhances both reflection and absorption losses, with absorption consistently dominating across all levels. Furthermore, this study not only provides insight into the synergistic contributions of rCFs and carbon nanofillers to shielding effectiveness but also paves the way for the design of sustainable, lightweight EMI shielding composite films for applications in electric vehicles, medical equipment, and portable electronics.

carbon nanofiller↗

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Effects of Field Line Curvature (FLC) Scattering on Ring Current Dynamics and Isotropic Boundary

Abstract In the ring current dynamics, various loss mechanisms contribute to the ring current decay, including losses to the upper atmosphere through particle precipitation. This study implements the field line curvature (FLC) scattering mechanism in a kinetic ring current model and investigates its role in precipitating ions into the ionosphere during the 17 March 2013 storm. Simulation results indicate that (1) the FLC scattering process exerts on energetic ions on the nightside where the magnetospheric configuration is more stretching. It is more effective on heavy ions (e.g., O + ). These ion losses thereafter lead to a faster recovery of the ring current. (2) The FLC‐associated ion precipitation mainly occurs in the outer region (L > 5 for protons and L > 4.5 for oxygen ions) on the nightside. The O + precipitation takes places in a wider region than protons although its intensity is much lower. Comparisons with POES observations suggest that more proton precipitation is needed in the inner region. This is probably caused by the less stretched configuration in the simulation that prevents more precipitation. It may also imply that other loss process is required in the model such as wave‐particle interactions. (3) The storm time precipitating proton flux of tens of keV due to the FLC scattering sometimes becomes comparable to that of electrons on the nightside, although electrons usually dominate the ionospheric energy deposition from the midnight eastward toward the dayside. (4) The FLC scattering process seems to be capable of explaining the formation of isotropic boundary in the ionosphere during the investigated event.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of the Langdon effect on crossed-beam energy transfer

The prediction that laser plasma heating distorts the electron distribution function away from Maxwellian and towards a super-Gaussian distribution dates back four decades. In conditions relevant to inertial confinement fusion, yet, no direct evidence of this so-called “Langdon effect” has previously been observed. Here, measurements of the spatially and temporally resolved Thomson scattering spectrum indicate the presence of super-Gaussian electron distribution functions that are consistent with existing theory. In such plasmas, ion acoustic wave frequencies increase monotonically with the super-Gaussian exponent. Our findings show that the measured power transfer between crossed laser beams mediated by ion acoustic waves requires a model that accounts for the non-Maxwellian electron distribution function, whereas the standard Maxwellian calculations overpredict power transfer over a wide region of parameter space. Including this effect is expected to improve the predictive capability of crossed-beam energy transfer modeling at the National Ignition Facility and may restore a larger operable design space for inertial confinement fusion experiments. Furthermore, this is expected to motivate further inquiry in other areas impacted by non-Maxwellian electron distribution functions, such as laser absorption, heat transport, and x-ray spectroscopy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Classical discrete time crystals

The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: the discrete time crystal (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator–prey models to Faraday waves and a.c.-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially long timescales. Power-law correlations are present along a first-order line, which terminates at a critical point. We analyse the transition by mapping it to the locked-to-sliding transition of a d.c.-driven charge density wave. To conclude, building upon results from the field of probabilistic cellular automata, we conjecture the existence of classical time crystals with true long-range order, where time-translation symmetry is broken out to infinite times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗