Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

On the Importance of Geolithological Features for the Estimate of the Site Response: The Case of Catania Metropolitan Area (Italy)

The evaluation of seismic site response in the urban area of Catania was tackled by selecting test areas having peculiar lithological and structural features, potentially favourable to large local amplifications of ground motion. The two selected areas are located in the historical downtown and in the northern part of Catania where the presence of a fault is evident. Site response was evaluated using spectral ratio technique taking the horizontal- to-vertical component ratio of ambient noise. Inferences from microtremor measurements are compared with results from synthetic accelerograms and response spectra computed at all drillings available for this area. Such method is particularly suitable in urban areas where the nature of the outcropping geological units is masked by city growth and anthropic intervention on the surface geology. The microtremor H/V spectral ratios evaluated at soft sites located within the downtown profile tend to be smaller than that usually reported in the literature for such soils. A tendency for amplifications to peaks near 2 Hz is observed only in some sites located on recent alluvial deposits. Evidences for amplifications of site effects (frequency range 4–8 Hz) were observed in the sampling sites located on the fault, with a rapid decrease of spectral amplitude just a few tenths of metres away from the discontinuity. Numerical simulations evidenced the importance of geolithological features at depth levels even greater than 20–30 m. Besides this, the results strongly confirm the importance of the subsurface geological conditions, in the estimate of seismic hazard at urban scale.     

Computer controlled cutting of multiple identical joints in real rock

Summary This novel method for reproducibly cutting joints in real rock will allow shearing of identical joints under varying stress paths. Shearing of such identical joints will greatly reduce the effect of sample variation and the need for scaling results.Joint surfaces are cut with a diamond band saw through a computer controlled feed table. As constructed, the system is capable of cutting samples 20.3 cm (8 in) thick and 30.5 cm (12 in) long. The 6.35 mm (0.25 in) blade used in these tests allows cuts with a maximum radius of curvature of 1.6 cm (0.625 in). A maximum asperity inclination (i) of 28 degrees is possible with the largest samples. These characteristics, in general, are not inherent in the method but rather imposed by the current dimensions of the device.The resulting surfaces, which have been produced in limestone, are analyzed in terms of surface finish at a micro scale, geometrical conformity to the desired surface, reproducibility of the surface and matchability of the two surfaces produced by one cut. These measures indicate a high level of accuracy which should produce repeatable test results when the joints are sheared under laboratory conditions.     

Argo floats revealing bimodality of large-scale mid-depth circulation in the North Atlantic

Analysis of Argo float trajectories at 1 000 m and temperature at 950 m in the North Atlantic between November 2003 and January 2005 demonstrates the existence of two different circulation modes with fast transition between them. Each mode has a pair of cyclonic - anticyclonic gyres. The difference is the location of the cyclonic gyre. The cyclonic gyre stretches from southeast to northwest in the first mode and from the southwest to the northeast in the second mode. The observed modes strongly affect the heat and salt transport in the North Atlantic. In particular, the second mode slows down the westward transport of the warm and saline water from the Mediterranean Sea.     

LU Vel (GJ 375): A M3.5Ve Flare and Double-Lined Spectroscopic Binary

High resolution echelle spectroscopic observations taken with the FEROS spectrograph at the 2.2 m telescope ESO confirm the binary nature of the flare M3.5V star LU Vel (GJ 375, RE J0958-462) previously reported by Christian and Mathioudakis (2002). Emission of similar intensity from both components is detected in the Balmer, Na i D₁&D₂, He i D₃, Ca ii H&K, and Ca ii IRT lines. We have determined precise radial velocities by cross correlation with radial velocity standard stars, which have allowed us to obtain for the first time the orbital solution of the system. The binary consists of two near-equal M3.5V components with an orbital period shorter than 2 days. We have analyzed the behaviour of the chromospheric activity indicators (variability and possible flares). In addition, we have determined its rotational velocity and kinematics.     

Optimization of SAGD Well Elevation

The application of steam-assisted gravity drainage (SAGD) to recover heavy oil sands is becoming increasingly important in the northern Alberta McMurray Formation because of the vast resources/reserves accessible with this mechanism. Selecting the stratigraphic elevations of SAGD well pairs is a vital decision for reservoir evaluation and planning. The inherent uncertainty in the distribution of geological variables significantly influences this decision. Geostatistical simulation is used to capture geological uncertainty, which is used can be transformed into a distribution of the best possible well pair elevations. A simple exhaustive calculation scheme is used to determine the optimum stratigraphic location of a SAGD well pair where the recovery R is maximized. There are three basic steps to the methodology: (1) model the uncertainty in the top continuous bitumen (TCB) and bottom continuous bitumen (BCB) surfaces, (2) calculate the recovery at all possible elevation increments within the TCB and BCB interval, and (3) identify the elevation that maximizes R. This is repeated for multiple TCB/BCB pairs of surfaces to assess uncertainty. The methodology is described and implemented on a subset of data from the Athabasca Oilsands in Fort McMurray, Alberta.     

The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

A One-dimensional Ensemble Forecast and Assimilation System for Fog Prediction

A probabilistic fog forecast system was designed based on two high resolution numerical 1-D models called COBEL and PAFOG. The 1-D models are coupled to several 3-D numerical weather prediction models and thus are able to consider the effects of advection. To deal with the large uncertainty inherent to fog forecasts, a whole ensemble of 1-D runs is computed using the two different numerical models and a set of different initial conditions in combination with distinct boundary conditions. Initial conditions are obtained from variational data assimilation, which optimally combines observations with a first guess taken from operational 3-D models. The design of the ensemble scheme computes members that should fairly well represent the uncertainty of the current meteorological regime. Verification for an entire fog season reveals the importance of advection in complex terrain. The skill of 1-D fog forecasts is significantly improved if advection is considered. Thus the probabilistic forecast system has the potential to support the forecaster and therefore to provide more accurate fog forecasts.     

MANAGEMENT AND SUBSURFACE EFFECTS ON RUNOFF AND SEDIMENT YIELD FROM SMALL WATERSHEDS

lINTR0DUCTIONSoiler0si0ncanbeaseriousproblemonsl0pinguplands0fthes0uthemUnitedStatesthatareintensivelycr0pped.Conservati0nprovisionsofthel985U.S.FarmBillmandateddevelopmentofaconservati0nplanonhighlyer0diblelands.Researchwasinitiatedinl9870na28-haareaattheA-E.Nels0nFarrninn0rthernMississippit0determineiftheconservati0nprovisionscouldbeachievedwhileprofitablygr0wingrowcropsincludingcott0n(GosSyPiumhirsutumL.),soybean(GlyciDemax(L.)Merr.),sorghum(SOrghumbicolor(L.)Moench),c0rn(Zeam…