Additional scaled solutions to Richards' equation for infiltration and drainage

Warrick and Hussen developed in the nineties of the last century a method to scale Richards' equation (RE) for similar soils. In this paper, new scaled solutions are added to the method of Warrick and Hussen considering a wider range of soils regardless of their dissimilarity. Gardner-Kozeny hydraulic functions are adopted instead of Brooks-Corey functions used originally by Warrick and Hussen. These functions allow to reduce the dependence of the scaled RE on the soil properties. To evaluate the proposed method (PM), the scaled RE was solved numerically using a finite difference method with a fully implicit scheme. Three cases were considered: constant-head infiltration, constant-flux infiltration, and drainage of an initially uniform wet soil. The results for five texturally different soils ranging from sand to clay (adopted from the literature) showed that the scaled solutions were invariant to a satisfactory degree. However, slight deviations were observed mainly for the sandy soil. Moreover, the scaled solutions deviated when the soil profile was initially wet in the infiltration case or when deeply wet in the drainage condition. Based on the PM, a Philip-type model was also developed to approximate RE solutions for the constant-head infiltration. The model showed a good agreement with the scaled RE for the same range of soils and conditions, however only for Gardner-Kozeny soils. Such a procedure reduces numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils. (C) 2011 Elsevier B.V. All rights reserved.

Noninvasive pressure difference mapping derived from 4D flow MRI in patients with unrepaired and repaired aortic coarctation.

PURPOSE To develop a method for computing and visualizing pressure differences derived from time-resolved velocity-encoded three-dimensional phase-contrast magnetic resonance imaging (4D flow MRI) and to compare pressure difference maps of patients with unrepaired and repaired aortic coarctation to young healthy volunteers. METHODS 4D flow MRI data of four patients with aortic coarctation either before or after repair (mean age 17 years, age range 3-28, one female, three males) and four young healthy volunteers without history of cardiovascular disease (mean age 24 years, age range 20-27, one female, three males) was acquired using a 1.5-T clinical MR scanner. Image analysis was performed with in-house developed image processing software. Relative pressures were computed based on the Navier-Stokes equation. RESULTS A standardized method for intuitive visualization of pressure difference maps was developed and successfully applied to all included patients and volunteers. Young healthy volunteers exhibited smooth and regular distribution of relative pressures in the thoracic aorta at mid systole with very similar distribution in all analyzed volunteers. Patients demonstrated disturbed pressures compared to volunteers. Changes included a pressure drop at the aortic isthmus in all patients, increased relative pressures in the aortic arch in patients with residual narrowing after repair, and increased relative pressures in the descending aorta in a patient after patch aortoplasty. CONCLUSIONS Pressure difference maps derived from 4D flow MRI can depict alterations of spatial pressure distribution in patients with repaired and unrepaired aortic coarctation. The technique might allow identifying pathophysiological conditions underlying complications after aortic coarctation repair.

Noninvasive 4D pressure difference mapping derived from 4D flow MRI in patients with repaired aortic coarctation: comparison with young healthy volunteers

To assess spatial and temporal pressure characteristics in patients with repaired aortic coarctation compared to young healthy volunteers using time-resolved velocity-encoded three-dimensional phase-contrast magnetic resonance imaging (4D flow MRI) and derived 4D pressure difference maps. After in vitro validation against invasive catheterization as gold standard, 4D flow MRI of the thoracic aorta was performed at 1.5T in 13 consecutive patients after aortic coarctation repair without recoarctation and 13 healthy volunteers. Using in-house developed processing software, 4D pressure difference maps were computed based on the Navier-Stokes equation. Pressure difference amplitudes, maximum slope of pressure amplitudes and spatial pressure range at mid systole were retrospectively measured by three readers, and twice by one reader to assess inter- and intraobserver agreement. In vitro, pressure differences derived from 4D flow MRI showed excellent agreement to invasive catheter measurements. In vivo, pressure difference amplitudes, maximum slope of pressure difference amplitudes and spatial pressure range at mid systole were significantly increased in patients compared to volunteers in the aortic arch, the proximal descending and the distal descending thoracic aorta (p < 0.05). Greatest differences occurred in the proximal descending aorta with values of the three parameters for patients versus volunteers being 19.7 ± 7.5 versus 10.0 ± 2.0 (p < 0.001), 10.9 ± 10.4 versus 1.9 ± 0.4 (p = 0.002), and 8.7 ± 6.3 versus 1.6 ± 0.9 (p < 0.001). Inter- and intraobserver agreements were excellent (p < 0.001). Noninvasive 4D pressure difference mapping derived from 4D flow MRI enables detection of altered intraluminal aortic pressures and showed significant spatial and temporal changes in patients with repaired aortic coarctation.

Extension of the characteristic equation to absorption chillers with adiabatic absorbers

Various researchers have developed models of conventional H2O–LiBr absorption machines with the aim of predicting their performance. In this paper, the methodology of characteristic equations developed by Hellmann et al. (1998) is applied. This model is able to represent the capacity of single effect absorption chillers and heat pumps by means of simple algebraic equations. An extended characteristic equation based on a characteristic temperature difference has been obtained, considering the facility features. As a result, it is concluded that for adiabatic absorbers a subcooling temperature must be specified. The effect of evaporator overflow has been characterized. Its influence on cooling capacity has been included in the extended characteristic equation. Taking into account the particular design and operation features, a good agreement between experimental performance data and those obtained through the extended characteristic equation has been achieved at off-design operation. This allows its use for simulation and control purposes.

