A class of analytic solutions of the magnetic force-free field equations

Abstract

Formation of electric current sheets in the corona is thought to play an important role in solar flares, prominences and coronal heating. It is therefore of great interest to identify magnetic field geometries whose evolution leads to variations in B over small length-scales. This paper considers a uniform field B ₀[zcirc], line-tied to rigid plates z = ±l, which are then subject to in-plane displacements modeling the effect of photospheric motion. The force-free field equations are formulated in terms of field-line displacements, and when the imposed plate motion is a linear function of position, these reduce to a 4 × 4 system of nonlinear, second-order ordinary differential equations. Simple analytic solutions are derived for the cases of plate rotation and shear, which both tend to form singularities in certain parameter limits. In the case of plate shear there are two solution branches—a simple example of non-uniqueness.     

Dynamic interaction between flexible strip foundations

This paper is concerned with the dynamic response of a system of flexible strip foundations resting in smooth contact with a homogeneous isotropic viscoelastic half space. An arbitrary number of foundations with different flexibilities and geometries subjected to time-harmonic distributed loadings are considered in the formulation. The response of each strip foundation is governed by the classical plate theory and its transverse deflection is represented by an admissible function containing a set of generalized co-ordinates. A coupled variational-Green's function scheme is employed to establish the equations of motion of the strip foundation system. The numerical stability and convergence of the solution scheme are established. The influence of the foundation flexibility, distance between adjacent foundations and frequency of motion on the response of the foundation system is investigated in the numerical study.     

Nonlinear sloshing in a floating‐roofed oil storage tank under long‐period seismic ground motion

A hybrid analytical and FEM is proposed to investigate the nonlinear sloshing in a floating‐roofed oil storage tank under long‐period seismic ground motion. The tank is composed of a rigid cylindrical wall and a flat bottom, whereas the floating roof is treated as an elastic plate undergoing large deflection. The contained liquid is assumed to be inviscid and incompressible, and the flow is assumed to be irrotational. The method of analysis is based on representation of the liquid motion by superposing the analytical modes that satisfy the Laplace equation and the rigid wall and bottom boundary conditions. The FEM is then applied to solve the remaining kinematic and dynamic boundary conditions at the moving liquid surface coupled with the nonlinear equation of motion of the floating roof. This requires only the discretization of the liquid surface and the floating roof into finite elements, thus leading to a computationally efficient and accurate method compared with full numerical analysis. As numerical examples to illustrate the applicability of the proposed method, two oil storage tanks with single‐deck type floating roofs damaged during the 2003 Tokachioki earthquake are studied. It is shown that the nonlinear oscillation modes with the circumferential wave numbers 0, 2 and 3 caused by the finite liquid surface elevation as well as the membrane action due to large deflection of the deck produce excessively large stresses in the pontoon, which may cause the catastrophic failure of pontoon followed by the submergence of the roof. Copyright © 2012 John Wiley & Sons, Ltd.     

An experimental assessment of the added mass of some plates vibrating in water

A reservoir of water is contained by a concrete valley block, a ferrocement wall and a steel plate. Both wall and plate contain an array of pressure transducer sockets (Figures 1 and 2). Using the M.A.M.A.¹ equipment pure modes of vibration are excited. Frequency and mode shape are measured with the reservoir empty. When the reservoir is full hydrodynamic pressure is also measured. These hydrodynamic pressures are compared with Chopra's² two-dimensional, series solution, which includes compressibility of water, and with two- and three-dimensional finite element solutions of Laplace's equation, which do not include compressibility. Chopra's solution is unsatisfactory for modes which contain a vertical node line. The best agreement between experimental and theoretical hydrodynamic pressure is obtained when the latter is obtained from three-dimensional solutions of Laplace's equations, indicating that compressibility does not play a significant rǒle. This conclusion is supported by agreement between experimental frequencies (reservoir full) and those calculated using added mass obtained from the Laplace solution. Similar conclusions were reached from tests on a floating steel plate, suspended in the surface of the reservoir by a soft spring. Here, dynamic pressure measurements were not made, reliance being placed on agreement between calculated and measured frequencies and mode shapes.     

Study of forward–backward solute dispersion profiles in a semi-infinite groundwater system

ABSTRACT

