Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse waves”

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,549 records · Page 86

Quantum Mechanical Study of Nanoscale MOSFET

The steady state characteristics of MOSFETS that are of practical Interest are the drive current, off-current, dope of drain current versus drain voltage, and threshold voltage. In this section, we show that quantum mechanical simulations yield significantly different results from drift-diffusion based methods. These differences arise because of the following quantum mechanical features: (I) polysilicon gate depletion in a manner opposite to the classical case (II) dependence of the resonant levels in the channel on the gate voltage, (III) tunneling of charge across the gate oxide and from source to drain, (IV) quasi-ballistic flow of electrons. Conclusions dI/dV versus V does not increase in a manner commensurate with the increase in number of subbands. - The increase in dI/dV with bias is much smaller then the increase in the number of subbands - a consequence of bragg reflection. Our calculations show an increase in transmission with length of contact, as seen in experiments. It is desirable for molecular electronics applications to have a small contact area, yet large coupling. In this case, the circumferential dependence of the nanotube wave function dictates: - Transmission in armchair tubes saturates around unity - Transmission in zigzag tubes saturates at two.

Svizhenko, Alexei↗

Aerosol Optical Tweezers Elucidate the Chemistry, Acidity, Phase Separations, and Morphology of Atmospheric Microdroplets

Aerosol particles represent unique chemical environments because of their high surface area-to-volume ratio that promotes the effects of interfacial chemistry in confined environments. Properties such as viscosity, diffusivity, water content, pH, and morphology—following liquid–liquid phase separation—can strongly alter how a particle interacts with condensable vapors and reactive trace gases, thus modifying its continual evolution and environmental effects. Our understanding of this chemical evolution of atmospheric particulate matter and its environmental impacts is largely limited by our ability to directly observe how these critical particle properties respond to the addition or reactive uptake of new chemical components. Aerosol optical tweezers (AOT) stably trap particles in focused laser beams, providing positional control and the retrieval of many of these critical properties required to understand and predict the chemistry of aerosolized microdroplets. The analytical power of the AOT stems from the retrieval of the cavity-enhanced Raman spectrum induced by the trapping laser. Analysis of the whispering gallery modes (WGMs) that resonate as a standing wave around the droplet’s interface, provide high accuracy measurements of the droplet’s size, refractive index (and thus a measurement of composition), and can distinguish between core–shell, partially engulfed, and homogeneous morphologies. In this study, we have advanced the ability to determine the properties of the core and shell phases in biphasic droplets, including obtaining high-accuracy pH measurements. These capabilities were applied to perform AOT physical chemistry experiments on authentic secondary organic aerosol (SOA) produced directly in the AOT chamber by ozonolysis of terpene vapors. The propensity of the SOA to phase separate as a shell from a wide range of nonpolar to polar core phases was observed, along with the discovery of a stable emulsified state of SOA particles in an aqueous salt droplet. Micron-thick SOA shells did not impede the gain or loss of water or squalane from the core to the surrounding air, indicating no significant diffusional limitations to condensational growth or partitioning even under dry conditions. These experiments formed the foundation of a new framework that predicts how the phase-separated morphology of complex aerosols containing organic carbon evolves during continual atmospheric oxidation processes. Increases in oxidation state will quickly drive conversion from a partially engulfed to core–shell morphology that has dramatically different chemical reactivity since the core phase is completely concealed by the shell. The recent advances in the experimental capabilities of the AOT technique such as presented here enable novel experimental methodologies that provide insights into the chemistry and multidimensional properties of aerosol microdroplets, and how these coevolve and respond to continual chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magneto-thermal reconnection processes, related mode momentum and formation of high energy particle populations

