Ion uptake from soils by plant roots, subject to the Epstein-Hagen relation

Soil Research ◽  
1968 ◽  
Vol 6 (2) ◽  
pp. 159 ◽  
Author(s):  
J Kautsky ◽  
KP Barley ◽  
DK Fiddaman
Keyword(s):  
Differential Equation ◽  
Linear Approximation ◽  
Numerical Solutions ◽  
Diffusion Theory ◽  
Nonlinear Function ◽  
Solution Concentration ◽  
Root Surface ◽  
Ion Uptake ◽  
Plant Roots

Diffusion theory has been used to describe the uptake of ions from soils by single roots. Numerical solutions of the relevant differential equation are given for the case where rate of uptake is a nonlinear function of the solution concentration at the root surface. The function employed is derived from a physiological relation described by Epstein and Hagen (1962). Results are compared with those obtained using a common linear approximation.

A Meshless Method for Solving the Mathematical Model Associated the Leakage Problem in Gas Pipeline

2014 ◽  
Vol 986-987 ◽  
pp. 1418-1421
Author(s):  
Jun Shan Li
Keyword(s):  
Differential Equation ◽  
Mathematical Model ◽  
Partial Differential Equation ◽  
Inverse Problem ◽  
Basis Function ◽  
Meshless Method ◽  
Numerical Solutions ◽  
Basis Functions ◽  
Gas Pipeline ◽  
The Mathematical Model

In this paper, we propose a meshless method for solving the mathematical model concerning the leakage problem when the pressure is tested in the gas pipeline. The method of radial basis function (RBF) can be used for solving partial differential equation by writing the solution in the form of linear combination of radius basis functions, that is, when integrating the definite conditions, one can find the combination coefficients and then the numerical solution. The leak problem is a kind of inverse problem that is focused by many engineers or mathematical researchers. The strength of the leak can find easily by the additional conditions and the numerical solutions.

Integral Solutions of Phase Change With Internal Heat Generation

2004 ◽  
Author(s):  
Ali Siahpush ◽  
John Crepeau
Keyword(s):  
Differential Equation ◽  
Phase Change ◽  
Heat Generation ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Melting Rate ◽  
Internal Heat Generation ◽  
Stefan Number ◽  
Temperature Boundary ◽  
Solid Liquid

This paper presents solutions to a one-dimensional solid-liquid phase change problem using the integral method for a semi-infinite material that generates internal heat. The analysis assumed a quadratic temperature profile and a constant temperature boundary condition on the exposed surface. We derived a differential equation for the solidification thickness as a function of the internal heat generation (IHG) and the Stefan number, which includes the temperature of the boundary. Plots of the numerical solutions for various values of the IHG and Stefan number show the time-dependant behavior of both the melting and solidification distances and rates. The IHG of the material opposes solidification and enhances melting. The differential equation shows that in steady-state, the thickness of the solidification band is inversely related to the square root of the IHG. The model also shows that the melting rate initially decreases and reaches a local minimum, then increases to an asymptotic value.

The Partial Transformational Decomposition Method for a Hybrid Analytical/Numerical Solution of the 3D Gas-Flow Problem in a Hydraulically Fractured Ultralow-Permeability Reservoir

SPE Journal ◽  
2021 ◽  
pp. 1-28
Author(s):  
George Moridis ◽  
Niwit Anantraksakul ◽  
Thomas A. Blasingame
Keyword(s):  
Differential Equation ◽  
Decomposition Method ◽  
Gas Flow ◽  
Numerical Solutions ◽  
Flow Equation ◽  
Time Variable ◽  
Strong Nonlinearity ◽  
Fully Implicit ◽  
Ultralow Permeability ◽  
Klinkenberg Effects

