Effects of Some Common Geological Features on Two-Dimensional Variably Saturated Flow

1992 ◽  
Vol 294 ◽  
Author(s):  
Amvrossios C. Bagtzoglou ◽  
Rachid Ababou ◽  
Budhi Sagar ◽  
M. Rashidul Islam
Keyword(s):  
Steady State ◽  
Unsaturated Flow ◽  
Fault Zones ◽  
Oscillatory Behavior ◽  
Stable Steady State ◽  
Two Dimensional ◽  
Time Step ◽  
Groundwater Flux ◽  
Geological Features ◽  
License Application

ABSTRACTThis paper presents results of unsaturated flow simulations undertaken as an auxiliary analysis for the Iterative Performance Assessment (IPA) project, one of the approaches adopted' by the U.S. NRC to develop repository license application review capabilities. The effects on flow of common geological features, such as nonhorizontal stratification and vertical or near-vertical fault zones intersecting the strata, in a two-dimensional (2D) domain are studied. Results indicate that the presence of layers and crosscutting fault zones tend to induce three-dimensional (3D) unstable flows in the unsaturated zone. The instability is manifested in our simulations by an oscillatory behavior of steady state. This numerical instability imposes extremely stringent criteria on the time step used in the simulation. Finally, once stable steady-state solutions are attained, the effect of the crossing point in the matrix-fault unsaturated hydraulic conductivity curve on groundwater flux vectors and moisture content distributions is studied.

An Experimental Model Investigation of the Opening of a Collapsed Untethered Pulmonary Airway

1995 ◽  
Vol 117 (3) ◽  
pp. 245-253 ◽  
Author(s):  
Matthew L. Perun ◽  
Donald P. Gaver
Keyword(s):  
Surface Tension ◽  
Steady State ◽  
Structural Characteristics ◽  
Stable Steady State ◽  
Effective Diameter ◽  
Fluid Viscosity ◽  
Two Dimensional ◽  
Fluid Properties ◽  
Wall Deformation ◽  
Fluid Film Thickness

We developed an essentially two-dimensional planar benchtop model of an untethered collapsed airway to investigate the influence of fluid properties (viscosity, μ and surface tension, γ) and the structural characteristics (effective diameter, D, longitudinal tension, T, and fluid film thickness, H) on airway reopening. This simplified model was used to quantify the relationship between wall deformation and meniscus curvature during reopening. We measured the pressure (P) required to move the meniscus at a constant velocity (U), and found the dimensionless post-startup pressure (PD/γ) increased monotonically with the capillary number (Ca = μU/γ). Startup pressures depend on the fluid viscosity and piston acceleration, and may significantly increase reopening pressures. Consistently stable steady-state pressures existed when Ca > 0.5. D was the most dominant structural characteristic, which caused an increase in the post-startup pressure (P) for a decrease in D. An increase in H caused a slight decrease in the reopening pressure, but a spatial variation in H resulted in only a transient increase in pressure. T did not significantly affect the reopening pressure. From our planar two-dimensional experiments an effective yield pressure of 3.69 γ/D was extrapolated from the steady-state pressures. Based on these results, we predicted airway pressures and reopening times for axisymmetrically collapsed airways under various disease states. These predictions indicate that increasing surface tension (as occurs in Respiratory Distress Syndrome) increases the yield pressure necessary to reopen the airways, and increasing viscosity (as in cystic fibrosis) increases the time to reopen once the yield pressure has been exceeded.

scholarly journals A two-dimensional glacier–fjord coupled model applied to estimate submarine melt rates and front position changes of Hansbreen, Svalbard

2018 ◽  
Vol 64 (247) ◽  
pp. 745-758 ◽  
Author(s):  
E. DE ANDRÉS ◽  
J. OTERO ◽  
F. NAVARRO ◽  
A. PROMIŃSKA ◽  
J. LAPAZARAN ◽  
...  
Keyword(s):  
Coupled Model ◽  
Small Scale ◽  
Two Dimensional ◽  
Time Step ◽  
Software Packages ◽  
Front Position ◽  
Circulation Patterns ◽  
Glacier Front ◽  
Order Of Magnitude ◽  
Frontal Ablation

ABSTRACTWe have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

scholarly journals Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP

2012 ◽  
Vol 6 (3) ◽  
pp. 573-588 ◽  
Author(s):  
F. Pattyn ◽  
C. Schoof ◽  
L. Perichon ◽  
R. C. A. Hindmarsh ◽  
E. Bueler ◽  
...  
Keyword(s):  
Steady State ◽  
Adaptive Mesh Refinement ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
Stable Steady State ◽  
Fixed Grid ◽  
Approximate Analytical Solutions ◽  
Grounding Line ◽  
Intercomparison Project ◽  
Line Positions

Abstract. Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients.

Numerical study of electric field effects on the deformation of two-dimensional liquid drops in simple shear flow at arbitrary Reynolds number

2009 ◽  
Vol 626 ◽  
pp. 367-393 ◽  
Author(s):  
STEFAN MÄHLMANN ◽  
DEMETRIOS T. PAPAGEORGIOU
Keyword(s):  
Steady State ◽  
Shear Flow ◽  
Electric Field ◽  
Electric Fields ◽  
Simple Shear ◽  
Simple Shear Flow ◽  
Physical Parameters ◽  
Shear Stresses ◽  
Two Dimensional ◽  
Liquid Drops