In the context of a two-fluid theory of magnetic reconnection, when the longitudinal electron thermal conductivity is relatively large, the perturbed electron temperature tends to become singular in the presence of a reconnected field component and an electron temperature gradient. A finite transverse thermal diffusivity removes this singularity while a finite ‘inductivity’ can remove the singularity of the relevant plasma displacement. Then (i) a new ‘magneto-thermal’ reconnection producing mode, is found with characteristic widths of the reconnection layer remaining significant even when the macroscopic distances involved are very large; (ii) the mode phase velocities can be both in the direction of the electron diamagnetic velocity as well in the opposite (ion) direction. A numerical solution of the complete set of equations has been carried out with a simplified analytical reformulation of the problem. A sequence of processes is analyzed to point out that high-energy particle populations can be produced as a result of reconnection events. These processes involve mode-particle resonances transferring energy of the reconnecting mode to a superthermal ion population and the excitation of lower hybrid waves that can lead to a significant superthermal electron population. The same modes excited in axisymmetric (e.g. toroidal) confinement configurations can extract angular momentum from the main body of the plasma column and thereby sustain a local ‘spontaneous rotation’ of it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Electrostatic noise bands associated with the electron gyrofrequency and plasma frequency in the outer magnetosphere

Naturally occurring noise bands near the electron plasma frequency are frequently detected by the University of Iowa plasma wave experiment on the Imp 6 satellite in the region from just inside the plasmapause to radial distances of about 10 earth radii in the outer magnetosphere. The electric field strength of these noise bands is usually small with typical broad band electric field strengths of about 2 microvolts per m. A magnetic field has been detected only in a few unusually intense cases, and in these cases the magnetic field energy density is several orders of magnitude smaller than the electric field energy density. The bands are observed at all magnetic latitudes covered by the Imp 6 orbit and appear to be a permanent feature of the outer magnetosphere. They are found at all local times and occur least frequently in the quadrant from 18 to 24 hours. The bands appear to consist of two distinct spectral types called diffuse and narrow band. In both types the center frequency of the noise band is bounded by consecutive harmonics of the electron gyrofrequency, and these noise bands occur most often between harmonics that are near the local electron plasma frequency.

Shaw, R. R.↗

Negligible magnetic losses at low temperatures in liquid phase epitaxy grown Y 3 Fe 5 O 12 films

Yttrium iron garnet (Y 3 Fe 5 O 12 ; YIG) has a unique combination of low magnetic damping, high spin-wave conductivity, and insulating properties that make it a highly attractive material for a variety of applications in the fields of magnetics and spintronics. While the room-temperature magnetization dynamics of YIG have been extensively studied, there are limited reports correlating the low-temperature magnetization dynamics to the material structure or growth method. Here, in this study, we investigate liquid phase epitaxy grown YIG films and their magnetization dynamics at temperatures down to 10 K. We show there is a negligible increase in the ferromagnetic resonance linewidth down to 10 K, which is unique when compared with YIG films grown by other deposition methods. From the broadband ferromagnetic resonance measurements, polarized neutron reflectivity, and scanning transmission electron microscopy, we conclude that these liquid phase epitaxy grown films have negligible rare-earth impurities present, specifically the suppression of Gd diffusion from the Gd 3 Ga 5 O 12 (GGG) substrate into the Y 3 Fe 5 O 12 film, and therefore negligible magnetic losses attributed to the slow-relaxation mechanism. Overall, liquid phase epitaxy YIG films have a YIG/GGG interface that is five times sharper and have ten times lower ferromagnetic resonance linewidths below 50 K than comparable YIG films by other deposition methods. Thus, liquid phase epitaxy grown YIG films are ideal for low-temperature experiments/applications that require low magnetic losses, such as quantum transduction and manipulation via magnon coupling.

36 MATERIALS SCIENCE↗

Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.

08 HYDROGEN↗

Recombination rate analysis in long minority carrier lifetime mid-wave infrared InGaAs/InAsSb superlattices

