Beach water table fluctuations due to spring-neap tides: moving boundary effects

Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection-diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring-neap tides (a bichromatic tidal system with the fundamental frequencies wt and wt) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency omega(1) - omega(2), i.e., the spring-neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations. (C) 2000 Elsevier Science Ltd. All rights reserved.

Consistent Dirichlet Boundary Conditions for Numerical Solution of Moving Boundary Problems

We consider the imposition of Dirichlet boundary conditions in the finite element modelling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass.

Moving boundary value problems for the wave equation

We study certain boundary value problems for the one-dimensional wave equation posed in a time-dependent domain. The approach we propose is based on a general transform method for solving boundary value problems for integrable nonlinear PDE in two variables, that has been applied extensively to the study of linear parabolic and elliptic equations. Here we analyse the wave equation as a simple illustrative example to discuss the particular features of this method in the context of linear hyperbolic PDEs, which have not been studied before in this framework.

Moving Boundary: A Virtual Museum of the Ethnic Groups

The virtual museum of ethnic groups “Moving boundary” is a proposition for dialogue on an important topic in today's world – ethnic coexistence. The project's goal is to make various resources (photos, documents, etc.) available for use by teachers, professors and students, thus forming and strengthening connections between educational institutions.

On the Moving Boundary Hitting Probability for the Brownian Motion

2000 Mathematics Subject Classification: 60J65.

On the solute coupling at the moving solid/liquid interface during equiaxed solidification

Integral mass conservation was widely accepted for the solute coupling to solve solute redistribution during equiaxed solidification so far. The present study revealed that the integral form was invalid for moving boundary problems as it could not represent the mass balance at the moving interface. Accordingly, differential mass conservation at the solid/liquid interface was used to solve solute diffusion for spherical geometry. The model was applied for hydrogen diffusion in solidification to validate that the hydrogen enrichment was significant and depended on the growth rate. (c) 2006 American Institute of Physics.

Factorization of elliptic boundary value problems by invariant embedding and application to overdetermined problems

Dissertação para obtenção do Grau de Doutor em Matemática

Boundary value problems for third-order linear PDEs in time-dependent domains

We present the extension of a methodology to solve moving boundary value problems from the second-order case to the case of the third-order linear evolution PDE qt + qxxx = 0. This extension is the crucial step needed to generalize this methodology to PDEs of arbitrary order. The methodology is based on the derivation of inversion formulae for a class of integral transforms that generalize the Fourier transform and on the analysis of the global relation associated with the PDE. The study of this relation and its inversion using the appropriate generalized transform are the main elements of the proof of our results.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

A moving mesh approach for modelling avascular tumour growth

A key step in many numerical schemes for time-dependent partial differential equations with moving boundaries is to rescale the problem to a fixed numerical mesh. An alternative approach is to use a moving mesh that can be adapted to focus on specific features of the model. In this paper we present and discuss two different velocity-based moving mesh methods applied to a two-phase model of avascular tumour growth formulated by Breward et al. (2002) J. Math. Biol. 45(2), 125-152. Each method has one moving node which tracks the moving boundary. The first moving mesh method uses a mesh velocity proportional to the boundary velocity. The second moving mesh method uses local conservation of volume fraction of cells (masses). Our results demonstrate that these moving mesh methods produce accurate results, offering higher resolution where desired whilst preserving the balance of fluxes and sources in the governing equations.

Ion flows and heating at a contracting polar-cap boundary: GISMOS observations indicating viscous-like interaction on the flanks of the magnetotail

This paper complements that in this issue by Clauer et al. concerning the international GISMOS campaign of 3–5 June 1987. From a detailed study of the EISCAT data, the polar-cap boundary, as defined by an almost shear east-west convection reversal, is found to contract across the EISCAT field of view between 04 and 07 MLT. An annulus of enhanced ion temperature and non-thermal plasma is observed immediately equatorward of the contracting boundary due to the lag in the response of the neutral-wind pattern to the change in ion flows. The ion flow inside the polar cap and at the boundary is shown to be relatively smooth, compared with that in the auroral oval, at 15-second resolution. The flow at the boundary is directed poleward, with velocities which exceed that of the boundary itself. The effect of velocity shears on the beamswinging technique used to derive the ion flows has been analysed in detail and it is found that spurious flows across a moving boundary can be generated. However, these are much smaller than the observed flows into the polar cap and cannot explain the 7 kV potential difference across the observed segment of the cap boundary between 04:30–06:30 UT. The ion temperature enhancements at the two observing azimuths is used to define the boundary orientation. The results are consistent with recent observations of slow anti-sunward flow of closed field lines on the flanks of the geomagnetic tail, which appears to be generated by some form of “viscous” coupling to the magnetosheath plasma.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.