Global air-sea flux of CO2: An estimate based on measurements of sea–air pCO2 difference

Approximately 250,000 measurements made for the pCO2 difference between surface water and the marine atmosphere, ΔpCO2, have been assembled for the global oceans. Observations made in the equatorial Pacific during El Nino events have been excluded from the data set. These observations are mapped on the global 4° × 5° grid for a single virtual calendar year (chosen arbitrarily to be 1990) representing a non-El Nino year. Monthly global distributions of ΔpCO2 have been constructed using an interpolation method based on a lateral advection–diffusion transport equation. The net flux of CO2 across the sea surface has been computed using ΔpCO2 distributions and CO2 gas transfer coefficients across sea surface. The annual net uptake flux of CO2 by the global oceans thus estimated ranges from 0.60 to 1.34 Gt-C⋅yr−1 depending on different formulations used for wind speed dependence on the gas transfer coefficient. These estimates are subject to an error of up to 75% resulting from the numerical interpolation method used to estimate the distribution of ΔpCO2 over the global oceans. Temperate and polar oceans of the both hemispheres are the major sinks for atmospheric CO2, whereas the equatorial oceans are the major sources for CO2. The Atlantic Ocean is the most important CO2 sink, providing about 60% of the global ocean uptake, while the Pacific Ocean is neutral because of its equatorial source flux being balanced by the sink flux of the temperate oceans. The Indian and Southern Oceans take up about 20% each.

The development of a numerical solution to transport equation /

"Prepared for U.S. Army Engineer District, Mobile."

On the Schrodinger equation involving a critical Sobolev exponent and magnetic field

We consider the semilinear Schrodinger equation -Delta(A)u + V(x)u = Q(x)vertical bar u vertical bar(2* -2) u. Assuming that V changes sign, we establish the existence of a solution u not equal 0 in the Sobolev space H-A,V(1) + (R-N). The solution is obtained by a min-max type argument based on a topological linking. We also establish certain regularity properties of solutions for a rather general class of equations involving the operator -Delta(A).

Computer methods for the heat conduction equation

The merits of various numerical methods for the solution of the one and two dimensional heat conduction equation with a radiation boundary condition have been examined from a practical standpoint in order to determine accuracies and efficiencies. It is found that the use of five increments to approximate the space derivatives gives sufficiently accurate results provided the time step is not too large; further, the implicit backward difference method of Liebmann (27) is found to be the most accurate method. On this basis, a new implicit method is proposed for the solution of the three-dimensional heat conduction equation with radiation boundary conditions. The accuracies of the integral and analogue computer methods are also investigated.

Optimal Control of a Second Order Parabolic Heat Equation

In this paper, we are concerned with the optimal control boundary control of a second order parabolic heat equation. Using the results in [Evtushenko, 1997] and spatial central finite difference with diagonally implicit Runge-Kutta method (DIRK) is applied to solve the parabolic heat equation. The conjugate gradient method (CGM) is applied to solve the distributed control problem. Numerical results are reported.

Nonstandard Finite Difference Schemes with Application to Finance: Option Pricing

2000 Mathematics Subject Classification: 65M06, 65M12.

Low Volatility Options and Numerical Diffusion of Finite Difference Schemes

2000 Mathematics Subject Classification: 65M06, 65M12.

Graded Meshes In Bio-thermal Problems With Transmission-line Modeling Method.

In this study, the transmission-line modeling (TLM) applied to bio-thermal problems was improved by incorporating several novel computational techniques, which include application of graded meshes which resulted in 9 times faster in computational time and uses only a fraction (16%) of the computational resources used by regular meshes in analyzing heat flow through heterogeneous media. Graded meshes, unlike regular meshes, allow heat sources to be modeled in all segments of the mesh. A new boundary condition that considers thermal properties and thus resulting in a more realistic modeling of complex problems is introduced. Also, a new way of calculating an error parameter is introduced. The calculated temperatures between nodes were compared against the results obtained from the literature and agreed within less than 1% difference. It is reasonable, therefore, to conclude that the improved TLM model described herein has great potential in heat transfer of biological systems.

Nasal Potential Difference In Cystic Fibrosis Considering Severe Cftr Mutations.

The gold standard for diagnosing cystic fibrosis (CF) is a sweat chloride value above 60 mEq/L. However, this historical and important tool has limitations; other techniques should be studied, including the nasal potential difference (NPD) test. CFTR gene sequencing can identify CFTR mutations, but this method is time-consuming and too expensive to be used in all CF centers. The present study compared CF patients with two classes I-III CFTR mutations (10 patients) (G1), CF patients with classes IV-VI CFTR mutations (five patients) (G2), and 21 healthy subjects (G3). The CF patients and healthy subjects also underwent the NPD test. A statistical analysis was performed using the Mann-Whitney, Kruskal-Wallis, χ(2), and Fisher's exact tests, α = 0.05. No differences were observed between the CF patients and healthy controls for the PDMax, Δamiloride, and Δchloride + free + amiloride markers from the NPD test. For the finger value, a difference between G2 and G3 was described. The Wilschanski index values were different between G1 and G3. In conclusion, our data showed that NPD is useful for CF diagnosis when classes I-III CFTR mutations are screened. However, if classes IV-VI are considered, the NPD test showed an overlap in values with healthy subjects.

Modelo para análise do ataque no futebol = : A model to analyze soccer offensive phase

Universidade Estadual de Campinas . Faculdade de Educação Física