Summary The analysis of gas production from fractured ultralow-permeability (ULP) reservoirs is most often accomplished using numerical simulation, which requires large 3D grids, many inputs, and typically long execution times. We propose a new hybrid analytical/numerical method that reduces the 3D equation of gas flow into either a simple ordinary-differential equation (ODE) in time or a 1D partial-differential equation (PDE) in space and time without compromising the strong nonlinearity of the gas-flow relation, thus vastly decreasing the size of the simulation problem and the execution time. We first expand the concept of pseudopressure of Al-Hussainy et al. (1966) to account for the pressure dependence of permeability and Klinkenberg effects, and we also expand the corresponding gas-flow equation to account for Langmuir sorption. In the proposed hybrid partial transformational decomposition method (TDM) (PTDM), successive finite cosine transforms (FCTs) are applied to the expanded, pseudopressure-based 3D diffusivity equation of gas flow, leading to the elimination of the corresponding physical dimensions. For production under a constant- or time-variable rate (q) regime, three levels of FCTs yield a first-order ODE in time. For production under a constant- or time-variable pressure (pwf) regime, two levels of FCTs lead to a 1D second-order PDE in space and time. The fully implicit numerical solutions for the FCT-based equations in the multitransformed spaces are inverted, providing solutions that are analytical in 2D or 3D and account for the nonlinearity of gas flow. The PTDM solution was coded in a FORTRAN95 program that used the Laplace-transform (LT) analytical solution for the q-problem and a finite-difference method for the pwf problem in their respective multitransformed spaces. Using a 3D stencil (the minimum repeatable element in the horizontal well and hydraulically fractured system), solutions over an extended production time and a substantial pressure drop were obtained for a range of isotropic and anisotropic matrix and fracture properties, constant and time-variableQ and pwf production schemes, combinations of stimulated-reservoir-volume (SRV) and non-SRV subdomains, sorbing and nonsorbing gases of different compositions and at different temperatures, Klinkenberg effects, and the dependence of matrix permeability on porosity. The limits of applicability of PTDM were also explored. The results were compared with the numerical solutions from a widely used, fully implicit 3D simulator that involved a finely discretized (high-definition) 3D domain involving 220,000 elements and show that the PTDM solutions can provide accurate results for long times for large well drawdowns even under challenging conditions. Of the two versions of PTDM, the PTD-1D was by far the better option and its solutions were shown to be in very good agreement with the full numerical solutions, while requiring a fraction of the memory and orders-of-magnitude lower execution times because these solutions require discretization of only the time domain and a single axis (instead of three). The PTD-0D method was slower than PTD-1D (but still much faster than the numerical solution), and although its solutions were accurate for t < 6 months, these solutions deteriorated beyond that point. The PTDM is an entirely new approach to the analysis of gas flow in hydraulically fractured ULP reservoirs. The PTDM solutions preserve the strong nonlinearity of the gas-flow equation and are analytical in 2D or 3D. This being a semianalytical approach, it needs very limited input data and requires computer storage and computational times that are orders-of-magnitude smaller than those in conventional (numerical) simulators because its discretization is limited to time and (possibly) a single spatial dimension.

Impact of a Mass on a Damped Elastically Supported Beam

1948 ◽  
Vol 15 (2) ◽  
pp. 125-136
Author(s):  
W. H. Hoppmann
Keyword(s):  
Differential Equation ◽  
Elastic Foundation ◽  
Forced Vibration ◽  
Numerical Solutions ◽  
Coefficient Of Restitution ◽  
Tabular Form ◽  
Single Degree Of Freedom ◽  
Simply Supported ◽  
Single Degree ◽  
Elastically Supported

Abstract In this paper a study is made of the problem of the central impact of a mass on a simply supported beam on an elastic foundation with considerations of internal and external damping. The differential equation for the forced vibration of the beam is developed. It is solved for the case in which the force is a function of time and is concentrated at the center of the beam. Formulas are obtained for the deflections. An expression is developed for the coefficient of restitution which is essential in determining the deflections and the strains. Criteria are devised for determining the cases in which the beam may be considered as a single-degree-of-freedom system when damping and an elastic foundation are considered. The importance of these criteria is discussed. A numerical example illustrating the theory developed in the paper is worked out in detail. Results of computations for several numerical solutions are given in tabular form.

Numerical analysis of time-dependent stagnation point flow of Oldroyd-B fluid subject to modified Fourier’s law

2021 ◽  
pp. 2150187
Author(s):  
Foukeea Qasim ◽  
Tian-Chuan Sun ◽  
S. Z. Abbas ◽  
W. A. Khan ◽  
M. Y. Malik
Keyword(s):  
Differential Equation ◽  
Fluid Flow ◽  
Stagnation Point ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Nonlinear Partial Differential Equation ◽  
Thermal Relaxation ◽  
Stagnation Point Flow ◽  
Time Dependent ◽  
Parameter Values