Forward–backward solute dispersion with an intermediate point source in one-dimensional semi-infinite homogeneous porous media is studied in this paper. Solute transport under sorption conditions, first-order decay and zero-order production terms are included. The first type of boundary condition is taken as a constant point source at an intermediate point from where forward and backward solute dispersion is examined. The Laplace transform method is adopted to solve the governing equation analytically. All the analytical results are obtained in graphical form to investigate the forward–backward solute transport in porous media for various hydrological input data. The graphical nature of the analytical solution is compared with numerical data taken from existing literature and similar results are obtained. Also, numerical solution of the governing equation is obtained by the Crank-Nicolson finite difference scheme and validated with the analytical solution, which demonstrates good agreement between them. Accuracy of the solution is also observed by using RMSE.     

Hydro-elastic wave proliferation over an impermeable seabed with bottom deformation

ABSTRACT

A hydro-elastic frame has been considered to investigate the proliferation of waves over small base deformation on an infinitely extended flexible seabed. The flexible base surface is assumed as a thin elastic plate of very small thickness and it depends on the Euler–Bernoulli beam equation. For any particular frequency, there are two different modes of time-harmonic propagating wave exists rather than one mode of propagating wave along the positive horizontal direction. The waves with smaller wavenumber spread along the free-surface of the sea (say, free-surface mode) and the waves with higher wavenumber spread along the flexible base surface (say, flexural mode). A simplified perturbation approach is utilised to bring down the entire equations which govern the original boundary value problem (bvp) to a less complex bvp for the first-order velocity potential function. The first-order potential function along with the first-order reflection and transmission coefficients for both modes are calculated by a procedure based upon Fourier transform approach. A shape of sinusoidal swells flexible base surface is taken as an example to approve the scientific results. It is observed that when the train of normal incident propagating wave spreads over base distortion because of either the free-surface unsettling influence or the flexural wave movement in the sea, the reflected and transmitted energy are always feasible to be exchanged from one particular wave mode to another wave mode. Furthermore, we notice that the realistic changes in the flexural rigidity behaviour on the flexible base surface of the sea have a significant effect on the problem of water wave proliferation over small base deformation. Moreover, the energy conservation equation is derived with the help of the Green's integral theorem. The results for the values of reflection and transmission coefficients obtained for both the free-surface unsettling influence as well as flexural wave movement in the fluid are found to satisfy the energy conservation equation almost accurately.     

A double-porosity model for a fractured aquifer with non-Darcian flow in fractures

Abstract

Non-Darcian flow in a finite fractured confined aquifer is studied. A stream bounds the aquifer at one side and an impervious stratum at the other. The aquifer consists of fractures capable of transmitting water rapidly, and porous blocks which mainly store water. Unsteady flow in the aquifer due to a sudden rise in the stream level is analysed by the double-porosity conceptual model. Governing equations for the flow in fractures and blocks are developed using the continuity equation. The fluid velocity in fractures is often too high for the linear Darcian flow so that the governing equation for fracture flow is modified by Forcheimer's equation, which incorporates a nonlinear term. Governing equations are coupled by an interaction term that controls the quasi-steady-state fracture—block interflow. Governing equations are solved numerically by the Crank-Nicolson implicit scheme. The numerical results are compared to the analytical results for the same problem which assumes Darcian flow in both fractures and blocks. Numerical and analytical solutions give the same results when the Reynolds number is less than 0.1. The effect of nonlinearity on the flow appears when the Reynolds number is greater than 0.1. The higher the rate of flow from the stream to the aquifer, the higher the degree of nonlinearity. The effect of aquifer parameters on the flow is also investigated. The proposed model and its numerical solution provide a useful application of nonlinear flow models to fractured aquifers. It is possible to extend the model to different types of aquifer, as well as boundary conditions at the stream side. Time-dependent flow rates in the analysis of recession hydrographs could also be evaluated by this model.     

Oblique wave diffraction by a flexible floating structure in the presence of a submerged flexible structure