The effect of an electric field on a periodic array of two-dimensional liquid drops suspended in simple shear flow is studied numerically. The shear is produced by moving the parallel walls of the channel containing the fluids at equal speeds but in opposite directions and an electric field is generated by imposing a constant voltage difference across the channel walls. The level set method is adapted to electrohydrodynamics problems that include a background flow in order to compute the effects of permittivity and conductivity differences between the two phases on the dynamics and drop configurations. The electric field introduces additional interfacial stresses at the drop interface and we perform extensive computations to assess the combined effects of electric fields, surface tension and inertia. Our computations for perfect dielectric systems indicate that the electric field increases the drop deformation to generate elongated drops at steady state, and at the same time alters the drop orientation by increasing alignment with the vertical, which is the direction of the underlying electric field. These phenomena are observed for a range of values of Reynolds and capillary numbers. Computations using the leaky dielectric model also indicate that for certain combinations of electric properties the drop can undergo enhanced alignment with the vertical or the horizontal, as compared to perfect dielectric systems. For cases of enhanced elongation and alignment with the vertical, the flow positions the droplets closer to the channel walls where they cause larger wall shear stresses. We also establish that a sufficiently strong electric field can be used to destabilize the flow in the sense that steady-state droplets that can exist in its absence for a set of physical parameters, become increasingly and indefinitely elongated until additional mechanisms can lead to rupture. It is suggested that electric fields can be used to enhance such phenomena.

Two-Dimensional Spectroscopy for Vacuum Arc in Steady State and Decaying State under Free-Recovery Condition

2021 ◽  
Author(s):  
Yuto Hatanaka ◽  
Takamitsu Miyazaki ◽  
Yusuke Nakano ◽  
Yasunori Tanaka ◽  
Yoshihiko Uesugi ◽  
...  
Keyword(s):  
Steady State ◽  
Vacuum Arc ◽  
Two Dimensional ◽  
Recovery Condition

Aeroacoustic Analysis of a UAV Propeller Operable at Various Altitudes

2021 ◽  
Author(s):  
Ji-Hun Song ◽  
Seungsoo Jang ◽  
Youn-Jea Kim
Keyword(s):  
Steady State ◽  
Turbulence Model ◽  
Technological Development ◽  
Basic Research ◽  
Noise Pollution ◽  
Maximum Error ◽  
Unmanned Aircraft ◽  
Time Step ◽  
Application Range ◽  
Noise Characteristics

Abstract With technological development and the wide application range of unmanned aerial vehicles (UAVs), the regulation of UAV altitude limits in many countries is further alleviated, and the problem of UAV noise pollution has emerged with the recent advent of urban air mobility (UAM) and personal air vehicle (PAV) markets. In this study, one typical propeller, the T-motor 15 × 5 propeller, was analyzed by use of the commercial CFD software, ANSYS FLUENT V19.3. The effects of gravity and convection were analyzed to determine the noise characteristics at altitude using the FW-H equation. A high-altitude drone, which operates at heights from 0 to 10 km with 1,000 to 5,000 revolutions per minute, was analyzed using the steady-state k-ω SST turbulence model. And using the steady-state data to initialize values, an unsteady analysis was performed with the LES turbulence model. The time step was divided based on the 1-degree rotational time, and the velocity residual on each axis was calculated until a value of 10−7 or less was achieved and there was no fluctuation of thrust, at which point it was considered converged. The CFD results were validated with the experimental results for thrust and their results show that the maximum error was 8.64%. The overall sound pressure level was calculated, and noise characteristics in the audible frequency range according to receiver points were also compared. Through this study, thrust and noise data according to altitude were provided. The aerodynamic and aeroacoustic characteristics at high-altitudes, which are generally difficult to measure by experiment, are also presented. Therefore, the appropriate operating altitudes and rotational speeds will be presented through the aeroacoustics analysis corresponding to operational altitude, and the basic research data can then be applied to upcoming unmanned aircraft system (UAS) market.

Enhanced Explicit Scheme to Solve Transient Heat Conduction Problem

2007 ◽  
Vol 2 (1) ◽  
Author(s):  
Ganesh Hegde ◽  
Madhu Gattumane
Keyword(s):  
Heat Conduction ◽  
Correction Factor ◽  
Truncation Error ◽  
Transient Heat Conduction ◽  
Explicit Scheme ◽  
Stability And Convergence ◽  
Two Dimensional ◽  
Time Step ◽  
One Dimensional ◽  
Partial Derivatives

Improvement in accuracy without sacrificing stability and convergence of the solution to unsteady diffusion heat transfer problems by computational method of enhanced explicit scheme (EES), has been achieved and demonstrated, through transient one dimensional and two dimensional heat conduction. The truncation error induced in the explicit scheme using finite difference technique is eliminated by optimization of partial derivatives in the Taylor series expansion, by application of interface theory developed by the authors. This theory, in its simple terms gives the optimum values to the decision vectors in a redundant linear equation. The time derivatives and the spatial partial derivatives in the transient heat conduction, take the values depending on the time step chosen and grid size assumed. The time correction factor and the space correction factor defined by step sizes govern the accuracy, stability and convergence of EES. The comparison of the results of EES with analytical results, show decreased error as compared to the result of explicit scheme. The paper has an objective of reducing error in the explicit scheme by elimination of truncation error introduced by neglecting the higher order terms in the expansion of the governing function. As the pilot examples of the exercise, the implementation is aimed at solving one-dimensional and two-dimensional problems of transient heat conduction and compared with the results cited in the referred literature.

Sign in / Sign up

Export Citation Format

Share Document