Gallium is incorporated into the strain-balanced In(Ga)As/InAsSb superlattice system to achieve the same mid-wave infrared cutoff tunability as conventional Ga-free InAs/InAsSb type-II superlattices, but with an additional degree of design freedom to enable optimization of absorption and transport properties. Time-resolved photoluminescence measurements of InGaAs/InAsSb superlattice characterization- and doped device structures are reported from 77 to 300 K and compared to InAs/InAsSb. The low-injection photoluminescence decay yields the minority carrier lifetime, which is analyzed with a recombination rate model, enabling the determination of the temperature-dependent Shockley–Read–Hall, radiative, and Auger recombination lifetimes and extraction of defect energy levels and capture cross section defect concentration products. The Shockley–Read–Hall-limited lifetime of undoped InGaAs/InAsSb is marginally reduced from 2.3 to 1.4 μs due to the inclusion of Ga; however, given that Ga improves the vertical hole mobility by a factor of >10×, a diffusion-limited InGaAs/InAsSb superlattice nBn could expect a lower bound of 2.5× improvement in diffusion length with significant impact on photodetector quantum efficiency and radiation hardness. At temperatures below 120 K, the doped device structures are Shockley–Read–Hall limited at 0.5 μs, which shows promise for detector applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Doping a bad metal: Origin of suppression of the metal-insulator transition in nonstoichiometric VO 2

Rutile ( R ) phase VO 2 is a quintessential example of a strongly correlated bad metal, which undergoes a metal-insulator transition (MIT) concomitant with a structural transition to a V-V dimerized monoclinic ( M 1 ) phase below T MIT ~ 340 K . It has been experimentally illustrated that one can control this transition by doping VO 2 . In particular, doping with oxygen vacancies ( V O ) has been shown to completely suppress this MIT without any structural transition. We explain this suppression by elucidating the influence of oxygen vacancies on the electronic structure of the metallic R phase VO 2 , explicitly treating strong electron-electron correlations using dynamical mean-field theory (DMFT) as well as diffusion Monte Carlo (DMC) flavor of quantum Monte Carlo (QMC) techniques. DMC calculations show a gap closure in the M 1 phase when vacancies are present, suggesting that when vacancies are introduced in the high-temperature rutile phase, the dimerized insulating phase cannot be reached when temperature is lowered. Both DMFT and DMC calculations of nonstoichiometric metallic rutile phase shows that this tendency not to dimerize in the presence of vacancies is because V O 's tend to change the V - 3 d filling away from its nominal half-filled value, with the e g π orbitals competing with the otherwise dominant a 1 g orbital. Loss of this near orbital polarization of the a 1 g orbital is associated with a weakening of electron correlations, especially along the V-V dimerization direction. This removes a charge-density wave (CDW) instability along this direction above a critical doping concentration, which further suppresses the metal-insulator transition. Our research further suggests that the MIT is predominantly driven by a correlation-induced CDW instability along the V-V dimerization direction.

36 MATERIALS SCIENCE↗

Optimization of a Deep Convective Cloud Technique in Evaluating the Long-Term Radiometric Stability of MODIS Reflective Solar Bands

MODIS reflective solar bands are calibrated on-orbit using a solar diffuser and near-monthly lunar observations. To monitor the performance and effectiveness of the on-orbit calibrations, pseudo-invariant targets such as deep convective clouds (DCCs), Libya-4, and Dome-C are used to track the long-term stability of MODIS Level 1B product. However, the current MODIS operational DCC technique (DCCT) simply uses the criteria set for the 0.65- m band. We optimize several critical DCCT parameters including the 11- micrometer IR-band Brightness Temperature (BT11) threshold for DCC identification, DCC core size and uniformity to help locate DCCs at convection centers, data collection time interval, and probability distribution function (PDF) bin increment for each channel. The mode reflectances corresponding to the PDF peaks are utilized as the DCC reflectances. Results show that the BT11 threshold and time interval are most critical for the Short Wave Infrared (SWIR) bands. The Bidirectional Reflectance Distribution Function model is most effective in reducing the DCC anisotropy for the visible channels. The uniformity filters and PDF bin size have minimal impacts on the visible channels and a larger impact on the SWIR bands. The newly optimized DCCT will be used for future evaluation of MODIS on-orbit calibration by MODIS Characterization Support Team.

Mu, Qiaozhen↗

Observations of Offshore Internal Boundary Layers