Expansion formulae associated with the interaction of oblique surface gravity waves with a floating flexible plate in the presence of a submerged horizontal flexible structure are derived using Green’s integral theorem in water of finite and infinite water depths. The associated Green’s functions are derived using the fundamental solution associated with the reduced wave equation. The integral forms of the Green’s functions and the velocity potentials are advantageous over the eigenfunction expansion method in situation when the roots of the dispersion relation coalesce. As an application of the expansion formulae, diffraction of oblique waves by a finite floating elastic plate in the presence of a submerged horizontal flexible membrane is investigated in water of finite depth. The accuracy of the numerical computation is demonstrated by analysing the convergence of the complex amplitude of the reflected waves and the energy relation. Effect of the submerged membrane on the diffraction of surface waves is studied by analysing the reflection and transmission coefficients for various parametric values. Further, the derivation of long wave equation under shallow water approximation is derived in a direct manner in the appendix. The concept and methodology can be easily extended to deal with acoustic wave interaction with flexible structures and related problems of mathematical physics and engineering.     

Coupled simulation of transient bed evolution in alluvial channels

Abstract

In dealing with the transient sediment transport problem, the commonly used uncoupled model may not be suitable. The uncoupling technique is intended to separate the physical coupling phenomenon of water flow and sediment transport into two independent processes. Very often, as a result, severe numerical oscillation and solution instability problems appear in the simulation of transient sediment transport in alluvial channels. The coupled model, which simultaneously solves water flow continuity, momentum and sediment continuity equations, gives fewer numerical oscillation and solution instability problems. In this article, a coupled model using a matrix double-sweep method to solve the system of nonlinear algebraic equations has been developed. Several test runs designed on the basis of a schematic model have been performed. The numerical oscillation and solution instability problems have been investigated through a comparison with those obtained from an uncoupled model. Based on the proposed case studies, it can be concluded that, for transient bed evolution, the performance of the coupled model is much better than that of the uncoupled model. The numerical oscillation is reduced and the solution is more stable. This newly developed coupled model was also applied to the Cho-Shui River in Taiwan. This application study implied that the effect of the peaky flood wave propagation on the bed evolution could be simulated better by the coupled model than by the uncoupled model.     

A note on the howard-Malkus-Whitehead floating heat sources

Abstract

Some deficiencies in a recent paper by Howard, Malkus and Whitehead are examined. The problem is reformulated in terms of an integro-differential equation, from which both asymptotic and numerical solutions are obtained.     

Non-local effects in the mean-field disc dynamo I. An asymptotic expansion

Abstract

We reconsider thin-disc global asymptotics for kinematic, axisymmetric mean-field dynamos with vacuum boundary conditions. Non-local terms arising from a small but finite radial field component at the disc surface are consistently taken into account for quadrupole modes. As in earlier approaches, the solution splits into a local part describing the field distribution along the vertical direction and a radial part describing the radial (global) variation of the eigenfunction. However, the radial part of the eigenfunction is now governed by an integro-differential equation whose kernel has a weak (logarithmic) singularity. The integral term arises from non-local interactions of magnetic fields at different radii through vacuum outside the disc. The non-local interaction can have a stronger effect on the solution than the local radial diffusion in a thin disc, however the effect of the integral term is still qualitatively similar to magnetic diffusion.     

Some results on the approach of a slowly-varyng gravity wave to a beach

Abstract

Results from the theory of slowly-varying solutions of the non-linear shallow water equations (Varley et al. (1971)) are used to estimate maximum wave height amplifications on a beach of arbitrary shape. Carrier's solution (1966) for a large disturbance out at sea is used as input. Generalizations of the slowly-varying results to non-plane situations are also presented.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Determination of Forchheimer equation coefficients a and b

This study focuses on the determination of the Forchheimer equation coefficients a and b for non‐Darcian flow in porous media. Original theoretical equations are evaluated and empirical relations are proposed based on an investigation of available data in the literature. The validity of these equations is checked using existing experimental data, and their accuracy versus existing approaches is studied. On the basis of this analysis, some insight into the physical background of the phenomenon is also provided. The dependence of the coefficients a and b on the Reynolds number is also detected, and potential future research areas, e.g. investigation of inertial effects for consolidated porous media, are pointed out. Copyright © 2006 John Wiley & Sons, Ltd.     

The unsteady plane flow of ice-sheets: A parabolic problem with two moving boundaries

Abstract