This paper aims to investigate the time-dependent stagnation point flow of an Oldroyd-B fluid subjected to the modified Fourier law. The flow into a vertically stretched cylinder at the stagnation point is discussed. The heat flux model of a non-Fourier is intended for the transfer of thermal energy in fluid flow. The study is carried out on the surface heating source, namely the surface temperature. The developed nonlinear partial differential equation for regulating fluid flow and heat transport is transformed via appropriate similarity variables into a nonlinear ordinary differential equation. The development and analysis of convergent series solutions were considered for velocity and temperature. Prandtl number numerical values are computed and investigated. This study’s findings are compared to the previous findings. By making use of the bvp4c Matlab method, numerical solutions are obtained. Besides, high buoyancy parameter values are found to increase the fluid velocity for the stimulating approach. By improving the thermal relaxation time parameter values, heat transfer in the fluid flow decreases. The temperature field effects are displayed graphically.

Solute Transport in the Soil near Root Surfaces

2000 ◽  
Author(s):  
Peter B. Tinker ◽  
Peter Nye
Keyword(s):  
Soil Solution ◽  
Dynamic Equilibrium ◽  
Pore Solution ◽  
Plant Root ◽  
Unit Length ◽  
Solution Concentration ◽  
Root Surface ◽  
Transport Processes ◽  
Soil Pores ◽  
Intact Root

We discussed in chapter 4 the movement of solute between small volumes of soil, and in chapter 5 some properties of plant roots and associated hairs, particularly the relation between the rate of uptake at the root surface and the concentration of solute in the ambient solution. In the chapters to follow, we consider the plant root in contact with the soil, and deal with their association in increasingly complex situations; first, when the root acts merely as a sink and, second, when it modifies its relations with the surrounding soil by changing its pH, excreting ions, stimulating microorganisms, or developing mycorrhizas. In this chapter, we take the simplest situation that can be studied in detail, namely, a single intact root alone in a volume of soil so large that it can be considered infinite. The essential transport processes occurring near the root surface are illustrated in figure 6.1. We have examined in chapter 3 the rapid dynamic equilibrium between solutes in the soil pore solution and those sorbed on the immediately adjacent solid surfaces. These sorbed solutes tend to buffer the soil solution against changes in concentration induced by root uptake. At the root surface, solutes are absorbed at a rate related to their concentration in the soil solution at the boundary (section 5.3.2); and the root demand coefficient, αa, is defined by the equation . . . I = 2παaCLa (6.1) . . . where I = inflow (rate of uptake per unit length), a = root radius, CLa = concentration in solution at the root surface. To calculate the inflow, we have to know CLa, and the main topic of this chapter is the relation between CLa, and the soil pore solution concentration CL. The root also absorbs water at its surface due to transpiration (chapter 2) so that the soil solution flows through the soil pores, thus carrying solutes to the root surface by mass flow (convection). Barber et al. (1962) calculated whether the nutrients in maize could be acquired solely by this process, by multiplying the composition of the soil solution by the amount of water the maize had transpired.

Fourier-Bessel Series Model for the Stefan Problem With Internal Heat Generation in Cylindrical Coordinates

2018 ◽  
Author(s):  
Lyudmyla Barannyk ◽  
John Crepeau ◽  
Patrick Paulus ◽  
Ali Siahpush
Keyword(s):  
Differential Equation ◽  
Steady State ◽  
Heat Generation ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Internal Heat Generation ◽  
Time Dependent ◽  
Cylindrical Coordinates ◽  
System Approaches ◽  
Melting And Solidification

A nonlinear, first-order ordinary differential equation that involves Fourier-Bessel series terms has been derived to model the time-dependent motion of the solid-liquid interface during melting and solidification of a material with constant internal heat generation in cylindrical coordinates. The model is valid for all Stefan numbers. One of the primary applications of this problem is for a nuclear fuel rod during meltdown. The numerical solutions to this differential equation are compared to the solutions of a previously derived model that was based on the quasi-steady approximation, which is valid only for Stefan numbers less than one. The model presented in this paper contains exponentially decaying terms in the form of Fourier-Bessel series for the temperature gradients in both the solid and liquid phases. The agreement between the two models is excellent in the low Stefan number regime. For higher Stefan numbers, where the quasi-steady model is not accurate, the new model differs from the approximate model since it incorporates the time-dependent terms for small times, and as the system approaches steady-state, the curves converge. At higher Stefan numbers, the system approaches steady-state faster than for lower Stefan numbers. During the transient process for both melting and solidification, the temperature profiles become parabolic.

Rough Contacts Modeled by Nonlinear Coatings

2007 ◽  
Author(s):  
I. I. Kudish
Keyword(s):  
Contact Problem ◽  
Numerical Solutions ◽  
Asymptotic Methods ◽  
Experimental Studies ◽  
Nonlinear Function ◽  
Normal Displacement ◽  
Experimental Fact ◽  
Axially Symmetric ◽  
Elastic Solids ◽  
Local Pressure