The growth of the marine internal boundary layer (MIBL, height h i ) with the shore-normal distance x, is a topic of continuing interest because of its applications in coastal dispersion, offshore wind farm siting, coastal air-sea fluxes and in evaporative ducting. Available data on MIBL are only scarce, given its difficult to measure the variability of coastal winds. During Coupled Air-Sea Processes and Electromagnetic Research (CASPER) campaigns an array of instrumentation was deployed to measure offshore spatial variability and its effect on electromagnetic wave propagation. Meteorological sensors (flux towers and remote sensing) were deployed along the coast of Point Mugu, California, on a research vessel and FLoating Instrument Platform (FLIP) provided surface layer and boundary layer observations. In this article, measurements from multiple remote sensing instruments, such as synchronized triple Doppler lidars, small boat operations with tethered lifting system, and radiosondes provide a holistic view of the MIBL growth and its spatial variability from the coast. A convective and stable MIBL observed during two intensive operating period days showed distinct growth characteristics off the coast of Point-Mugu. During stable stratified atmospheric conditions, an MIBL was observed at least up to 47 km from the coast. The growth of the MIBL within the nearshore adjustment zone was influenced by surrounding atmospheric, oceanographic, and topographic conditions. A new parameterization scheme is developed based on advection-diffusion balance equations accounting for upstream turbulence and compared with h i observations from Doppler lidar and profiles from small boat. An evaluation of existing IBL theories was also conducted.

54 ENVIRONMENTAL SCIENCES↗

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗

Power-Law Entanglement and Hilbert Space Fragmentation in Nonreciprocal Quantum Circuits

Quantum circuits utilizing measurement to evolve a quantum wave function offer a new and rich playground to engineer unconventional entanglement dynamics. Here, in this work, we introduce a hybrid, nonreciprocal setup featuring a quantum circuit, whose updates are conditioned on the state of a classical dynamical agent. In our example the circuit is represented by a Majorana quantum chain controlled by a classical N-state Potts chain undergoing pair flips. The local orientation of the classical spins controls whether randomly drawn local measurements on the quantum chain are allowed or not. This imposes a dynamical kinetic constraint on the entanglement growth, described by the transfer matrix of an N-colored loop model. It yields an equivalent description of the circuit by an SU(N)-symmetric Temperley-Lieb Hamiltonian or by a kinetically constrained surface growth model for an N-component height field. For N = 2, we find a diffusive growth of the half-chain entanglement toward a stationary profile S(L) ~ L 1/2 for L sites. For N ≥ 3, the kinetic constraints impose Hilbert space fragmentation, yielding subdiffusive growth toward S(L) ~ L 0.57 . This showcases how the control by a classical dynamical agent can enrich the entanglement dynamics in quantum circuits, paving a route toward novel entanglement dynamics in nonreciprocal hybrid circuit architectures.

1-dimensional spin chains↗

Revealing the Drag Instability in One-fluid Nonideal Magnetohydrodynamic Simulations of a 1D Isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in nonideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the post-shock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Aqua MODIS: 20 Years of On-orbit Calibration and Performance

Since its launch in May 2002, Aqua MODIS has successfully operated for more than 20 years and has continuously generated a wide range of data products that have enabled and supported the remote sensing community and users worldwide for their studies of the Earth’s system by monitoring changes in its key environmental parameters. Although Aqua MODIS, designed with a lifetime requirement of 6 years, is currently operated in its extended mission phase, it continues to make high quality global observations of the Earth’s surface via its 36 spectral bands that cover wavelengths from visible to long-wave infrared. To date, all instrument on-board calibrators (OBC) remain capable of performing their design functions, providing various calibration data sets to help monitor on-orbit changes in sensor responses and performance characteristics. In addition to the OBC, regularly scheduled lunar observations and select Earth-view targets are used extensively to support sensor on-orbit calibration, especially for the calibration of the visible channels (or bands). In this paper, we provide an overview of Aqua MODIS on-orbit calibration activities and methodologies for both reflective solar bands (RSB) and thermal emissive bands (TEB), illustrate its on-orbit performance over the past 20 years using examples derived from OBC measurements, lunar observations, and Earth-view response trends, and describe various calibration improvements made over its entire mission. We focus on key issues identified since launch, such as solar diffuser degradation, electronic crosstalk, and on-orbit changes in sensor response versus scan-angle, along with approaches and strategies developed to mitigate their impact on sensor calibration quality. Also discussed in this paper are some of the key calibration enhancements incorporated recently in the Collection 6.1 and the upcoming Collection 7 Level-1B algorithms.