Finite difference algorithms have been developed to solve a one-dimensional non-linear parabolic equation with one or two moving boundaries and to analyse the unsteady plane flow of ice-sheets. They are designed to investigate the response of an ice-sheet to changes in climate, and to reconstruct climatic changes implied by past ice-sheet variations inferred from glacial geological data. Two algorithms are presented and compared. The first, a fixed domain method, replaces time as an independent variable with span. The grid interval in real space is kept constant, and thus the number of grid points changes with span. The second, a moving mesh method, retains time as one of the independent variables, but normalises the spatial variable relative to the span, which now enters the diffusion and advection coeficients in the parabolic equation for the surface profile.

Crank-Nicholson schemes for the solution of the equations are constructed, and iterative schemes for the solution of the resulting non-linear equations are considered.

Boundary (margin) motion is governed by the surface slope at the margin. Differentiation of the evolution equations results in an evolution equation for the margin slopes. It is shown that incorporation of this evolution equation, while not formally increasing the accuracy of the finite difference schemes, in practice increases accuracy of the solution.     

Periodic and aperiodic dynamo waves

Abstract

In order to show that aperiodic magnetic cycles, with Maunder minima, can occur naturally in nonlinear hydromagnetic dynamos, we have investigated a simple nonlinear model of an oscillatory stellar dynamo. The parametrized mean field equations in plane geometry have a Hopf bifurcation when the dynamo number D=1, leading to Parker's dynamo waves. Including the nonlinear interaction between the magnetic field and the velocity shear results in a system of seven coupled nonlinear differential equations. For D>1 there is an exact nonlinear solution, corresponding to periodic dynamo waves. In the regime described by a fifth order system of equations this solution remains stable for all D and the velocity shear is progressively reduced by the Lorentz force. In a regime described by a sixth order system, the solution becomes unstable and successive transitions lead to chaotic behaviour. Oscillations are aperiodic and modulated to give episodes of reduced activity.     

Dynamic response of a flexible plate on saturated soil layer

An analytical approach is developed to study the dynamic response of a flexible plate on single-layered saturated soil. The analysis is based on Biot's two-phased theory of poroelasticity and also on the classical thin-plate theory. First, the governing differential equations for saturated soil are solved by the use of Hankel transform. The general solutions of the skeleton displacements, stresses, and pore pressures, derived in the transformed domain, are subsequently incorporated into the imposed boundary conditions, which leads to a set of dual integral equations describing the corresponding mixed boundary value problem. These governing integral equations are finally reduced to the Fredholm integral equations of the second kind and solved by standard numerical procedures. The accuracy of the present solution is validated via comparisons with existing solutions for an ideal elastic half-space. Furthermore, some numerical results are presented to show the influences of the layer depth, the plate flexibility, and the soil porosity on the dynamic compliances.     

Inferring groundwater system dynamics from hydrological time-series data

Abstract

The problem of identifying and reproducing the hydrological behaviour of groundwater systems can often be set in terms of ordinary differential equations relating the inputs and outputs of their physical components under simplifying assumptions. Conceptual linear and nonlinear models described as ordinary differential equations are widely used in hydrology and can be found in several studies. Groundwater systems can be described conceptually as an interlinked reservoir model structured as a series of nonlinear tanks, so that the groundwater table can be schematized as the water level in one of the interconnected tanks. In this work, we propose a methodology for inferring the dynamics of a groundwater system response to rainfall, based on recorded time series data. The use of evolutionary techniques to infer differential equations from data in order to obtain their intrinsic phenomenological dynamics has been investigated recently by a few authors and is referred to as evolutionary modelling. A strategy named Evolutionary Polynomial Regression (EPR) has been applied to a real hydrogeological system, the shallow unconfined aquifer of Brindisi, southern Italy, for which 528 recorded monthly data over a 44-year period are available. The EPR returns a set of non-dominated models, as ordinary differential equations, reproducing the system dynamics. The choice of the representative model can be made both on the basis of its performance against a test data set and based on its incorporation of terms that actually entail physical meaning with respect to the conceptualization of the system.

Citation Doglioni, A., Mancarella, D., Simeone, V. & Giustolisi, O. (2010) Inferring groundwater system dynamics from hydrological time-series data. Hydrol. Sci. J. 55(4), 593–608.