Water table fluctuation with a random initial condition

The nonlinear Boussinesq equation is used to understand water table fluctuations in various ditch drainage problems. An approximate solution of this equation with a random initial condition and deterministic boundary conditions, recharge rate and aquifer parameters has been developed to predict a transient water table in a ditch-drainage system. The effects of uncertainty in the initial condition on the water table are illustrated with the help of a synthetic example. These results would find applications in ditch-drainage design.Notation A / tanh t - a lower value of the random variable representing the initial water table height at the mid point - a+b Upper value of the random variable representing the initial water table height at the midpoint - B tanh t - C 4/ - h variable water table height - h mean of the variable water table height - h _m variable water table height at the mid point - h _m mean of the variable water table height at the mid point - K hydraulic conductivity - L half spacing between the ditches - m ₀ initial water table height at the mid point - N Uniform rate of recharge - S specific yield - t time of observation - x distance measured from the ditch boundary - (4/SL)(NK)^1/2 - (L/4)(N/K)^1/2 - dummy integral variable     

Dynamics of fluid circulation in coupled free and heterogeneous porous domains

Fluid flow in coupled free and porous domain, particularly when the porous medium is heterogeneous, is encountered in many hydro-environmental conditions, e.g., leakage from underground pipe, combined groundwater lake-subsurface interactions. One of the most difficult problems in the study of coupled flow behaviour has been the development of a universally applicable modelling scheme for combining the flow regimes. This is because the free/porous interfacial properties (e.g., shear-stress; velocity slip) that govern the coupled flow behaviour are difficult to determine experimentally under hydro-environmental conditions. On the other hand, the implications of various forms of heterogeneity in the porous media properties can be very different on the fluid-flow behaviour. Difficulties may also arise in direct coupling of the model equations that govern the fluid flow in the individual regions (e.g., Navier-Stokes for free-flow region and the Darcy's equation for the porous flow region). Consequently, models of coupled free and porous flow for hydro-environmental conditions are not very well developed at the moment. While there are some indications that fluids in coupled free and porous domains may circulate (i.e., development of flow cells), there is a lack of appropriate 3D analysis on how heterogeneities in porous media may affect such flow patterns. In this paper, we aim to analyse how porous media heterogeneity affects the dynamics of flow circulation in the porous side of a coupled free and porous domain. For this purpose, we analyse flow patterns in several model domains made up of two porous layers with differing permeabilities. The governing model equations are discretised and solved using the standard finite volume method on a staggered cell-centred mesh. The temporal discretisation is done using the explicit method. An in-house graphical user interface (GUI) has been created specifically to aid in the visualisation of otherwise complex flow patterns. The GUI contains many post-processing options and provides a comprehensive tool for the analysis of hydrodynamics and contaminant motion (not discussed in this paper) in coupled free and porous flow domains. This GUI is described in this paper briefly. The effects of altering the aspect ratio (i.e., multi-scale) of the domain on the coupled flow pattern have also been discussed.     

Methodology to resolve the transport equation with the discrete ordinates code TORT into the IPEN/MB-01 reactor

Resolution of the steady-state Neutron Transport Equation in a nuclear pool reactor is usually achieved by means of two different numerical methods: Monte Carlo (stochastic) and Discrete Ordinates (deterministic). The Discrete Ordinates method solves the Neutron Transport Equation for a set of selected directions, obtaining a set of directional equations and solutions for each equation which are the angular flux. In order to deal with the energy dependence, an energy multi-group approximation is commonly performed, obtaining a set of equations depending on the number of energy groups. In addition, spatial discretization is also required and the problem is solved by sweeping the geometry mesh. However, special cross-sections are required due to the energy and directional discretization, thus a methodology based on NJOY99 code capabilities has been used. Finally, in order to demonstrate the capability of this method, the 3D discrete ordinates code TORT has been applied to resolve the IPEN/MB-01 reactor.     

Bubble temperature data for phenol-alcohol systems: Phenol-I-Amyl alcohol and phenol-n-Amyl alcohol systems

In order to supplement the literature data on VLE, isobaric bubble temperature data were collected for the binary systems phenol-i-amyl alcohol and phenol-n-amyl alcohol over the entire composition range by indirect method using a standard Swietoslawski type ebulliometer. The measurements were made at three sub-atmospheric pressures. The experimental t-x data were correlated with both equation of state method and gamma-phi methods. In the former approach, the Soave equation of state is used, whereas in the latter case Wilson, NRTL and UNIQUAC equations are used.     

GENERALIZED DERIVATIONS FOR MIXTURES USING THE VIRIAL EQUATION: APPLICATION TO FUGACITY

Calculation of mixture properties using the virial equation of state (VEOS) becomes tedious when the expression contains more than the second virial term. This paper describes a derivation procedure which permits relatively simple extension to any number of virial terms. We have used the fugacity derivation as an example of using the new technique and provide an equation for the fugacity coefficient using any number of terms and any number of components.     

PHASE-SPACE DIFFUSION EQUATIONS FOR SINGLE BROWNIAN PARTICLE MOTION NEAR SURFACES

In numerous physical processes involving the motion of micron and submicron sized particles near surfaces, such as the filtration of hydrosols and aerosols, the particle motion is the net result of the combined effects of fluid convection, external forces, particle inertia, Brownian particle motion, and particle-surface fluid dynamic interactions. The most general method of describing particle motion under the combined action of these effects is through the so-called Fokker-Planck equation. In the absence of particle-surface fluid dynamic interactions, the Fokker-Planck equation is well-known, and it has been applied in a general way to problems involving the adsorption or deposition of Brownian particles onto surfaces through a solution technique known as the Brownian dynamics simulation method.

In this study, the Fokker-Planck equation for Brownian particle motion near surfaces is generalized to include particle-surface fluid dynamic interactions. The Fokker-Planck equation is shown to follow from the Liouville equation for the Brownian particle and n-fluid molecules present in the system, thus, establishing a firm theoretical foundation for the Fokker-Planck equation and the various other phase-space diffusion equations that follow from it.

Based on diagonalization of the Fokker-Planck equation, its short-time behavior is also derived here which enables a generalization of the Brownian dynamics method for the study of particle motion near surfaces including fluid dynamic interactions. Additionally, a perturbation solution of the Fokker-Planck equation under the conditions of small, but finite particle Stokes number is also derived. These solutions are shown to agree with previously given representations of the Smoluchowski or convective-diffusion equation for Brownian particle motion near surfaces, as well as with inertial corrections to the Smoluchowski equation available in the literature. This latter equation is also generalized here to include particle-surface fluid dynamic interactions.     

电磁领域中复超越方程的遗传算法求解方法

将遗传算法与参数跟踪策略有效结合,跟踪过程中进行搜索域压缩与位移操作,形成了一种功能强大的新算法,可成功应用于电磁领域中各种各样复超越方程的高精度求解问题.在算法实现过程中,使用参数跟踪策略有效地缩小了搜索区域,保证了解的单一性,提高了运算速度;使用动态搜索域提高了解的精度;应用三阶差商公式预估新的搜索中心,使运算速度得以进一步提高.应用本文发展的算法详细求解了终端短路法测量材料的复介电常数时得到的复超越方程和部分填充矩形波导的特征方程,计算结果表明该算法能够轻松地解决复超越方程中的多值问题,解集完备性好,算法鲁棒性强.     

广义KDV-MKDV方程的精确解

利用直接积分方法将广义KDV-MKDV方程化为一阶变系数非线性常微分方程组,然后用待定系数法确定相应的常数获得了广义KDV-MKDV方程新的精确解;利用先作假设变换后选取试探函数的方法来直接构造广义KDV-MKDV方程新的精确解.     

脉冲微分方程理论在药物动力学中的应用研究

为了揭示药物在机体内的ADME过程,将脉冲微分方程理论应用于研究药物动力学若干问题中,通过对系统动力学行为的研究进而设计最佳用药方案,丰富了脉冲微分方程理论并对指导临床实践具有积极意义.     

小波变换在求解重力测定反问题中的应用

针对重力测定反问题,提出将其离散为线性不适定问题,利用小波变换方法进行数值求解.该方法将小波变换和正则化方法相结合,选取小波函数作为一组基底,将原不适定问题转化为粗子空间上的适定问题,并给出选取粗子空间基的方法.通过数值模拟已有方法和小波变换方法求解结果的比较,表明了小波变换方法的可行性和有效性.