A number of experimental studies [1–3] revealed that the normal displacement in a contact of rough surfaces due to asperities presence is a nonlinear function of local pressure and it can be approximated by a power function of pressure. Originally, a linear mathematical model accounting for surface roughness of elastic solids in contact was introduced by I. Shtaerman [4]. He assumed that the effect of asperities present in a contact of elastic solids can be essentially replaced by the presence of a thin coating simulated by an additional normal displacement of solids’ surfaces proportional to a local pressure. Later, a similar but nonlinear problem formulation that accounted for the above mentioned experimental fact was proposed by L. Galin. In a series of papers this problem was studied by numerical and asymptotic methods [5–9]. The present paper has a dual purpose: to analyze the problem analytically and to provide some asymptotic and numerical solutions. The results presented below provide an overview of the results obtained on the topic and published by the author earlier in the journals hardly accessible to the international tribological community (such as Russian and mathematical journals) and, therefore, mostly unknown by tribologists. A number of recent publications on contacts of rough elastic solids supports the view that these results are still of value to the specialists involved in nanotribology. The existence and uniqueness of a solution of a contact problem for elastic bodies with rough (coated) surfaces is established based on the variational inequalities approach. Four different equivalent formulations of the problem including three variational ones were considered. A comparative analysis of solutions of the contact problem for different values of initial parameters (such as the indenter shape, parameters characterizing roughness, elastic parameters of the substrate material) is done with the help of calculus of variations and the Zaremba-Giraud principle of maximum for harmonic functions [10,11]. The results include the relations between the pressure and displacement distributions for rough and smooth solids as well as the relationships for solutions of the problems for rough solids with fixed and free contact boundaries. For plane and axially symmetric cases some asymptotic and numerical solutions are presented.

scholarly journals CRITICAL LENGTH FOR THE SPREADING–VANISHING DICHOTOMY IN HIGHER DIMENSIONS

2020 ◽  
Vol 62 (1) ◽  
pp. 3-17 ◽  
Author(s):  
MATTHEW J. SIMPSON
Keyword(s):  
Differential Equation ◽  
Partial Differential Equation ◽  
Moving Boundary ◽  
Numerical Solutions ◽  
Critical Length ◽  
Higher Dimensions ◽  
Moving Boundary Problem ◽  
Kolmogorov Equation ◽  
One Dimension ◽  
Spatial Dimensions

We consider an extension of the classical Fisher–Kolmogorov equation, called the “Fisher–Stefan” model, which is a moving boundary problem on $0<x<L(t)$. A key property of the Fisher–Stefan model is the “spreading–vanishing dichotomy”, where solutions with $L(t)>L_{\text{c}}$ will eventually spread as $t\rightarrow \infty$, whereas solutions where $L(t)\ngtr L_{\text{c}}$ will vanish as $t\rightarrow \infty$. In one dimension it is well known that the critical length is $L_{\text{c}}=\unicode[STIX]{x1D70B}/2$. In this work, we re-formulate the Fisher–Stefan model in higher dimensions and calculate $L_{\text{c}}$ as a function of spatial dimensions in a radially symmetric coordinate system. Our results show how $L_{\text{c}}$ depends upon the dimension of the problem, and numerical solutions of the governing partial differential equation are consistent with our calculations.

Invariant solutions of hyperbolic fuzzy fractional differential equations

2019 ◽  
Vol 34 (01) ◽  
pp. 2050015 ◽  
Author(s):  
C. Vinothkumar ◽  
J. J. Nieto ◽  
A. Deiveegan ◽  
P. Prakash
Keyword(s):  
Differential Equation ◽  
Fractional Differential Equation ◽  
Numerical Solutions ◽  
Banach Fixed Point Theorem ◽  
Invariant Solutions ◽  
Existence And Uniqueness Results ◽  
Uniqueness Results ◽  
Fractional Differential ◽  
Fuzzy Fractional Differential Equations ◽  
Fuzzy Fractional Differential Equation

We consider the hyperbolic type fuzzy fractional differential equation and derive the second-order fuzzy fractional differential equation using scaling transformation. We present a theoretical and a numerical method to find the invariant solutions of such equations. Also, we prove the existence and uniqueness results using Banach fixed point theorem. Numerical solutions are approximated using finite difference method. Finally, numerical examples are given to illustrate the obtained results.

Sign in / Sign up

Export Citation Format

Share Document