On the evolution and stability of interfacial ripples on and off the critical frequency

Abstract

Under consideration are interfaces between two media of different densities and which arise from the interaction between the Mth and Nth harmonics of the motion where 1 ≤ N < M. By means of the method of multiple scales in both space and time a pair of nonlinear coupled partial differential equations is derived which model the progression of the interface. The equations contain a detuning parameter [sgrave] which allow imperfections in the resonance to be taken into account. Stokes-type sinusoidal solutions to the equations were sought. It was found that solutions exist for all values of the interaction ratio M/N. In some situations interfaces exist at both exact and near resonance; while in others they are destroyed by amplifications in the detuning. In yet others, a quantity of detuning is actually necessary for the profiles to exist. In all cases, even when the parameters are fixed, a very large class of interface profiles is possible. Finally, the stability of the profiles is studied. It is found that some are quite stable, even to perturbations with wavenumbers close to the main flow.     

Non-linear channel–shoal dynamics in long tidal embayments

The dynamics of finite-amplitude bed forms in a tidal channel is studied with the use of an idealized morphodynamic model. The latter is based on depth-averaged equations for the tidal flow over a sandy bottom. The model considers phenomena on spatial scales of the order of the tidal excursion length. Transport of sediment mainly takes place as suspended load. The reference state of this model is characterized by a spatially uniform M₂ tidal current over a fixed horizontal bed. The temporal evolution of deviations from this reference state is governed by amplitude equations: these are a set of non-linear equations that describe the temporal evolution of bed forms. These equations are used to obtain new morphodynamic equilibria which may be either static or time-periodic. Several of these bottom profiles show strong similarity with the tidal bars that are observed in natural estuaries. The dependence of the equilibrium solutions on the value of bottom friction and channel width is investigated systematically. For narrow channels (width small compared to the tidal excursion length) stable static equilibria exist if bottom friction is slightly larger than r_cr. For channel widths more comparable to the tidal excursion length, multiple stable steady states may exist for bottom friction parameter values below r_cr. Regardless of channel width, stable time-periodic equilibria seem to emerge as the bottom friction is increased.Responsible Editor: Jens Kappenberg     

Symmetries of time-dependent magnetoconvection

Abstract

In the presence of a magnetic field, convection may set in at a stationary or an oscillatory bifurcation, giving rise to branches of steady, standing wave and travelling wave solutions. Numerical experiments provide examples of nonlinear solutions with a variety of different spatiotemporal symmetries, which can be classified by establishing an appropriate group structure. For the idealized problem of two-dimensional convection in a stratified layer the system has left-right spatial symmetry and a continuous symmetry with respect to translations in time. For solutions of period P the latter can be reduced to Z ₂ symmetry by sampling solutions at intervals of ½P. Then the fundamental steady solution has the spatiotemporal symmetry D ₂ = Z ₂ ? Z ₂ and symmetry-breaking yields solutions with Z ₂ symmetry corresponding to travelling waves, standing waves and pulsating waves. A further loss of symmetry leads to modulated waves. Interactions between the fundamental and its first harmonic are described by the group D _2h = D ₂ ? Z ₂ and its invariant subgroups, which describe solutions that are either steady or periodic in a uniformly moving frame. For a Boussinesq fluid in a layer with identical top and bottom boundary conditions there is also an up-down symmetry. With fixed lateral boundaries the spatiotemporal symmetries, again described by D _2h and its invariant subgroups, can be related to results obtained in numerical experiments and analysed by Nagata et al. (1990). With periodic boundary conditions, the full symmetry group, D _2h ?Z ₂, is of order 16. Its invariant subgroups describe pure and mixed-mode solutions, which may be steady states, standing waves, travelling waves, pulsating waves or modulated waves.     

Nonlinear waves on a coupled density front

Abstract

It is shown that the inclusion of the nonlinear terms in the equations of motion of a coupled density front of zero potential vorticity results in wave solutions which merely propagate with time. The linear theory, on the other hand, predicts an exponential temporal growth. The nonlinear equation admits steady solutions representing standing waves whereas if the nonlinear terms are omitted no steady solutions exist. The general initial value problem is difficult to solve numerically since the linear problem is ill posed.

In addition we prove that the general similarity solution of the nonlinear equation tends to zero for large times, at any point in space, regardless of the initial condition.     

Analytical models of steady-state plumes undergoing sequential first-order degradation

