scholarly journals Analytical solutions for linearized Richards Equation with arbitrary time-dependent surface fluxes

2001 ◽  
Vol 37 (4) ◽  
pp. 1091-1093 ◽  
Author(s):  
Jiann-Mou Chen ◽  
Yih-Chi Tan ◽  
Chu-Hui Chen ◽  
J.-Y. Parlange
2021 ◽  
Author(s):  
Timo Koch ◽  
Hanchuan Wu ◽  
Kent-André Mardal ◽  
Rainer Helmig ◽  
Martin Schneider

<p>1D-3D methods are used to describe root water and nutrient uptake in complex root networks. Root systems are described as networks of line segments embedded in a three-dimensional soil domain. Particularly for dry soils, local water pressure and nutrient concentration gradients can be become very large in the vicinity of roots. Commonly used discretization lengths (for example 1cm) in root-soil interaction models do not allow to capture these gradients accurately. We present a new numerical scheme for approximating root-soil interface fluxes. The scheme is formulated in the continuous PDE setting so that is it formally independent of the spatial discretization scheme (e.g. FVM, FD, FEM). The interface flux approximation is based on a reconstruction of interface quantities using local analytical solutions of the steady-rate Richards equation. The local mass exchange is numerically distributed in the vicinity of the root. The distribution results in a regularization of the soil pressure solution which is easier to approximate numerically. This technique allows for coarser grid resolutions while maintaining approximation accuracy. The new scheme is verified numerically against analytical solutions for simplified cases. We also explore limitations and possible errors in the flux approximation with numerical test cases. Finally, we present the results of a recently published benchmark case using this new method.</p>


2019 ◽  
Vol 2019 (1) ◽  
pp. 21-47 ◽  
Author(s):  
Shuzhan Ren

Abstract A solution to the 3D transport equation for passive tracers in the atmospheric boundary layer (ABL), formulated in terms of Green’s function (GF), is derived to show the connection between the concentration and surface fluxes of passive tracers through GF. Analytical solutions to the 1D vertical diffusion equation are derived to reveal the nonlinear dependence of the concentration and flux on the diffusivity, time, and height, and are employed to examine the impact of the diffusivity on the diurnal variations of CO2 in the ABL. The properties of transport operator H and their implications in inverse modeling are discussed. It is found that H has a significant contribution to the rectifier effect in the diurnal variations of CO2. Since H is the integral of GF in time, the narrow distribution of GF in time justifies the reduction of the size of H in inverse modeling. The exponential decay of GF with height suggests that the estimated surface fluxes in inverse modeling are more sensitive to the observations in the lower ABL. The solutions and first mean value theorem are employed to discuss the uncertainties associated with the concentration–mean surface flux equation used to link the concentrations and mean surface flux. Both analytical and numerical results show that the equation can introduce big errors, particularly when surface flux is sign indefinite. Numerical results show that the conclusions about the evolution properties of passive tracers based on the analytical solutions also hold in the cases with a more complicated diffusion coefficient and time-varying ABL height.


Author(s):  
M. Fakoor-Pakdaman ◽  
M. Andisheh-Tadbir ◽  
Majid Bahrami

A new all-time model is developed to predict transient laminar forced convection heat transfer inside a circular tube under arbitrary time-dependent heat flux. Slug flow condition is assumed for the velocity profile inside the tube. The solution to the time-dependent energy equation for a step heat flux boundary condition is generalized for arbitrary time variations in surface heat flux using a Duhamel’s integral technique. A cyclic time-dependent heat flux is considered and new compact closed-form relationships are proposed to predict: i) fluid temperature distribution inside the tube ii) fluid bulk temperature and iii) the Nusselt number. A new definition, cyclic fully-developed Nusselt number, is introduced and it is shown that in the thermally fully-developed region the Nusselt number is not a function of axial location, but it varies with time and the angular frequency of the imposed heat flux. Optimum conditions are found which maximize the heat transfer rate of the unsteady laminar forced-convective tube flow. We also performed an independent numerical simulation using ANSYS to validate the present analytical model. The comparison between the numerical and the present analytical model shows great agreement; a maximum relative difference less than 5.3%.


1982 ◽  
Vol 117 ◽  
pp. 211-231 ◽  
Author(s):  
Mostafa A. Foda

An analytical theory is developed to describe how negative pressure, (or ‘mud suction’, as it is sometimes referred to) develops underneath a body as it detaches itself from the ocean bottom. Biot's quasistatic equations of poro-elasticity are used to model the ocean bottom, and a general three-dimensional time-dependent analysis of the problem is worked out first using the boundary-layer approximation recently proposed by Mei and Foda. Then, explicit leading-order analytical solutions are presented for the problems of extrication of slender bodies as well as axisymmetric bodies from the ocean bottom.


Sign in / Sign up

Export Citation Format

Share Document