Polarization↗

Causes of fuel–ablator mix inferred from modeling of monochromatic time-gated radiography of OMEGA cryogenic implosions

Here, we present evidence, in the context of OMEGA cryogenic target implosions, that laser imprint, known to be capable of degrading laser-direct-drive target performance, plays a major role in generating fuel–ablator mix. OMEGA cryogenic target implosions show a performance boundary correlated with acceleration-phase shell stability; for sufficiently low adiabats (where the adiabat is the ratio of the pressure to the Fermi pressure) and high in-flight aspect ratios (IFAR's), the neutron-weighted shell areal density and neutron yield relative to the clean simulated values sharply decline. Direct evidence of Rayleigh–Taylor fuel–ablator mixing was previously obtained using a Si Heα backlighter driven by an ~20-ps short pulse generated by OMEGA EP. The shadow cast by the shell shortly prior to stagnation, as diagnosed using backlit radiographs, shows a softening near the limb, which is evidence of an ablator–fuel mix region for a low-adiabat implosion (α ~ 1.9, IFAR = 14) but not for a moderate adiabat implosion (α ~ 2.5, IFAR = 10). We find good agreement between experimental and synthetic radiographs in simulations that model laser imprint and account for uncertainty in the initial ablator thickness. Furthermore, we further explore the role of other mechanisms such as classical instability growth at the fuel–ablator interface, species concentration diffusion, and long-wavelength drive and target asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Effects of Ground and Space Processing on the Properties of Organic, Polymeric, and Colloidal Materials

In recent years, a great deal of interest has been directed toward the use of organic materials in the development of high-efficiency optoelectronic and phototonic devices. There is a myriad of possibilities among organic materials which allow flexibility in the design of unique structures with a variety of functional groups. The use of nonlinear optical (NLO) organic materials as thin film wave-guides allows full exploitation of their desirable qualifies by permitting long interaction lengths and large susceptibilities allowing modest power input. There are several methods in use to prepare thin films such as Langmuir-Blodgett (LB) and self-assembly techniques, vapor deposition, growth from sheared solution or melt, and melt growth between glass plates. Organic-based materials have many features that make them desirable for use in optical devices, such as high second-and third-order nonlinearity, flexibility of molecular design, and damage resistance to optical radiation. However, their use in devices has been hindered by processing difficulties for crystals and thin films. We discuss the potential role of microgravity processing of a few organic and polymeric materials. It is of interest to note how materials with second-and third-order NLO behavior may be improved in a diffusion-limited environment and ways in which convection may be detrimental to these materials. We focus our discussion on third-order materials for all-optical switching, and second-order materials for frequency conversion and electrooptics. The goal of minimizing optical loss obviously depends on processing methods. For solution-based processes, such as solution crystal growth and solution photopolymerization, it is well known that thermal and solutal density gradients can initiate buoyancy-driven convection. Resultant fluid flows can affect transport of material to and from growth interfaces and become manifest in the morphology and homogeneity of the growing film or crystal. Likewise, buoyancy-driven convection can hinder production of defect-free, high-quality crystals or films during crystal and film growth by vapor deposition.

Frazier, Donald O.↗

A residual-based shock capturing scheme for the continuous/discontinuous spectral element solution of the 2D shallow water equations

The high-order numerical solution of the non-linear shallow water equations is susceptible to Gibbs oscillations in the proximity of strong gradients. In this paper, we tackle this issue by presenting a shock capturing model based on the numerical residual of the solution. Via numerical tests, we demonstrate that the model removes the spurious oscillations in the proximity of strong wave fronts while preserving their strength. Furthermore, for coarse grids, it prevents energy from building up at small wave-numbers. When applied to the continuity equation to stabilize the water surface, the addition of the shock capturing scheme does not affect mass conservation. We found that our model improves the continuous and discontinuous Galerkin solutions alike in the proximity of sharp fronts propagating on wet surfaces. Furthermore, in the presence of wet/dry interfaces, however, the model needs to be enhanced with the addition of an inundation scheme which, however, we do not address in this paper.

58 GEOSCIENCES↗