An exact, closed‐form analytical solution is derived for one‐dimensional (1D), coupled, steady‐state advection‐dispersion equations with sequential first‐order degradation of three dissolved species in groundwater. Dimensionless and mathematical analyses are used to examine the sensitivity of longitudinal dispersivity in the parent and daughter analytical solutions. The results indicate that the relative error decreases to less than 15% for the 1D advection‐dominated and advection‐dispersion analytical solutions of the parent and daughter when the Damköhler number of the parent decreases to less than 1 (slow degradation rate) and the Peclet number increases to greater than 6 (advection‐dominated). To estimate first‐order daughter product rate constants in advection‐dominated zones, 1D, two‐dimensional (2D), and three‐dimensional (3D) steady‐state analytical solutions with zero longitudinal dispersivity are also derived for three first‐order sequentially degrading compounds. The closed form of these exact analytical solutions has the advantage of having (1) no numerical integration or evaluation of complex‐valued error function arguments, (2) computational efficiency compared to problems with long times to reach steady state, and (3) minimal effort for incorporation into spreadsheets. These multispecies analytical solutions indicate that BIOCHLOR produces accurate results for 1D steady‐state, applications with longitudinal dispersion. Although BIOCHLOR is inaccurate in multidimensional applications with longitudinal dispersion, these multidimensional multispecies analytical solutions indicate that BIOCHLOR produces accurate steady‐state results when the longitudinal dispersion is zero. As an application, the 1D advection‐dominated analytical solution is applied to estimate field‐scale rate constants of 0.81, 0.74, and 0.69/year for trichloroethene, cis‐1,2‐dichloroethene, and vinyl chloride, respectively, at the Harris Palm Bay, FL, CERCLA site.     

A review of: “Problems of Solar and Stellar Oscillations. Iau Colloquium No. 66. Edited by D. O. Gough. D. Reidel Publishing Company, 493 pp. Df1195.00/U.S. $84.50 (ISBN 90-277-1554-8). 1983.”

Abstract

The model equations describing two-dimensional thermohaline convection of a Boussinesq fluid in a rotating horizontal layer are known to support multiple instabilities, depending on the values of certain control parameters (Arneodo et al., 1985). Most of these multiple instabilities have already been studied for double or triple diffusive convection, where behaviours ranging from simple steady to irregular motions have been found. Here we consider the one remaining bifurcation mentioned by Arneodo et al. (1985): the interaction between a steady and an oscillatory convection roll when the linear spectrum for a single wavenumber comprises one zero and one pair of purely imaginary eigenvalues. The method of centre manifolds and normal forms is used to derive evolution equations for the amplitudes of the convection rolls close to bifurcation and the behaviours associated with the equations is discussed.     

Calculation of magnetic field turbulent diffusivities and α-effect

Abstract

First the exact numerical solutions of DIA system of equations describing the transportation of magnetic field in an infinite medium are presented. It is assumed that the turbulence is stationary, homogeneous, isotropic and incompressible. The spectra of turbulence of δ-type and Kolmogorov's type were used. The steady state values of magnetic field diffusivity D_T and the α-effect coefficient α _T were calculated for various values of space-scale and lifetimes of these spectra and the spectra of helicity. Also investigated is the dependence of D_T and α _T on the degree of helicity. The corrections to the α _T -coefficient due to the contribution of four-order velocity correlators are given. The results are compared with those due to the self-consistent technique.     

Two-dimensional direct current (DC) resistivity inversion: Data space Occam's approach

A data space Occam's inversion algorithm for 2D DC resistivity data has been developed to seek the smoothest structure subject to an appropriate fit to the data. For traditional model space Gauss–Newton (GN) type inversion, the system of equations has the dimensions of M × M, where M is the number of model parameter, resulting in extensive computing time and memory storage. However, the system of equations can be mathematically transformed to the data space, resulting in a dramatic drop in its dimensions to N × N, where N is the number of data parameter, which is usually less than M. The transformation has helped to significantly reduce both computing time and memory storage. Numerical experiments with synthetic data and field data show that applying the data space technique to 2D DC resistivity data for various configurations is robust and accurate when compared with the results from the model space method and the commercial software RES2DINV.     

Analytical solutions to steady state unsaturated flow in layered,randomly heterogeneous soils via Kirchhoff transformation

In this study, we derive analytical solutions of the first two moments (mean and variance) of pressure head for one-dimensional steady state unsaturated flow in a randomly heterogeneous layered soil column under random boundary conditions. We first linearize the steady state unsaturated flow equations by Kirchhoff transformation and solve the moments of the transformed variable up to second order in terms of σ_Y and σ_β, the standard deviations of log hydraulic conductivity Y=ln(K_s) and of the log pore size distribution parameter β=ln(α). In addition, we also give solutions for the mean and variance of the unsaturated hydraulic conductivity. The analytical solutions of moment equations are validated via Monte Carlo simulations.     

