Wave Simulations (wave + simulation)

Distribution by Scientific Domains


Selected Abstracts


A finite volume solver for 1D shallow-water equations applied to an actual river

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2002
N. Gouta
Abstract This paper describes the numerical solution of the 1D shallow-water equations by a finite volume scheme based on the Roe solver. In the first part, the 1D shallow-water equations are presented. These equations model the free-surface flows in a river. This set of equations is widely used for applications: dam-break waves, reservoir emptying, flooding, etc. The main feature of these equations is the presence of a non-conservative term in the momentum equation in the case of an actual river. In order to apply schemes well adapted to conservative equations, this term is split in two terms: a conservative one which is kept on the left-hand side of the equation of momentum and the non-conservative part is introduced as a source term on the right-hand side. In the second section, we describe the scheme based on a Roe Solver for the homogeneous problem. Next, the numerical treatment of the source term which is the essential point of the numerical modelisation is described. The source term is split in two components: one is upwinded and the other is treated according to a centred discretization. By using this method for the discretization of the source term, one gets the right behaviour for steady flow. Finally, in the last part, the problem of validation is tackled. Most of the numerical tests have been defined for a working group about dam-break wave simulation. A real dam-break wave simulation will be shown. Copyright © 2002 John Wiley & Sons, Ltd. [source]


Entropy sources in a dynamical core atmosphere model

THE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 614 2006
Tim Woollings
Abstract Numerical atmosphere models are not generally constructed to ensure accurate treatment of entropy, but little is known about the significance of the resulting errors. This paper examines the entropy changes during a baroclinic wave simulation in a typical dynamical core model, specifically a ,-coordinate spectral model, which includes scale-selective dissipation terms in the form of a numerical hyperdiffusion. Lagrangian entropy conservation is found to be badly represented, with numerical transport errors resulting in cross-isentrope mass fluxes which are of the same size as those associated with some real diabatic processes. In a global average, the total entropy increases at a rate of just 0.5 mW m,2K,1. This, however, is seen to be the residual of two opposing numerical effects which are several times larger, namely the destruction of entropy by dispersion and Gibbs errors, and its creation by diffusion. The entropy generated by diffusion is shown to be remarkably insensitive to the details of the diffusion scheme. This leads us to hypothesize that the entropy source from diffusion is determined by the rate at which small scales are generated by the deformation field of the large-scale flow so that, while the diffusion mechanism is clearly unrealistic, the magnitude of the entropy source is, we argue, representative of that generated by physical dissipative processes in the real atmosphere. Even in this simple model it is not possible to quantify precisely the different entropy sources and sinks which combine to give the overall entropy change. However, we can say that if there is a systematic spurious entropy source in this model, then it is small, i.e. of size 0.5 mW m,2K,1 or smaller. Copyright © 2006 Royal Meteorological Society [source]


Calculation of the wave propagation angle in complex media: application to turning wave simulations

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2009
Xiaofeng Jia
SUMMARY The wave propagation angle is one of the key factors in seismic processing methods. For the dual-domain propagators, it is sometimes necessary to acquire the wave propagation angle in the space-frequency domain instead of the wavenumber domain or the angle domain. We propose a method dealing with this problem, in which the wavefield gradient is used for the calculation of the wave propagation angle. The wavefield gradient can be directly obtained by either the finite difference approximation or the marching expression of the propagator. This method is not applicable in the case of extremely low frequency due to the comparability between the wavelength and the grid interval. Combined with the superwide-angle one-way propagator, this approach is instrumental in simulating the turning wave, which is hard to be handled by the traditional one-way propagator. Numerical examples show the good performance of the superwide-angle one-way propagator with our approach involved. The turning wave is modelled accurately; as a result, a high-quality image of the overhanging salt flank can be obtained. [source]


Reducing numerical diffusion in interfacial gravity wave simulations

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2005
O. B. Fringer
Abstract We demonstrate how the background potential energy is an excellent measure of the effective numerical diffusion or antidiffusion of an advection scheme by applying several advection schemes to a standing interfacial gravity wave. All existing advection schemes do not maintain the background potential energy because they are either diffusive, antidiffusive, or oscillatory. By taking advantage of the compressive nature of some schemes, which causes a decrease in the background potential energy, and the diffusive nature of others, which causes an increase in the background potential energy, we develop two background potential energy preserving advection schemes that are well-suited to study interfacial gravity waves at a density interface between two miscible fluids in closed domains such as lakes. The schemes employ total variation diminishing limiters and universal limiters in which the limiter is a function of both the upwind and local gradients as well as the background potential energy. The effectiveness of the schemes is validated by computing a sloshing interfacial gravity wave with a nonstaggered-grid Boussinesq solver, in which QUICK is employed for momentum and the pressure correction method is used, which is second-order accurate in time. For scalar advection, the present background potential energy preserving schemes are employed and compared to other TVD and non-TVD schemes, and we demonstrate that the schemes can control the change in the background potential energy due to numerical effects. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Multiple-Well, multiple-path unimolecular reaction systems.

INTERNATIONAL JOURNAL OF CHEMICAL KINETICS, Issue 4 2001

Vibrationally excited 2-methylhexyl radicals formed by shock wave activation or by chemical activation can isomerize by multiple pathways to form any of six stable isomers, can fragment by multiple CH and CC bond fission pathways, and can be collisionally stabilized. Master equation simulations of chemical activation and of shock wave activation are used to explore the generic behavior of this complicated coupled system. Selecting the argon pressure in chemical activation systems that produce the 2-methyl-1-hexyl radical isomer (1) can control the yield of specific isomers. Shock heating of 1 also shows a highly regular sequence of isomer formation. This regular behavior is because the first isomerization steps are faster than subsequent steps. Other radical isomers, such as 2-methyl-3-hexyl (3), do not show such regular behavior, because the first isomerization step is slower than subsequent steps. Incubation and unimolecular rate-constant fall-off are observed in the shock wave simulations. The unimolecular rate-constant fall-off for the coupled system produces low-pressure limiting rate constants proportional to [M]n, where n can be greater than unity. The fact that n can be greater than unity is a natural feature of multichannel coupled unimolecular reaction systems, but detection of the effect in experiments may be very demanding. © 2001 John Wiley & Sons, Inc. Int J Chem Kinet 33: 246,261, 2001 [source]