An exact peak capturing and essentially oscillation-free (EPCOF) algorithm, consisting of advection-dispersion decoupling, backward method of characteristics, forward node tracking, and adaptive local grid refinement, is developed to solve transport equations. This algorithm repr...
Fractional Advective-Dispersive Equation as a Model of Solute Transport in Porous Media
Technology Transfer Automated Retrieval System (TEKTRAN)
Understanding and modeling transport of solutes in porous media is a critical issue in the environmental protection. The common model is the advective-dispersive equation (ADE) describing the superposition of the advective transport and the Brownian motion in water-filled pore space. Deviations from...
Analytical solution for the advection-dispersion transport equation in layered media
Technology Transfer Automated Retrieval System (TEKTRAN)
The advection-dispersion transport equation with first-order decay was solved analytically for multi-layered media using the classic integral transform technique (CITT). The solution procedure used an associated non-self-adjoint advection-diffusion eigenvalue problem that had the same form and coef...
Technology Transfer Automated Retrieval System (TEKTRAN)
Analytical solutions of the advection-dispersion equation and related models are indispensable for predicting or analyzing contaminant transport processes in streams and rivers, as well as in other surface water bodies. Many useful analytical solutions originated in disciplines other than surface-w...
İbiş, Birol
2014-01-01
This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662
Ibiş, Birol; Bayram, Mustafa
2014-01-01
This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662
Technology Transfer Automated Retrieval System (TEKTRAN)
Analytical solutions of the advection-dispersion solute transport equation remain useful for a large number of applications in science and engineering. In this paper we extend the Duhamel theorem, originally established for diffusion type problems, to the case of advective-dispersive transport subj...
Technology Transfer Automated Retrieval System (TEKTRAN)
The classical model to describe solute transport in soil is based on the advective-dispersive equation where Fick’s law is used to explain dispersion. From the microscopic point of view this is equivalent to consider that the motion of the particles of solute may be simulated by the Brownian motion....
Solution of the advection-dispersion equation: Continuous load of finite duration
Runkel, R.L.
1996-01-01
Field studies of solute fate and transport in streams and rivers often involve an. experimental release of solutes at an upstream boundary for a finite period of time. A review of several standard references on surface-water-quality modeling indicates that the analytical solution to the constant-parameter advection-dispersion equation for this type of boundary condition has been generally overlooked. Here an exact analytical solution that considers a continuous load of unite duration is compared to an approximate analytical solution presented elsewhere. Results indicate that the exact analytical solution should be used for verification of numerical solutions and other solute-transport problems wherein a high level of accuracy is required. ?? ASCE.
Bad behavior of Godunov mixed methods for strongly anisotropic advection-dispersion equations
NASA Astrophysics Data System (ADS)
Mazzia, Annamaria; Manzini, Gianmarco; Putti, Mario
2011-09-01
We study the performance of Godunov mixed methods, which combine a mixed-hybrid finite element solver and a Godunov-like shock-capturing solver, for the numerical treatment of the advection-dispersion equation with strong anisotropic tensor coefficients. It turns out that a mesh locking phenomenon may cause ill-conditioning and reduce the accuracy of the numerical approximation especially on coarse meshes. This problem may be partially alleviated by substituting the mixed-hybrid finite element solver used in the discretization of the dispersive (diffusive) term with a linear Galerkin finite element solver, which does not display such a strong ill conditioning. To illustrate the different mechanisms that come into play, we investigate the spectral properties of such numerical discretizations when applied to a strongly anisotropic diffusive term on a small regular mesh. A thorough comparison of the stiffness matrix eigenvalues reveals that the accuracy loss of the Godunov mixed method is a structural feature of the mixed-hybrid method. In fact, the varied response of the two methods is due to the different way the smallest and largest eigenvalues of the dispersion (diffusion) tensor influence the diagonal and off-diagonal terms of the final stiffness matrix. One and two dimensional test cases support our findings.
NASA Astrophysics Data System (ADS)
Benson, D. A.; Zhang, Y.
2006-12-01
Conservative solute transport through natural media is typically "anomalous" or non-Fickian. The anomalous transport may be characterized by faster than linear growth of the centered second moment, or non-Gaussian leading or trailing edges of a plume emanating from a point source. These characteristics develop because of non-local dependence on either past (time) or far upstream (space) concentrations. Non-local equations developed to describe anomalous dispersion usually focus on constant transport parameters and/or independence of the transport on space dimension. These simplifications have been useful for fitting simple transport processes, such as laboratory column tests or 1-D projections of field data. However, they may be insufficient for real field settings, where direction-dependent depositional processes and nonstationary heterogeneity can occur. We develop a generalized, multi-dimensional, spatiotemporal fractional advection- dispersion equation (fADE) with variable parameters to characterize regional-scale anomalous dispersion processes including trapping in immobile zones and/or super-Fickian rapid transport. A Lagrangian numerical model of the space-time fractional transport equation is developed in which solute particles can disperse in both space and time, depending on the medium heterogeneity properties, such as the connectivity and statistical distributions of high versus low-permeability deposits. In the generalized fADE, the range of the order of fractional time derivative is (0 2], representing a wide range of possible trapping behavior. The extension of the order to the range (1 2] is novel to transport theory. We apply the numerical model in 1-D and 2-D to the MADE site tritium plumes, and results indicate that this method can capture the main behaviors of realistic plumes, including local variations of spreading, direction-dependent scaling rates, and arbitrary rapid transport along preferential flow paths. Since the governing equation
NASA Astrophysics Data System (ADS)
Dean, A. M.; Benson, D. A.; Major, E.
2010-12-01
By adding a fractional-in-time term to the traditional advection dispersion equation, a model is able to simulate a late-time heavy-tailed contaminant breakthrough curve. This heavy-tailed breakthrough curve is observed in data collected during a conservative tracer “push-pull” test at the Macrodispersion Experiment (MADE) site. A time fractional advection dispersion equation (fADE) is able to predict power law tailing of conservative solutes by accounting for solutes transferring between the mobile and relatively immobile phases. Solutes can become trapped in a low permeability zone where the transport is controlled by diffusion instead of advection. It has been observed that the late-time heavy-tailed breakthrough curve may follow a power law due to the movement into these low flow zones. By solving the time fADE in a particle tracking program (SLIM-FAST) the model accounts for mass transfer between various phases and produces the same power law tail as observed in field data. For the implementation of the time fADE, in SLIM-FAST, the particles move based on a random-walk motion but have the ability to transition into a relatively immobile phase after (exponentially) random mobile times. Following a period in the immobile phase, the particle re-enters the mobile phase to be moved by advection and Fickian dispersion. To test the fADE approach, a recent single-well push-pull tracer test at the MADE site is reproduced using a groundwater flow code (ParFlow) and a particle tracking code (SLIM-FAST) using various immobile residence-time distributions.
Knopman, Debra S.; Voss, Clifford I.
1987-01-01
The spatial and temporal variability of sensitivities has a significant impact on parameter estimation and sampling design for studies of solute transport in porous media. Physical insight into the behavior of sensitivities is offered through an analysis of analytically derived sensitivities for the one-dimensional form of the advection-dispersion equation. When parameters are estimated in regression models of one-dimensional transport, the spatial and temporal variability in sensitivities influences variance and covariance of parameter estimates. Several principles account for the observed influence of sensitivities on parameter uncertainty. (1) Information about a physical parameter may be most accurately gained at points in space and time. (2) As the distance of observation points from the upstream boundary increases, maximum sensitivity to velocity during passage of the solute front increases. (3) The frequency of sampling must be 'in phase' with the S shape of the dispersion sensitivity curve to yield the most information on dispersion. (4) The sensitivity to the dispersion coefficient is usually at least an order of magnitude less than the sensitivity to velocity. (5) The assumed probability distribution of random error in observations of solute concentration determines the form of the sensitivities. (6) If variance in random error in observations is large, trends in sensitivities of observation points may be obscured by noise. (7) Designs that minimize the variance of one parameter may not necessarily minimize the variance of other parameters.
NASA Astrophysics Data System (ADS)
Pelosi, A.; Schumer, R.; Parker, G.; Ferguson, R. I.
2016-03-01
Tracer pebbles are often used to study bed load transport processes in gravel bed rivers. Models have been proposed for their downstream dispersion, and also for vertical dispersion, but not for the combined effects of downstream and vertical movement. Here we use the Exner-Based Master Equation to characterize the transient coevolution of streamwise and vertical advection-diffusion of tracer pebbles under equilibrium transport conditions (no net aggradation or degradation). The coevolution of streamwise and vertical dispersion gives rise to behavior that can differ markedly from that associated with purely streamwise processes with no vertical exchange. One example is streamwise advective slowdown. Particles that are advected downward into zones where the probability of reentrainment becomes asymptotically small are essentially trapped and can no longer participate in streamwise advection. As a result, the mean streamwise velocity of the tracer plume declines in time. Qualitative and quantitative comparisons with two field experiments show encouraging agreement despite the simplified boundary conditions in the model.
Healy, R.W.; Russell, T.F.
1992-01-01
A finite-volume Eulerian-Lagrangian local adjoint method for solution of the advection-dispersion equation is developed and discussed. The method is mass conservative and can solve advection-dominated ground-water solute-transport problems accurately and efficiently. An integrated finite-difference approach is used in the method. A key component of the method is that the integral representing the mass-storage term is evaluated numerically at the current time level. Integration points, and the mass associated with these points, are then forward tracked up to the next time level. The number of integration points required to reach a specified level of accuracy is problem dependent and increases as the sharpness of the simulated solute front increases. Integration points are generally equally spaced within each grid cell. For problems involving variable coefficients it has been found to be advantageous to include additional integration points at strategic locations in each well. These locations are determined by backtracking. Forward tracking of boundary fluxes by the method alleviates problems that are encountered in the backtracking approaches of most characteristic methods. A test problem is used to illustrate that the new method offers substantial advantages over other numerical methods for a wide range of problems.
Embry, Irucka; Roland, Victor; Agbaje, Oluropo; Watson, Valetta; Martin, Marquan; Painter, Roger; Byl, Tom; Sharpe, Lonnie
2013-01-01
A new residence-time distribution (RTD) function has been developed and applied to quantitative dye studies as an alternative to the traditional advection-dispersion equation (AdDE). The new method is based on a jointly combined four-parameter gamma probability density function (PDF). The gamma residence-time distribution (RTD) function and its first and second moments are derived from the individual two-parameter gamma distributions of randomly distributed variables, tracer travel distance, and linear velocity, which are based on their relationship with time. The gamma RTD function was used on a steady-state, nonideal system modeled as a plug-flow reactor (PFR) in the laboratory to validate themore » effectiveness of the model. The normalized forms of the gamma RTD and the advection-dispersion equation RTD were compared with the normalized tracer RTD. The normalized gamma RTD had a lower mean-absolute deviation (MAD) (0.16) than the normalized form of the advection-dispersion equation (0.26) when compared to the normalized tracer RTD. The gamma RTD function is tied back to the actual physical site due to its randomly distributed variables. The results validate using the gamma RTD as a suitable alternative to the advection-dispersion equation for quantitative tracer studies of non-ideal flow systems.« less
Technology Transfer Automated Retrieval System (TEKTRAN)
It has been reported that this model cannot take into account several important features of solute movement through soil. Recently, a new model has been suggested that results in a solute transport equation with fractional spatial derivatives, or FADE. We have assembled a database on published solu...
Technology Transfer Automated Retrieval System (TEKTRAN)
Understanding and modeling transport of solutes in porous media is a critical issue in the environmental protection. Contaminants from various industrial and agricultural sources can travel in soil and ground water and eventually affect human and animal health. The parabolic advective-dispersive equ...
Parashar, R.; Cushman, J.H.
2008-06-20
Microbial motility is often characterized by 'run and tumble' behavior which consists of bacteria making sequences of runs followed by tumbles (random changes in direction). As a superset of Brownian motion, Levy motion seems to describe such a motility pattern. The Eulerian (Fokker-Planck) equation describing these motions is similar to the classical advection-diffusion equation except that the order of highest derivative is fractional, {alpha} element of (0, 2]. The Lagrangian equation, driven by a Levy measure with drift, is stochastic and employed to numerically explore the dynamics of microbes in a flow cell with sticky boundaries. The Eulerian equation is used to non-dimensionalize parameters. The amount of sorbed time on the boundaries is modeled as a random variable that can vary over a wide range of values. Salient features of first passage time are studied with respect to scaled parameters.
Healy, R.W.; Russell, T.F.
1993-01-01
Test results demonstrate that the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) outperforms standard finite-difference methods for solute transport problems that are dominated by advection. FVELLAM systematically conserves mass globally with all types of boundary conditions. Integrated finite differences, instead of finite elements, are used to approximate the governing equation. This approach, in conjunction with a forward tracking scheme, greatly facilitates mass conservation. The mass storage integral is numerically evaluated at the current time level, and quadrature points are then tracked forward in time to the next level. Forward tracking permits straightforward treatment of inflow boundaries, thus avoiding the inherent problem in backtracking of characteristic lines intersecting inflow boundaries. FVELLAM extends previous results by obtaining mass conservation locally on Lagrangian space-time elements. -from Authors
Healy, R.W.; Russell, T.F.
1998-01-01
We extend the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) for solution of the advection-dispersion equation to two dimensions. The method can conserve mass globally and is not limited by restrictions on the size of the grid Peclet or Courant number. Therefore, it is well suited for solution of advection-dominated ground-water solute transport problems. In test problem comparisons with standard finite differences, FVELLAM is able to attain accurate solutions on much coarser space and time grids. On fine grids, the accuracy of the two methods is comparable. A critical aspect of FVELLAM (and all other ELLAMs) is evaluation of the mass storage integral from the preceding time level. In FVELLAM this may be accomplished with either a forward or backtracking approach. The forward tracking approach conserves mass globally and is the preferred approach. The backtracking approach is less computationally intensive, but not globally mass conservative. Boundary terms are systematically represented as integrals in space and time which are evaluated by a common integration scheme in conjunction with forward tracking through time. Unlike the one-dimensional case, local mass conservation cannot be guaranteed, so slight oscillations in concentration can develop, particularly in the vicinity of inflow or outflow boundaries. Published by Elsevier Science Ltd.
NASA Astrophysics Data System (ADS)
Parker, Jack C.; Kim, Ungtae
2015-11-01
The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of
Parker, Jack C; Kim, Ungtae
2015-11-01
The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of
Backward fractional advection dispersion model for contaminant source prediction
NASA Astrophysics Data System (ADS)
Zhang, Yong; Meerschaert, Mark M.; Neupauer, Roseanna M.
2016-04-01
The forward Fractional Advection Dispersion Equation (FADE) provides a useful model for non-Fickian transport in heterogeneous porous media. The space FADE captures the long leading tail, skewness, and fast spreading typically seen in concentration profiles from field data. This paper develops the corresponding backward FADE model, to identify source location and release time. The backward method is developed from the theory of inverse problems, and then explained from a stochastic point of view. The resultant backward FADE differs significantly from the traditional backward Advection Dispersion Equation (ADE) because the fractional derivative is not self-adjoint and the probability density function for backward locations is highly skewed. Finally, the method is validated using tracer data from a well-known field experiment, where the peak of the backward FADE curve predicts source release time, while the median or a range of percentiles can be used to determine the most likely source location for the observed plume. The backward ADE cannot reliably identify the source in this application, since the forward ADE does not provide an adequate fit to the concentration data.
A Lagrangian-Eulerian finite element method with adaptive gridding for advection-dispersion problems
Ijiri, Y.; Karasaki, K.
1994-02-01
In the present paper, a Lagrangian-Eulerian finite element method with adaptive gridding for solving advection-dispersion equations is described. The code creates new grid points in the vicinity of sharp fronts at every time step in order to reduce numerical dispersion. The code yields quite accurate solutions for a wide range of mesh Peclet numbers and for mesh Courant numbers well in excess of 1.
Garges, J.A.; Baehr, A.L.
1998-01-01
The relative importance of advection and dispersion for both solute and vapor transport can be determined from type curves or concentration, flux, or cumulative flux. The dimensionless form of the type curves provides a means to directly evaluate the importance of mass transport by advection relative to that of mass transport by diffusion and dispersion. Type curves based on an analytical solution to the advection-dispersion equation are plotted in terms of dimensionless time and Peclet number. Flux and cumulative flux type curves provide additional rationale for transport regime determination in addition to the traditional concentration type curves. The extension of type curves to include vapor transport with phase partitioning in the unsaturated zone is a new development. Type curves for negative Peclet numbers also are presented. A negative Peclet number characterizes a problem in which one direction of flow is toward the contamination source, and thereby diffusion and advection can act in opposite directions. Examples are the diffusion of solutes away from the downgradient edge of a pump-and-treat capture zone, the upward diffusion of vapors through the unsaturated zone with recharge, and the diffusion of solutes through a low hydraulic conductivity cutoff wall with an inward advective gradient.
NASA Astrophysics Data System (ADS)
Ancey, Christophe; Bohorquez, Patricio; Heyman, Joris
2016-04-01
The advection-diffusion equation arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Stochastic models can also be used to derive this equation, with the significant advantage that they provide information on the statistical properties of particle activity. Stochastic models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. We develop an approach based on birth-death Markov processes, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received little attention. We show that particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due to velocity fluctuations), with the important consequence that local measurements depend on both the intrinsic properties of particle displacement and the dimensions of the measurement system.
Lattice Boltzmann method for the fractional advection-diffusion equation
NASA Astrophysics Data System (ADS)
Zhou, J. G.; Haygarth, P. M.; Withers, P. J. A.; Macleod, C. J. A.; Falloon, P. D.; Beven, K. J.; Ockenden, M. C.; Forber, K. J.; Hollaway, M. J.; Evans, R.; Collins, A. L.; Hiscock, K. M.; Wearing, C.; Kahana, R.; Villamizar Velez, M. L.
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β , the fractional order α , and the single relaxation time τ , the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering.
Lattice Boltzmann method for the fractional advection-diffusion equation.
Zhou, J G; Haygarth, P M; Withers, P J A; Macleod, C J A; Falloon, P D; Beven, K J; Ockenden, M C; Forber, K J; Hollaway, M J; Evans, R; Collins, A L; Hiscock, K M; Wearing, C; Kahana, R; Villamizar Velez, M L
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β, the fractional order α, and the single relaxation time τ, the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering. PMID:27176431
Purely Lagrangian Simulation of Advection, Dispersion, Precipitation, and Dissolution
NASA Astrophysics Data System (ADS)
Benson, D.; Zhang, Y.; Reeves, D. M.
2008-05-01
We extend the advantages of Lagrangian random walk particle tracking (RWPT) methods that have long been used to simulate advection and dispersion in highly heterogeneous media. By formulating dissolution as a random, independent decay process, the classical continuum rate law is recovered. Formulating the random precipitation process requires a consideration of the probability that two nearby particles will coincide in a given time period. This depends on local mixing (as by diffusion) and the total domain particle number density, which are fixed and therefore easy to calculate. The result is that the classical law of mass action for equilibrium reactions can be reproduced in an ensemble sense. The same number of parameters for A+B ⇌ C are needed in a probabilistic versus continuum reaction simulation-- —one each for forward and backward probabilities that correspond to rates. The random nature of the simulations allows for significant disequilibrium in any given region at any time that is independent of the numerical details such as time stepping or particle density. This is exemplified by nearby or intermingled groups of reactants and little or no product--—a result that is often noted in the field that is difficult to reconcile with continuum methods or coarse-grained Eulerian models. Our results support recent results of perturbed advection-dispersion-reaction continuum models (Luo et al., WRR 44, 2008), and suggest that many different kinds of reactions can be easily added to existing RWPT codes.
Nikolaevskiy equation with dispersion.
Simbawa, Eman; Matthews, Paul C; Cox, Stephen M
2010-03-01
The Nikolaevskiy equation was originally proposed as a model for seismic waves and is also a model for a wide variety of systems incorporating a neutral "Goldstone" mode, including electroconvection and reaction-diffusion systems. It is known to exhibit chaotic dynamics at the onset of pattern formation, at least when the dispersive terms in the equation are suppressed, as is commonly the practice in previous analyses. In this paper, the effects of reinstating the dispersive terms are examined. It is shown that such terms can stabilize some of the spatially periodic traveling waves; this allows us to study the loss of stability and transition to chaos of the waves. The secondary stability diagram ("Busse balloon") for the traveling waves can be remarkably complicated. PMID:20365845
Multirate Runge-Kutta schemes for advection equations
NASA Astrophysics Data System (ADS)
Schlegel, Martin; Knoth, Oswald; Arnold, Martin; Wolke, Ralf
2009-04-01
Explicit time integration methods can be employed to simulate a broad spectrum of physical phenomena. The wide range of scales encountered lead to the problem that the fastest cell of the simulation dictates the global time step. Multirate time integration methods can be employed to alter the time step locally so that slower components take longer and fewer time steps, resulting in a moderate to substantial reduction of the computational cost, depending on the scenario to simulate [S. Osher, R. Sanders, Numerical approximations to nonlinear conservation laws with locally varying time and space grids, Math. Comput. 41 (1983) 321-336; H. Tang, G. Warnecke, A class of high resolution schemes for hyperbolic conservation laws and convection-diffusion equations with varying time and pace grids, SIAM J. Sci. Comput. 26 (4) (2005) 1415-1431; E. Constantinescu, A. Sandu, Multirate timestepping methods for hyperbolic conservation laws, SIAM J. Sci. Comput. 33 (3) (2007) 239-278]. In air pollution modeling the advection part is usually integrated explicitly in time, where the time step is constrained by a locally varying Courant-Friedrichs-Lewy (CFL) number. Multirate schemes are a useful tool to decouple different physical regions so that this constraint becomes a local instead of a global restriction. Therefore it is of major interest to apply multirate schemes to the advection equation. We introduce a generic recursive multirate Runge-Kutta scheme that can be easily adapted to an arbitrary number of refinement levels. It preserves the linear invariants of the system and is of third order accuracy when applied to certain explicit Runge-Kutta methods as base method.
NASA Astrophysics Data System (ADS)
Ancey, C.; Bohorquez, P.; Heyman, J.
2015-12-01
The advection-diffusion equation is one of the most widespread equations in physics. It arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Phenomenological laws are usually sufficient to derive this equation and interpret its terms. Stochastic models can also be used to derive it, with the significant advantage that they provide information on the statistical properties of particle activity. These models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. Among these stochastic models, the most common approach consists of random walk models. For instance, they have been used to model the random displacement of tracers in rivers. Here we explore an alternative approach, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. Birth-death Markov processes are well suited to this objective. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received no attention. We therefore look into the possibility of deriving the advection-diffusion equation (with a source term) within the framework of birth-death Markov processes. We show that in the continuum limit (when the cell size becomes vanishingly small), we can derive an advection-diffusion equation for particle activity. Yet while this derivation is formally valid in the continuum limit, it runs into difficulty in practical applications involving cells or meshes of finite length. Indeed, within our stochastic framework, particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due
Probabilistic exposure risk assessment with advective-dispersive well vulnerability criteria
NASA Astrophysics Data System (ADS)
Enzenhoefer, Rainer; Nowak, Wolfgang; Helmig, Rainer
2012-02-01
Time-related advection-based well-head protection zones are commonly used to manage the contamination risk of drinking water wells. According to current water safety plans advanced risk management schemes are needed to better control and monitor all possible hazards within catchments. The goal of this work is to cast the four advective-dispersive intrinsic well vulnerability criteria by Frind et al. [1] into a framework of probabilistic risk assessment framework. These criteria are: (i) arrival time, (ii) level of peak concentration, (iii) time until first arrival of critical concentrations and (iv) exposure time. Our probabilistic framework yields catchment-wide maps of probabilities to not comply with these criteria. This provides indispensable information for catchment managers to perform probabilistic exposure risk assessment and thus improves the basis for risk-informed well-head management. We resolve heterogeneity with high-resolution Monte Carlo simulations and use a new reverse formulation of temporal moment transport equations to keep computational costs low. Our method is independent of dimensionality and boundary conditions, and can account for arbitrary sources of uncertainty. It can be coupled with any method for conditioning on available data. For simplicity, we demonstrate the concept on a 2D example that includes conditioning on synthetic data.
Lewis, F.M.; Voss, C.I.; Rubin, Jacob
1986-01-01
A model was developed that can simulate the effect of certain chemical and sorption reactions simultaneously among solutes involved in advective-dispersive transport through porous media. The model is based on a methodology that utilizes physical-chemical relationships in the development of the basic solute mass-balance equations; however, the form of these equations allows their solution to be obtained by methods that do not depend on the chemical processes. The chemical environment is governed by the condition of local chemical equilibrium, and may be defined either by the linear sorption of a single species and two soluble complexation reactions which also involve that species, or binary ion exchange and one complexation reaction involving a common ion. Partial differential equations that describe solute mass balance entirely in the liquid phase are developed for each tenad (a chemical entity whose total mass is independent of the reaction process) in terms of their total dissolved concentration. These equations are solved numerically in two dimensions through the modification of an existing groundwater flow/transport computer code. (Author 's abstract)
NASA Astrophysics Data System (ADS)
Appadu, A. R.
2016-06-01
An unconditionally positive definite scheme has been derived in [1] to approximate a linear advection-diffusion-reaction equation which models exponential travelling waves and the coefficients of advective, diffusive and reactive terms have been chosen as one. The scheme has been baptised as Unconditionally Positive Finite Difference (UPFD). In this work, we use the UPFD scheme to solve the advection-diffusion-reaction problem in [1] and we also extend our study to three other important regimes involved in this model. The temporal step size is varied while fixing the spatial step size. We compute some errors namely; L1 error, dispersion, dissipation errors. We also study the variation of the modulus of the exact amplification factor, modulus of amplification factor of the scheme and relative phase error, all vs the phase angle for the four different regimes.
General solution of a fractional diffusion-advection equation for solar cosmic-ray transport
NASA Astrophysics Data System (ADS)
Rocca, M. C.; Plastino, A. R.; Plastino, A.; Ferri, G. L.; de Paoli, A.
2016-04-01
In this effort we exactly solve the fractional diffusion-advection equation for solar cosmic-ray transport and give its general solution in terms of hypergeometric distributions. Numerical analysis of this equation shows that its solutions resemble power-laws.
NASA Astrophysics Data System (ADS)
Sofiev, M.; Vira, J.; Kouznetsov, R.; Prank, M.; Soares, J.; Genikhovich, E.
2015-11-01
The paper presents the transport module of the System for Integrated modeLling of Atmospheric coMposition SILAM v.5 based on the advection algorithm of Michael Galperin. This advection routine, so far weakly presented in the international literature, is positively defined, stable at any Courant number, and efficient computationally. We present the rigorous description of its original version, along with several updates that improve its monotonicity and shape preservation, allowing for applications to long-living species in conditions of complex atmospheric flows. The scheme is connected with other parts of the model in a way that preserves the sub-grid mass distribution information that is a cornerstone of the advection algorithm. The other parts include the previously developed vertical diffusion algorithm combined with dry deposition, a meteorological pre-processor, and chemical transformation modules. The quality of the advection routine is evaluated using a large set of tests. The original approach has been previously compared with several classic algorithms widely used in operational dispersion models. The basic tests were repeated for the updated scheme and extended with real-wind simulations and demanding global 2-D tests recently suggested in the literature, which allowed one to position the scheme with regard to sophisticated state-of-the-art approaches. The advection scheme performance was fully comparable with other algorithms, with a modest computational cost. This work was the last project of Dr. Sci. Michael Galperin, who passed away on 18 March 2008.
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents a formal exact solution of the linear advection-diffusion transport equation with constant coefficients for both transient and steady-state regimes. A classical mathematical substitution transforms the original advection-diffusion equation into an exclusively diffusive equation. ...
A global spectral element model for poisson equations and advective flow over a sphere
NASA Astrophysics Data System (ADS)
Mei, Huan; Wang, Faming; Zeng, Zhong; Qiu, Zhouhua; Yin, Linmao; Li, Liang
2016-03-01
A global spherical Fourier-Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre-Gauss-Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity, with only one element. Then, the partial differential equations are solved on the longitude-latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. Highprecision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.
NASA Astrophysics Data System (ADS)
Bernabé, Y.; Wang, Y.; Qi, T.; Li, M.
2016-02-01
The main purpose of this work is to investigate the relationship between passive advection-dispersion and permeability in porous materials presumed to be statistically homogeneous at scales larger than the pore scale but smaller than the reservoir scale. We simulated fluid flow through pipe network realizations with different pipe radius distributions and different levels of connectivity. The flow simulations used periodic boundary conditions, allowing monitoring of the advective motion of solute particles in a large periodic array of identical network realizations. In order to simulate dispersion, we assumed that the solute particles obeyed Taylor dispersion in individual pipes. When a particle entered a pipe, a residence time consistent with local Taylor dispersion was randomly assigned to it. When exiting the pipe, the particle randomly proceeded into one of the pipes connected to the original one according to probabilities proportional to the outgoing volumetric flow in each pipe. For each simulation we tracked the motion of at least 6000 solute particles. The mean fluid velocity was 10-3 ms-1, and the distance traveled was on the order of 10 m. Macroscopic dispersion was quantified using the method of moments. Despite differences arising from using different types of lattices (simple cubic, body-centered cubic, and face-centered cubic), a number of general observations were made. Longitudinal dispersion was at least 1 order of magnitude greater than transverse dispersion, and both strongly increased with decreasing pore connectivity and/or pore size variability. In conditions of variable hydraulic radius and fixed pore connectivity and pore size variability, the simulated dispersivities increased as power laws of the hydraulic radius and, consequently, of permeability, in agreement with previously published experimental results. Based on these observations, we were able to resolve some of the complexity of the relationship between dispersivity and permeability.
It is well known that the fate and transport of contaminants in the subsurface are controlled by complex processes including advection, dispersion-diffusion, and chemical reactions. However, the interplay between the physical transport processes and chemical reactions, and their...
Technology Transfer Automated Retrieval System (TEKTRAN)
The convective-dispersive, or advective-dispersive, or CDE, equation has long been the model of choice for solute transport in soils. Using the total mass of soluble salts in soil profile to evaluate changes in salinity due to irrigation can be beneficial when the spatial variability of soil salini...
CDF Solutions of Advection-Reaction equations with uncertain parameters (Invited)
NASA Astrophysics Data System (ADS)
Boso, F.; Tartakovsky, D. M.
2013-12-01
Flow and transport models are affected by parametric uncertainty. Quantitative forecasting of such processes in natural porous media are especially prone to uncertainty because of the inaccessibility and multi-scale nature of the subsurface. We consider a reduced-complexity stochastic transport system which takes into account advection and nonlinear reactions in advection-reaction equations (AREs) with uncertain (random) velocity and reaction parameters. We derive a deterministic equation that governs the evolution of cumulative distribution function (CDF) of a solution of the underlying ARE. Although requiring closure, this differential equation benefits from uniquely defined boundary and initial conditions and can be solved with classic techniques. Here we analyze the accuracy and robustness of the large-eddy-diffusivity closure by comparison with Monte Carlo simulations for different correlation structures and parameters.
NASA Astrophysics Data System (ADS)
Cartwright, Ian
Advection-dispersion fluid flow models implicitly assume that the infiltrating fluid flows through an already fluid-saturated medium. However, whether rocks contain a fluid depends on their reaction history, and whether any initial fluid escapes. The behaviour of different rocks may be illustrated using hypothetical marble compositions. Marbles with diverse chemistries (e.g. calcite + dolomite + quartz) are relatively reactive, and will generally produce a fluid during heating. By contrast, marbles with more restricted chemistries (e.g. calcite + quartz or calcite-only) may not. If the rock is not fluid bearing when fluid infiltration commences, mineralogical reactions may produce a reaction-enhanced permeability in calcite + dolomite + quartz or calcite + quartz, but not in calcite-only marbles. The permeability production controls the pattern of mineralogical, isotopic, and geochemical resetting during fluid flow. Tracers retarded behind the mineralogical fronts will probably be reset as predicted by the advection-dispersion models; however, tracers that are expected to be reset ahead of the mineralogical fronts cannot progress beyond the permeability generating reaction. In the case of very unreactive lithologies (e.g. pure calcite marbles, cherts, and quartzites), the first reaction to affect the rocks may be a metasomatic one ahead of which there is little pervasive resetting of any tracer. Centimetre-scale layering may lead to the formation of self-perpetuating fluid channels in rocks that are not fluid saturated due to the juxtaposition of reactants. Such layered rocks may show patterns of mineralogical resetting that are not predicted by advection-dispersion models. Patterns of mineralogical and isotopic resetting in marbles from a number of terrains, for example: Chillagoe, Marulan South, Reynolds Range (Australia); Adirondack Mountains, Old Woman Mountains, Notch Peak (USA); and Stephen Cross Quarry (Canada) vary as predicted by these models.
Cauchy's dispersion equation reconsidered : dispersion in silicate glasses.
Smith, D. Y.; Inokuti, M.; Karstens, W.; Physics; Univ. of Vermont; St. Michael's College
2002-01-01
We formulate a novel method of characterizing optically transparent substances using dispersion theory. The refractive index is given by a generalized Cauchy dispersion equation with coefficients that are moments of the uv and ir absorptions. Mean dispersion, Abbe number, and partial dispersion are combinations of these moments. The empirical relation between index and dispersion for families of glasses appears as a consequence of Beer's law applied to the uv spectra.
Lichtner, P.C.; Helgeson, H.C.
1986-06-20
A general formulation of multi-phase fluid flow coupled to chemical reactions was developed based on a continuum description of porous media. A preliminary version of the computer code MCCTM was constructed which implemented the general equations for a single phase fluid. The computer code MCCTM incorporates mass transport by advection-diffusion/dispersion in a one-dimensional porous medium coupled to reversible and irreversible, homogeneous and heterogeneous chemical reactions. These reactions include aqueous complexing, oxidation/reduction reactions, ion exchange, and hydrolysis reactions of stoichiometric minerals. The code MCCTM uses a fully implicit finite difference algorithm. The code was tested against analytical calculations. Applications of the code included investigation of the propagation of sharp chemical reaction fronts, metasomatic alteration of microcline at elevated temperatures and pressures, and ion-exchange in a porous column. Finally numerical calculations describing fluid flow in crystalline rock in the presence of a temperature gradient were compared with experimental results for quartzite.
Dynamic typology of hydrothermal systems: competing effects of advection, dispersion and reactivity
NASA Astrophysics Data System (ADS)
Dolejs, David
2016-04-01
Genetic interpretation hydrothermal systems relies on recognition of (i) hydrothermal fluid source, (ii) fluid migration pathways, and (iii) deposition site identified by hydrothermal alteration and/or mineralization. Frequently, only the last object is of interest or accessible to direct observation, but constraints on the fluid source (volume) and pathways can be obtained from evaluation of the time-integrated fluid flux during hydrothermal event. Successful interpretation of the petrological record, that is, progress of alteration reactions, relies on identification of individual contributions arising from solute advection (to the deposition site), its lateral dispersion, and reaction efficiency. Although these terms are all applicable in a mass-conservation relationship within the framework of the transport theory, they are rarely considered simultaneously and their relative magnitudes evaluated. These phenomena operate on variable length and time scales, and may in turn provide insight into the system dynamics such as flow, diffusion and reaction rates, or continuous vs. episodic behavior of hydrothermal events. In addition, here we demonstrate that they also affect estimate of the net fluid flux, frequently by several orders of magnitude. The extent of alteration and mineralization reactions between the hydrothermal fluid and the host environment is determined by: (i) temperature, pressure or any other gradients across the mineralization site, (ii) magnitude of disequilibrium at inflow to the mineralization site, which is related to physico-chemical gradient between the fluid source and the mineralization site, and (iii) chemical redistribution (dispersion) within the mineralization site. We introduce quantitative mass-transport descriptors - Péclet and Damköhler II numbers - to introduce division into dispersion-dominated, advection-dominated and reaction-constrained systems. Dispersive systems are characterized by lateral solute redistribution, driven by
NASA Astrophysics Data System (ADS)
Cornaton, F.; Perrochet, P.
2006-09-01
Groundwater age and life expectancy probability density functions (pdf) have been defined, and solved in a general three-dimensional context by means of forward and backward advection-dispersion equations [Cornaton F, Perrochet P. Groundwater age, life expectancy and transit time distributions in advective-dispersive systems; 1. Generalized reservoir theory. Adv Water Res (xxxx)]. The discharge and recharge zones transit time pdfs were then derived by applying the reservoir theory (RT) to the global system, thus considering as ensemble the union of all inlet boundaries on one hand, and the union of all outlet boundaries on the other hand. The main advantages in using the RT to calculate the transit time pdf is that the outlet boundary geometry does not represent a computational limiting factor (e.g. outlets of small sizes), since the methodology is based on the integration over the entire domain of each age, or life expectancy, occurrence. In the present paper, we extend the applicability of the RT to sub-drainage basins of groundwater reservoirs by treating the reservoir flow systems as compartments which transfer the water fluxes to a particular discharge zone, and inside which mixing and dispersion processes can take place. Drainage basins are defined by the field of probability of exit at outlet. In this way, we make the RT applicable to each sub-drainage system of an aquifer of arbitrary complexity and configuration. The case of the well-head protection problem is taken as illustrative example, and sensitivity analysis of the effect of pore velocity variations on the simulated ages is carried out.
Cox, T.J.; Runkel, R.L.
2008-01-01
Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme's limitations. ?? 2008 ASCE.
Huang, Y.H.; Saiers, J.E.; Harvey, J.W.; Noe, G.B.; Mylon, S.
2008-01-01
The movement of particulate matter within wetland surface waters affects nutrient cycling, contaminant mobility, and the evolution of the wetland landscape. Despite the importance of particle transport in influencing wetland form and function, there are few data sets that illuminate, in a quantitative way, the transport behavior of particulate matter within surface waters containing emergent vegetation. We report observations from experiments on the transport of 1 ??m latex microspheres at a wetland field site located in Water Conservation Area 3A of the Florida Everglades. The experiments involved line source injections of particles inside two 4.8-m-long surface water flumes constructed within a transition zone between an Eleocharis slough and Cladium jamaicense ridge and within a Cladium jamaicense ridge. We compared the measurements of particle transport to calculations of two-dimensional advection-dispersion model that accounted for a linear increase in water velocities with elevation above the ground surface. The results of this analysis revealed that particle spreading by longitudinal and vertical dispersion was substantially greater in the ridge than within the transition zone and that particle capture by aquatic vegetation lowered surface water particle concentrations and, at least for the timescale of our experiments, could be represented as an irreversible, first-order kinetics process. We found generally good agreement between our field-based estimates of particle dispersion and water velocity and estimates determined from published theory, suggesting that the advective-dispersive transport of particulate matter within complex wetland environments can be approximated on the basis of measurable properties of the flow and aquatic vegetation. Copyright 2008 by the American Geophysical Union.
NASA Astrophysics Data System (ADS)
Huber, Markus; Tailleux, Remi; Ferreira, David; Kuhlbrodt, Till; Gregory, Jonathan
2015-04-01
The classic vertical advection-diffusion (VAD) balance is a central concept in studying the ocean heat budget, in particular in simple climate models (SCMs). Here we present a new framework to calibrate the parameters of the VAD equation to the vertical ocean heat balance of two fully-coupled climate models that is traceable to the models' circulation as well as to vertical mixing and diffusion processes. Based on temperature diagnostics, we derive an effective vertical velocity w∗ and turbulent diffusivity kν∗ for each individual physical process. In steady state, we find that the residual vertical velocity and diffusivity change sign in middepth, highlighting the different regional contributions of isopycnal and diapycnal diffusion in balancing the models' residual advection and vertical mixing. We quantify the impacts of the time evolution of the effective quantities under a transient 1% CO2 simulation and make the link to the parameters of currently employed SCMs.
NASA Astrophysics Data System (ADS)
Jiang, Tian; Zhang, Yong-Tao
2016-04-01
Implicit integration factor (IIF) methods were developed in the literature for solving time-dependent stiff partial differential equations (PDEs). Recently, IIF methods were combined with weighted essentially non-oscillatory (WENO) schemes in Jiang and Zhang (2013) [19] to efficiently solve stiff nonlinear advection-diffusion-reaction equations. The methods can be designed for arbitrary order of accuracy. The stiffness of the system is resolved well and the methods are stable by using time step sizes which are just determined by the non-stiff hyperbolic part of the system. To efficiently calculate large matrix exponentials, Krylov subspace approximation is directly applied to the implicit integration factor (IIF) methods. So far, the IIF methods developed in the literature are multistep methods. In this paper, we develop Krylov single-step IIF-WENO methods for solving stiff advection-diffusion-reaction equations. The methods are designed carefully to avoid generating positive exponentials in the matrix exponentials, which is necessary for the stability of the schemes. We analyze the stability and truncation errors of the single-step IIF schemes. Numerical examples of both scalar equations and systems are shown to demonstrate the accuracy, efficiency and robustness of the new methods.
Looney, B.B.; Scott, M.T.
1988-12-31
Recent field and laboratory data have confirmed that apparent dispersivity is a function of the flow distance of the measurement. This scale effect is not consistent with classical advection dispersion modeling often used to describe the transport of solutes in saturated porous media. Many investigators attribute this anomalous behavior to the fact that the spreading of solute is actually the result of the heterogeneity of subsurface materials and the wide distribution of flow paths and velocities available in such systems. An analysis using straightforward analytical equations confirms this hypothesis. An analytical equation based on a flow variance approach matches available field data when a variance description of approximately 0.4 is employed. Also, current field data provide a basis for statistical selection of the variance parameter based on the level of concern related to the resulting calculated concentration. While the advection dispersion approach often yielded reasonable predictions, continued development of statistical and stochastic techniques will provide more defendable and mechanistically descriptive models.
NASA Technical Reports Server (NTRS)
Leonard, B. P.
1988-01-01
A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.
Barth, Andrea Lang, Annika
2012-12-15
In this paper, the strong approximation of a stochastic partial differential equation, whose differential operator is of advection-diffusion type and which is driven by a multiplicative, infinite dimensional, cadlag, square integrable martingale, is presented. A finite dimensional projection of the infinite dimensional equation, for example a Galerkin projection, with nonequidistant time stepping is used. Error estimates for the discretized equation are derived in L{sup 2} and almost sure senses. Besides space and time discretizations, noise approximations are also provided, where the Milstein double stochastic integral is approximated in such a way that the overall complexity is not increased compared to an Euler-Maruyama approximation. Finally, simulations complete the paper.
Exact PDF equations and closure approximations for advective-reactive transport
Venturi, D.; Tartakovsky, Daniel M.; Tartakovsky, Alexandre M.; Karniadakis, George E.
2013-06-01
Mathematical models of advection–reaction phenomena rely on advective flow velocity and (bio) chemical reaction rates that are notoriously random. By using functional integral methods, we derive exact evolution equations for the probability density function (PDF) of the state variables of the advection–reaction system in the presence of random transport velocity and random reaction rates with rather arbitrary distributions. These PDF equations are solved analytically for transport with deterministic flow velocity and a linear reaction rate represented mathematically by a heterog eneous and strongly-correlated random field. Our analytical solution is then used to investigate the accuracy and robustness of the recently proposed large-eddy diffusivity (LED) closure approximation [1]. We find that the solution to the LED-based PDF equation, which is exact for uncorrelated reaction rates, is accurate even in the presence of strong correlations and it provides an upper bound of predictive uncertainty.
Variational integration for ideal magnetohydrodynamics with built-in advection equations
Zhou, Yao; Burby, J. W.; Bhattacharjee, A.; Qin, Hong
2014-10-15
Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.
Preconditioned iterative methods for space-time fractional advection-diffusion equations
NASA Astrophysics Data System (ADS)
Zhao, Zhi; Jin, Xiao-Qing; Lin, Matthew M.
2016-08-01
In this paper, we propose practical numerical methods for solving a class of initial-boundary value problems of space-time fractional advection-diffusion equations. First, we propose an implicit method based on two-sided Grünwald formulae and discuss its stability and consistency. Then, we develop the preconditioned generalized minimal residual (preconditioned GMRES) method and preconditioned conjugate gradient normal residual (preconditioned CGNR) method with easily constructed preconditioners. Importantly, because resulting systems are Toeplitz-like, fast Fourier transform can be applied to significantly reduce the computational cost. We perform numerical experiments to demonstrate the efficiency of our preconditioners, even in cases with variable coefficients.
Variational Integration for Ideal MHD with Built-in Advection Equations
Zhou, Yao; Qin, Hong; Burby, J. W.; Bhattacharjee, A.
2014-08-05
Newcomb's Lagrangian for ideal MHD in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.
Horizontal advection and dispersion in a stratified shelf sea: The role of inertial oscillations
NASA Astrophysics Data System (ADS)
Inall, Mark E.; Aleynik, Dmitry; Neil, Clare
2013-10-01
The role played by inertial motions in horizontal dispersion within the thermocline of a broad, mid-latitude shelf sea is examined through the analysis of a deliberately released dye tracer. Our analysis is of the horizontal and vertical evolution over 40 h of a dye tracer injected into the seasonally stratified thermocline of the Celtic Sea on the NW European Shelf. The inferred diapycnal diffusivity was 1.3-1.5 × 10-5 m2 s-1, and the radial horizontal diffusivities of the depth integrated dye patch ranged from 1.9 to 4.0 m2 s-1. The inferred vertical diffusivity is in agreement with microstructure based estimates, and the depth integrated horizontal diffusivity is broadly in agreement with previous dye release derived estimates made over similar scales and time periods. Asymmetry in the horizontal evolution of the dye patch was evident. We argue that mean shear dispersion was responsible for lateral elongation of the dye patch, particularly between hours 23 and 35 after release, during which time horizontal diffusivity along the major axis, Ka, exceeded that along the minor axis, Kb, by more than a factor of 10. We further show that along-patch shear was predominantly a result of differential advection between a deep residual flow to the south-east and an oscillating wind-driven surface Ekman layer. In this region of strong low frequency (inertial) shear a time dependent model of shear dispersion (Young et al., 1982) was able to account for the observed rate of horizontal dispersion calculated on the target isopycnal surface.
Application of a Particle Method to the Advection-Diffusion-Reaction Equation
NASA Astrophysics Data System (ADS)
Paster, A.; Bolster, D.; Benson, D. A.
2012-12-01
A reaction between two chemical species can only happen if molecules collide and react. Thus, the mixing of a system can become a limiting factor in the onset of reaction. Solving for reaction rate in a well-mixed system is typically a straightforward task. However, when incomplete mixing kicks in, obtaining a solution becomes more challenging. Since reaction can only happen in regions where both reactants co-exist, the incomplete mixing may slow down the reaction rate, when compared to a well-mixed system. The effect of incomplete mixing upon reaction is a highly important aspect of various processes in natural and engineered systems, ranging from mineral precipitation in geological formations to groundwater remediation in aquifers. We study a relatively simple system with a bi-molecular irreversible kinetic reaction A+B → Ø where the underlying transport of reactants is governed by an advection-diffusion equation, and the initial concentrations are given in terms of an average and a perturbation. Such a system does not have an analytical solution to date, even for the zero advection case. We model the system by a Monte Carlo particle tracking method, where particles represent some reactant mass. In this method, diffusion is modeled by a random walk of the particles, and reaction is modeled by annihilation of particles. The probability of the annihilation is proportional to the reaction rate constant and the probability density associated with particle co-location. We study the numerical method in depth, characterizing typical numerical errors and time step restrictions. In particular, we show that the numerical method converges to the advection-diffusion-reaction equation at the limit Δt →0. We also rigorously derive the relationship between the initial number of particles in the system and the initial concentrations perturbations represented by that number. We then use the particle simulations of zero-advection system to demonstrate the well
Xu, Bruce S; Lollar, Barbara Sherwood; Passeport, Elodie; Sleep, Brent E
2016-04-15
Aqueous phase diffusion-related isotope fractionation (DRIF) for carbon isotopes was investigated for common groundwater contaminants in systems in which transport could be considered to be one-dimensional. This paper focuses not only on theoretically observable DRIF effects in these systems but introduces the important concept of constraining "observable" DRIF based on constraints imposed by the scale of measurements in the field, and on standard limits of detection and analytical uncertainty. Specifically, constraints for the detection of DRIF were determined in terms of the diffusive fractionation factor, the initial concentration of contaminants (C0), the method detection limit (MDL) for isotopic analysis, the transport time, and the ratio of the longitudinal mechanical dispersion coefficient to effective molecular diffusion coefficient (Dmech/Deff). The results allow a determination of field conditions under which DRIF may be an important factor in the use of stable carbon isotope measurements for evaluation of contaminant transport and transformation for one-dimensional advective-dispersive transport. This study demonstrates that for diffusion-dominated transport of BTEX, MTBE, and chlorinated ethenes, DRIF effects are only detectable for the smaller molar mass compounds such as vinyl chloride for C0/MDL ratios of 50 or higher. Much larger C0/MDL ratios, corresponding to higher source concentrations or lower detection limits, are necessary for DRIF to be detectable for the higher molar mass compounds. The distance over which DRIF is observable for VC is small (less than 1m) for a relatively young diffusive plume (<100years), and DRIF will not easily be detected by using the conventional sampling approach with "typical" well spacing (at least several meters). With contaminant transport by advection, mechanical dispersion, and molecular diffusion this study suggests that in field sites where Dmech/Deff is larger than 10, DRIF effects will likely not be
Kordilla, Jannes; Pan, Wenxiao; Tartakovsky, Alexandre M.
2014-12-14
We propose a novel Smoothed Particle Hydrodynamics (SPH) discretization of the fully-coupled Landau-Lifshitz-Navier-Stokes (LLNS) and advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations are found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for the coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study the formation of the so-called giant fluctuations of the front between light and heavy fluids with and without gravity, where the light fluid lays on the top of the heavy fluid. We find that the power spectra of the simulated concentration field is in good agreement with the experiments and analytical solutions. In the absence of gravity the the power spectra decays as the power -4 of the wave number except for small wave numbers which diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations resulting in the much weaker dependence of the power spectra on the wave number. Finally the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlying a light fluid. The front dynamics is shown to agree well with the analytical solutions.
NASA Astrophysics Data System (ADS)
Kordilla, Jannes; Pan, Wenxiao; Tartakovsky, Alexandre
2014-12-01
We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called "giant fluctuations" of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power -4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.
Kordilla, Jannes; Pan, Wenxiao Tartakovsky, Alexandre
2014-12-14
We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.
Analytical solutions of a fractional diffusion-advection equation for solar cosmic-ray transport
Litvinenko, Yuri E.; Effenberger, Frederic
2014-12-01
Motivated by recent applications of superdiffusive transport models to shock-accelerated particle distributions in the heliosphere, we analytically solve a one-dimensional fractional diffusion-advection equation for the particle density. We derive an exact Fourier transform solution, simplify it in a weak diffusion approximation, and compare the new solution with previously available analytical results and with a semi-numerical solution based on a Fourier series expansion. We apply the results to the problem of describing the transport of energetic particles, accelerated at a traveling heliospheric shock. Our analysis shows that significant errors may result from assuming an infinite initial distance between the shock and the observer. We argue that the shock travel time should be a parameter of a realistic superdiffusive transport model.
Analytical Solutions of a Fractional Diffusion-advection Equation for Solar Cosmic-Ray Transport
NASA Astrophysics Data System (ADS)
Litvinenko, Yuri E.; Effenberger, Frederic
2014-12-01
Motivated by recent applications of superdiffusive transport models to shock-accelerated particle distributions in the heliosphere, we analytically solve a one-dimensional fractional diffusion-advection equation for the particle density. We derive an exact Fourier transform solution, simplify it in a weak diffusion approximation, and compare the new solution with previously available analytical results and with a semi-numerical solution based on a Fourier series expansion. We apply the results to the problem of describing the transport of energetic particles, accelerated at a traveling heliospheric shock. Our analysis shows that significant errors may result from assuming an infinite initial distance between the shock and the observer. We argue that the shock travel time should be a parameter of a realistic superdiffusive transport model.
Preconditioned time-difference methods for advection-diffusion-reaction equations
Aro, C.; Rodrigue, G.; Wolitzer, D.
1994-12-31
Explicit time differencing methods for solving differential equations are advantageous in that they are easy to implement on a computer and are intrinsically very parallel. The disadvantage of explicit methods is the severe restrictions placed on stepsize due to stability. Stability bounds for explicit time differencing methods on advection-diffusion-reaction problems are generally quite severe and implicit methods are used instead. The linear systems arising from these implicit methods are large and sparse so that iterative methods must be used to solve them. In this paper the authors develop a methodology for increasing the stability bounds of standard explicit finite differencing methods by combining explicit methods, implicit methods, and iterative methods in a novel way to generate new time-difference schemes, called preconditioned time-difference methods.
Dispersion in tidally averaged transport equation
Cheng, R.T.; Casulli, V.
1992-01-01
A general governing inter-tidal transport equation for conservative solutes has been derived without invoking the weakly nonlinear approximation. The governing inter-tidal transport equation is a convection-dispersion equation in which the convective velocity is a mean Lagrangian residual current, and the inter-tidal dispersion coefficient is defined by a dispersion patch. When the weakly nonlinear condition is violated, the physical significance of the Stokes' drift, as used in tidal dynamics, becomes questionable. For nonlinear problems, analytical solutions for the mean Lagrangian residual current and for the inter-tidal dispersion coefficient do not exist, they must be determined numerically. A rectangular tidal inlet with a constriction is used in the first example. The solutions of the residual currents and the computed properties of the inter-tidal dispersion coefficient are used to illuminate the mechanisms of the inter-tidal transport processes. Then, the present formulation is tested in a geometrically complex tidal estuary – San Francisco Bay, California. The computed inter-tidal dispersion coefficients are in the range between 5×104 and 5×106 cm2/sec., which are consistent with the values reported in the literature
NASA Astrophysics Data System (ADS)
Moiseev, N. Ya.; Silant'eva, I. Yu.
2008-07-01
An approach to the construction of second-and higher order accurate difference schemes in time and space is described for solving the linear one-and multidimensional advection equations with constant coefficients by the Godunov method with antidiffusion. The differential approximations for schemes of up to the fifth order are constructed and written. For multidimensional advection equations with constant coefficients, it is shown that Godunov schemes with splitting over spatial variables are preferable, since they have a smaller truncation error than schemes without splitting. The high resolution and efficiency of the difference schemes are demonstrated using test computations.
NASA Astrophysics Data System (ADS)
Bouchelaghem, F.; Vulliet, L.
2001-10-01
The development of a predictive model of behaviour of porous media during injection of miscible grout, taking into account convection, dilution and filtration of grout solution with interstitial water, as well as consolidation aspects, is presented. Model assumptions are reviewed and discussed first. During the establishment of the model, we insist on surface terms and their physical relevance in expressing adsorption effects. Constitutive laws such as Fick's law for diffusive mass transport, hydrodynamic dispersion tensor dealing with miscibility, are modified by taking into account filtration effects. A new surface term appears in mass balance equations as a consequence of filtration. According to the filtration laws used, an initial filtration rate is estimated on the basis of a one-dimensional experimental campaign. The field equations are discretized by using Galerkin finite element and -scheme standard method. For transport equation, Streamline Upwind Petrov Galerkin method is employed to prevent numerical oscillations. Lastly, confrontation of numerical results with laboratory experiments constitutes a first step to validate the model on a realistic basis.
Brine heterogeneity and dispersed interstitial advective flow underneath the sea of galilee, Israel
NASA Astrophysics Data System (ADS)
Weinstein, Y.; Katz, A.; Kastner, M.; Nishri, A.; Jannasch, H.
2003-04-01
Saline groundwater from submerged sources is today the main source of salts to the Sea of Galilee, supplying annually 72,000 tons chloride to the lake. MOSQUITO flux-meters were deployed during April to September 2001 at seven shallow Kinneret sites in order to study the dispersed interstitial flow. Each instrument carried 3-5 Osmo-Samplers that continuously sampled pore fluids at various depths in the sediment (0-45 cm). Samples were analyzed for their chemistry and for concentration of Na-fluorescein that was previously injected into the sediment. In general, oncentrations of conservative elements (e.g. Cl, Na, B) increase with depth into the sediment at all sites. However, concentrations vary significantly from one site to another. For instance, in sub-lacustrine brines next to the Tiberias Hot Springs Cl concentration reaches more than 14,500 mg/l at 35 cm below lake floor (not much less than the 18,000 mg/l in the nearby on-shore springs), while at other stations, brines are diluted by meteoric water to less than 1,500 mg/l Cl. Ion ratios in pore water indicate that the shallow parts of Lake Kinneret are underlain by several, separate, brine pockets, that are sometimes located very close to each other and discharge to the same area. Pore water at the east and northwest of the lake have Na/Cl ionic ratios between 0.7 and 0.8, similar to that of the overlying lake water, while at the west and south, ratios are significantly lower (<0.6 and <0.5, respectively), indicating larger degree of evaporation of the original brine end-member. Brines next to Tiberias Hot Springs have significantly higher Br/Cl and lower Mg/Cl ratios than pore water from other sites. Eastern shore sub-lacusrine brines (next to Gofra) are distinguished by their very high Sr/Cl and B/Cl ratios. Advective flow rates were derived from temporal patterns of Na-fluorescein concentration in pore water. Flow rates were between 0 and 80 cm/yr. Annual fluxes through the shallow part of the lake are
AN EULERIAN-LAGRANGIAN LOCALIZED ADJOINT METHOD FOR THE ADVECTION-DIFFUSION EQUATION
Many numerical methods use characteristic analysis to accommodate the advective component of transport. Such characteristic methods include Eulerian-Lagrangian methods (ELM), modified method of characteristics (MMOC), and operator splitting methods. A generalization of characteri...
EULERIAN-LAGRANGIAN LOCALIZED ADJOINT METHOD FOR THE ADVECTION-DIFFUSION EQUATION
Many numerical methods use characteristic analysis to accommodate the advective component of transport. uch characteristic methods include Eulerian-Lagrangian methods (ELM), modified method of characteristics (MMOC), and operator splitting methods. eneralization of characteristic...
Downes, Barbara J; Lancaster, Jill
2010-01-01
1. In advection-dominated systems (both freshwater and marine), population dynamics are usually presumed to be dominated by the effects of migrants dispersing by advection, especially over the small spatial scales at which populations can be studied, but few studies have tested this presumption. We tested the hypothesis that benthic densities are controlled by densities of dispersers for two aquatic insects in upland streams. 2. Our study animals were two species of caddisflies (Hydropsychidae), which become sedentary filter-feeders following settlement onto substrata. Densities of dispersers in the drift (advective dispersal) were quantified using nets placed along the upstream edges of riffles, where the latter abruptly abutted a slower, upstream run. Settlement was estimated at each site using brick pavers, half of which had been fenced to prevent colonization of their top surfaces by walking hydropsychids, thus allowing us to distinguish also the mode of movement during settlement. 3. First through fifth instars of two species, Smicrophylax sp. AV2 and Asmicridea sp. AV1, were abundant and showed disparate results. Drift and settlement were relatively strongly related for Smicrophylax. The best fit lines were shown by second and third instars settling on plain bricks, suggesting that drift played a strong role in settlement, but that some drifters dropped to the bottom and located substrata by walking. Quantile regression suggested that drift sets limits to settlement in this species and that settlement success was highly variable. In contrast, settlement by Asmicridea was poorly related to drift; settlers were mainly individuals re-dispersing within sites. 4. Smicrophylax densities appear to be controlled by dispersal from upstream, but benthic density of Asmicridea is more likely linked to local demography. Our data demonstrate the dangers of assuming that supposedly drift-prone species can all be modelled in the same way. Alternative models emphasizing
NASA Astrophysics Data System (ADS)
Rubbab, Qammar; Mirza, Itrat Abbas; Qureshi, M. Zubair Akbar
2016-07-01
The time-fractional advection-diffusion equation with Caputo-Fabrizio fractional derivatives (fractional derivatives without singular kernel) is considered under the time-dependent emissions on the boundary and the first order chemical reaction. The non-dimensional problem is formulated by using suitable dimensionless variables and the fundamental solutions to the Dirichlet problem for the fractional advection-diffusion equation are determined using the integral transforms technique. The fundamental solutions for the ordinary advection-diffusion equation, fractional and ordinary diffusion equation are obtained as limiting cases of the previous model. Using Duhamel's principle, the analytical solutions to the Dirichlet problem with time-dependent boundary pulses have been obtained. The influence of the fractional parameter and of the drift parameter on the solute concentration in various spatial positions was analyzed by numerical calculations. It is found that the variation of the fractional parameter has a significant effect on the solute concentration, namely, the memory effects lead to the retardation of the mass transport.
A dispersive wave equation using nonlocal elasticity
NASA Astrophysics Data System (ADS)
Challamel, Noël; Rakotomanana, Lalaonirina; Le Marrec, Loïc
2009-08-01
Nonlocal continuum mechanics allows one to account for the small length scale effect that becomes significant when dealing with micro- or nano-structures. This Note investigates a model of wave propagation in a nonlocal elastic material. We show that a dispersive wave equation is obtained from a nonlocal elastic constitutive law, based on a mixture of a local and a nonlocal strain. This model comprises both the classical gradient model and the Eringen's integral model. The dynamic properties of the model are discussed, and corroborate well some recent theoretical studies published to unify both static and dynamics gradient elasticity theories. Moreover, an excellent matching of the dispersive curve of the Born-Kármán model of lattice dynamics is obtained with such nonlocal model. To cite this article: N. Challamel et al., C. R. Mecanique 337 (2009).
NASA Astrophysics Data System (ADS)
Mudunuru, M. K.; Nakshatrala, K.
2012-12-01
Advection-Diffusion-Reaction (ADR) equations naturally arises in many physical phenomena, which include seepage of contaminants in heterogeneous porous media, transport of injected tracers due to the flow of oil in a petroleum reservoir, and degradation of a deformable solid due to diffusing chemical species. Vast literature exists on how to solve this equation in the cases when the medium is isotropic, velocity field being divergence free, and for advection-dominated problems. However, it is well know that many popular finite element formulations (e.g., the standard Galerkin formulation, stabilized methods, variational multi-scale methods, subgrid-scale methods, and primitive least-squares formulations) do not satisfy element-by-element mass/species balance and do not produce non-negative solutions on general computational grids. Various post-processing based methods were developed in order to recover some properties of computed numerical solutions. Most of these post-processing techniques are ad hoc, and are not variationally consistent. In this poster, we shall present a novel numerical methodology for ADR equations that satisfy discrete maximum principles, the non-negative constraint, and element-by-element mass/species balance. The methodology can handle general computational grids, no additional restrictions on time-step, and for heterogeneous anisotropic media. Several numerical results pertinent to advection-dominated ADR problems will be presented to illustrate the performance of the proposed numerical formulation.
Dispersion equation of gravito-MHD waves
NASA Astrophysics Data System (ADS)
Jovanović, Gordana
2016-03-01
We derive the dispersion equation for gravito-MHD waves in an isothermal, gravitationally stratified plasma with a horizontal inhomogeneous magnetic field. In the present model the sound and the Alfvén speeds are constant. It is known that in this model analytical solutions can be obtained for linearized perturbations. There are three modes propagating in the considered plasma: the fast, the slow and the Alfvén mode, all modified by gravity. In the extreme short wavelength limit, these waves propagate in a locally uniform plasma. The waves with larger wavelengths will be affected by the nonuniformity of the medium resulting from the action of gravitational force ρg. In the case without magnetic field these waves become gravito-acoustic waves.
Lancaster, Jill; Downes, Barbara J
2014-12-01
Many communities comprise species that select resources that are patchily distributed in an environment that is otherwise unsuitable or suboptimal. Effects of this patchiness can depend on the characteristics of patch arrays and animal movements, and produce non-intuitive outcomes in which population densities are unrelated to resource abundance. Resource mosaics are predicted to have only weak effects, however, where patches are ephemeral or organisms are transported advectively. The running waters of streams and benthic invertebrates epitomize such systems, but empirical tests of resource mosaics are scarce. We sampled 15 common macroinvertebrates inhabiting distinct detritus patches at four sites within a sand-bed stream, where detritus formed a major resource of food and living space. At each site, environmental variables were measured for 100 leaf packs; invertebrates were counted in 50 leaf packs. Sites differed in total abundance of detritus, leaf pack sizes and invertebrate densities. Multivariate analysis indicated that patch size was the dominant environmental variable, but invertebrate densities differed significantly between sites even after accounting for patch size. Leaf specialists showed positive and strong density-area relationships, except where the patch size range was small and patches were aggregated. In contrast, generalist species had weaker and variable responses to patch sizes. Population densities were not associated with total resource abundance, with the highest densities of leaf specialists in sites with the least detritus. Our results demonstrate that patchy resources can affect species even in communities where species are mobile, have advective dispersal, and patches are relatively ephemeral. PMID:25190216
NASA Technical Reports Server (NTRS)
Hung, R. J.; Liaw, G. S.
1980-01-01
The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.
Stochastic differential equations and turbulent dispersion
NASA Technical Reports Server (NTRS)
Durbin, P. A.
1983-01-01
Aspects of the theory of continuous stochastic processes that seem to contribute to an understanding of turbulent dispersion are introduced and the theory and philosophy of modelling turbulent transport is emphasized. Examples of eddy diffusion examined include shear dispersion, the surface layer, and channel flow. Modeling dispersion with finite-time scale is considered including the Langevin model for homogeneous turbulence, dispersion in nonhomogeneous turbulence, and the asymptotic behavior of the Langevin model for nonhomogeneous turbulence.
NASA Astrophysics Data System (ADS)
Zaramella, M.; Marion, A.; Lewandowski, J.; Nützmann, G.
2016-07-01
Solute transport in rivers is controlled by surface flow hydrodynamics and by transient storage in dead zones, pockets of vegetation and hyporheic sediments where mass exchange and retention are governed by complex mechanisms. The physics of these processes are generally investigated by optimization of transient storage models (TSMs) to experimental data often yielding inconsistent and equifinal parameter sets. Uncertainty on parameters estimation is found to depend not only on the rates of exchange between the stream and storage zones, the stream-water velocity and the stream reach length according to the experimental Damkohler number (DaI), but also on the relative significance between transient storage and longitudinal dispersion on breakthrough curves (BTCs). An optimization strategy was developed and applied to an experimental dataset obtained from tracer tests in a small lowland river, analyzing BTCs generated through tracer injections under different conditions. The method supplies a tool to estimate model parameters from observed data through the analysis of the relative parameter significance. To analyze model performance a double compartment TSM was optimized by a regular fit procedure based on simple root mean square error minimization and by a fit based on a relative significance analysis of mechanism signatures. As a result consistent longitudinal dispersion and transient storage parameters were obtained when the signature targeted optimization was used.
NASA Astrophysics Data System (ADS)
Ho, D. T.; Schlosser, P.; Schlosser, P.; Schlosser, P.; Caplow, T.; Garrison, M. R.
2001-12-01
In recent years, deliberate tracer release experiments have been used in the ocean, rivers, lakes, and groundwater flow systems to study advection, mixing, air-water gas exchange, and exchanges between subsystems. Here we report results from a recent deliberate tracer release experiment conducted in the tidal Hudson River. On July 25, 2001, ca. 3.3 moles of the inert gas SF6 were injected into the Hudson River near Newburgh, NY at a depth of about 6 m. Subsequently, the SF6 was monitored from a boat (Riverkeeper) by pumping water (from 2 m depth) via a submersed pump mounted on the front of the boat, through a gas extraction unit, followed by measurement using an onboard gas chromatograph. The measurement interval was about 2 minutes and the maximum speed of the boat was about 15 km h-1. This allowed us to obtain detailed surveys of the temporal evolution of the tracer plume for 14 days. Initial results from the experiment show that during July/August 2001, there was virtually no net downward advection of the water body originally tagged with SF6. Instead, we observed rapid mixing of the tracer-tagged water up- and down-river. After one week, the tracer-tagged water could be detected over a stretch of 70 km along the axis of the river channel. At this time, the stretch of river labeled with concentrations >50% of the peak value was about 14 km. After two weeks, the tracer-tagged water had extended to over 90 km, while 28 km had SF6 concentrations >50% of the peak value. Vertical mixing into depressions on the bottom of the river reaching more than 175 feet seemed to be rapid. Dispersion coefficients and vertical turbulent exchange coefficients will be discussed.
NASA Astrophysics Data System (ADS)
Bouchelaghem, F.; Vulliet, L.; Leroy, D.; Laloui, L.; Descoeudres, F.
2001-10-01
A model was developed, to describe miscible grout propagation in a saturated deformable porous medium, based on Bear's statistical model with spatial volume averaging. In a previous paper, the model was first successfully confronted to one-dimensional laboratory experiments.In the present paper, the numerical model is used to simulate practical grouting operation in a cylindrical injection model. The cylindrical injection model lends itself to study main flow and propagation character istics for a dispersed suspension-type grout, under axisymmetric conditions close to real scale conditions.Comparison between numerical solutions and experimental results is essential to confirm the validity and accuracy of the proposed model from a phenomenological standpoint. The numerical model performances show that the underlying mathematical model constitutes a realistic predictive model reproducing most prominent features during injection of a suspension-type grout into a deformable porous medium. The basic mechanism by which injected miscible grout permeates a soil mass is discussed in detail. Such a tool leads to quality control criteria for grouting on a theoretical basis, which complements existing criteria acquired through engineering practice.
Webb, S.W.
1996-05-01
Two models for gas-phase diffusion and advection in porous media, the Advective-Dispersive Model (ADM) and the Dusty-Gas Model (DGM), are reviewed. The ADM, which is more widely used, is based on a linear addition of advection calculated by Darcy`s Law and ordinary diffusion using Fick`s Law. Knudsen diffusion is often included through the use of a Klinkenberg factor for advection, while the effect of a porous medium on the diffusion process is through a porosity-tortuosity-gas saturation multiplier. Another, more comprehensive approach for gas-phase transport in porous media has been formulated by Evans and Mason, and is referred to as the Dusty- Gas Model (DGM). This model applies the kinetic theory of gases to the gaseous components and the porous media (or ``dust``) to develop an approach for combined transport due to ordinary and Knudsen diffusion and advection including porous medium effects. While these two models both consider advection and diffusion, the formulations are considerably different, especially for ordinary diffusion. The various components of flow (advection and diffusion) are compared for both models. Results from these two models are compared to isothermal experimental data for He-Ar gas diffusion in a low-permeability graphite. Air-water vapor comparisons have also been performed, although data are not available, for the low-permeability graphite system used for the helium-argon data. Radial and linear air-water heat pipes involving heat, advection, capillary transport, and diffusion under nonisothermal conditions have also been considered.
NASA Astrophysics Data System (ADS)
Vikas, Kumar; K. Gupta, R.; Ram, Jiwari
2014-03-01
In this paper, the variable-coefficient diffusion—advection (DA) equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by determining the complete sets of point symmetries of this equation, and then exact and numerical solutions are reported for the reduced second-order nonlinear ordinary differential equations. Further, an extended (G'/G)-expansion method is applied to the DA equation to construct some new non-traveling wave solutions.
Power-law spatial dispersion from fractional Liouville equation
Tarasov, Vasily E.
2013-10-15
A microscopic model in the framework of fractional kinetics to describe spatial dispersion of power-law type is suggested. The Liouville equation with the Caputo fractional derivatives is used to obtain the power-law dependence of the absolute permittivity on the wave vector. The fractional differential equations for electrostatic potential in the media with power-law spatial dispersion are derived. The particular solutions of these equations for the electric potential of point charge in this media are considered.
Symmetry analysis for a class of nonlinear dispersive equations
NASA Astrophysics Data System (ADS)
Charalambous, K.; Sophocleous, C.
2015-05-01
A class of dispersive equations is studied within the framework of group analysis of differential equations. The enhanced Lie group classification is achieved. The complete list of equivalence transformations is presented. It is shown that certain equations from the class admit nonclassical reductions. Potential and potential nonclassical symmetries are also considered.
Alcaraz, Mar; García-Gil, Alejandro; Vázquez-Suñé, Enric; Velasco, Violeta
2016-02-01
Borehole Heat Exchangers (BHEs) are increasingly being used to exploit shallow geothermal energy. This paper presents a new methodology to provide a response to the need for a regional quantification of the geothermal potential that can be extracted by BHEs and the associated environmental impacts. A set of analytical solutions facilitates accurate calculation of the heat exchange of BHEs with the ground and its environmental impacts. For the first time, advection and dispersion heat transport mechanisms and the temporal evolution from the start of operation of the BHE are taken into account in the regional estimation of shallow geothermal resources. This methodology is integrated in a GIS environment, which facilitates the management of input and output data at a regional scale. An example of the methodology's application is presented for Barcelona, in Spain. As a result of the application, it is possible to show the strengths and improvements of this methodology in the development of potential maps of low temperature geothermal energy as well as maps of environmental impacts. The minimum and maximum energy potential values for the study site are 50 and 1800 W/m(2) for a drilled depth of 100 m, proportionally to Darcy velocity. Regarding to thermal impacts, the higher the groundwater velocity and the energy potential, the higher the size of the thermal plume after 6 months of exploitation, whose length ranges from 10 to 27 m long. A sensitivity analysis was carried out in the calculation of heat exchange rate and its impacts for different scenarios and for a wide range of Darcy velocities. The results of this analysis lead to the conclusion that the consideration of dispersion effects and temporal evolution of the exploitation prevent significant differences up to a factor 2.5 in the heat exchange rate accuracy and up to several orders of magnitude in the impacts generated. PMID:26605833
Realistic dispersion kernels applied to cohabitation reaction dispersion equations
NASA Astrophysics Data System (ADS)
Isern, Neus; Fort, Joaquim; Pérez-Losada, Joaquim
2008-10-01
We develop front spreading models for several jump distance probability distributions (dispersion kernels). We derive expressions for a cohabitation model (cohabitation of parents and children) and a non-cohabitation model, and apply them to the Neolithic using data from real human populations. The speeds that we obtain are consistent with observations of the Neolithic transition. The correction due to the cohabitation effect is up to 38%.
An integrable shallow water equation with linear and nonlinear dispersion.
Dullin, H R; Gottwald, G A; Holm, D D
2001-11-01
We use asymptotic analysis and a near-identity normal form transformation from water wave theory to derive a 1+1 unidirectional nonlinear wave equation that combines the linear dispersion of the Korteweg-deVries (KdV) equation with the nonlinear/nonlocal dispersion of the Camassa-Holm (CH) equation. This equation is one order more accurate in asymptotic approximation beyond KdV, yet it still preserves complete integrability via the inverse scattering transform method. Its traveling wave solutions contain both the KdV solitons and the CH peakons as limiting cases. PMID:11690414
Erratum: A Comparison of Closures for Stochastic Advection-Diffusion Equations
Jarman, Kenneth D.; Tartakovsky, Alexandre M.
2015-01-01
This note corrects an error in the authors' article [SIAM/ASA J. Uncertain. Quantif., 1 (2013), pp. 319 347] in which the cited work [Neuman, Water Resour. Res., 29(3) (1993), pp. 633 645] was incorrectly represented and attributed. Concentration covariance equations presented in our article as new were in fact previously derived in the latter work. In the original abstract, the phrase " . . .we propose a closed-form approximation to two-point covariance as a measure of uncertainty. . ." should be replaced by the phrase " . . .we study a closed-form approximation to two-point covariance, previously derived in [Neuman 1993], as a measure of uncertainty." The primary results in our article--the analytical and numerical comparison of existing closure methods for specific example problems are not changed by this correction.
Wang, Wei; Shu, Chi-Wang; Yee, H.C.; Sjögreen, Björn
2012-01-01
A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.
Wang, H.; Man, S.; Ewing, R.E.; Qin, G.; Lyons, S.L.; Al-Lawatia, M.
1999-06-10
Many difficult problems arise in the numerical simulation of fluid flow processes within porous media in petroleum reservoir simulation and in subsurface contaminant transport and remediation. The authors develop a family of Eulerian-Lagrangian localized adjoint methods for the solution of the initial-boundary value problems for first-order advection-reaction equations on general multi-dimensional domains. Different tracking algorithms, including the Euler and Runge-Kutta algorithms, are used. The derived schemes, which are full mass conservative, naturally incorporate inflow boundary conditions into their formulations and do not need any artificial outflow boundary conditions. Moreover, they have regularly structured, well-conditioned, symmetric, and positive-definite coefficient matrices, which can be efficiently solved by the conjugate gradient method in an optimal order number of iterations without any preconditioning needed. Numerical results are presented to compare the performance of the ELLAM schemes with many well studied and widely used methods, including the upwind finite difference method, the Galerkin and the Petrov-Galerkin finite element methods with backward-Euler or Crank-Nicolson temporal discretization, the streamline diffusion finite element methods, the monotonic upstream-centered scheme for conservation laws (MUSCL), and the Minmod scheme.
Drift-free kinetic equations for turbulent dispersion.
Bragg, A; Swailes, D C; Skartlien, R
2012-11-01
The dispersion of passive scalars and inertial particles in a turbulent flow can be described in terms of probability density functions (PDFs) defining the statistical distribution of relevant scalar or particle variables. The construction of transport equations governing the evolution of such PDFs has been the subject of numerous studies, and various authors have presented formulations for this type of equation, usually referred to as a kinetic equation. In the literature it is often stated, and widely assumed, that these PDF kinetic equation formulations are equivalent. In this paper it is shown that this is not the case, and the significance of differences among the various forms is considered. In particular, consideration is given to which form of equation is most appropriate for modeling dispersion in inhomogeneous turbulence and most consistent with the underlying particle equation of motion. In this regard the PDF equations for inertial particles are considered in the limit of zero particle Stokes number and assessed against the fully mixed (zero-drift) condition for fluid points. A long-standing question regarding the validity of kinetic equations in the fluid-point limit is answered; it is demonstrated formally that one version of the kinetic equation (derived using the Furutsu-Novikov method) provides a model that satisfies this zero-drift condition exactly in both homogeneous and inhomogeneous systems. In contrast, other forms of the kinetic equation do not satisfy this limit or apply only in a limited regime. PMID:23214875
Dispersion equation for ballooning modes in two-component plasma
NASA Astrophysics Data System (ADS)
Kozlov, D. A.; Mazur, N. G.; Pilipenko, V. A.; Fedorov, E. N.; Fedorov
2014-06-01
The ballooning magnetohydrodynamic (MHD) modes have been often suggested as a possible instability trigger of the substorm onset, and a mechanism of compressional waves in the outer magnetosphere and magnetotail. Commonly, these disturbances are characterized by the local dispersion equation that is widely applied for the description of ultra-low-frequency oscillatory disturbances and instabilities in the nightside magnetosphere. In realistic situations, especially in the inner magnetosphere, the magnetospheric plasma is composed of two components: background `cold' plasma and `hot' particles. The ballooning disturbances in a two-component plasma immersed into a curved magnetic field are described with the system of coupled equations for the Alfvén and slow magnetosonic (SMS) modes. We have reduced the basic system of MHD equations to the dispersion equation for the small-scale in transverse direction disturbances, and applied WKB approximation along a field line. As a result, we have derived a dispersion equation that can be used for geophysical applications. In particular, from this relationship the dispersion, instability threshold, and stop-bands of the Alfvén and SMS modes in two-component plasma have been examined.
Guyonvarch, Estelle; Ramin, Elham; Kulahci, Murat; Plósz, Benedek Gy
2015-10-15
The present study aims at using statistically designed computational fluid dynamics (CFD) simulations as numerical experiments for the identification of one-dimensional (1-D) advection-dispersion models - computationally light tools, used e.g., as sub-models in systems analysis. The objective is to develop a new 1-D framework, referred to as interpreted CFD (iCFD) models, in which statistical meta-models are used to calculate the pseudo-dispersion coefficient (D) as a function of design and flow boundary conditions. The method - presented in a straightforward and transparent way - is illustrated using the example of a circular secondary settling tank (SST). First, the significant design and flow factors are screened out by applying the statistical method of two-level fractional factorial design of experiments. Second, based on the number of significant factors identified through the factor screening study and system understanding, 50 different sets of design and flow conditions are selected using Latin Hypercube Sampling (LHS). The boundary condition sets are imposed on a 2-D axi-symmetrical CFD simulation model of the SST. In the framework, to degenerate the 2-D model structure, CFD model outputs are approximated by the 1-D model through the calibration of three different model structures for D. Correlation equations for the D parameter then are identified as a function of the selected design and flow boundary conditions (meta-models), and their accuracy is evaluated against D values estimated in each numerical experiment. The evaluation and validation of the iCFD model structure is carried out using scenario simulation results obtained with parameters sampled from the corners of the LHS experimental region. For the studied SST, additional iCFD model development was carried out in terms of (i) assessing different density current sub-models; (ii) implementation of a combined flocculation, hindered, transient and compression settling velocity function; and (iii
Global well-posedness for the 2D dispersive SQG equation and inviscid Boussinesq equations
NASA Astrophysics Data System (ADS)
Wan, Renhui; Chen, Jiecheng
2016-08-01
In this paper, we obtain global well-posedness for the 2D dispersive SQG equation and inviscid Boussinesq equations. Our works are consistent with the corresponding works by Elgindi-Widmayer (SIAM J Math Anal 47:4672-4684, 2015) in the special case {A=κ=1}. In addition, our result concerning the SQG equation can be regarded as the borderline case of the work by Cannone et al. (Proc Lond Math Soc 106:650-674, 2013).
Romero-González, J; Walton, J C; Peralta-Videa, J R; Rodríguez, E; Romero, J; Gardea-Torresdey, J L
2009-01-15
The biosorption of Cr(III) onto packed columns of Agave lechuguilla was analyzed using an advective-dispersive (AD) model and its analytical solution. Characteristic parameters such as axial dispersion coefficients, retardation factors, and distribution coefficients were predicted as functions of inlet ion metal concentration, time, flow rate, bed density, cross-sectional column area, and bed length. The root-mean-square-error (RMSE) values 0.122, 0.232, and 0.285 corresponding to the flow rates of 1, 2, and 3 (10(-3))dm3min(-1), respectively, indicated that the AD model provides an excellent approximation of the simulation of lumped breakthrough curves for the adsorption of Cr(III) by lechuguilla biomass. Therefore, the model can be used for design purposes to predict the effect of varying operational conditions. PMID:18462882
Merritt, M.L.
1993-01-01
The simulation of the transport of injected freshwater in a thin brackish aquifer, overlain and underlain by confining layers containing more saline water, is shown to be influenced by the choice of the finite-difference approximation method, the algorithm for representing vertical advective and dispersive fluxes, and the values assigned to parametric coefficients that specify the degree of vertical dispersion and molecular diffusion that occurs. Computed potable water recovery efficiencies will differ depending upon the choice of algorithm and approximation method, as will dispersion coefficients estimated based on the calibration of simulations to match measured data. A comparison of centered and backward finite-difference approximation methods shows that substantially different transition zones between injected and native waters are depicted by the different methods, and computed recovery efficiencies vary greatly. Standard and experimental algorithms and a variety of values for molecular diffusivity, transverse dispersivity, and vertical scaling factor were compared in simulations of freshwater storage in a thin brackish aquifer. Computed recovery efficiencies vary considerably, and appreciable differences are observed in the distribution of injected freshwater in the various cases tested. The results demonstrate both a qualitatively different description of transport using the experimental algorithms and the interrelated influences of molecular diffusion and transverse dispersion on simulated recovery efficiency. When simulating natural aquifer flow in cross-section, flushing of the aquifer occurred for all tested coefficient choices using both standard and experimental algorithms. ?? 1993.
Langevin equation model of dispersion in the convective boundary layer
Nasstrom, J S
1998-08-01
This dissertation presents the development and evaluation of a Lagrangian stochastic model of vertical dispersion of trace material in the convective boundary layer (CBL). This model is based on a Langevin equation of motion for a fluid particle, and assumes the fluid vertical velocity probability distribution is skewed and spatially homogeneous. This approach can account for the effect of large-scale, long-lived turbulent structures and skewed vertical velocity distributions found in the CBL. The form of the Langevin equation used has a linear (in velocity) deterministic acceleration and a skewed randomacceleration. For the case of homogeneous fluid velocity statistics, this ""linear-skewed" Langevin equation can be integrated explicitly, resulting in a relatively efficient numerical simulation method. It is shown that this approach is more efficient than an alternative using a "nonlinear-Gaussian" Langevin equation (with a nonlinear deterministic acceleration and a Gaussian random acceleration) assuming homogeneous turbulence, and much more efficient than alternative approaches using Langevin equation models assuming inhomogeneous turbulence. "Reflection" boundary conditions for selecting a new velocity for a particle that encounters a boundary at the top or bottom of the CBL were investigated. These include one method using the standard assumption that the magnitudes of the particle incident and reflected velocities are positively correlated, and two alternatives in which the magnitudes of these velocities are negatively correlated and uncorrelated. The constraint that spatial and velocity distributions of a well-mixed tracer must be the same as those of the fluid, was used to develop the Langevin equation models and the reflection boundary conditions. The two Langevin equation models and three reflection methods were successfully tested using cases for which exact, analytic statistical properties of particle velocity and position are known, including well
The dispersion equation of the induced Smith-Purcell instability
NASA Astrophysics Data System (ADS)
Klochkov, Dmitry N.; Artemyev, Alexander I.; Oganesyan, Koryun B.; Rostovtsev, Yuri V.; Scully, Marlan O.; Hu, Chin-Kun
2010-09-01
The simplest model of a magnetized infinitely thin electron beam is considered. For the grating, where the depth of grooves is a small parameter, the dispersion equation of the induced Smith-Purcell instability was obtained. It was found that the condition of the Thompson or the Raman regimes of excitation does not depend on beam current but depends on the height of the beam above the grating surface. The growth rate of instability in both cases is proportional to the square root of the electron beam current.
NASA Astrophysics Data System (ADS)
Vukadinovic, J.; Dedits, E.; Poje, A. C.; Schäfer, T.
2015-08-01
We consider the two-dimensional advection-diffusion equation (ADE) on a bounded domain subject to Dirichlet or von Neumann boundary conditions involving a Liouville integrable Hamiltonian. Transformation to action-angle coordinates permits averaging in time and angle, resulting in an equation that allows for separation of variables. The Fourier transform in the angle coordinate transforms the equation into an effective diffusive equation and a countable family of non-self-adjoint Schrödinger equations. For the corresponding Liouville-Sturm problem, we apply complex-plane WKB methods to study the spectrum in the semi-classical limit for vanishing diffusivity. The spectral limit graph is found to consist of analytic curves (branches) related to Stokes graphs forming a tree-structure. Eigenvalues in the neighborhood of branches emanating from the imaginary axis are subject to various sublinear power laws with respect to diffusivity, leading to convection-enhanced rates of dissipation of the corresponding modes. The solution of the ADE converges in the limit of vanishing diffusivity to the solution of the effective diffusion equation on convective time scales that are sublinear with respect to the diffusive time scales.
Wagner, Brian J.; Gorelick, Steven M.
1986-01-01
A simulation nonlinear multiple-regression methodology for estimating parameters that characterize the transport of contaminants is developed and demonstrated. Finite difference containment transport simulation is combined with a nonlinear weighted least squares multiple-regression procedure. The technique provides optimal parameter estimates and gives statistics for assessing the reliability of these estimates under certain general assumptions about the distributions of the random measurement errors. Monte Carlo analysis is used to estimate parameter reliability for a hypothetical homogeneous soil column for which concentration data contain large random measurement errors. The value of data collected spatially versus data collected temporally was investigated for estimation of velocity, dispersion coefficient, effective porosity, first-order decay rate, and zero-order production. The use of spatial data gave estimates that were 2-3 times more reliable than estimates based on temporal data for all parameters except velocity. (Estimated author abstract) Refs.
NASA Astrophysics Data System (ADS)
Furbish, David Jon; Childs, Elise M.; Haff, Peter K.; Schmeeckle, Mark W.
2009-09-01
We formulate soil grain transport by rain splash as a stochastic advection-dispersion process. By taking into account the intermittency of grain motions activated by raindrop impacts, the formulation indicates that gradients in raindrop intensity, and thus grain activity (the volume of grains in motion per unit area) can be as important as gradients in grain concentration and surface slope in effecting transport. This idea is confirmed by rain splash experiments and manifest in topographic roughening via mound growth beneath desert shrubs. The formulation provides a framework for describing transport and dispersal of any soil material moveable by rain splash, including soil grains, soil-borne pathogens and nutrients, seeds, or debitage. As such it shows how classic models of topographic "diffusion" reflect effects of slope-dependent grain drift, not diffusion, and it highlights the role of rain splash in the ecological behavior of desert shrubs as "resource islands." Specifically, the growth of mounds beneath shrub canopies, where differential rain splash initially causes more grains to be splashed inward beneath the protective canopy than outward, involves the "harvesting" of nearby soil material, including nutrients. Mounds thus represent temporary storage of soil derived from areas surrounding the shrubs. As the inward grain flux associated with differential rain splash is sustained over the shrub lifetime, mound material is effectively sequestered from erosional processes that might otherwise move this material downslope. With shrub death and loss of the protective canopy, differential rain splash vanishes and the mound material is dispersed to the surrounding area, again subject to downslope movement.
The zero dispersion limits of nonlinear wave equations
Tso, T.
1992-01-01
In chapter 2 the author uses functional analytic methods and conservation laws to solve the initial-value problem for the Korteweg-de Vries equation, the Benjamin-Bona-Mahony equation, and the nonlinear Schroedinger equation for initial data that satisfy some suitable conditions. In chapter 3 the energy estimates are used to show that the strong convergence of the family of the solutions of the KdV equation obtained in chapter 2 in H[sup 3](R) as [epsilon] [yields] 0; also, it is shown that the strong L[sup 2](R)-limit of the solutions of the BBM equation as [epsilon] [yields] 0 before a critical time. In chapter 4 the author uses the Whitham modulation theory and averaging method to find the 2[pi]-periodic solutions and the modulation equations of the KdV equation, the BBM equation, the Klein-Gordon equation, the NLS equation, the mKdV equation, and the P-system. It is shown that the modulation equations of the KdV equation, the K-G equation, the NLS equation, and the mKdV equation are hyperbolic but those of the BBM equation and the P-system are not hyperbolic. Also, the relations are studied of the KdV equation and the mKdV equation. Finally, the author studies the complex mKdV equation to compare with the NLS equation, and then study the complex gKdV equation.
NASA Astrophysics Data System (ADS)
Yan, Zhenya; Bluman, George
2002-11-01
The special exact solutions of nonlinearly dispersive Boussinesq equations (called B( m, n) equations), utt- uxx- a( un) xx+ b( um) xxxx=0, is investigated by using four direct ansatze. As a result, abundant new compactons: solitons with the absence of infinite wings, solitary patterns solutions having infinite slopes or cups, solitary waves and singular periodic wave solutions of these two equations are obtained. The variant is extended to include linear dispersion to support compactons and solitary patterns in the linearly dispersive Boussinesq equations with m=1. Moreover, another new compacton solution of the special case, B(2,2) equation, is also found.
Local volume-time averaged equations of motion for dispersed, turbulent, multiphase flows
Sha, W.T.; Slattery, J.C.
1980-11-01
In most flows of liquids and their vapors, the phases are dispersed randomly in both space and time. These dispersed flows can be described only statistically or in terms of averages. Local volume-time averaging is used here to derive a self-consistent set of equations governing momentum and energy transfer in dispersed, turbulent, multiphase flows. The empiricisms required for use with these equations are the subject of current research.
CONTAMINANT TRANSPORT IN SEDIMENT UNDER THE INFLUENCE OF ADVECTIVE FLUX
Chemical flux across the sediment/water interface is controlled by a combination of diffusive, dispersive and advective processes. The advective process is a function of submarine groundwater discharge and tidal effects. In areas where surface water interacts with groundwater, ...
Propagation estimates for dispersive wave equations: Application to the stratified wave equation
NASA Astrophysics Data System (ADS)
Pravica, David W.
1999-01-01
The plane-stratified wave equation (∂t2+H)ψ=0 with H=-c(y)2∇z2 is studied, where z=x⊕y, x∈Rk, y∈R1 and |c(y)-c∞|→0 as |y|→∞. Solutions to such an equation are solved for the propagation of waves through a layered medium and can include waves which propagate in the x-directions only (i.e., trapped modes). This leads to a consideration of the pseudo-differential wave equation (∂t2+ω(-Δx))ψ=0 such that the dispersion relation ω(ξ2) is analytic and satisfies c1⩽ω'(ξ2)⩽c2 for c*>0. Uniform propagation estimates like ∫|x|⩽|t|αE(UtP±φ0)dkx⩽Cα,β(1+|t|)-β∫E(φ0)dkx are obtained where Ut is the evolution group, P± are projection operators onto the Hilbert space of initial conditions φ∈H and E(ṡ) is the local energy density. In special cases scattering of trapped modes off a local perturbation satisfies the causality estimate ||P+ρΛjSP-ρΛk||⩽Cνρ-ν for each ν<1/2. Here P+ρΛj (P-ρΛk) are remote outgoing/detector (incoming/transmitter) projections for the jth (kth) trapped mode. Also Λ⋐R+ is compact, so the projections localize onto formally-incoming (eventually-outgoing) states.
Integrodifference equations in patchy landscapes : I. Dispersal Kernels.
Musgrave, Jeffrey; Lutscher, Frithjof
2014-09-01
What is the effect of individual movement behavior in patchy landscapes on redistribution kernels? To answer this question, we derive a number of redistribution kernels from a random walk model with patch dependent diffusion, settling, and mortality rates. At the interface of two patch types, we integrate recent results on individual behavior at the interface. In general, these interface conditions result in the probability density function of the random walker being discontinuous at an interface. We show that the dispersal kernel can be characterized as the Green's function of a second-order differential operator. Using this characterization, we illustrate the kind of (discontinuous) dispersal kernels that result from our approach, using three scenarios. First, we assume that dispersal distance is small compared to patch size, so that a typical disperser crosses at most one interface during the dispersal phase. Then we consider a single bounded patch and generate kernels that will be useful to study the critical patch size problem in our sequel paper. Finally, we explore dispersal kernels in a periodic landscape and study the dependence of certain dispersal characteristics on model parameters. PMID:23907527
This manual describes the next generation of the modular three-dimensional transport model, MT3D, with significantly expanded capabilities, including the addition of (a) a third-order total-variation-diminishing (TVD) scheme for solving the advection term that is mass conservativ...
Concentration through large advection
NASA Astrophysics Data System (ADS)
Aleja, D.; López-Gómez, J.
2014-11-01
In this paper we extend the elegant results of Chen, Lam and Lou [6, Section 2], where a concentration phenomenon was established as the advection blows up, to a general class of adventive-diffusive generalized logistic equations of degenerate type. Our improvements are really sharp as we allow the carrying capacity of the species to vanish in some subdomain with non-empty interior. The main technical devices used in the derivation of the concentration phenomenon are Proposition 3.2 of Cano-Casanova and López-Gómez [5], Theorem 2.4 of Amann and López-Gómez [1] and the classical Harnack inequality. By the relevance of these results in spatial ecology, complete technical details seem imperative, because the proof of Theorem 2.2 of [6] contains some gaps originated by an “optimistic” use of Proposition 3.2 of [5]. Some of the general assumptions of [6] are substantially relaxed.
Finite-Difference Time-Domain solution of Maxwell's equations for the dispersive ionosphere
NASA Astrophysics Data System (ADS)
Nickisch, L. J.; Franke, P. M.
1992-10-01
The Finite-Difference Time-Domain (FDTD) technique is a conceptually simple, yet powerful, method for obtaining numerical solutions to electromagnetic propagation problems. However, the application of FDTD methods to problems in ionospheric radiowave propagation is complicated by the dispersive nature of the ionospheric plasma. In the time domain, the electric displacement is the convolution of the dielectric tensor with the electric field, and thus requires information from the entire signal history. This difficulty can be avoided by returning to the dynamical equations from which the dielectric tensor is derived. By integrating these differential equations simultaneously with the Maxwell equations, temporal dispersion is fully incorporated.
On a hierarchy of nonlinearly dispersive generalized Korteweg - de Vries evolution equations
Christov, Ivan C.
2015-08-20
We propose a hierarchy of nonlinearly dispersive generalized Korteweg–de Vries (KdV) evolution equations based on a modification of the Lagrangian density whose induced action functional the KdV equation extremizes. Two recent nonlinear evolution equations describing wave propagation in certain generalized continua with an inherent material length scale are members of the proposed hierarchy. Like KdV, the equations from the proposed hierarchy possess Hamiltonian structure. Unlike KdV, the solutions to these equations can be compact (i.e., they vanish outside of some open interval) and, in addition, peaked. Implicit solutions for these peaked, compact traveling waves (“peakompactons”) are presented.
The zero dispersion limit for the Korteweg-deVries KdV equation.
Lax, P D; Levermore, C D
1979-08-01
We use the inverse scattering method to determine the weak limit of solutions of the Korteweg-deVries equation as dispersion tends to zero. The limit, valid for all time, is characterized in terms of a quadratic programming problem which can be solved with the aid of function theoretic methods. For large t, the solutions satisfy Whitham's averaged equations at some times and the equations found by Flaschka et al. at other times. PMID:16592690
A two-sided fractional conservation of mass equation
NASA Astrophysics Data System (ADS)
Olsen, Jeffrey S.; Mortensen, Jeff; Telyakovskiy, Aleksey S.
2016-05-01
A two-sided fractional conservation of mass equation is derived by using left and right fractional Mean Value Theorems. This equation extends the one-sided fractional conservation of mass equation of Wheatcraft and Meerschaert. Also, a two-sided fractional advection-dispersion equation is derived. The derivations are based on Caputo fractional derivatives.
NASA Astrophysics Data System (ADS)
Mudunuru, M. K.; Nakshatrala, K. B.
2016-01-01
We present a robust computational framework for advective-diffusive-reactive systems that satisfies maximum principles, the non-negative constraint, and element-wise species balance property. The proposed methodology is valid on general computational grids, can handle heterogeneous anisotropic media, and provides accurate numerical solutions even for very high Péclet numbers. The significant contribution of this paper is to incorporate advection (which makes the spatial part of the differential operator non-self-adjoint) into the non-negative computational framework, and overcome numerical challenges associated with advection. We employ low-order mixed finite element formulations based on least-squares formalism, and enforce explicit constraints on the discrete problem to meet the desired properties. The resulting constrained discrete problem belongs to convex quadratic programming for which a unique solution exists. Maximum principles and the non-negative constraint give rise to bound constraints while element-wise species balance gives rise to equality constraints. The resulting convex quadratic programming problems are solved using an interior-point algorithm. Several numerical results pertaining to advection-dominated problems are presented to illustrate the robustness, convergence, and the overall performance of the proposed computational framework.
NASA Astrophysics Data System (ADS)
Mohammed, K. Elboree
2015-10-01
In this paper, we investigate the traveling wave solutions for the nonlinear dispersive equation, Korteweg-de Vries Zakharov-Kuznetsov (KdV-ZK) equation and complex coupled KdV system by using extended simplest equation method, and then derive the hyperbolic function solutions include soliton solutions, trigonometric function solutions include periodic solutions with special values for double parameters and rational solutions. The properties of such solutions are shown by figures. The results show that this method is an effective and a powerful tool for handling the solutions of nonlinear partial differential equations (NLEEs) in mathematical physics.
Eulerian-Lagrangian solution of the convection-dispersion equation in natural co-ordinates.
Cheng, R.T.; Casulli, V.; Milford, S.N.
1984-01-01
The vast majority of numerical investigations of transport phenomena use an Eulerian formulation for the convenience that the computational grids are fixed in space. An Eulerian-Lagrangian method (ELM) of solution for the convection-dispersion equation is discussed and analyzed. The ELM uses the Lagrangian concept in an Eulerian computational grid system.-from Authors
Dispersion Interactions between Rare Gas Atoms: Testing the London Equation Using ab Initio Methods
ERIC Educational Resources Information Center
Halpern, Arthur M.
2011-01-01
A computational chemistry experiment is described in which students can use advanced ab initio quantum mechanical methods to test the ability of the London equation to account quantitatively for the attractive (dispersion) interactions between rare gas atoms. Using readily available electronic structure applications, students can calculate the…
Symmetries of the TDNLS equations for weakly nonlinear dispersive MHD waves
NASA Technical Reports Server (NTRS)
Webb, G. M.; Brio, M.; Zank, G. P.
1995-01-01
In this paper we consider the symmetries and conservation laws for the TDNLS equations derived by Hada (1993) and Brio, Hunter and Johnson, to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. The equations describe the interaction of the Alfven and magnetoacoustic modes near the triple umbilic, where the fast magnetosonic, slow magnetosonic and Alfven speeds coincide and a(g)(exp 2) = V(A)(exp 2) where a(g) is the gas sound speed and V(A) is the Alfven speed. We discuss Lagrangian and Hamiltonian formulations, and similarity solutions for the equations.
Monger, Gregg R.; Duncan, Candice Morrison; Brusseau, Mark L.
2015-01-01
A gas-phase tracer test (GTT) was conducted at a landfill in Tucson, AZ, to help elucidate the impact of landfill gas generation on the transport and fate of chlorinated aliphatic volatile organic contaminants (VOCs). Sulfur hexafluoride (SF6) was used as the non-reactive gas tracer. Gas samples were collected from a multiport monitoring well located 15.2 m from the injection well, and analyzed for SF6, CH4, CO2, and VOCs. The travel times determined for SF6 from the tracer test are approximately two to ten times smaller than estimated travel times that incorporate transport by only gas-phase diffusion. In addition, significant concentrations of CH4 and CO2 were measured, indicating production of landfill gas. Based on these results, it is hypothesized that the enhanced rates of transport observed for SF6 are caused by advective transport associated with landfill gas generation. The rates of transport varied vertically, which is attributed to multiple factors including spatial variability of water content, refuse mass, refuse permeability, and gas generation. PMID:26380532
Wave-packet formation at the zero-dispersion point in the Gardner-Ostrovsky equation.
Whitfield, A J; Johnson, E R
2015-05-01
The long-time effect of weak rotation on an internal solitary wave is the decay into inertia-gravity waves and the eventual emergence of a coherent, steadily propagating, nonlinear wave packet. There is currently no entirely satisfactory explanation as to why these wave packets form. Here the initial value problem is considered within the context of the Gardner-Ostrovsky, or rotation-modified extended Korteweg-de Vries, equation. The linear Gardner-Ostrovsky equation has maximum group velocity at a critical wave number, often called the zero-dispersion point. It is found here that a nonlinear splitting of the wave-number spectrum at the zero-dispersion point, where energy is shifted into the modulationally unstable regime of the Gardner-Ostrovsky equation, is responsible for the wave-packet formation. Numerical comparisons of the decay of a solitary wave in the Gardner-Ostrovsky equation and a derived nonlinear Schrödinger equation at the zero-dispersion point are used to confirm the spectral splitting. PMID:26066112
Nonlinear modulation of periodic waves in the small dispersion limit of the Benjamin-Ono equation
NASA Astrophysics Data System (ADS)
Matsuno, Y.
1998-12-01
The Whitham modulation theory is used to construct large time asymptotic solutions of the Benjamin-Ono (BO) equation in the small dispersion limit. For a wide class of initial data, asymptotic solutions are represented by a single-phase periodic solution of the BO equation with slowly varying amplitude and wave number. The Whitham system of modulation equations for these wave parameters has a very simple structure, and it can be solved exactly under appropriate boundary conditions. It is found that the oscillating zone expands with time, and eventually evolves into a train of solitary waves. In the case of localized initial data, the number density function of solitary waves is derived in a closed form. The resulting expression coincides with the corresponding formula obtained from the asymptotic theory based on the conservation laws of the BO equation. For steplike initial data, the total number of created solitary waves increases without limit in proportion to time.
NASA Astrophysics Data System (ADS)
Kim, Bong-Sik
Three dimensional (3D) Navier-Stokes-alpha equations are considered for uniformly rotating geophysical fluid flows (large Coriolis parameter f = 2O). The Navier-Stokes-alpha equations are a nonlinear dispersive regularization of usual Navier-Stokes equations obtained by Lagrangian averaging. The focus is on the existence and global regularity of solutions of the 3D rotating Navier-Stokes-alpha equations and the uniform convergence of these solutions to those of the original 3D rotating Navier-Stokes equations for large Coriolis parameters f as alpha → 0. Methods are based on fast singular oscillating limits and results are obtained for periodic boundary conditions for all domain aspect ratios, including the case of three wave resonances which yields nonlinear "2½-dimensional" limit resonant equations for f → 0. The existence and global regularity of solutions of limit resonant equations is established, uniformly in alpha. Bootstrapping from global regularity of the limit equations, the existence of a regular solution of the full 3D rotating Navier-Stokes-alpha equations for large f for an infinite time is established. Then, the uniform convergence of a regular solution of the 3D rotating Navier-Stokes-alpha equations (alpha ≠ 0) to the one of the original 3D rotating NavierStokes equations (alpha = 0) for f large but fixed as alpha → 0 follows; this implies "shadowing" of trajectories of the limit dynamical systems by those of the perturbed alpha-dynamical systems. All the estimates are uniform in alpha, in contrast with previous estimates in the literature which blow up as alpha → 0. Finally, the existence of global attractors as well as exponential attractors is established for large f and the estimates are uniform in alpha.
Antidiffusive velocities for multipass donor cell advection
Margolin, L.; Smolarkiewicz, P.K.
1999-01-01
Multidimensional positive definite advection transport algorithm (MPDATA) is an iterative process for approximating the advection equation, which uses a donor cell approximation to compensate for the truncation error of the originally specified donor cell scheme. This step may be repeated an arbitrary number of times, leading to successfully more accurate solutions to the advection equation. In this paper, the authors show how to sum the successive approximations analytically to find a single antidiffusive velocity that represents the effects of an arbitrary number of passes. The analysis is first done in one dimension to illustrate the method and then is repeated in two dimensions. The existence of cross terms in the truncation analysis of the two-dimensional equations introduces an extra complication into the calculation. The authors discuss the implementation of the antidiffusive velocities and provide some examples of applications, including a third-order accurate scheme.
NASA Astrophysics Data System (ADS)
Arqub, Omar Abu; El-Ajou, Ahmad; Momani, Shaher
2015-07-01
Building fractional mathematical models for specific phenomena and developing numerical or analytical solutions for these fractional mathematical models are crucial issues in mathematics, physics, and engineering. In this work, a new analytical technique for constructing and predicting solitary pattern solutions of time-fractional dispersive partial differential equations is proposed based on the generalized Taylor series formula and residual error function. The new approach provides solutions in the form of a rapidly convergent series with easily computable components using symbolic computation software. For method evaluation and validation, the proposed technique was applied to three different models and compared with some of the well-known methods. The resultant simulations clearly demonstrate the superiority and potentiality of the proposed technique in terms of the quality performance and accuracy of substructure preservation in the construct, as well as the prediction of solitary pattern solutions for time-fractional dispersive partial differential equations.
The modulational instability for the TDNLS equations for weakly nonlinear dispersive MHD waves
NASA Technical Reports Server (NTRS)
Webb, G. M.; Brio, M.; Zank, G. P.
1995-01-01
In this paper we study the modulational instability for the TDNLS equations derived by Hada (1993) and Brio, Hunter, and Johnson to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. We employ Whitham's averaged Lagrangian method to study the modulational instability. This complements studies of the modulational instability by Hada (1993) and Hollweg (1994), who did not use the averaged Lagrangian approach.
An optimally blended finite-spectral element scheme with minimal dispersion for Maxwell equations
NASA Astrophysics Data System (ADS)
Wajid, Hafiz Abdul; Ayub, Sobia
2012-10-01
We study the dispersive properties of the time harmonic Maxwell equations for optimally blended finite-spectral element scheme using tensor product elements defined on rectangular grid in d-dimensions. We prove and give analytical expressions for the discrete dispersion relations for this scheme. We find that for a rectangular grid (a) the analytical expressions for the discrete dispersion error in higher dimensions can be obtained using one dimensional discrete dispersion error expressions; (b) the optimum value of the blending parameter is p/(p+1) for all p∈N and for any number of spatial dimensions; (c) analytical expressions for the discrete dispersion relations for finite element and spectral element schemes can be obtained when the value of blending parameter is chosen to be 0 and 1 respectively; (d) the optimally blended scheme guarantees two additional orders of accuracy compared with standard finite element and spectral element schemes; and (e) the absolute accuracy of the optimally blended scheme is O(p-2) and O(p-1) times better than that of the pure finite element and spectral element schemes respectively.
NASA Astrophysics Data System (ADS)
Bekki, Naoaki; Shintani, Seine A.
2015-12-01
We consider the Rayleigh-Lamb-type equation for propagating pulsive waves excited by aortic-valve closure at end-systole in the human heart wall. We theoretically investigate the transcendental dispersion equation of pulsive waves for the asymmetrical zero-order mode of the Lamb wave. We analytically find a simple dispersion equation with a universal constant for a small Lamb wavenumber. We show that the simple dispersion equation can qualitatively explain the myocardial noninvasive measurements in vivo of pulsive waves in the human heart wall. We can also consistently estimate the viscoelastic constant of the myocardium in the human heart wall using the simple dispersion equation for a small Lamb wavenumber instead of using a complex nonlinear optimization.
Construction of the wave operator for non-linear dispersive equations
NASA Astrophysics Data System (ADS)
Tsuruta, Kai Erik
In this thesis, we will study non-linear dispersive equations. The primary focus will be on the construction of the positive-time wave operator for such equations. The positive-time wave operator problem arises in the study of the asymptotics of a partial differential equation. It is a map from a space of initial data X into itself, and is loosely defined as follows: Suppose that for a solution ψlin to the dispersive equation with no non-linearity and initial data ψ +, there exists a unique solution ψ to the non-linear equation with initial data ψ0 such that ψ behaves as ψ lin as t → infinity. Then the wave operator is the map W+ that takes ψ + to ψ0. By its definition, W+ is injective. An important additional question is whether or not the map is also surjective. If so, then every non-linear solution emanating from X behaves, in some sense, linearly as it evolves (this is known as asymptotic completeness). Thus, there is some justification for treating these solutions as their much simpler linear counterparts. The main results presented in this thesis revolve around the construction of the wave operator(s) at critical non-linearities. We will study the "semi-relativistic" Schrodinger equation as well as the Klein-Gordon-Schrodinger system on R2 . In both cases, we will impose fairly general quadratic non-linearities for which conservation laws cannot be relied upon. These non-linearities fall below the scaling required to employ such tools as the Strichartz estimates. We instead adapt the "first iteration method" of Jang, Li, and Zhang to our setting which depends crucially on the critical decay of the non-linear interaction of the linear evolution. To see the critical decay in our problem, careful analysis is needed to treat the regime where one has spatial and/or time resonance.
MAGNETIC ADVECTION DUE TO DIFFUSIVITY GRADIENTS
NASA Astrophysics Data System (ADS)
Zita, E. J.
2009-12-01
We derive and discuss an important source of advection of magnetic fields in plasmas, for a completely general case. Magnetic diffusivity is proportional to electrical resistivity: where the value this parameter is high, it is well known that magnetic fields can leak (or diffuse) rapidly into (or out) of the plasma. Magnetohydrodynamic lore has it that where gradients, or changes in space, of the value of the diffusivity are high, magnetic fields can have enhanced flow (or advection). We derive this phenomenon rigorously, compare our results to other work in the literature, and discuss its implications, especially for kinematic dynamos. As an extra mathematical bonus, we find that the magnetic advection due to diffusivity gradients can be expressed in terms of a diffusion equation within the induction equation, making its computational implementation especially simple.
NASA Astrophysics Data System (ADS)
Hamdi, Adel
2009-11-01
This paper deals with the identification of a point source (localization of its position and recovering the history of its time-varying intensity function) that constitutes the right-hand side of the first equation in a system of two coupled 1D linear transport equations. Assuming that the source intensity function vanishes before reaching the final control time, we prove the identifiability of the sought point source from recording the state relative to the second coupled transport equation at two observation points framing the source region. Note that at least one of the two observation points should be strategic. We establish an identification method that uses these records to identify the source position as the root of a continuous and strictly monotonic function. Whereas the source intensity function is recovered using a recursive formula without any need of an iterative process. Some numerical experiments on a variant of the surface water pollution BOD-OD coupled model are presented.
Tracer dispersion simulation in low wind speed conditions with a new 2D Langevin equation system
NASA Astrophysics Data System (ADS)
Anfossi, D.; Alessandrini, S.; Trini Castelli, S.; Ferrero, E.; Oettl, D.; Degrazia, G.
The simulation of atmospheric dispersion in low wind speed conditions (LW) is still recognised as a challenge for modellers. Recently, a new system of two coupled Langevin equations that explicitly accounts for meandering has been proposed. It is based on the study of turbulence and dispersion properties in LW. The new system was implemented in the Lagrangian stochastic particle models LAMBDA and GRAL. In this paper we present simulations with this new approach applying it to the tracer experiments carried out in LW by Idaho National Engineering Laboratory (INEL, USA) in 1974 and by the Graz University of Technology and CNR-Torino near Graz in 2003. To assess the improvement obtained with the present model with respect to previous models not taking into account the meandering effect, the simulations for the INEL experiments were also performed with the old version of LAMBDA. The results of the comparisons clearly indicate that the new approach improves the simulation results.
NASA Astrophysics Data System (ADS)
Cuperman, S.; Heristchi, D.
1992-08-01
The transcendental dispersion equation for electromagnetic waves propagating in the slow mode in sheared non-neutral relativistic cylindrical electron beams in strong applied magnetic fields is solved exactly. Thus, rather than truncated power series for the modified Bessel functions involved, use is made of modern algorithms able to compute such functions up to 18-digit accuracy. Consequently, new and significantly more important branches of the velocity shear instability are found. When the shear-factor and/or the geometrical parameter a/b (pipe-to-beam radius ratio) are increased, the unstable branches join, and the higher-frequency, larger-wavenumber modes are significantly enhanced. Since analytical solutions of the exact dispersion relation cannot be obtained, it is suggested that in all similar cases the methods proposed and demonstrated here should be used to carry out a rigorous stability analysis.
Adaptive domain decomposition methods for advection-diffusion problems
Carlenzoli, C.; Quarteroni, A.
1995-12-31
Domain decomposition methods can perform poorly on advection-diffusion equations if diffusion is dominated by advection. Indeed, the hyperpolic part of the equations could affect the behavior of iterative schemes among subdomains slowing down dramatically their rate of convergence. Taking into account the direction of the characteristic lines we introduce suitable adaptive algorithms which are stable with respect to the magnitude of the convective field in the equations and very effective on bear boundary value problems.
Modelling drug degradation in a spray dried polymer dispersion using a modified Arrhenius equation.
Patterson, Adele; Ferreira, Ana P; Banks, Elizabeth; Skeene, Kirsty; Clarke, Graham; Nicholson, Sarah; Rawlinson-Malone, Clare
2015-01-15
The Pharmaceutical industry is increasingly utilizing amorphous technologies to overcome solubility challenges. A common approach is the use of drug in polymer dispersions to prevent recrystallization of the amorphous drug. Understanding the factors affecting chemical and physical degradation of the drug within these complex systems, e.g., temperature and relative humidity, is an important step in the selection of a lead formulation, and development of appropriate packaging/storage control strategies. The Arrhenius equation has been used as the basis of a number of models to predict the chemical stability of formulated product. In this work, we investigate the increase in chemical degradation seen for one particular spray dried dispersion formulation using hydroxypropyl methylcellulose acetate succinate (HPMC-AS). Samples, prepared using polymers with different substitution levels, were placed on storage for 6 months under a range of different temperature and relative humidity conditions and the degradant level monitored using high-performance liquid chromatography (HPLC). While the data clearly illustrates the impact of temperature and relative humidity on the degradant levels detected, it also highlighted that these terms do not account for all the variability in the data. An extension of the Arrhenius equation to include a term for the polymer chemistry, specifically the degree of succinoyl substitution on the polymer backbone, was shown to improve the fit of the model to the data. PMID:25450477
NASA Astrophysics Data System (ADS)
Tóth, Gábor; Keppens, Rony
2012-07-01
The Versatile Advection Code (VAC) is a freely available general hydrodynamic and magnetohydrodynamic simulation software that works in 1, 2 or 3 dimensions on Cartesian and logically Cartesian grids. VAC runs on any Unix/Linux system with a Fortran 90 (or 77) compiler and Perl interpreter. VAC can run on parallel machines using either the Message Passing Interface (MPI) library or a High Performance Fortran (HPF) compiler.
Classical non-Markovian Boltzmann equation
Alexanian, Moorad
2014-08-01
The modeling of particle transport involves anomalous diffusion, (x²(t) ) ∝ t{sup α} with α ≠ 1, with subdiffusive transport corresponding to 0 < α < 1 and superdiffusive transport to α > 1. These anomalies give rise to fractional advection-dispersion equations with memory in space and time. The usual Boltzmann equation, with only isolated binary collisions, is Markovian and, in particular, the contributions of the three-particle distribution function are neglected. We show that the inclusion of higher-order distribution functions give rise to an exact, non-Markovian Boltzmann equation with resulting transport equations for mass, momentum, and kinetic energy with memory in both time and space. The two- and the three-particle distribution functions are considered under the assumption that the two- and the three-particle correlation functions are translationally invariant that allows us to obtain advection-dispersion equations for modeling transport in terms of spatial and temporal fractional derivatives.
Optimal dispersion with minimized Poisson equations for non-hydrostatic free surface flows
NASA Astrophysics Data System (ADS)
Cui, Haiyang; Pietrzak, J. D.; Stelling, G. S.
2014-09-01
A non-hydrostatic shallow-water model is proposed to simulate the wave propagation in situations where the ratio of the wave length to the water depth is small. It exploits the reduced-size stencil in the Poisson pressure solver to make the model less expensive in terms of memory and CPU time. We refer to this new technique as the minimized Poisson equations formulation. In the simplest case when the method applied to a two-layer model, the new model requires the same computational effort as depth-integrated non-hydrostatic models, but can provide a much better description of dispersive waves. To allow an easy implementation of the new method in depth-integrated models, the governing equations are transformed into a depth-integrated system, in which the velocity difference serves as an extra variable. The non-hydrostatic shallow-water model with minimized Poisson equations formulation produces good results in a series of numerical experiments, including a standing wave in a basin, a non-linear wave test, solitary wave propagation in a channel and a wave propagation over a submerged bar.
Antidiffusive velocities for multipass donor cell advection
Margolin, L.G. ); Smolarkiewicz, P.K. )
1989-12-01
Smolarkiewicz describes an iterative process for approximating the advection equation. Basically, he uses a donor cell approximation to correct for the truncation error of the originally specified donor cell scheme. This step may be repeated an arbitrary number of times leading to successively more accurate solutions to the advection equation. In this report, we show how to sum the successive approximations analytically to find a single antidiffusive velocity that represents the effects of an arbitrary number of passes. The analysis is first done dimension to illustrate the method. The analysis is then repeated in two dimensions. The existence of cross terms in the truncation analysis of the two-dimensional equations introduces an extra complication into the calculation. We discuss the implementation of our new antidiffusive velocities and provide some examples of applications. 6 refs., 5 figs., 4 tabs.
Anomalous scaling of a scalar field advected by turbulence
Kraichnan, R.H.
1995-12-31
Recent work leading to deduction of anomalous scaling exponents for the inertial range of an advected passive field from the equations of motion is reviewed. Implications for other turbulence problems are discussed.
NASA Astrophysics Data System (ADS)
Li, Wan-Tong; Wang, Jia-Bing; Zhang, Li
2016-08-01
This paper is concerned with the new types of entire solutions other than traveling wave solutions of nonlocal dispersal equations with monostable nonlinearity in space periodic habitats. We first establish the existence and properties of spatially periodic solutions connecting two steady states. Then new types of entire solutions are constructed by combining the rightward and leftward pulsating traveling fronts with different speeds and a spatially periodic solution. Finally, for a class of special heterogeneous reaction, we further establish the uniqueness of entire solutions and the continuous dependence of such an entire solution on parameters, such as wave speeds and the shifted variables. In other words, we build a five-dimensional manifold of solutions and the traveling wave solutions are on the boundary of the manifold.
Equations of state for crystalline zirconium iodide: The role of dispersion
NASA Astrophysics Data System (ADS)
Rossi, Matthew L.; Taylor, Christopher D.
2013-02-01
We present the first-principle equations of state of several zirconium iodides, ZrI2, ZrI3, and ZrI4, computed using density functional theory methods that apply various methods for introducing the dispersion correction. Iodides formed due to reaction of molecular or atomic iodine with zirconium and zircaloys are of particular interest due to their application to the cladding material used in the fabrication of nuclear fuel rods. Stress corrosion cracking (SCC), associated with fission product chemistry with the clad material, is a major concern in the life cycle of nuclear fuels, as many of the observed rod failures have occurred due to pellet-cladding chemical interactions (PCCI) [A. Atrens, G. Dannhäuser, G. Bäro, Stress-corrosion-cracking of zircaloy-4 cladding tubes, Journal of Nuclear Materials 126 (1984) 91-102; P. Rudling, R. Adamson, B. Cox, F. Garzarolli, A. Strasser, High burn-up fuel issues, Nuclear Engineering and Technology 40 (2008) 1-8]. A proper understanding of the physical properties of the corrosion products is, therefore, required for the development of a comprehensive SCC model. In this particular work, we emphasize that, while existing modeling techniques include methods to compute crystal structures and associated properties, it is important to capture intermolecular forces not traditionally included, such as van der Waals (dispersion) correction. Furthermore, crystal structures with stoichiometries favoring a high I:Zr ratio are found to be particularly sensitive, such that traditional density functional theory approaches that do not incorporate dispersion incorrectly predict significantly larger volumes of the lattice. This latter point is related to the diffuse nature of the iodide electron cloud.
Hydrodynamic dispersion within porous biofilms.
Davit, Y; Byrne, H; Osborne, J; Pitt-Francis, J; Gavaghan, D; Quintard, M
2013-01-01
Many microorganisms live within surface-associated consortia, termed biofilms, that can form intricate porous structures interspersed with a network of fluid channels. In such systems, transport phenomena, including flow and advection, regulate various aspects of cell behavior by controlling nutrient supply, evacuation of waste products, and permeation of antimicrobial agents. This study presents multiscale analysis of solute transport in these porous biofilms. We start our analysis with a channel-scale description of mass transport and use the method of volume averaging to derive a set of homogenized equations at the biofilm-scale in the case where the width of the channels is significantly smaller than the thickness of the biofilm. We show that solute transport may be described via two coupled partial differential equations or telegrapher's equations for the averaged concentrations. These models are particularly relevant for chemicals, such as some antimicrobial agents, that penetrate cell clusters very slowly. In most cases, especially for nutrients, solute penetration is faster, and transport can be described via an advection-dispersion equation. In this simpler case, the effective diffusion is characterized by a second-order tensor whose components depend on (1) the topology of the channels' network; (2) the solute's diffusion coefficients in the fluid and the cell clusters; (3) hydrodynamic dispersion effects; and (4) an additional dispersion term intrinsic to the two-phase configuration. Although solute transport in biofilms is commonly thought to be diffusion dominated, this analysis shows that hydrodynamic dispersion effects may significantly contribute to transport. PMID:23410370
Multi-moment advection scheme in three dimension for Vlasov simulations of magnetized plasma
Minoshima, Takashi; Matsumoto, Yosuke; Amano, Takanobu
2013-03-01
We present an extension of the multi-moment advection scheme [T. Minoshima, Y. Matsumoto, T. Amano, Multi-moment advection scheme for Vlasov simulations, Journal of Computational Physics 230 (2011) 6800–6823] to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.
A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory
NASA Astrophysics Data System (ADS)
Stolk, Christiaan C.
2016-06-01
We develop a new dispersion minimizing compact finite difference scheme for the Helmholtz equation in 2 and 3 dimensions. The scheme is based on a newly developed ray theory for difference equations. A discrete Helmholtz operator and a discrete operator to be applied to the source and the wavefields are constructed. Their coefficients are piecewise polynomial functions of hk, chosen such that phase and amplitude errors are minimal. The phase errors of the scheme are very small, approximately as small as those of the 2-D quasi-stabilized FEM method and substantially smaller than those of alternatives in 3-D, assuming the same number of gridpoints per wavelength is used. In numerical experiments, accurate solutions are obtained in constant and smoothly varying media using meshes with only five to six points per wavelength and wave propagation over hundreds of wavelengths. When used as a coarse level discretization in a multigrid method the scheme can even be used with down to three points per wavelength. Tests on 3-D examples with up to 108 degrees of freedom show that with a recently developed hybrid solver, the use of coarser meshes can lead to corresponding savings in computation time, resulting in good simulation times compared to the literature.
NASA Astrophysics Data System (ADS)
Zhang, Lijun; Chen, Li-Qun; Zhang, Jianming
2013-10-01
Bifurcation and exact solutions of the modified nonlinearly dispersive mK (m,n,k) equation with nonlinear dispersion um-1ut+a(un)x+b(uk)xxx = 0,nk≠0 are investigated in this paper. As a result, under different parameter conditions, abundant compactons, peakons and solitary solutions including not only some known results but also some new ones are obtained. We also point out the original reason of the existence of the non-smooth traveling wave solutions. The approach we used here is also suitable for the study of traveling wave solutions of some other nonlinear equations.
2001-01-01
The ability of species to migrate and disperse is a trait that has interested ecologists for many years. Now that so many species and ecosystems face major environmental threats from habitat fragmentation and global climate change, the ability of species to adapt to these changes by dispersing, migrating, or moving between patches of habitat can be crucial to ensuring their survival. This book provides a timely and wide-ranging overview of the study of dispersal and incorporates much of the latest research. The causes, mechanisms, and consequences of dispersal at the individual, population, species and community levels are considered. The potential of new techniques and models for studying dispersal, drawn from molecular biology and demography, is also explored. Perspectives and insights are offered from the fields of evolution, conservation biology and genetics. Throughout the book, theoretical approaches are combined with empirical data, and care has been taken to include examples from as wide a range of species as possible.
NASA Technical Reports Server (NTRS)
Wang, Xiao-Yen; Chow, Chuen-Yen; Chang, Sin-Chung
1999-01-01
Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.
Numerical dispersion analysis for three-dimensional Laplace-Fourier-domain scalar wave equation
NASA Astrophysics Data System (ADS)
Chen, Jing-Bo
2016-06-01
Based on the phase velocity and attenuation propagation velocity, a method for performing numerical dispersion analysis of three-dimensional Laplace-Fourier-domain scalar wave equation is presented. This method is applied to a 27-point average-derivative optimal scheme and a 27-point finite-element scheme. Within the relative error of 1%, the 27-point average-derivative optimal scheme requires seven grid points per wavelength and pseudo-wavelength while the 27-point finite-element scheme requires 23 grid points per wavelength and pseudo-wavelength for equal and unequal directional sampling intervals. Numerical examples show that the 27-point Laplace-Fourier-domain average-derivative optimal scheme is more accurate than the 27-point Laplace-Fourier-domain finite-element scheme for the same computational cost. By using larger directional sampling intervals while maintaining accuracy, the 27-point Laplace-Fourier-domain average-derivative optimal scheme can greatly reduce the computational cost of three-dimensional Laplace-Fourier-domain modelling.
Non-dispersive carrier transport in molecularly doped polymers and the convection-diffusion equation
NASA Astrophysics Data System (ADS)
Tyutnev, A. P.; Parris, P. E.; Saenko, V. S.
2015-08-01
We reinvestigate the applicability of the concept of trap-free carrier transport in molecularly doped polymers and the possibility of realistically describing time-of-flight (TOF) current transients in these materials using the classical convection-diffusion equation (CDE). The problem is treated as rigorously as possible using boundary conditions appropriate to conventional time of flight experiments. Two types of pulsed carrier generation are considered. In addition to the traditional case of surface excitation, we also consider the case where carrier generation is spatially uniform. In our analysis, the front electrode is treated as a reflecting boundary, while the counter electrode is assumed to act either as a neutral contact (not disturbing the current flow) or as an absorbing boundary at which the carrier concentration vanishes. As expected, at low fields transient currents exhibit unusual behavior, as diffusion currents overwhelm drift currents to such an extent that it becomes impossible to determine transit times (and hence, carrier mobilities). At high fields, computed transients are more like those typically observed, with well-defined plateaus and sharp transit times. Careful analysis, however, reveals that the non-dispersive picture, and predictions of the CDE contradict both experiment and existing disorder-based theories in important ways, and that the CDE should be applied rather cautiously, and even then only for engineering purposes.
LAYER DEPENDENT ADVECTION IN CMAQ
The advection methods used in CMAQ require that the Courant-Friedrichs-Lewy (CFL) condition be satisfied for numerical stability and accuracy. In CMAQ prior to version 4.3, the ADVSTEP algorithm established CFL-safe synchronization and advection timesteps that were uniform throu...
Theory of advection-driven long range biotic transport
Technology Transfer Automated Retrieval System (TEKTRAN)
We propose a simple mechanistic model to examine the effects of advective flow on the spread of fungal diseases spread by wind-blown spores. The model is defined by a set of two coupled non-linear partial differential equations for spore densities. One equation describes the long-distance advectiv...
NASA Astrophysics Data System (ADS)
Remiddi, Ettore; Tancredi, Lorenzo
2016-06-01
It is shown that the study of the imaginary part and of the corresponding dispersion relations of Feynman graph amplitudes within the differential equations method can provide a powerful tool for the solution of the equations, especially in the massive case. The main features of the approach are illustrated by discussing the simple cases of the 1-loop self-mass and of a particular vertex amplitude, and then used for the evaluation of the two-loop massive sunrise and the QED kite graph (the problem studied by Sabry in 1962), up to first order in the (d - 4) expansion.
Technology Transfer Automated Retrieval System (TEKTRAN)
Different parts of soil solution move with different velocities, and therefore chemicals are leached gradually from soil with infiltrating water. Solute dispersivity is the soil parameter characterizing this phenomenon. To characterize the dispersivity of soil profile at field scale, it is desirable...
Andrea Prosperetti
2006-03-24
The report briefly describes the activities carried out in the course of the project. A first line of research was the development of systematic closure relations for averaged equations for disperse multiphase flow. A second line was the development of efficient numerical methods for the simulation of Navier-Stokes flows with many suspended particles. The report also lists the 21 journal articles in which this work is more fully decsribed.
NASA Astrophysics Data System (ADS)
Künzli, Pierre; Tsunematsu, Kae; Albuquerque, Paul; Falcone, Jean-Luc; Chopard, Bastien; Bonadonna, Costanza
2016-04-01
Volcanic ash transport and dispersal models typically describe particle motion via a turbulent velocity field. Particles are advected inside this field from the moment they leave the vent of the volcano until they deposit on the ground. Several techniques exist to simulate particles in an advection field such as finite difference Eulerian, Lagrangian-puff or pure Lagrangian techniques. In this paper, we present a new flexible simulation tool called TETRAS (TEphra TRAnsport Simulator) based on a hybrid Eulerian-Lagrangian model. This scheme offers the advantages of being numerically stable with no numerical diffusion and easily parallelizable. It also allows us to output particle atmospheric concentration or ground mass load at any given time. The model is validated using the advection-diffusion analytical equation. We also obtained a good agreement with field observations of the tephra deposit associated with the 2450 BP Pululagua (Ecuador) and the 1996 Ruapehu (New Zealand) eruptions. As this kind of model can lead to computationally intensive simulations, a parallelization on a distributed memory architecture was developed. A related performance model, taking into account load imbalance, is proposed and its accuracy tested.
Anisotropic Turbulent Advection of a Passive Vector Field: Effects of the Finite Correlation Time
NASA Astrophysics Data System (ADS)
Antonov, N. V.; Gulitskiy, N. M.
2016-02-01
The turbulent passive advection under the environment (velocity) field with finite correlation time is studied. Inertial-range asymptotic behavior of a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow, is investigated by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and prescribed pair correlation function. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to nontrivial fixed points of the RG equations and depend on the relation between the exponents in the energy energy spectrum ɛ ∝ k⊥1-ξ and the dispersion law ω ∝ k⊥2-η . The corresponding anomalous exponents are associated with the critical dimensions of tensor composite operators built solely of the passive vector field itself. In contrast to the well-known isotropic Kraichnan model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the dependence on the integral turbulence scale L has a logarithmic behavior: instead of power-like corrections to ordinary scaling, determined by naive (canonical) dimensions, the anomalies manifest themselves as polynomials of logarithms of L. Due to the presence of the anisotropy in the model, all multiloop diagrams are equal to zero, thus this result is exact.
NASA Astrophysics Data System (ADS)
Klein, C.; Peter, R.
2015-06-01
We present a detailed numerical study of solutions to general Korteweg-de Vries equations with critical and supercritical nonlinearity, both in the context of dispersive shocks and blow-up. We study the stability of solitons and show that they are unstable against being radiated away and blow-up. In the L2 critical case, the blow-up mechanism by Martel, Merle and Raphaël can be numerically identified. In the limit of small dispersion, it is shown that a dispersive shock always appears before an eventual blow-up. In the latter case, always the first soliton to appear will blow up. It is shown that the same type of blow-up as for the perturbations of the soliton can be observed which indicates that the theory by Martel, Merle and Raphaël is also applicable to initial data with a mass much larger than the soliton mass. We study the scaling of the blow-up time t∗ in dependence of the small dispersion parameter ɛ and find an exponential dependence t∗(ɛ) and that there is a minimal blow-up time t0∗ greater than the critical time of the corresponding Hopf solution for ɛ → 0. To study the cases with blow-up in detail, we apply the first dynamic rescaling for generalized Korteweg-de Vries equations. This allows to identify the type of the singularity.
High Order Semi-Lagrangian Advection Scheme
NASA Astrophysics Data System (ADS)
Malaga, Carlos; Mandujano, Francisco; Becerra, Julian
2014-11-01
In most fluid phenomena, advection plays an important roll. A numerical scheme capable of making quantitative predictions and simulations must compute correctly the advection terms appearing in the equations governing fluid flow. Here we present a high order forward semi-Lagrangian numerical scheme specifically tailored to compute material derivatives. The scheme relies on the geometrical interpretation of material derivatives to compute the time evolution of fields on grids that deform with the material fluid domain, an interpolating procedure of arbitrary order that preserves the moments of the interpolated distributions, and a nonlinear mapping strategy to perform interpolations between undeformed and deformed grids. Additionally, a discontinuity criterion was implemented to deal with discontinuous fields and shocks. Tests of pure advection, shock formation and nonlinear phenomena are presented to show performance and convergence of the scheme. The high computational cost is considerably reduced when implemented on massively parallel architectures found in graphic cards. The authors acknowledge funding from Fondo Sectorial CONACYT-SENER Grant Number 42536 (DGAJ-SPI-34-170412-217).
NASA Astrophysics Data System (ADS)
Matsuno, Yoshimasa
2014-03-01
We develop a direct method for solving a modified Camassa-Holm equation with cubic nonlinearity and linear dispersion under the rapidly decreasing boundary condition. We obtain a compact parametric representation for the multisoliton solutions and investigate their properties. We show that the introduction of a linear dispersive term exhibits various new features in the structure of solutions. In particular, we find the smooth solitons whose characteristics are different from those of the Camassa-Holm equation, as well as the novel types of singular solitons. A remarkable feature of the soliton solutions is that the underlying structure of the associated tau-functions is the same as that of a model equation for shallow-water waves introduced by Ablowitz et al (1974 Stud. Appl. Math. 53 249-315). Finally, we demonstrate that the short-wave limit of the soliton solutions recovers the soliton solutions of the short pulse equation which describes the propagation of ultra-short optical pulses in nonlinear media.
Feynman-Kac Equation for Convection-Dispersion with Mobile and Immobile Walkers
NASA Astrophysics Data System (ADS)
Choquet, C.; Neel, M. C.
2012-07-01
Some experiments for determining transport properties are based on distributions data of paths integrals of particles undergoing displacement. Unfortunately occurrence of abnormally long immobile periods renders the celebrated Feynman-Kac equation inappropriate to rule the joint density of positions and path integrals of walkers. Nevertheless, a natural partition of the distribution of particles (mobile and immobile) may then be reunified by a fractional integral operator. An adequate Feynman-Kac type's equation follows from mass conservation principle.
Advection around ventilated U-shaped burrows: A model study
NASA Astrophysics Data System (ADS)
Brand, Andreas; Lewandowski, JöRg; Hamann, Enrico; Nützmann, Gunnar
2013-05-01
Advective transport in the porous matrix of sediments surrounding burrows formed by fauna such as Chironomus plumosus has been generally neglected. A positron emission tomography study recently revealed that the pumping activity of the midge larvae can indeed induce fluid flow in the sediment. We present a numerical model study which explores the conditions at which advective transport in the sediment becomes relevant. A 0.15 m deep U-shaped burrow with a diameter of 0.002 m within the sediment was represented in a 3-D domain. Fluid flow in the burrow was calculated using the Navier-Stokes equation for incompressible laminar flow in the burrow, and flow in the sediment was described by Darcy's law. Nonreactive and reactive transport scenarios were simulated considering diffusion and advection. The pumping activity of the model larva results in considerable advective flow in the sediment at reasonable high permeabilities with flow velocities of up to 7.0 × 10-6 m s-1 close to the larva for a permeability of 3 × 10-12 m2. At permeabilities below 7 × 10-13 m2 advection is negligible compared to diffusion. Reactive transport simulations using first-order kinetics for oxygen revealed that advective flux into the sediment downstream of the pumping larva enhances sedimentary uptake, while the advective flux into the burrow upstream of the larvae inhibits diffusive sedimentary uptake. Despite the fact that both effects cancel each other with respect to total solute uptake, the advection-induced asymmetry in concentration distribution can lead to a heterogeneous solute and redox distribution in the sediment relevant to complex reaction networks.
FRACTIONAL SOLUTE TRANSPORT EQUATION EVALUATED WITH THE MISCIBLE DISPLACEMENT EXPERIMENTAL DATA
Technology Transfer Automated Retrieval System (TEKTRAN)
A new solute transport model has been recently developed assuming that the movements of solute particles in hierarchically-structured porous media belongs to the family of Lévy motions rather than to the Brownian motion. The one-dimensional fractional advective-dispersive transport equation, or FADE...
Classification of dispersion equations for homogeneous, dielectric-magnetic, uniaxial materials
NASA Astrophysics Data System (ADS)
Depine, Ricardo A.; Inchaussandague, Marina E.; Lakhtakia, Akhlesh
2006-04-01
The geometric representation at a fixed frequency of the wave vector (or dispersion) surface ω(k) for lossless, homogeneous, dielectric-magnetic uniaxial materials is explored for the case when the elements of the relative permittivity and permeability tensors of the material can have any sign. Electromagnetic plane waves propagating inside the material can exhibit dispersion surfaces in the form of ellipsoids of revolution, hyperboloids of one sheet, or hyperboloids of two sheets. Furthermore, depending on the relative orientation of the optic axis, the intersections of these surfaces with fixed planes of propagation can be circles, ellipses, hyperbolas, or straight lines. The understanding obtained is used to study the reflection and refraction of electromagnetic plane waves due to a planar interface with an isotropic medium.
Classification of dispersion equations for homogeneous, dielectric-magnetic, uniaxial materials.
Depine, Ricardo A; Inchaussandague, Marina E; Lakhtakia, Akhlesh
2006-04-01
The geometric representation at a fixed frequency of the wave vector (or dispersion) surface omega(k) for lossless, homogeneous, dielectric-magnetic uniaxial materials is explored for the case when the elements of the relative permittivity and permeability tensors of the material can have any sign. Electromagnetic plane waves propagating inside the material can exhibit dispersion surfaces in the form of ellipsoids of revolution, hyperboloids of one sheet, or hyperboloids of two sheets. Furthermore, depending on the relative orientation of the optic axis, the intersections of these surfaces with fixed planes of propagation can be circles, ellipses, hyperbolas, or straight lines. The understanding obtained is used to study the reflection and refraction of electromagnetic plane waves due to a planar interface with an isotropic medium. PMID:16604780
Dark solitons in the Lugiato-Lefever equation with normal dispersion
NASA Astrophysics Data System (ADS)
Parra-Rivas, P.; Knobloch, E.; Gomila, D.; Gelens, L.
2016-06-01
The regions of existence and stability of dark solitons in the Lugiato-Lefever model with normal chromatic dispersion are described. These localized states are shown to be organized in a bifurcation structure known as collapsed snaking implying the presence of a region in parameter space with a finite multiplicity of dark solitons. For some parameter values dynamical instabilities are responsible for the appearance of oscillations and temporal chaos. The importance of the results for understanding frequency comb generation in microresonators is emphasized.
NASA Astrophysics Data System (ADS)
Charalampidis, E. G.; Kevrekidis, P. G.; Frantzeskakis, D. J.; Malomed, B. A.
2016-08-01
We consider a two-component, two-dimensional nonlinear Schrödinger system with unequal dispersion coefficients and self-defocusing nonlinearities, chiefly with equal strengths of the self- and cross-interactions. In this setting, a natural waveform with a nonvanishing background in one component is a vortex, which induces an effective potential well in the second component, via the nonlinear coupling of the two components. We show that the potential well may support not only the fundamental bound state, but also multiring excited radial state complexes for suitable ranges of values of the dispersion coefficient of the second component. We systematically explore the existence, stability, and nonlinear dynamics of these states. The complexes involving the excited radial states are weakly unstable, with a growth rate depending on the dispersion of the second component. Their evolution leads to transformation of the multiring complexes into stable vortex-bright solitons ones with the fundamental state in the second component. The excited states may be stabilized by a harmonic-oscillator trapping potential, as well as by unequal strengths of the self- and cross-repulsive nonlinearities.
Markel, Vadim A; Schotland, John C
2012-06-01
We consider the problem of homogenizing the Maxwell equations for periodic composites. The analysis is based on Bloch-Floquet theory. We calculate explicitly the reflection coefficient for a half space and derive and implement a computationally efficient continued-fraction expansion for the effective permittivity. Our results are illustrated by numerical computations for the case of two-dimensional systems. The homogenization theory of this paper is designed to predict various physically measurable quantities rather than to simply approximate certain coefficients in a partial differential equation. PMID:23005233
Spectral approximation to advection-diffusion problems by the fictitious interface method
NASA Astrophysics Data System (ADS)
Frati, A.; Pasquarelli, F.; Quarteroni, A.
1993-08-01
The algorithmic aspects of the 'fictitious interface' method used in numerical approximations of convection-dominated flows are discussed. The solution algorithm presented alternates the advection-equation solution with that of the advection-diffusion equation within complementary subdomains. For the problems presently considered, spatial discretization is obtained by the spectral collocation method via Legendre-Gaussian modes. Attention is given to the the fictitious interface method's application to the Burgers equation.
NASA Astrophysics Data System (ADS)
Bona, G.; Chen, J. A.; Saut, Jing Ping
2002-08-01
Considered herein are a number of variants of the classical Boussinesq system and their higher-order generalizations. Such equations were first derived by Boussinesq to describe the two-way propagation of small-amplitude, long wavelength, gravity waves on the surface of water in a canal. These systems arise also when modeling the propagation of long-crested waves on large lakes or the ocean and in other contexts. Depending on the modeling of dispersion, the resulting system may or may not have a linearization about the rest state which is well posed. Even when well posed, the linearized system may exhibit a lack of conservation of energy that is at odds with its status as an approximation to the Euler equations. In the present script, we derive a four-parameter family of Boussinesq systems from the two-dimensional Euler equations for free-surface flow and formulate criteria to help decide which of these equations one might choose in a given modeling situation. The analysis of the systems according to these criteria is initiated.
Experiments in Advective and Turbulent Hyporheic Pumping
NASA Astrophysics Data System (ADS)
Mccluskey, A. H.; Grant, S.; Stewardson, M. J.
2014-12-01
Hyporheic exchange (HE) is the mixing of stream and subsurface waters beneath the sediment-water interface (SWI). At the patch and reach scales, HE is dominated by periodic upwelling and downwelling zones, induced by pressure variation and processes within the turbulent boundary layer (TBL). This can be caused by (1) the geometry of the stream, imposing a stationary wave at the SWI or (2) by a travelling wave associated with the propagation of turbulent pressure waves generated from the TBL. Case (1) has generally been the favoured model of hyporheic exchange and has been referred to as hyporheic 'pumping' by Elliott and Brooks, and subsequently others. Case (2) can be termed turbulent pumping, and has been proposed as a mechanism to model the combined effects of turbulent dispersion alongside steady-state advection. While this has been represented numerically and analytically, conjecture remains about the physical representation of these combined processes. We present initial results from experiments undertaken to classify the spatial and temporal characteristics of pressure variation at and beneath the SWI, with a periodic sinusoidal geometry of wavelength 0.28m and height 0.02m. As an initial characterisation, the advective flow profile has been examined using time-lapse photography of dyes released across the span of a periodic downwelling zone. These tracer tests confirmed delineation of isolated upwelling and downwelling cells as noted by previous authors in modelling studies. However, their distribution deviates from the typical pumping pattern with increased discharge and stream gradient. Empirical orthogonal function (EOF) analysis of high frequency (250Hz) pressure measurements, sampled at an array along the centroid of the flume underneath one wavelength gave further insight into the spatial distribution of turbulent signatures arising from roughness-generated turbulence. A turbulent frequency of 6-10Hz dominates, however the penetration depth appears to
NASA Astrophysics Data System (ADS)
Haghi, M.; Anand, L.
1990-01-01
An experimental study of the constitutive response of the oxide dispersion-strengthened (ODS) superalloy MA 956, which consists of an Fe-Cr-Al matrix dispersion strengthened with yttria, has been performed. Single-crystal specimens of MA 956 having remarkably simple initial microstructures have been tested in compression in the temperature range of 900 °C to 1200 °C and in the axial strain-rate range of 1.8 x 10-4 s-1 to 10-2 s-1. The deformation response of the material has been examined by performing constant true strain-rate tests, strain-rate jump tests, and stress relaxation tests. The orientation dependence of the stress-strain response of the single crystals has been compensated for by determining the operative slip systems and resolving the stresses and strains accordingly. These experiments, together with electron-microscopic observations of deformed and quenched specimens, allow a number of conclusions to be drawn about the physics of particle strengthening in this simple ODS alloy at high temperatures. Further, drawing on this physical understanding, a set of phenomenological internal variable constitutive equations which model the high-temperature deformation behavior of this alloy is also developed. These equations reasonably well model not only the temperature and strain-rate sensitivity of the flow stress but also the strain-hardening behavior of the material.
Feng, Xin; Ye, Xingyou; Park, Jun-Bom; Lu, Wenli; Morott, Joe; Beissner, Brad; Lian, Zhuoyang John; Pinto, Elanor; Bi, Vivian; Porter, Stu; Durig, Tom; Majumdar, Soumyajit; Repka, Michael A
2015-01-01
The recrystallization of an amorphous drug in a solid dispersion system could lead to a loss in the drug solubility and bioavailability. The primary objective of the current research was to use an improved kinetic model to evaluate the recrystallization kinetics of amorphous structures and to further understand the factors influencing the physical stability of amorphous solid dispersions. Amorphous solid dispersions of fenofibrate with different molecular weights of hydroxypropylcellulose, HPC (Klucel™ LF, EF, ELF) were prepared utilizing hot-melt extrusion technology. Differential scanning calorimetry was utilized to quantitatively analyze the extent of recrystallization in the samples stored at different temperatures and relative humidity (RH) conditions. The experimental data were fitted into the improved kinetics model of a modified Avrami equation to calculate the recrystallization rate constants. Klucel LF, the largest molecular weight among the HPCs used, demonstrated the greatest inhibition of fenofibrate recrystallization. Additionally, the recrystallization rate (k) decreased with increasing polymer content, however exponentially increased with higher temperature. Also k increased linearly rather than exponentially over the range of RH studied. PMID:25224341
Waves, advection, and cloud patterns on Venus
NASA Technical Reports Server (NTRS)
Schinder, Paul J.; Gierasch, Peter J.; Leroy, Stephen S.; Smith, Michael D.
1990-01-01
The stable layers adjacent to the nearly neutral layer within the Venus clouds are found to be capable of supporting vertically trapped, horizontally propagating waves with horizontal wavelengths of about 10 km and speeds of a few meters per second relative to the mean wind in the neutral layer. These waves may possibly be excited by turbulence within the neutral layer. Here, the properties of the waves, and the patterns which they might produce within the visible clouds if excited near the subsolar point are examined. The patterns can be in agreement with many features in images. The waves are capable of transferring momentum latitudinally to help maintain the general atmospheric spin, but at present we are not able to evaluate wave amplitudes. We also examine an alternative possibility that the cloud patterns are produced by advection and shearing by the mean zonal and meridional flow of blobs formed near the equator. It is concluded that advection and shearing by the mean flow is the most likely explanation for the general pattern of small scale striations.
Dispersive transport of charge carriers in disordered nanostructured materials
NASA Astrophysics Data System (ADS)
Sibatov, R. T.; Uchaikin, V. V.
2015-07-01
Dispersive transport of charge carriers in disordered nanostructured semiconductors is described in terms of integral diffusion equations nonlocal in time. Transient photocurrent kinetics is analyzed for different situations. Relation to the fractional differential approach is demonstrated. Using this relation provides specifications in interpretation of the time-of-flight data. Joint influence of morphology and energy distribution of localized states is described in frames of the trap-limited advection-diffusion on a comb structure modeling a percolation cluster.
Arditi, Roger; Lobry, Claude; Sari, Tewfik
2015-12-01
The standard model for the dynamics of a fragmented density-dependent population is built from several local logistic models coupled by migrations. First introduced in the 1970s and used in innumerable articles, this standard model applied to a two-patch situation has never been completely analysed. Here, we complete this analysis and we delineate the conditions under which fragmentation associated to dispersal is either beneficial or detrimental to total population abundance. Therefore, this is a contribution to the SLOSS question. Importantly, we also show that, depending on the underlying mechanism, there is no unique way to generalize the logistic model to a patchy situation. In many cases, the standard model is not the correct generalization. We analyse several alternative models and compare their predictions. Finally, we emphasize the shortcomings of the logistic model when written in the r-K parameterization and we explain why Verhulst's original polynomial expression is to be preferred. PMID:26472064
NASA Astrophysics Data System (ADS)
Cardenas, M. Bayani
2009-12-01
The transition from non-Fickian to Fickian macroscale transport is explicitly demonstrated for an increasing array of three-dimensional pores with vortices in between a lattice of cubic packed spheres by microscale finite element Navier-Stokes flow and transport simulations. Solute residence time distribution begins with a power law for one pore but gradually and eventually transforms to an exponential distribution typical of classic dispersive transport after about ten pores. Parameter fitting of an analytical solution to the 1-D advection-dispersion equation using the simulated breakthrough curves leads to fitted pore velocities within 1% of actual values and an asymptotic fitted dispersion coefficient after a few pores. Therefore, after dozens of pores, bulk transport can be described by the advection-dispersion equation. Persistent vortices in similarly structured porous media subjected to similar grain-scale Reynolds and Peclet numbers may have minimal contribution to anomalous transport observed at larger scales.
NASA Astrophysics Data System (ADS)
Zhang, Guo-Bao; Ma, Ruyun
2014-10-01
This paper is concerned with the traveling wave solutions and the spreading speeds for a nonlocal dispersal equation with convolution-type crossing-monostable nonlinearity, which is motivated by an age-structured population model with time delay. We first prove the existence of traveling wave solution with critical wave speed c = c*. By introducing two auxiliary monotone birth functions and using a fluctuation method, we further show that the number c = c* is also the spreading speed of the corresponding initial value problem with compact support. Then, the nonexistence of traveling wave solutions for c < c* is established. Finally, by means of the (technical) weighted energy method, we prove that the traveling wave with large speed is exponentially stable, when the initial perturbation around the wave is relatively small in a weighted norm.
Giraldo, F.X.; Neta, B.
1995-04-21
An Eulerian and semi-Lagrangian finite element methods for the solution of the two dimensional advection equation were developed. Bilinear rectangular elements were used. Linear stability analysis of the method is given.
Technology Transfer Automated Retrieval System (TEKTRAN)
Solute transport in soils and sediments is commonly simulated with the parabolic advective-dispersive equation, or ADE. In the last decades, it has been reported that this model cannot take in account several important features of solute movement through soil. Recently, a new model base on the assu...
NASA Astrophysics Data System (ADS)
Bona, J. L.; Chen, M.; Saut, J.-C.
2004-05-01
In part I of this work (Bona J L, Chen M and Saut J-C 2002 Boussinesq equations and other systems for small-amplitude long waves in nonlinear dispersive media I: Derivation and the linear theory J. Nonlinear Sci. 12 283-318), a four-parameter family of Boussinesq systems was derived to describe the propagation of surface water waves. Similar systems are expected to arise in other physical settings where the dominant aspects of propagation are a balance between the nonlinear effects of convection and the linear effects of frequency dispersion. In addition to deriving these systems, we determined in part I exactly which of them are linearly well posed in various natural function classes. It was argued that linear well-posedness is a natural necessary requirement for the possible physical relevance of the model in question. In this paper, it is shown that the first-order correct models that are linearly well posed are in fact locally nonlinearly well posed. Moreover, in certain specific cases, global well-posedness is established for physically relevant initial data. In part I, higher-order correct models were also derived. A preliminary analysis of a promising subclass of these models shows them to be well posed.
Consistency problem with tracer advection in the Atmospheric Model GAMIL
NASA Astrophysics Data System (ADS)
Zhang, Kai; Wan, Hui; Wang, Bin; Zhang, Meigen
2008-03-01
The radon transport test, which is a widely used test case for atmospheric transport models, is carried out to evaluate the tracer advection schemes in the Grid-Point Atmospheric Model of IAP-LASG (GAMIL). Two of the three available schemes in the model are found to be associated with significant biases in the polar regions and in the upper part of the atmosphere, which implies potentially large errors in the simulation of ozone-like tracers. Theoretical analyses show that inconsistency exists between the advection schemes and the discrete continuity equation in the dynamical core of GAMIL and consequently leads to spurious sources and sinks in the tracer transport equation. The impact of this type of inconsistency is demonstrated by idealized tests and identified as the cause of the aforementioned biases. Other potential effects of this inconsistency are also discussed. Results of this study provide some hints for choosing suitable advection schemes in the GAMIL model. At least for the polar-region-concentrated atmospheric components and the closely correlated chemical species, the Flux-Form Semi-Lagrangian advection scheme produces more reasonable simulations of the large-scale transport processes without significantly increasing the computational expense.
NASA Astrophysics Data System (ADS)
Ramadan, Omar
2015-09-01
In this paper, systematic wave-equation finite difference time domain (WE-FDTD) formulations are presented for modeling electromagnetic wave-propagation in linear and nonlinear dispersive materials. In the proposed formulations, the complex conjugate pole residue (CCPR) pairs model is adopted in deriving a unified dispersive WE-FDTD algorithm that allows modeling different dispersive materials, such as Debye, Drude and Lorentz, in the same manner with the minimal additional auxiliary variables. Moreover, the proposed formulations are incorporated with the wave-equation perfectly matched layer (WE-PML) to construct a material independent mesh truncating technique that can be used for modeling general frequency-dependent open region problems. Several numerical examples involving linear and nonlinear dispersive materials are included to show the validity of the proposed formulations.
Numerical simulation of life cycles of advection warm fog
NASA Technical Reports Server (NTRS)
Hung, R. J.; Vaughan, O. H.
1977-01-01
The formation, development and dissipation of advection warm fog is investigated. The equations employed in the model include the equation of continuity, momentum and energy for the descriptions of density, wind component and potential temperature, respectively, together with two diffusion equations for the modification of water-vapor mixing ratio and liquid-water mixing ratios. A description of the vertical turbulent transfer of heat, moisture and momentum has been taken into consideration. The turbulent exchange coefficients adopted in the model are based on empirical flux-gradient relations.
Spiral defect chaos in an advection-reaction-diffusion system
NASA Astrophysics Data System (ADS)
Affan, H.; Friedrich, R.
2014-06-01
This paper comprises numerical and theoretical studies of spatiotemporal patterns in advection-reaction-diffusion systems in which the chemical species interact with the hydrodynamic fluid. Due to the interplay between the two, we obtained the spiral defect chaos in the activator-inhibitor-type model. We formulated the generalized Swift-Hohenberg-type model for this system. Then the evolution of fractal boundaries due to the effect of the strong nonlinearity at the interface of the two chemical species is studied numerically. The purpose of the present paper is to point out that spiral defect chaos, observed in model equations of the extended Swift-Hohenberg equation for low Prandtl number convection, may actually be obtained also in certain advection-reaction-diffusion systems.
Laser speckle contrast imaging is sensitive to advective flux
NASA Astrophysics Data System (ADS)
Khaksari, Kosar; Kirkpatrick, Sean J.
2016-07-01
Unlike laser Doppler flowmetry, there has yet to be presented a clear description of the physical variables that laser speckle contrast imaging (LSCI) is sensitive to. Herein, we present a theoretical basis for demonstrating that LSCI is sensitive to total flux and, in particular, the summation of diffusive flux and advective flux. We view LSCI from the perspective of mass transport and briefly derive the diffusion with drift equation in terms of an LSCI experiment. This equation reveals the relative sensitivity of LSCI to both diffusive flux and advective flux and, thereby, to both concentration and the ordered velocity of the scattering particles. We demonstrate this dependence through a short series of flow experiments that yield relationships between the calculated speckle contrast and the concentration of the scatterers (manifesting as changes in scattering coefficient), between speckle contrast and the velocity of the scattering fluid, and ultimately between speckle contrast and advective flux. Finally, we argue that the diffusion with drift equation can be used to support both Lorentzian and Gaussian correlation models that relate observed contrast to the movement of the scattering particles and that a weighted linear combination of these two models is likely the most appropriate model for relating speckle contrast to particle motion.
NASA Astrophysics Data System (ADS)
Beech, Robert; Osman, Frederick
2005-10-01
This paper will present the nonlinearity and dispersion effects involved in propagation of optical solitons, which can be understood by using a numerical routine to solve the nonlinear Schrödinger equation (NLSE). Here, Mathematica v5© (Wolfram, 2003) is used to explore in depth several features of optical solitons formation and propagation. These numerical routines were implemented through the use of Mathematica v5© and the results give a very clear idea of this interesting and important practical phenomenon. It is hoped that this work will open up an important new approach to the cause, effect, and correction of interference from secondary radiation found in the uses of soliton waves in lasers and in optical fiber telecommunication. It is believed that these results will be of considerable use in any work or research in this field and in self-focusing properties of the soliton (Osman et al., 2004a, 2004b; Hora, 1991). In a previous paper on this topic (Beech & Osman, 2004), it was shown that solitons of NLSE radiate. This paper goes on from there to show that these radiations only occur in solitons derived from cubic, or odd-numbered higher orders of NLSE, and that there are no such radiations from solitons of quadratic, or even-numbered higher order of NLSE. It is anticipated that this will stimulate research into practical means to control or eliminate such radiations.
Comparison of thermal advection measurements by clear-air radar and radiosonde techniques
Crochet, M.; Rougier, G.; Bazile, G. Meteorologie Nationale, Trappes )
1990-10-01
Vertical profiles of the horizontal wind have been measured every 4 min by a clear-air radar (stratospheric-troposphere radar), and vertical profiles of temperature have been obtained every 2 hours by three radiosonde soundings in the same zone in Brittany during the Mesoscale Frontal Dynamics Project FRONTS 87 campaign. Radar thermal advection is deduced from the thermal wind equation using the measured real horizontal wind instead of the geostrophic wind. Radiosonde thermal advection is determined directly from the sounding station data sets of temperature gradients and also approximately from the thermodynamic equation by the temperature tendency. These approximations, applied during a frontal passage, show the same general features and magnitude of the thermal advection, giving a preliminary but encouraging conclusion for a possible real-time utilization of clear-air radars to monitor thermal advection and to identify its characteristic features. 6 refs.
NASA Astrophysics Data System (ADS)
Moura, R. C.; Sherwin, S. J.; Peiró, J.
2016-02-01
This study addresses linear dispersion-diffusion analysis for the spectral/hp continuous Galerkin (CG) formulation in one dimension. First, numerical dispersion and diffusion curves are obtained for the advection-diffusion problem and the role of multiple eigencurves peculiar to spectral/hp methods is discussed. From the eigencurves' behaviour, we observe that CG might feature potentially undesirable non-smooth dispersion/diffusion characteristics for under-resolved simulations of problems strongly dominated by either convection or diffusion. Subsequently, the linear advection equation augmented with spectral vanishing viscosity (SVV) is analysed. Dispersion and diffusion characteristics of CG with SVV-based stabilization are verified to display similar non-smooth features in flow regions where convection is much stronger than dissipation or vice-versa, owing to a dependency of the standard SVV operator on a local Péclet number. First a modification is proposed to the traditional SVV scaling that enforces a globally constant Péclet number so as to avoid the previous issues. In addition, a new SVV kernel function is suggested and shown to provide a more regular behaviour for the eigencurves along with a consistent increase in resolution power for higher-order discretizations, as measured by the extent of the wavenumber range where numerical errors are negligible. The dissipation characteristics of CG with the SVV modifications suggested are then verified to be broadly equivalent to those obtained through upwinding in the discontinuous Galerkin (DG) scheme. Nevertheless, for the kernel function proposed, the full upwind DG scheme is found to have a slightly higher resolution power for the same dissipation levels. These results show that improved CG-SVV characteristics can be pursued via different kernel functions with the aid of optimization algorithms.
Vertical Structure of Advection-dominated Accretion Flows
NASA Astrophysics Data System (ADS)
Zahra Zeraatgari, Fateme; Abbassi, Shahram
2015-08-01
We solve the set of hydrodynamic equations for optically thin advection-dominated accretion flows by assuming a radially self-similar spherical coordinate system (r,θ ,φ ). The disk is considered to be in steady state and axisymmetric. We define the boundary conditions at the pole and the equator of the disk and, to avoid singularity at the rotation axis, the disk is taken to be symmetric with respect to this axis. Moreover, only the {τ }rφ component of the viscous stress tensor is assumed, and we have set {v}θ =0. The main purpose of this study is to investigate the variation of dynamical quantities of the flow in the vertical direction by finding an analytical solution. As a consequence, we found that the advection parameter, {f}{adv}, varies along the θ direction and reaches its maximum near the rotation axis. Our results also show that, in terms of the no-outflow solution, thermal equilibrium still exists and consequently advection cooling can balance viscous heating.
NASA Astrophysics Data System (ADS)
Hornby, P. G.
2005-12-01
Understanding chemical and thermal processes taking place in hydrothermal mineral deposition systems could well be a key to unlocking new mineral reserves through improved targeting of exploration efforts. To aid in this understanding it is very helpful to be able to model such processes with sufficient fidelity to test process hypotheses. To gain understanding, it is often sufficient to obtain semi-quantitative results that model the broad aspects of the complex set of thermal and chemical effects taking place in hydrothermal systems. For example, it is often sufficient to gain an understanding of where thermal, geometric and chemical factors converge to precipitate gold (say) without being perfectly precise about how much gold is precipitated. The traditional approach is to use incompressible Darcy flow together with the Boussinesq approximation. From the flow field, the heat equation is used to advect-conduct the heat. The flow field is also used to transport solutes by solving an advection-dispersion-diffusion equation. The reactions in the fluid and between fluid and rock act as source terms for these advection-dispersion equations. Many existing modelling systems that are used for simulating such systems use explicit time marching schemes and finite differences. The disadvantage of this approach is the need to work on rectilinear grids and the number of time steps required by the Courant condition in the solute transport step. The second factor can be particularly significant if the chemical system is complex, requiring (at a minimum) an equilibrium calculation at each grid point at each time step. In the approach we describe, we use finite elements rather than finite differences, and the pressure, heat and advection-dispersion equations are solved implicitly. The general idea is to put unconditional numerical stability of the time integration first, and let accuracy assume a secondary role. It is in this sense that the method is semi-quantiative. However
A fully implicit method for 3D quasi-steady state magnetic advection-diffusion.
Siefert, Christopher; Robinson, Allen Conrad
2009-09-01
We describe the implementation of a prototype fully implicit method for solving three-dimensional quasi-steady state magnetic advection-diffusion problems. This method allows us to solve the magnetic advection diffusion equations in an Eulerian frame with a fixed, user-prescribed velocity field. We have verified the correctness of method and implementation on two standard verification problems, the Solberg-White magnetic shear problem and the Perry-Jones-White rotating cylinder problem.
Advection-Induced Spectrotemporal Defects in a Free-Electron Laser
Bielawski, S.; Szwaj, C.; Bruni, C.; Garzella, D.; Orlandi, G.L.; Couprie, M.E.
2005-07-15
We evidence numerically and experimentally that advection can induce spectrotemporal defects in a system presenting a localized structure. Those defects in the spectrum are associated with the breakings induced by the drift of the localized solution. The results are based on simulations and experiments performed on the super-ACO free-electron laser. However, we show that this instability can be generalized using a real Ginzburg-Landau equation with (i) advection and (ii) a finite-size supercritical region.
NASA Technical Reports Server (NTRS)
Mazaheri, Alireza; Ricchiuto, Mario; Nishikawa, Hiroaki
2016-01-01
In this paper, we introduce a new hyperbolic first-order system for general dispersive partial differential equations (PDEs). We then extend the proposed system to general advection-diffusion-dispersion PDEs. We apply the fourth-order RD scheme of Ref. 1 to the proposed hyperbolic system, and solve time-dependent dispersive equations, including the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted results, including the gradient and Hessian (second derivative), are in a very good agreement with the exact solutions. We then show that the RD scheme applied to the proposed system accurately captures dispersive shocks without numerical oscillations. We also verify that the solution, gradient and Hessian are predicted with equal order of accuracy.
Dispersive processes in models of regional radionuclide migration. Technical memorandum
Evenson, D.E.; Dettinger, M.D.
1980-05-01
Three broad areas of concern in the development of aquifer scale transport models will be local scale diffusion and dispersion processes, regional scale dispersion processes, and numerical problems associated with the advection-dispersion equation. Local scale dispersion processes are fairly well understood and accessible to observation. These processes will generally be dominated in large scale systems by regional processes, or macro-dispersion. Macro-dispersion is primarily the result of large scale heterogeneities in aquifer properties. In addition, the effects of many modeling approximations are often included in the process. Because difficulties arise in parameterization of this large scale phenomenon, parameterization should be based on field measurements made at the same scale as the transport process of interest or else partially circumvented through the application of a probabilistic advection model. Other problems associated with numerical transport models include difficulties with conservation of mass, stability, numerical dissipation, overshoot, flexibility, and efficiency. We recommend the random-walk model formulation for Lawrence Livermore Laboratory's purposes as the most flexible, accurate and relatively efficient modeling approach that overcomes these difficulties.
Super-diffusion versus competitive advection: a simulation
NASA Astrophysics Data System (ADS)
Del Moro, D.; Giannattasio, F.; Berrilli, F.; Consolini, G.; Lepreti, F.; Gošić, M.
2015-04-01
Context. Magnetic element tracking is often used to study the transport and diffusion of the magnetic field on the solar photosphere. From the analysis of the displacement spectrum of these tracers, it has recently been agreed that a regime of super-diffusivity dominates the solar surface. Quite habitually this result is discussed in the framework of fully developed turbulence. Aims: However, the debate whether the super-diffusivity is generated by a turbulent dispersion process, by the advection due to the convective pattern, or even by another process is still open, as is the question of the amount of diffusivity at the scales relevant to the local dynamo process. Methods: To understand how such peculiar diffusion in the solar atmosphere takes place, we compared the results from two different data sets (ground-based and space-borne) and developed a simulation of passive tracers advection by the deformation of a Voronoi network. Results: The displacement spectra of the magnetic elements obtained by the data sets are consistent in retrieving a super-diffusive regime for the solar photosphere, but the simulation also shows a super-diffusive displacement spectrum: its competitive advection process can reproduce the signature of super-diffusion. Conclusions: Therefore, it is not necessary to hypothesize a totally developed turbulence regime to explain the motion of the magnetic elements on the solar surface.
NASA Astrophysics Data System (ADS)
Chen, Yong; Yan, Zhenya
2016-03-01
Solitons are of the important significant in many fields of nonlinear science such as nonlinear optics, Bose-Einstein condensates, plamas physics, biology, fluid mechanics, and etc. The stable solitons have been captured not only theoretically and experimentally in both linear and nonlinear Schrödinger (NLS) equations in the presence of non-Hermitian potentials since the concept of the parity-time -symmetry was introduced in 1998. In this paper, we present novel bright solitons of the NLS equation with third-order dispersion in some complex -symmetric potentials (e.g., physically relevant -symmetric Scarff-II-like and harmonic-Gaussian potentials). We find stable nonlinear modes even if the respective linear -symmetric phases are broken. Moreover, we also use the adiabatic changes of the control parameters to excite the initial modes related to exact solitons to reach stable nonlinear modes. The elastic interactions of two solitons are exhibited in the third-order NLS equation with -symmetric potentials. Our results predict the dynamical phenomena of soliton equations in the presence of third-order dispersion and -symmetric potentials arising in nonlinear fiber optics and other physically relevant fields.
Chen, Yong; Yan, Zhenya
2016-01-01
Solitons are of the important significant in many fields of nonlinear science such as nonlinear optics, Bose-Einstein condensates, plamas physics, biology, fluid mechanics, and etc. The stable solitons have been captured not only theoretically and experimentally in both linear and nonlinear Schrödinger (NLS) equations in the presence of non-Hermitian potentials since the concept of the parity-time -symmetry was introduced in 1998. In this paper, we present novel bright solitons of the NLS equation with third-order dispersion in some complex -symmetric potentials (e.g., physically relevant -symmetric Scarff-II-like and harmonic-Gaussian potentials). We find stable nonlinear modes even if the respective linear -symmetric phases are broken. Moreover, we also use the adiabatic changes of the control parameters to excite the initial modes related to exact solitons to reach stable nonlinear modes. The elastic interactions of two solitons are exhibited in the third-order NLS equation with -symmetric potentials. Our results predict the dynamical phenomena of soliton equations in the presence of third-order dispersion and -symmetric potentials arising in nonlinear fiber optics and other physically relevant fields. PMID:27002543
Chen, Yong; Yan, Zhenya
2016-01-01
Solitons are of the important significant in many fields of nonlinear science such as nonlinear optics, Bose-Einstein condensates, plamas physics, biology, fluid mechanics, and etc. The stable solitons have been captured not only theoretically and experimentally in both linear and nonlinear Schrödinger (NLS) equations in the presence of non-Hermitian potentials since the concept of the parity-time -symmetry was introduced in 1998. In this paper, we present novel bright solitons of the NLS equation with third-order dispersion in some complex -symmetric potentials (e.g., physically relevant -symmetric Scarff-II-like and harmonic-Gaussian potentials). We find stable nonlinear modes even if the respective linear -symmetric phases are broken. Moreover, we also use the adiabatic changes of the control parameters to excite the initial modes related to exact solitons to reach stable nonlinear modes. The elastic interactions of two solitons are exhibited in the third-order NLS equation with -symmetric potentials. Our results predict the dynamical phenomena of soliton equations in the presence of third-order dispersion and -symmetric potentials arising in nonlinear fiber optics and other physically relevant fields. PMID:27002543
A Quasi-Conservative Adaptive Semi-Lagrangian Advection-Diffusion Scheme
NASA Astrophysics Data System (ADS)
Behrens, Joern
2014-05-01
Many processes in atmospheric or oceanic tracer transport are conveniently represented by advection-diffusion type equations. Depending on the magnitudes of both components, the mathematical representation and consequently the discretization is a non-trivial problem. We will focus on advection-dominated situations and will introduce a semi-Lagrangian scheme with adaptive mesh refinement for high local resolution. This scheme is well suited for pollutant transport from point sources, or transport processes featuring fine filamentation with corresponding local concentration maxima. In order to achieve stability, accuracy and conservation, we combine an adaptive mesh refinement quasi-conservative semi-Lagrangian scheme, based on an integral formulation of the underlying advective conservation law (Behrens, 2006), with an advection diffusion scheme as described by Spiegelman and Katz (2006). The resulting scheme proves to be conservative and stable, while maintaining high computational efficiency and accuracy.
NASA Astrophysics Data System (ADS)
Bodin, Jacques
2015-03-01
In this study, new multi-dimensional time-domain random walk (TDRW) algorithms are derived from approximate one-dimensional (1-D), two-dimensional (2-D), and three-dimensional (3-D) analytical solutions of the advection-dispersion equation and from exact 1-D, 2-D, and 3-D analytical solutions of the pure-diffusion equation. These algorithms enable the calculation of both the time required for a particle to travel a specified distance in a homogeneous medium and the mass recovery at the observation point, which may be incomplete due to 2-D or 3-D transverse dispersion or diffusion. The method is extended to heterogeneous media, represented as a piecewise collection of homogeneous media. The particle motion is then decomposed along a series of intermediate checkpoints located on the medium interface boundaries. The accuracy of the multi-dimensional TDRW method is verified against (i) exact analytical solutions of solute transport in homogeneous media and (ii) finite-difference simulations in a synthetic 2-D heterogeneous medium of simple geometry. The results demonstrate that the method is ideally suited to purely diffusive transport and to advection-dispersion transport problems dominated by advection. Conversely, the method is not recommended for highly dispersive transport problems because the accuracy of the advection-dispersion TDRW algorithms degrades rapidly for a low Péclet number, consistent with the accuracy limit of the approximate analytical solutions. The proposed approach provides a unified methodology for deriving multi-dimensional time-domain particle equations and may be applicable to other mathematical transport models, provided that appropriate analytical solutions are available.
Advective and diffusive cosmic ray transport in galactic haloes
NASA Astrophysics Data System (ADS)
Heesen, Volker; Dettmar, Ralf-Jürgen; Krause, Marita; Beck, Rainer; Stein, Yelena
2016-05-01
We present 1D cosmic ray transport models, numerically solving equations of pure advection and diffusion for the electrons and calculating synchrotron emission spectra. We find that for exponential halo magnetic field distributions advection leads to approximately exponential radio continuum intensity profiles, whereas diffusion leads to profiles that can be better approximated by a Gaussian function. Accordingly, the vertical radio spectral profiles for advection are approximately linear, whereas for diffusion they are of `parabolic' shape. We compare our models with deep Australia Telescope Compact Array observations of two edge-on galaxies, NGC 7090 and 7462, at λλ 22 and 6 cm. Our result is that the cosmic ray transport in NGC 7090 is advection dominated with V=150^{+80}_{-30} km s^{-1}, and that the one in NGC 7462 is diffusion dominated with D=3.0± 1.0 × 10^{28}E_GeV^{0.5} cm^2 s^{-1}. NGC 7090 has both a thin and thick radio disc with respective magnetic field scale heights of hB1 = 0.8 ± 0.1 kpc and hB2 = 4.7 ± 1.0 kpc. NGC 7462 has only a thick radio disc with hB2 = 3.8 ± 1.0 kpc. In both galaxies, the magnetic field scale heights are significantly smaller than what estimates from energy equipartition would suggest. A non-negligible fraction of cosmic ray electrons can escape from NGC 7090, so that this galaxy is not an electron calorimeter.
Shear dispersion in a capillary tube with a porous wall
NASA Astrophysics Data System (ADS)
Dejam, Morteza; Hassanzadeh, Hassan; Chen, Zhangxin
2016-02-01
An analytical expression is presented for the shear dispersion during solute transport in a coupled system comprised of a capillary tube and a porous medium. The dispersion coefficient is derived in a capillary tube with a porous wall by considering an accurate boundary condition, which is the continuity of concentration and mass flux, at the interface between the capillary tube and porous medium. A comparison of the obtained results with that in a non-coupled system identifies three regimes including: diffusion-dominated, transition, and advection-dominated. The results reveal that it is essential to include the exchange of solute between the capillary tube and porous medium in development of the shear dispersion coefficient for the last two regimes. The resulting equivalent transport equation revealed that due to mass transfer between the capillary tube and the porous medium, the dispersion coefficient is decreased while the effective velocity in the capillary tube increases. However, a larger effective advection term leads to faster breakthrough of a solute and enhances mass delivery to the porous medium as compared with the classical double-porosity model with a non-coupled dispersion coefficient. The obtained results also indicate that the finite porous medium gives faster breakthrough of a solute as compared with the infinite one. These results find applications in solute transport in porous capillaries and membranes.
Shear dispersion in a capillary tube with a porous wall.
Dejam, Morteza; Hassanzadeh, Hassan; Chen, Zhangxin
2016-01-01
An analytical expression is presented for the shear dispersion during solute transport in a coupled system comprised of a capillary tube and a porous medium. The dispersion coefficient is derived in a capillary tube with a porous wall by considering an accurate boundary condition, which is the continuity of concentration and mass flux, at the interface between the capillary tube and porous medium. A comparison of the obtained results with that in a non-coupled system identifies three regimes including: diffusion-dominated, transition, and advection-dominated. The results reveal that it is essential to include the exchange of solute between the capillary tube and porous medium in development of the shear dispersion coefficient for the last two regimes. The resulting equivalent transport equation revealed that due to mass transfer between the capillary tube and the porous medium, the dispersion coefficient is decreased while the effective velocity in the capillary tube increases. However, a larger effective advection term leads to faster breakthrough of a solute and enhances mass delivery to the porous medium as compared with the classical double-porosity model with a non-coupled dispersion coefficient. The obtained results also indicate that the finite porous medium gives faster breakthrough of a solute as compared with the infinite one. These results find applications in solute transport in porous capillaries and membranes. PMID:26845232
Antonov, N V; Gulitskiy, N M
2015-10-01
In this work we study the generalization of the problem considered in [Phys. Rev. E 91, 013002 (2015)] to the case of finite correlation time of the environment (velocity) field. The model describes a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow. Inertial-range asymptotic behavior is studied by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and preassigned pair correlation function. Due to the presence of distinguished direction n, all the multiloop diagrams in this model vanish, so that the results obtained are exact. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to the two nontrivial fixed points of the RG equations. Their stability depends on the relation between the exponents in the energy spectrum E∝k(⊥)(1-ξ) and the dispersion law ω∝k(⊥)(2-η). In contrast to the well-known isotropic Kraichnan's model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the corrections to ordinary scaling are polynomials of logarithms of the integral turbulence scale L. PMID:26565343
NASA Astrophysics Data System (ADS)
Antonov, N. V.; Gulitskiy, N. M.
2015-10-01
In this work we study the generalization of the problem considered in [Phys. Rev. E 91, 013002 (2015), 10.1103/PhysRevE.91.013002] to the case of finite correlation time of the environment (velocity) field. The model describes a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow. Inertial-range asymptotic behavior is studied by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and preassigned pair correlation function. Due to the presence of distinguished direction n , all the multiloop diagrams in this model vanish, so that the results obtained are exact. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to the two nontrivial fixed points of the RG equations. Their stability depends on the relation between the exponents in the energy spectrum E ∝k⊥1 -ξ and the dispersion law ω ∝k⊥2 -η . In contrast to the well-known isotropic Kraichnan's model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the corrections to ordinary scaling are polynomials of logarithms of the integral turbulence scale L .
NASA Astrophysics Data System (ADS)
Lissy, Pierre
2015-11-01
In this paper, we prove explicit lower bounds for the cost of fast boundary controls for a class of linear equations of parabolic or dispersive type involving the spectral fractional Laplace operator. We notably deduce the following striking result: in the case of the heat equation controlled on the boundary, Miller's conjecture formulated in Miller (2004) [16] is not verified. Moreover, we also give a new lower bound for the minimal time needed to ensure the uniform controllability of the one-dimensional convection-diffusion equation with negative speed controlled on the left boundary, proving that the conjecture formulated in Coron and Guerrero (2005) [2] concerning this problem is also not verified at least for negative speeds. The proof is based on complex analysis, and more precisely on a representation formula for entire functions of exponential type, and is quite related to the moment method.
Burgers turbulence and passive random advection
NASA Astrophysics Data System (ADS)
Boldyrev, Stanislav Anatolievich
1999-10-01
The thesis is devoted to development of new methods in the theory of strong turbulence. These methods are illustrated with the so-called Burgers' model of turbulence, i.e., the Navier-Stokes equation without pressure, supplemented by a Gaussian, short-time correlated external force. The main goal of the theory is to describe the statistics of the velocity field. Since the Navier-Stokes equation is nonlinear, the problem is highly nontrivial; it is sometimes referred to as the ``Ising model'' of strong turbulence. The importance of the problem for plasma physics, astrophysics, physics of self-organized criticality, disordered systems, etc., is discussed in Chapter 1. In this thesis a new self-consistent theoretical approach to the problem is developed. The problem is treated from the field-theoretical point of view, and, therefore, appropriate methods such as regularization, operator product expansion, and an assumption about scaling invariance are employed. The scheme ``from particular to general'' is adopted. The main ideas of the approach are first developed in detail for the one-dimensional Burgers model in Chapter 2 and then generalized to the multidimensional case in Chapter 3. In all of the cases the velocity- difference and velocity-gradient probability density functions are obtained. Their derivation is based on the self-consistent conjecture about the operator product expansion for the dissipative term, introduced by Polyakov [1995]. Comparison of the obtained results with the available direct numerical simulations shows a very good agreement. The practically important longitudinal velocity-difference PDF and div v PDF in the multidimensional case are discussed within the approach. In Chapter 4 the statistics of passive quantities (such as temperature, concentration, magnetic field) ``frozen'' into the turbulent fluid are obtained by using the methods developed in Chapters 2 and 3. The velocity field is assumed to be Gaussian, and short-time correlated
Global Solutions of the Boltzmann Equation Over {{R}^D} Near Global Maxwellians with Small Mass
NASA Astrophysics Data System (ADS)
Bardos, Claude; Gamba, Irene M.; Golse, François; Levermore, C. David
2016-07-01
We study the dynamics defined by the Boltzmann equation set in the Euclidean space {{R}^D} in the vicinity of global Maxwellians with finite mass. A global Maxwellian is a special solution of the Boltzmann equation for which the collision integral vanishes identically. In this setting, the dispersion due to the advection operator quenches the dissipative effect of the Boltzmann collision integral. As a result, the large time limit of solutions of the Boltzmann equation in this regime is given by noninteracting, freely transported states and can be described with the tools of scattering theory.
Surfzone alongshore advective accelerations: observations and modeling
NASA Astrophysics Data System (ADS)
Hansen, J.; Raubenheimer, B.; Elgar, S.
2014-12-01
The sources, magnitudes, and impacts of non-linear advective accelerations on alongshore surfzone currents are investigated with observations and a numerical model. Previous numerical modeling results have indicated that advective accelerations are an important contribution to the alongshore force balance, and are required to understand spatial variations in alongshore currents (which may result in spatially variable morphological change). However, most prior observational studies have neglected advective accelerations in the alongshore force balance. Using a numerical model (Delft3D) to predict optimal sensor locations, a dense array of 26 colocated current meters and pressure sensors was deployed between the shoreline and 3-m water depth over a 200 by 115 m region near Duck, NC in fall 2013. The array included 7 cross- and 3 alongshore transects. Here, observational and numerical estimates of the dominant forcing terms in the alongshore balance (pressure and radiation-stress gradients) and the advective acceleration terms will be compared with each other. In addition, the numerical model will be used to examine the force balance, including sources of velocity gradients, at a higher spatial resolution than possible with the instrument array. Preliminary numerical results indicate that at O(10-100 m) alongshore scales, bathymetric variations and the ensuing alongshore variations in the wave field and subsequent forcing are the dominant sources of the modeled velocity gradients and advective accelerations. Additional simulations and analysis of the observations will be presented. Funded by NSF and ASDR&E.
NASA Astrophysics Data System (ADS)
Hammer, René; Pötz, Walter; Arnold, Anton
2014-01-01
A finite difference scheme is presented for the Dirac equation in (1+1)D. It can handle space- and time-dependent mass and potential terms and utilizes exact discrete transparent boundary conditions (DTBCs). Based on a space- and time-staggered leap-frog scheme it avoids fermion doubling and preserves the dispersion relation of the continuum problem for mass zero (Weyl equation) exactly. Considering boundary regions, each with a constant mass and potential term, the associated DTBCs are derived by first applying this finite difference scheme and then using the Z-transform in the discrete time variable. The resulting constant coefficient difference equation in space can be solved exactly on each of the two semi-infinite exterior domains. Admitting only solutions in l2 which vanish at infinity is equivalent to imposing outgoing boundary conditions. An inverse Z-transformation leads to exact DTBCs in form of a convolution in discrete time which suppress spurious reflections at the boundaries and enforce stability of the whole space-time scheme. An exactly preserved functional for the norm of the Dirac spinor on the staggered grid is presented. Simulations of Gaussian wave packets, leaving the computational domain without reflection, demonstrate the quality of the DTBCs numerically, as well as the importance of a faithful representation of the energy-momentum dispersion relation on a grid.
NASA Technical Reports Server (NTRS)
Goorjian, Peter M.; Taflove, Allen
1992-01-01
The initial results for femtosecond electromagnetic soliton propagation and collision obtained from first principles, i.e., by a direct time integration of Maxwell's equations are reported. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit the modeling of 2D and 3D optical soliton propagation, scattering, and switching from the full-vector Maxwell's equations.
Advection fog formation in a polluted atmosphere
Hung, R.J.; Liaw, G.S.
1981-01-01
Large quantities of atmospheric aerosols with composition SO/sub 4//sup 2 -/, NO/sub 3//sup -/ and NH/sub 4//sup +/ have been detected in highly industrialized areas. The major portions of aerosol products are the results of energy related fuel combustion. Both microphysical and macrophysical processes are considered in investigating the time dependent evolution of the saturation spectra of condensation nuclei associated with both polluted and clean atmospheres during the time periods of advection fog formation. The results show that the condensation nuclei associated with a polluted atmosphere provide more favorable conditions than condensation nuclei associated with a clean atmosphere to produce dense advection fog, and that attaining a certain degree of supersaturation is not necessarily required for the formation of advection fog with condensation nuclei associated with a polluted atmosphere for monodisperse distribution.
Efficient mass transport by optical advection
NASA Astrophysics Data System (ADS)
Kajorndejnukul, Veerachart; Sukhov, Sergey; Dogariu, Aristide
2015-10-01
Advection is critical for efficient mass transport. For instance, bare diffusion cannot explain the spatial and temporal scales of some of the cellular processes. The regulation of intracellular functions is strongly influenced by the transport of mass at low Reynolds numbers where viscous drag dominates inertia. Mimicking the efficacy and specificity of the cellular machinery has been a long time pursuit and, due to inherent flexibility, optical manipulation is of particular interest. However, optical forces are relatively small and cannot significantly modify diffusion properties. Here we show that the effectiveness of microparticle transport can be dramatically enhanced by recycling the optical energy through an effective optical advection process. We demonstrate theoretically and experimentally that this new advection mechanism permits an efficient control of collective and directional mass transport in colloidal systems. The cooperative long-range interaction between large numbers of particles can be optically manipulated to create complex flow patterns, enabling efficient and tunable transport in microfluidic lab-on-chip platforms.
First-Order Hyperbolic System Method for Time-Dependent Advection-Diffusion Problems
NASA Technical Reports Server (NTRS)
Mazaheri, Alireza; Nishikawa, Hiroaki
2014-01-01
A time-dependent extension of the first-order hyperbolic system method for advection-diffusion problems is introduced. Diffusive/viscous terms are written and discretized as a hyperbolic system, which recovers the original equation in the steady state. The resulting scheme offers advantages over traditional schemes: a dramatic simplification in the discretization, high-order accuracy in the solution gradients, and orders-of-magnitude convergence acceleration. The hyperbolic advection-diffusion system is discretized by the second-order upwind residual-distribution scheme in a unified manner, and the system of implicit-residual-equations is solved by Newton's method over every physical time step. The numerical results are presented for linear and nonlinear advection-diffusion problems, demonstrating solutions and gradients produced to the same order of accuracy, with rapid convergence over each physical time step, typically less than five Newton iterations.
Deift, P; Venakides, S; Zhou, X
1998-01-20
This paper extends the steepest descent method for Riemann-Hilbert problems introduced by Deift and Zhou in a critical new way. We present, in particular, an algorithm, to obtain the support of the Riemann-Hilbert problem for leading asymptotics. Applying this extended method to small dispersion KdV (Korteweg-de Vries) equation, we (i) recover the variational formulation of P. D. Lax and C. D. Levermore [(1979) Proc. Natl. Acad. Sci. USA76, 3602-3606] for the weak limit of the solution, (ii) derive, without using an ansatz, the hyperelliptic asymptotic solution of S. Venakides that describes the oscillations; and (iii) are now able to compute the phase shifts, integrating the modulation equations exactly. The procedure of this paper is a version of fully nonlinear geometrical optics for integrable systems. With some additional analysis the theory can provide rigorous error estimates between the solution and its computed asymptotic expression. PMID:11038618
NASA Astrophysics Data System (ADS)
Molinari, Vincenzo; Mostacci, Domiziano
2015-10-01
He-4 is known to become superfluid at very low temperatures. This effect is now generally accepted to be connected with BEC (Bose-Einstein Condensation). The dispersion relation of pressure waves in superfluid He-4 has been determined at 1.1 °K by Yarnell et al., and exhibits a non monotonic behavior-with a maximum and a minimum-usually explained in terms of excitations called rotons, introduced by Landau. In the present work an attempt is made to describe the phenomenon within the Bohmian interpretation of QM. To this end, the effects of the intermolecular potential, taken to be essentially of the Lennard-Jones type modified to account for molecule finiteness, are included as a Vlasov-type self-consistent field. A dispersion relation is found, that is in quite good agreement with Yarnell's curve.
Scalora, Michael; Syrchin, Maxim S; Akozbek, Neset; Poliakov, Evgeni Y; D'Aguanno, Giuseppe; Mattiucci, Nadia; Bloemer, Mark J; Zheltikov, Aleksei M
2005-07-01
A new generalized nonlinear Schrödinger equation describing the propagation of ultrashort pulses in bulk media exhibiting frequency dependent dielectric susceptibility and magnetic permeability is derived and used to characterize wave propagation in a negative index material. The equation has new features that are distinct from ordinary materials (mu=1): the linear and nonlinear coefficients can be tailored through the linear properties of the medium to attain any combination of signs unachievable in ordinary matter, with significant potential to realize a wide class of solitary waves. PMID:16090616
NASA Technical Reports Server (NTRS)
Joseph, Rose M.; Hagness, Susan C.; Taflove, Allen
1991-01-01
The initial results for femtosecond pulse propagation and scattering interactions for a Lorentz medium obtained by a direct time integration of Maxwell's equations are reported. The computational approach provides reflection coefficients accurate to better than 6 parts in 10,000 over the frequency range of dc to 3 x 10 to the 16th Hz for a single 0.2-fs Gaussian pulse incident upon a Lorentz-medium half-space. New results for Sommerfeld and Brillouin precursors are shown and compared with previous analyses. The present approach is robust and permits 2D and 3D electromagnetic pulse propagation directly from the full-vector Maxwell's equations.
Fujioka, J; Espinosa, A
2015-11-01
In this article, we show that if the nonlinear Schrödinger (NLS) equation is generalized by simultaneously taking into account higher-order dispersion, a quintic nonlinearity, and self-steepening terms, the resulting equation is interesting as it has exact soliton solutions which may be (depending on the values of the coefficients) stable or unstable, standard or "embedded," fixed or "moving" (i.e., solitons which advance along the retarded-time axis). We investigate the stability of these solitons by means of a modified version of the Vakhitov-Kolokolov criterion, and numerical tests are carried out to corroborate that these solitons respond differently to perturbations. It is shown that this generalized NLS equation can be derived from a Lagrangian density which contains an auxiliary variable, and Noether's theorem is then used to show that the invariance of the action integral under infinitesimal gauge transformations generates a whole family of conserved quantities. Finally, we study if this equation has the Painlevé property. PMID:26627574
Estimation of the advection effects induced by surface heterogeneities in the surface energy budget
NASA Astrophysics Data System (ADS)
Cuxart, Joan; Wrenger, Burkhard; Martínez-Villagrasa, Daniel; Reuder, Joachim; Jonassen, Marius O.; Jiménez, Maria A.; Lothon, Marie; Lohou, Fabienne; Hartogensis, Oscar; Dünnermann, Jens; Conangla, Laura; Garai, Anirban
2016-07-01
The effect of terrain heterogeneities in one-point measurements is a continuous subject of discussion. Here we focus on the order of magnitude of the advection term in the equation of the evolution of temperature as generated by documented terrain heterogeneities and we estimate its importance as a term in the surface energy budget (SEB), for which the turbulent fluxes are computed using the eddy-correlation method. The heterogeneities are estimated from satellite and model fields for scales near 1 km or broader, while the smaller scales are estimated through direct measurements with remotely piloted aircraft and thermal cameras and also by high-resolution modelling. The variability of the surface temperature fields is not found to decrease clearly with increasing resolution, and consequently the advection term becomes more important as the scales become finer. The advection term provides non-significant values to the SEB at scales larger than a few kilometres. In contrast, surface heterogeneities at the metre scale yield large values of the advection, which are probably only significant in the first centimetres above the ground. The motions that seem to contribute significantly to the advection term in the SEB equation in our case are roughly those around the hectometre scales.
Technology Transfer Automated Retrieval System (TEKTRAN)
Identification of flow and transport processes under natural boundary conditions in field soils is a complex task since most model parameters are not measurable at that scale. However, combining a numerical solution method of the governing flow and transport equations with an inverse optimization al...
SEPs Dropout Events Associated with Advected Interplanetary Magnetic Structures
NASA Astrophysics Data System (ADS)
Bruno, R.; Trenchi, L.; Telloni, D.; D'Amicis, R.; Marcucci, F.; Zurbuchen, T.; Weberg, M. J.
2013-05-01
The intensity profile of energetic particles from impulsive solar flares (SEP) often shows abrupt dropouts affecting all energies simultaneously, without time-dispersion. Part of the community thinks that these modulations are directly related to the presence of magnetic structures with a different magnetic topology advected by the wind, a sort of magnetic flux tubes. During the expansion, following the dynamical interaction between plasma regions travelling at different speed, these structures would be partially tangled up in a sort of spaghetti-like bundle. These flux tubes would be alternatively connected or not connected with the flare site and, consequently, they would be filled or devoid of SEPs. When the observer passes through them, he would observe clear particles dropout signatures. We will report about results from a detailed analysis of SEP events which showed several signatures in the local magnetic field and/or plasma parameters associated with SEP modulations. These findings corroborate the idea of a possible link between these particles events observed at the Earth's orbit and magnetic connection or disconnection of the ambient magnetic field with the flare region at the Sun. We will also discuss the advantages represented by future Solar Orbiter in-situ observations. As a matter of fact, Solar Orbiter, from its orbital vantage point during the quasi corotation phase, will be a priviledged observer of this kind of phenomenon since it will observe the advected structure of the solar wind not yet reprocessed by dynamical interaction due to wind expansion.
3D Flow Visualization Using Texture Advection
NASA Technical Reports Server (NTRS)
Kao, David; Zhang, Bing; Kim, Kwansik; Pang, Alex; Moran, Pat (Technical Monitor)
2001-01-01
Texture advection is an effective tool for animating and investigating 2D flows. In this paper, we discuss how this technique can be extended to 3D flows. In particular, we examine the use of 3D and 4D textures on 3D synthetic and computational fluid dynamics flow fields.
NASA Astrophysics Data System (ADS)
Dubina, Sean Hyun; Wedgewood, Lewis Edward
2016-07-01
Ferrofluids are often favored for their ability to be remotely positioned via external magnetic fields. The behavior of particles in ferromagnetic clusters under uniformly applied magnetic fields has been computationally simulated using the Brownian dynamics, Stokesian dynamics, and Monte Carlo methods. However, few methods have been established that effectively handle the basic principles of magnetic materials, namely, Maxwell's equations. An iterative constraint method was developed to satisfy Maxwell's equations when a uniform magnetic field is imposed on ferrofluids in a heterogeneous Brownian dynamics simulation that examines the impact of ferromagnetic clusters in a mesoscale particle collection. This was accomplished by allowing a particulate system in a simple shear flow to advance by a time step under a uniformly applied magnetic field, then adjusting the ferroparticles via an iterative constraint method applied over sub-volume length scales until Maxwell's equations were satisfied. The resultant ferrofluid model with constraints demonstrates that the magnetoviscosity contribution is not as substantial when compared to homogeneous simulations that assume the material's magnetism is a direct response to the external magnetic field. This was detected across varying intensities of particle-particle interaction, Brownian motion, and shear flow. Ferroparticle aggregation was still extensively present but less so than typically observed.
Fluid dispersion effects on density-driven thermohaline flow and transport in porous media
NASA Astrophysics Data System (ADS)
Jamshidzadeh, Zahra; Tsai, Frank T.-C.; Mirbagheri, Seyed Ahmad; Ghasemzadeh, Hasan
2013-11-01
This study introduces the dispersive fluid flux of total fluid mass to the density-driven flow equation to improve thermohaline modeling of salt and heat transports in porous media. The dispersive fluid flux in the flow equation is derived to account for an additional fluid flux driven by the density gradient and mechanical dispersion. The coupled flow, salt transport and heat transport governing equations are numerically solved by a fully implicit finite difference method to investigate solution changes due to the dispersive fluid flux. The numerical solutions are verified by the Henry problem and the thermal Elder problem under a moderate density effect and by the brine Elder problem under a strong density effect. It is found that increment of the maximum ratio of the dispersive fluid flux to the advective fluid flux results in increasing dispersivity for the Henry problem and the brine Elder problem. The effects of the dispersive fluid flux on salt and heat transports under high density differences and high dispersivities are more noticeable than under low density differences and low dispersivities. Values of quantitative indicators such as the Nusselt number, mass flux, salt mass stored and maximum penetration depth in the brine Elder problem show noticeable changes by the dispersive fluid flux. In the thermohaline Elder problem, the dispersive fluid flux shows a considerable effect on the shape and the number of developed fingers and makes either an upwelling or a downwelling flow in the center of the domain. In conclusion, for the general case that involves strong density-driven flow and transport modeling in porous media, the dispersive fluid flux should be considered in the flow equation.
NASA Astrophysics Data System (ADS)
Jaunet, V.; Collin, E.; Delville, J.
2016-05-01
This article describes a model obtained by applying proper orthogonal decomposition to the advection equation. The resulting set of equations links the POD modes, their temporal and spatial derivatives and the flow convection velocity. It provides a technique to calculate the convection velocity of coherent structures. It follows, from the model, that a priori knowledge of the convection velocity suffices to construct a dynamical model of the flow. This is demonstrated using experimental data.
On diagonalization of coupled hydrologic transport and geochemical reaction equations
Yeh, Gour-Tsyh; Cheng, Hwai-Ping
1996-12-31
Two basic ingredients present in modeling the transport of reactive multi-components: the transport is described by a set of advection-dispersion-reactive partial differential equations (PDEs) based on the principle of mass balance; the chemical reactions, under the assumptions of local equilibrium, are described by a set of highly nonlinear algebraic equations (AEs) base on the principles of mole balance and mass action. For a typical application, the complete set of nonlinear PDEs and AEs consist of more than one hundred simultaneous equations. Thus, it is impractical to solve this set of equations simultaneously. General practice is to divide this set of equations into two subsets: one is the primary governing equations (PGEs) consisting of mainly the transport equations and the other one is the secondary governing equations consisting of mainly the geochemical reaction equations. The PGEs are solved for the chosen primary dependent variables (PDVs) and the SGEs are used to compute for the secondary dependent variables (SDVs). The major difficulties in simulating the reactive transport is the numerical solution of PGEs. From the computational point of view, the solution of the set of highly nonlinear PDEs are solved either with the direct substitution approach (DSA) or with the sequential iteration approach (SIA). For DSA, geochemical equilibrium reaction equations are substituted into the hydrologic transport equations to results in a set of nonlinear partial differential equations.
NASA Astrophysics Data System (ADS)
Boybeyi, Zafer
A mesoscale atmospheric dispersion modeling system has been developed to investigate mesoscale circulations and associated air pollution dispersion, including effects of terrain topography, large water bodies and urban areas. The system is based on a three-dimensional mesoscale meteorological model coupled with two dispersion models (an Eulerian dispersion model and a Lagrangian particle dispersion model). The mesoscale model is hydrostatic and based on primitive equations formulated in a terrain-following coordinate system with a E-varepsilon turbulence closure scheme. The Eulerian dispersion model is based on numerical solution of the advection-diffusion equation to allow one to simulate releases of non-buoyant pollutants (especially from area and volume sources). The Lagrangian particle dispersion model allows one to simulate releases of buoyant pollutants from arbitrary sources (particularly from point and line sources). The air pollution dispersion models included in the system are driven by the meteorological information provided by the mesoscale model. Mesoscale atmospheric circulations associated with sea and lake breezes have been examined using the mesoscale model. A series of model sensitivity studies were performed to investigate the effects of different environmental parameters on these circulations. It was found that the spatial and temporal variation of the sea and lake breeze convergence zones and the associated convective activities depend to a large extent on the direction and the magnitude of the ambient wind. Dispersion of methyl isocyanate gas from the Bhopal accident was investigated using the mesoscale atmospheric dispersion modeling system. A series of numerical experiments were performed to investigate the possible role of the mesoscale circulations on this industrial gas episode. The temporal and spatial variations of the wind and turbulence fields were simulated with the mesoscale model. The dispersion characteristics of the accidental
Evolution of Advection Upstream Splitting Method Schemes
NASA Technical Reports Server (NTRS)
Liou, Meng-Sing
2010-01-01
This paper focuses on the evolution of advection upstream splitting method(AUSM) schemes. The main ingredients that have led to the development of modern computational fluid dynamics (CFD) methods have been reviewed, thus the ideas behind AUSM. First and foremost is the concept of upwinding. Second, the use of Riemann problem in constructing the numerical flux in the finite-volume setting. Third, the necessity of including all physical processes, as characterised by the linear (convection) and nonlinear (acoustic) fields. Fourth, the realisation of separating the flux into convection and pressure fluxes. The rest of this review briefly outlines the technical evolution of AUSM and more details can be found in the cited references. Keywords: Computational fluid dynamics methods, hyperbolic systems, advection upstream splitting method, conservation laws, upwinding, CFD
Passive advection in a collisionless plasma
NASA Astrophysics Data System (ADS)
Kanekar, Anjor; Schekochihin, Alexander; Hammett, Greg; Dorland, William; Loureiro, Nuno
2014-10-01
We consider a simple kinetic model for the evolution of the particle distribution function in a magnetized turbulent plasma that includes both phase mixing (Landau damping) and advection by a stochastic velocity field: a ``kinetic passive scalar'' in the Batchelor regime. The advection due to stochastic velocity field allows for a stochastic version of the plasma echo by coupling the ``phase-mixing'' and the ``un-phase-mixing'' components of the free energy. We have developed a new analytical framework to diagnose the efficiency of such coupling. We have also developed a new GPU code named Gandalf that solves this kinetic model numerically. In this poster, we shall present numerical and analytical results related to this work.
Advection of methane in the hydrate zone: model, analysis and examples
NASA Astrophysics Data System (ADS)
Peszynska, Malgorzata; Showalter, Ralph E.; Webster, Justin T.
2015-12-01
A two-phase two-component model is formulated for the advective-diffusive transport of methane in liquid phase through sediment with the accompanying formation and dissolution of methane hydrate. This free-boundary problem has a unique generalized solution in $L^1$; the proof combines analysis of the stationary semilinear elliptic Dirichlet problem with the nonlinear semigroup theory in Banach space for an m-accretive multi-valued operator. Additional estimates of maximum principle type are obtained, and these permit appropriate maximal extensions of the phase-change relations. An example with pure advection indicates the limitations of these estimates and of the model developed here. We also consider and analyze the coupled pressure equation that determines the advective flux in the transport model.
Distributed Parallel Particle Advection using Work Requesting
Muller, Cornelius; Camp, David; Hentschel, Bernd; Garth, Christoph
2013-09-30
Particle advection is an important vector field visualization technique that is difficult to apply to very large data sets in a distributed setting due to scalability limitations in existing algorithms. In this paper, we report on several experiments using work requesting dynamic scheduling which achieves balanced work distribution on arbitrary problems with minimal communication overhead. We present a corresponding prototype implementation, provide and analyze benchmark results, and compare our results to an existing algorithm.
Efficient mass transport by optical advection
Kajorndejnukul, Veerachart; Sukhov, Sergey; Dogariu, Aristide
2015-01-01
Advection is critical for efficient mass transport. For instance, bare diffusion cannot explain the spatial and temporal scales of some of the cellular processes. The regulation of intracellular functions is strongly influenced by the transport of mass at low Reynolds numbers where viscous drag dominates inertia. Mimicking the efficacy and specificity of the cellular machinery has been a long time pursuit and, due to inherent flexibility, optical manipulation is of particular interest. However, optical forces are relatively small and cannot significantly modify diffusion properties. Here we show that the effectiveness of microparticle transport can be dramatically enhanced by recycling the optical energy through an effective optical advection process. We demonstrate theoretically and experimentally that this new advection mechanism permits an efficient control of collective and directional mass transport in colloidal systems. The cooperative long-range interaction between large numbers of particles can be optically manipulated to create complex flow patterns, enabling efficient and tunable transport in microfluidic lab-on-chip platforms. PMID:26440069
An implicit dispersive transport algorithm for the US Geological Survey MOC3D solute-transport model
Kipp, K.L., Jr.; Konikow, L.F.; Hornberger, G.Z.
1998-01-01
This report documents an extension to the U.S. Geological Survey MOC3D transport model that incorporates an implicit-in-time difference approximation for the dispersive transport equation, including source/sink terms. The original MOC3D transport model (Version 1) uses the method of characteristics to solve the transport equation on the basis of the velocity field. The original MOC3D solution algorithm incorporates particle tracking to represent advective processes and an explicit finite-difference formulation to calculate dispersive fluxes. The new implicit procedure eliminates several stability criteria required for the previous explicit formulation. This allows much larger transport time increments to be used in dispersion-dominated problems. The decoupling of advective and dispersive transport in MOC3D, however, is unchanged. With the implicit extension, the MOC3D model is upgraded to Version 2. A description of the numerical method of the implicit dispersion calculation, the data-input requirements and output options, and the results of simulator testing and evaluation are presented. Version 2 of MOC3D was evaluated for the same set of problems used for verification of Version 1. These test results indicate that the implicit calculation of Version 2 matches the accuracy of Version 1, yet is more efficient than the explicit calculation for transport problems that are characterized by a grid Peclet number less than about 1.0.
Kubo formulas for dispersion in heterogeneous periodic nonequilibrium systems
NASA Astrophysics Data System (ADS)
Guérin, T.; Dean, D. S.
2015-12-01
We consider the dispersion properties of tracer particles moving in nonequilibrium heterogeneous periodic media. The tracer motion is described by a Fokker-Planck equation with arbitrary spatially periodic (but constant in time) local diffusion tensors and drifts, eventually with the presence of obstacles. We derive a Kubo-like formula for the time-dependent effective diffusion tensor valid in any dimension. From this general formula, we derive expressions for the late time effective diffusion tensor and drift in these systems. In addition, we find an explicit formula for the late finite-time corrections to these transport coefficients. In one dimension, we give a closed analytical formula for the transport coefficients. The formulas derived here are very general and provide a straightforward method to compute the dispersion properties in arbitrary nonequilibrium periodic advection-diffusion systems.
Kubo formulas for dispersion in heterogeneous periodic nonequilibrium systems.
Guérin, T; Dean, D S
2015-12-01
We consider the dispersion properties of tracer particles moving in nonequilibrium heterogeneous periodic media. The tracer motion is described by a Fokker-Planck equation with arbitrary spatially periodic (but constant in time) local diffusion tensors and drifts, eventually with the presence of obstacles. We derive a Kubo-like formula for the time-dependent effective diffusion tensor valid in any dimension. From this general formula, we derive expressions for the late time effective diffusion tensor and drift in these systems. In addition, we find an explicit formula for the late finite-time corrections to these transport coefficients. In one dimension, we give a closed analytical formula for the transport coefficients. The formulas derived here are very general and provide a straightforward method to compute the dispersion properties in arbitrary nonequilibrium periodic advection-diffusion systems. PMID:26764628
Advection and diffusion in shoreline change prediction
NASA Astrophysics Data System (ADS)
Anderson, T. R.; Frazer, L. N.
2010-12-01
We added longshore advection and diffusion to the simple cross-shore rate calculation method, as used widely by the USGS and others, to model historic shorelines and to predict future shoreline positions; and applied this to Hawaiian Island beach data. Aerial photographs, sporadically taken throughout the past century, yield usable, albeit limited, historic shoreline data. These photographs provide excellent spatial coverage, but poor temporal resolution, of the shoreline. Due to the sparse historic shoreline data, and the many natural and anthropogenic events influencing coastlines, we constructed a simplistic shoreline change model that can identify long-term behavior of a beach. Our new, two-dimensional model combines the simple rate method to accommodate for cross-shore sediment transport with the classic Pelnard-Considère model for diffusion, as well as a longshore advection speed term. Inverse methods identify cross-shore rate, longshore advection speed, and longshore diffusivity down a sandy coastline. A spatial averaging technique then identifies shoreline segments where one parameter can reasonably account for the cross-shore and longshore transport rates in that area. This produces model results with spatial resolution more appropriate to the temporal spacing of the data. Because changes in historic data can be accounted for by varying degrees of cross-shore and longshore sediment transport - for example, beach erosion can equally be explained by sand moving either off-shore or laterally - we tested several different model scenarios on the data: allowing only cross-shore sediment movement, only longshore movement, and a combination of the two. We used statistical information criteria to determine both the optimal spatial resolution and best-fitting scenario. Finally, we employed a voting method predicting the relaxed shoreline position over time.
Coupled ensemble flow line advection and analysis.
Guo, Hanqi; Yuan, Xiaoru; Huang, Jian; Zhu, Xiaomin
2013-12-01
Ensemble run simulations are becoming increasingly widespread. In this work, we couple particle advection with pathline analysis to visualize and reveal the differences among the flow fields of ensemble runs. Our method first constructs a variation field using a Lagrangian-based distance metric. The variation field characterizes the variation between vector fields of the ensemble runs, by extracting and visualizing the variation of pathlines within ensemble. Parallelism in a MapReduce style is leveraged to handle data processing and computing at scale. Using our prototype system, we demonstrate how scientists can effectively explore and investigate differences within ensemble simulations. PMID:24051840
Electrokinetic induced solute dispersion in porous media; pore network modeling
NASA Astrophysics Data System (ADS)
Li, Shuai; Schotting, Ruud; Raoof, Amir
2013-04-01
Electrokinetic flow plays an important role in remediation process, separation technique, and chromatography. The solute dispersion is a key parameter to determine transport efficiency. In this study, we present the electrokinetic effects on solute dispersion in porous media at the pore scale, using a pore network model. The analytical solution of the electrokinetic coupling coefficient was obtained to quantity the fluid flow velocity in a cylinder capillary. The effect of electrical double layer on the electrokinetic coupling coefficient was investigated by applying different ionic concentration. By averaging the velocity over cross section within a single pore, the average flux was obtained. Applying such single pore relationships, in the thin electrical double layer limit, to each and every pore within the pore network, potential distribution and the induced fluid flow was calculated for the whole domain. The resulting pore velocities were used to simulate solute transport within the pore network. By averaging the results, we obtained the breakthrough curve (BTC) of the average concentration at the outlet of the pore network. Optimizing the solution of continuum scale advection-dispersion equation to such a BTC, solute dispersion coefficient was estimated. We have compared the dispersion caused by electrokinetic flow and pure pressure driven flow under different Peclet number values. In addition, the effect of microstructure and topological properties of porous media on fluid flow and solute dispersion is presented, mainly based on different pore coordination numbers.
Advective turbulent transport in the fluid plasma
NASA Astrophysics Data System (ADS)
Min, Byung-Hoon; An, Chan-Yong; Kim, Chang-Bae
2013-10-01
The Hasegawa-Wakatani model (HWM) has been employed in pedagogical analyses of the physics behind the behavior of the tokamak plasmas. In addition to the geometric simplicity HWM has an appealing feature of sustaining autonomous quasi-steady state, unstable modes providing the power that is being transported by the nonlinear interactions and is eventually dissipated by the collisional damping at small scales. Emergence of the zonal flow out of the turbulence is a main candidate to cause the transition from the low plasma confinement to the high mode. In the study of such LH transition with the HWM, the adiabaticity parameter has been shown to play an important role in forcing the zonal flow that results in the regulation of the drift-wave turbulence. Instead of concentrating on the physics of the feedback loop between the turbulence and the zonal flow the present study focuses on the presence of the advective transport of the energy. Numerical simulations of HWM are performed and the connections between the advective transport and the zonal flow will be presented. This work was supported by the Supercpmputing Center/Korea Institute of Science and Technology Information with supercomputing resources including technical support (KSC-2013-C1-009).
Advection, diffusion and delivery over a network
Heaton, Luke L.M.; López, Eduardo; Maini, Philip K.; Fricker, Mark D.; Jones, Nick S.
2014-01-01
Many biological, geophysical and technological systems involve the transport of resource over a network. In this paper we present an algorithm for calculating the exact concentration of resource at any point in space or time, given that the resource in the network is lost or delivered out of the network at a given rate, while being subject to advection and diffusion. We consider the implications of advection, diffusion and delivery for simple models of glucose delivery through a vascular network, and conclude that in certain circumstances, increasing the volume of blood and the number of glucose transporters can actually decrease the total rate of glucose delivery. We also consider the case of empirically determined fungal networks, and analyze the distribution of resource that emerges as such networks grow over time. Fungal growth involves the expansion of fluid filled vessels, which necessarily involves the movement of fluid. In three empirically determined fungal networks we found that the minimum currents consistent with the observed growth would effectively transport resource throughout the network over the time-scale of growth. This suggests that in foraging fungi, the active transport mechanisms observed in the growing tips may not be required for long range transport. PMID:23005783
Normal and Anomalous Dispersion in Fluvial Sediment Transport
NASA Astrophysics Data System (ADS)
Bradley, D. N.; Tucker, G. E.
2005-12-01
Understanding the rate of motion and pattern of dispersion in fluvial sediment transport is essential for a variety of applications, including predicting the fate and transport of solid-phase contaminants and modeling the cosmogenic-nuclide inheritance of water-borne sediment. In order to create a probabilistic model of sediment particle motion, it is necessary to characterize the statistical properties of fluvial sediment dispersion. In general, two modes of behavior have been observed in advective-diffusive transport systems: normal and anomalous dispersion. Normal dispersion is characterized by a well-defined mean position and spatial variance and the time evolution of particle concentration is described by a simple advection-diffusion equation. In contrast, a transport system that exhibits anomalous dispersion will tend to have a heavy-tailed spatial distribution, a mean position that is different from the peak concentration, and a large variance. The fundamental difference lies in the probability distribution of individual particle velocities. When the distribution is sufficiently heavy-tailed, the resulting dispersion pattern will be anomalous. Anomalous dispersion has been observed in geophysical systems ranging from turbulent flow to transport in heterogeneous porous media. Several lines of evidence from the sediment transport literature suggest that fluvial sediment may undergo anomalous dispersion. Tracer experiments show a preference for right-skewed travel distance distributions, a characteristic of anomalous diffusion. Studies suggest that large inputs of sediment to rivers (such as a landslide) tend to disperse in place rather than translate downstream. In addition, the fact that sediment grains can become trapped in flood plains and bars for long periods of time and then move long distances in rare, short duration events such as floods suggests a potential for anomalous dispersion due to a broad distribution of particle residence times. We develop a
Kukačka, Libor; Nosek, Štĕpán; Kellnerová, Radka; Jurčáková, Klára; Jaňour, Zbyněk
2012-01-01
The objective of this study is to determine processes of pollution ventilation in the X-shaped street intersection in an idealized symmetric urban area for the changing approach flow direction. A unique experimental setup for simultaneous wind tunnel measurement of the flow velocity and the tracer gas concentration in a high temporal resolution is assembled. Advective horizontal and vertical scalar fluxes are computed from averaged measured velocity and concentration data within the street intersection. Vertical advective and turbulent scalar fluxes are computed from synchronized velocity and concentration signals measured in the plane above the intersection. All the results are obtained for five approach flow directions. The influence of the approach flow on the advective and turbulent fluxes is determined. The contribution of the advective and turbulent flux to the ventilation is discussed. Wind direction with the best dispersive conditions in the area is found. The quadrant analysis is applied to the synchronized signals of velocity and concentration fluctuation to determine events with the dominant contribution to the momentum flux and turbulent scalar flux. PMID:22649290
A balancing domain decomposition method by constraints for advection-diffusion problems
Tu, Xuemin; Li, Jing
2008-12-10
The balancing domain decomposition methods by constraints are extended to solving nonsymmetric, positive definite linear systems resulting from the finite element discretization of advection-diffusion equations. A pre-conditioned GMRES iteration is used to solve a Schur complement system of equations for the subdomain interface variables. In the preconditioning step of each iteration, a partially sub-assembled finite element problem is solved. A convergence rate estimate for the GMRES iteration is established, under the condition that the diameters of subdomains are small enough. It is independent of the number of subdomains and grows only slowly with the subdomain problem size. Numerical experiments for several two-dimensional advection-diffusion problems illustrate the fast convergence of the proposed algorithm.