Kinematic dynamo action in a sphere: Effects of periodic time-dependent flows on solutions with axial dipole symmetry

Choosing a simple class of flows, with characteristics that may be present in the Earth's core, we study the ability to generate a magnetic field when the flow is permitted to oscillate periodically in time. The flow characteristics are parameterised by D, representing a differential rotation, M, a meridional circulation, and C, a component characterising convective rolls. The dynamo action of all solutions with fixed parameters (steady flows) is known from earlier studies. Dynamo action is sensitive to these flow parameters and fails spectacularly for much of the parameter space where magnetic flux is concentrated into small regions, leading to high diffusion. In addition, steady flows generate only steady or regularly reversing oscillatory fields and cannot therefore reproduce irregular geomagnetic-type reversal behaviour. Oscillations of the flow are introduced by varying the flow parameters in time, defining a closed orbit in the space ( D,?M ). When the frequency of the oscillation is small, the net growth rate of the magnetic field over one period approaches the average of the growth rates for steady flows along the orbit. At increased frequency time-dependence appears to smooth out flux concentrations, often enhancing dynamo action. Dynamo action can be impaired, however, when flux concentrations of opposite signs occur close together as smoothing destroys the flux by cancellation. It is possible to produce geomagnetic-type reversals by making the orbit stray into a region where the steady flows generate oscillatory fields. In this case, however, dynamo action was not found to be enhanced by the time-dependence. A novel approach is being taken to solve the time-dependent eigenvalue problem where, by combining Floquet theory with a matrix-free Krylov-subspace method, we can avoid large memory requirements for storing the matrix required by the standard approach.     

Steady flows in solar magnetic solutions structures-a class of exact MHD solutions

Abstract

A class of exact solutions to the steady, two-dimensional magnetohydrodynamic equations ina cylindrical geometry is presented. These may model both closed and open magnetic structures found in the solar atmosphere. For closed structures, it is found that increasing the flow speed causes the summit of the arcade of closed magnetic fieldlines to rise. Parameter ranges also exist where the solution has regions of open and closed field, and so the solutions may be relevant for modelling flows in solar magnetic structures such as coronal streamers, X-ray bright points coronal plumes and coronal holes.     

High-order finite volume WENO schemes for the shallow water equations with dry states

The shallow water equations are used to model flows in rivers and coastal areas, and have wide applications in ocean, hydraulic engineering, and atmospheric modeling. These equations have still water steady state solutions in which the flux gradients are balanced by the source term. It is desirable to develop numerical methods which preserve exactly these steady state solutions. Another main difficulty usually arising from the simulation of dam breaks and flood waves flows is the appearance of dry areas where no water is present. If no special attention is paid, standard numerical methods may fail near dry/wet front and produce non-physical negative water height. A high-order accurate finite volume weighted essentially non-oscillatory (WENO) scheme is proposed in this paper to address these difficulties and to provide an efficient and robust method for solving the shallow water equations. A simple, easy-to-implement positivity-preserving limiter is introduced. One- and two-dimensional numerical examples are provided to verify the positivity-preserving property, well-balanced property, high-order accuracy, and good resolution for smooth and discontinuous solutions.     

Convection in a rotating cylinder with its sidewall having a finite thermal conductance

Abstract

An investigation is made of steady thermal convection of a Boussinesq fluid confined in a vertically-mounted rotating cylinder. The top and bottom endwall disks are thermal conductors at temperatures T_t and T_b with δT = T_t ? T_b >0. The vertical sidewall has a finite thermal conductance. A Newtonian heat flux condition is adopted at the sidewall. The Rayleigh number of the fluid system is large to render a boundary layer-type flow. Finite-difference numerical solutions to the full Navier-Stokes equations are obtained. The vertical motions within the buoyancy layer along the sidewall induce weak meridional flows in the interior. Because of the Coriolis acceleration, the meridional flows give rise to azimuthal flows relative to the rotating container. Strong vertical gradients of azimuthal flows exist in the regions near the endwalls. As the stratification effect increases, concentration of flow gradients in thin endwall boundary layers becomes more pronounced. The azimuthal flow field exhibits considerable horizontal gradients. The temperature field develops horizontal variations superposed on the dominant vertical distribution. As either the sidewall thermal conductance or the stratification effect decreases, the temperature distribution tends to the profile varying linearly with height. Comparisons of the sizes of the dynamic effects demonstrate that, in the bulk of flow field, the vertical shear of azimuthal velocity is supported by the horizontal temperature gradient, resulting in a thermal-wind relation.     

Empirical shear wave velocity equations in terms of characteristic soil indexes

An investigation to systematize empirical equations for the shear wave velocity of soils was made in terms of four characteristic indexes. The adopted indexes are the N-value of the Standard Penetration Test, depth where the soil is situated, geological epoch and soil type. As some of these indexes are variates belonging to interval scales while others belong to nominal or ordinal scales, the technique known as a multivariate analysis cannot be employed. A new approach to the theory of quantification, after C. Hayashi, was introduced and developed for solving this difficulty. Fifteen sets of empirical equations to estimate low strain shear wave velocity theoretically may be obtained by combining the above four indexes. All of these sets were derived by use of about 300 data, and their accuracies were evaluated by means of correlation coefficients between the measured and estimated shear wave velocities. The best equation was found to be the one which included all the indexes, and its correlation coefficient was 0.86. The empirical equation relating the standard penetration N-value to the shear wave velocity provided a correlation of only 0.72, and is one of the lowest ranking among the 15 sets of equations.     

The progression of flow rates in variable declining rate filter systems

It has been demonstrated that the mathematical model of variable declining rate filters developed by Di Bernardo may be described by (z + 1) non‐linear equations, where z is the number of filters in a bank. Three approximate solutions to this system of equations have been developed and then verified by comparison with numerical solution and published experimental data. Two of these solutions appeared to be very accurate, while the third showed higher, but still acceptable errors of calculation. According to this approximation, flow rates through filters are elements of a geometric progression.     

A methodology for locating strong motion arrays

This paper investigates a methodology for locating strong motion accelerographs in a seismically active region. Starting with the probability density of earthquakes in a given region, the paper attempts, within the framework of optimization theory, to formulate the following two questions: (1) given N accelerographs, where should they be located in a seismically active zone, and (2) having fixed these N accelerographs, where should the next M be located? Three different cost functions are presented. Some closed form solutions are illustrated for problems when N and M are small. For larger arrays, numerical optimization is resorted to. To demonstrate the methodology, a region with J faults, with given spatial locations is selected. An efficient algorithm for optimization is utilized and the technique illustrated. Good agreement with closed form solutions obtained in some simple cases is indicated. Specific application of the method to the placing of twenty strong motion instruments in a seismically active area has been carried out, and the patterns of sensor location, for each of the cost functions, illustrated.     

Nonlinear dynamics of a stellar dynamo in a spherical shell

Abstract

A simple mean-field model of a nonlinear stellar dynamo is considered, in which dynamo action is supposed to occur in a spherical shell, and where the only nonlinearity retained is the influence of the Lorentz forces on the zonal flow field. The equations are simplified by truncating in the radial direction, while full latitudinal dependence is retained. The resulting nonlinear p.d.e.'s in latitude and time are solved numerically, and it is found that while regular dynamo wave type solutions are stable when the dynamo number D is sufficiently close to its critical value, there is a wide variety of stable solutions at larger values of D. Furthermore, two different types of dynamo can coexist at the same parameter values. Implications for fields in late-type stars are discussed.     

A theoretical investigation of bed-form shapes

The deformation of movable boundaries under the action of an applied turbulent shear stress is well known. The resulting bed forms often are highly organized and nearly two-dimensional, which makes them an intriguing focus of study considering that they are generated in both steady and oscillatory turbulent flows. Many past studies share a common approach in which an infinitesimal perturbation is prescribed and the resulting growth or decay patterns are examined. In this approach, the bed forms are usually sinusoidal and the perturbation analysis does not provide a theoretical prediction of equilibrium bed-form geometry. An alternative approach is suggested here in which the forcing terms (pressure and stress) are prescribed parametrically and the governing equations are solved for the flow velocity and the associated boundary deformation. Using a multilayered approach, in which the bottom boundary layer is divided into a discrete, yet, arbitrary number of finite layers, analytical solutions for the horizontal current and bed profile are derived. The derivations identify two nondimensional parameters, p₀/u₀² and ₀/kh₀u₀², which modulate the amplitude of the velocity fluctuations and boundary deformation. For the case of combined pressure and stress divergence anomalies, the magnitude of the front face and lee slopes exhibit an asymmetry that is consistent with observed bed forms in steady two-dimensional flows.Responsible Editor: Jens Kappenberg     

Linear theory of rotating fluids using spherical harmonics part I: Steady flows

Abstract

It is shown that a systematic development of physical quantities using spherical harmonics provides analytical solutions to a whole class of linear problems of rotating fluids.

These solutions are regular throughout the whole domain of the fluid and are not much affected by the equatorial singularity of steady boundary layers in spherical geometries.

A comparison between this method and the one based on boundary layer theory is carried out in the case of the steady spin-up of a fluid inside a sphere.