Sensitivity Analysis of Dimensionless Parameters for Physical Simulation of Water-Flooding Reservoir

A numerical approach to optimize dimensionless parameters of water-flooding porous media flows is proposed based on the analysis of the sensitivity factor defined as the variation ration of a target function with respect to the variation of dimensionless parameters. A complete set of scaling criteria for water-flooding reservoir of five-spot well pattern case is derived from the 3-D governing equations, involving the gravitational force, the capillary force and the compressibility of water, oil and rock. By using this approach, we have estimated the influences of each dimensionless parameter on experimental results and thus sorted out the dominant ones with larger sensitivity factors ranging from10-4to10-0 .

Long-Range Stress Redistribution Resulting From Damage in Heterogeneous Media

It has been shown in CA simulations and data analysis of earthquakes that declustered or characteristic large earthquakes may occur with long-range stress redistribution. In order to understand long-range stress redistribution, we propose a linear-elastic but heterogeneous-brittle model. The stress redistribution in the heterogeneous-brittle medium implies a longer-range interaction than that in an elastic medium. Therefore, it is surmised that the longer-range stress redistribution resulting from damage in heterogeneous media may be a plausible mechanism governing main shocks.

A Nonlinear Threshold Model Applied to Spallation Analysis

Spallation in heterogeneous media is a complex, dynamic process. Generally speaking, the spallation process is relevant to multiple scales and the diversity and coupling of physics at different scales present two fundamental difficulties for spallation modeling and simulation. More importantly, these difficulties can be greatly enhanced by the disordered heterogeneity on multi-scales. In this paper, a driven nonlinear threshold model for damage evolution in heterogeneous materials is presented and a trans-scale formulation of damage evolution is obtained. The damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that some dimensionless numbers govern the whole process of deformation and damage evolution. The effects of heterogeneity in terms of Weibull modulus on damage evolution in spallation process are also investigated.

Stability Analysis Of The Moving Interface In Piston- And Non-Piston-Like Displacements

From the macroscopic point of view, expressions involving reservoir and operational parameters are established for investigating the stability of moving interface in piston- and non-piston-like displacements. In the case of axi-symmetrical piston-like displacement, the stability is related to the moving interface position and water to oil mobility ratio. The capillary effect on the stability of moving interface depends on whether or not the moving interface is already stable and correlates with the wettability of the reservoir rock. In the case of non-piston-like displacement, the stability of the front is governed by both the relative permeability and the mobility ratio.

Sensitivity Analysis Of The Dimensionless Parameters In Scaling A Polymer Flooding Reservoir

A set of scaling criteria of a polymer flooding reservoir is derived from the governing equations, which involve gravity and capillary force, compressibility of water, oil, and rock, non-Newtonian behavior of the polymer solution, absorption, dispersion, and diffusion, etc. A numerical approach to quantify the dominance degree of each dimensionless parameter is proposed. With this approach, the sensitivity factor of each dimensionless parameter is evaluated. The results show that in polymer flooding, the order of the sensitivity factor ranges from 10(-5) to 10(0) and the dominant dimensionless parameters are generally the ratio of the oil permeability under the condition of the irreducible water saturation to water permeability under the condition of residual oil saturation, density, and viscosity ratios between water and oil, the reduced initial oleic phase saturation and the shear rate exponent of the polymer solution. It is also revealed that the dominant dimensionless parameters may be different from case to case. The effect of some physical variables, such as oil viscosity, injection rate, and permeability, on the dominance degree of the dimensionless parameters is analyzed and the dominant ones are determined for different cases.

Fractal dimensions for unsaturated porous media

The analytical expressions of the fractal dimensions for wetting and non-wetting phases for unsaturated porous media are derived and are found to be a function of porosity, maximum and minimum pore sizes as well as saturation. There is no empirical constant in the proposed fractal dimensions. It is also found that the fractal dimensions increase with porosity of a medium and are meaningful only in a certain range of saturation S-w, i.e. S-w > S-min for wetting phase and S-w < S-max for non-wetting phase at a given porosity, based on real porous media for requirements from both fractal theory and experimental observations. The present analysis of the fractal dimensions is verified to be consistent with the existing experimental observations and it makes possible to analyze the transport properties such as permeability, thermal dispersion in unsaturated porous media by fractal theory and technique.

El estupro. Sexualidad delictiva en la Baja Edad Media y primera Edad Moderna

[ES]Análisis del delito de estupro en los siglos XIV y XVII en los reinos de Castilla y de Navarra a través de la legislación y las fuentes judiciales.

Propagation of Explosion Wave in Hard-Soft-Hard Layered Media

The paper presents a reasonable analysis for dynamic response and failure process of a plane multi-layered media, which are subjected to a blast loading. This blast loading is induced by a cylindric explosive put on the center of top surface of the layered media. With the help of numerical simulation technique provided by LS-DYNA software, the whole process of explosion wave propagation and attenuation can be revealed. The feature of local failure around the blasting site is also discussed in some detail. Our focus will be on the explosion wave attenuation for the hard-soft-hard sandwich layers. As seen in the paper, the computational results are delivered in a feasible way by comparing with experimental data.

A quantitative analysis of the phytoplankton of the Gulf of Panama: I. Results of the regional phytoplankton surveys during July and November 1957 and March 1958

ENGLISH: 1. Quantitative phytoplankton samples were collected by the Inter-American Tropical Tuna Commission at the surface and ten meters in the Gulf of Panama, as follows: a) 18-21 March, 1958 (31 stations)-during the height of the upwelling season, b) 10-12 July, 1957 (10 stations)-during the transition to the rainy season at a time when mild upwelling winds reappear, c) 7-8 November, 1957 (15 stations)-during the height of the rainy season. 2. Maximum phytoplankton populations occurred during the upwelling season, followed by a considerable decline during July, and a further Subsidence during November. 3. A remarkable regional uniformity in species composition was observed during the surveys despite regional differences in growth conditions. Diatoms overwhelmingly dominated the communities. 4. During all surveys, the innermost regions, generally north of 8°30'N, were the most productive. The least productive areas were in the offing of San Miguel Bay and Parita Bay, suggesting that nutrient accretion via runoff is inadequate to sustain sizeable autotrophic plant populations in those regions. 5. During all surveys, phytoplankton growth appeared to be limited by nutrient availability. 6. During all surveys, phytoplankton growth appeared to be related to depth of the water column. 7. Although below average rainfall contributed to unusually favorable growth conditions (reduced stability, increased transparency and, presumably, nutrient reserves) during the November survey relative to November 1955 and 1956 at 8°45'N, 79°23'W, the anticipated heightened phytoplankton response was not observed. 8. During the November survey, the local diatom responses and their regional fluctuations could be satisfactorily related to the accompanying surface salinity conditions. However, this correspondence is undoubtedly attributable to factors associated with the observed salinity levels, probably nutrients, rather than salinity directly. 9. Unusually warm conditions occurred during the March survey, attributable to considerably weaker upwelling winds than normally occurring then, which contributed to a considerably lower standing crop and a retardation in succession of three to five weeks relative to that observed during 1955-1957 at 8°45'N, 79°23'W in the Gulf of Panama. 10. During the March survey, a well defined inverse relationship existed between mean temperature and mean diatom abundance in the upper ten meters, and between transparency and mean diatom abundance. A direct relationship occurred between surface salinity and mean diatom abundance in the upper ten meters. These relationships are interpreted to indicate that diatom abundance primarily reflected the nutrient concentrations associated with a given upwelling intensity, rather than describing casual relationships. 11. The survey results indicate that the phytoplankton dynamics observed at 8°45'N, 79°23'W from November, 1954 through May, 1957 are generally representative of the Gulf of Panama. 12. The following new forms, to be described in a later publication, were observed during the surveys: Actinoptychus undulatus f. catenata n.f., Asterionella japonica f. tropicum n.f., Leptocylindrus maximus n. sp., Skeletonema costatum f. tropicum n.f. SPANISH: 1. La Comisión Interamericana del Atun Tropical recolectó en el Golfo de Panama muestras cuantitativas de fitoplancton en la superficie y a los diez metros, como sigue: a) Del 18 al 21 de marzo de 1958 (31 estaciones)-durante el maximum de la estación de afloramiento. b) Del 10 al 12 de julio de 1957 (10 estaciones)-durante la epóca de transición a la estación lluviosa cuando reaparecen los vientos ligeros que causan el afloramiento. c) Del 7 al 8 de noviembre de 1957 (15 estaciones)-durante el maximum de la estación lluviosa. 2. Las poblaciones maximas de fitoplancton aparecieron durante la estación de afloramiento, seguido por una considerable disminución durante el mes de julio y una calma durante noviembre. 3. Durante la investigación se observó una remarcable uniformidad regional en la composición de las especies a pesar de las diferencias regionales en las condiciones de crecimiento. Las diatomeas predominaban en gran numero en las comunidades. 4. Durante todas las investigaciones, las regiones mas cerca de la costa, generalmente al norte de los 8°30'N, eran las mas productivas. Las areas menos productivas fueron las mar afuera de las Bahias de San Miguel y Parita, lo que sugiere que el aumento en las sales nutritivas causado por las escorrentias es inadecuado para sostener poblaciones grandes de plantas autotróficas en estas regiones. 5. Durante todas las investigaciones, el crecimiento del fitoplancton parecio estar limitado por la disponibilidad de las. sales nutritivas. 6. Durante todas las investigaciones el crecimiento del fitoplancton parecio estar relacionado con la profundidad de la columna de agua. 7. Aunque las precipitacion por debajo del promedio normal contribuyo a condiciones desusadamente favorables de crecimiento (estabilidad reducida, aumento de la transparencia y, presumiblemente, de la reserva de sales nutritivas) durante la investigación de noviembre en relación a noviembre de 1955 y de 1956 en los 8°45'N, 79°23'W, no se observo-la alta reacción de fitoplancton que se esperaba. 8. Durante la investigación de noviembre, las reacciones locales de las diatomeas y sus fluctuaciones regionales pudieron relacionarse en forma satisfactoria con condiciones asociadas con la salinidad de la superficie. Sin embargo, esta correspondencia puede atribuirse sin duda a factores asociados con los niveles observados de salinidad, probablemente con las sales nutritivas, en lugar de directamente con la salinidad. 9. Condiciones calurosas no comunes ocurrieron durante la investigación de marzo, las que pueden atribuirse a que los vientos que ocasionan el afloramiento fueran mas debiles que los normales, lo que contribuyó a que la cosecha estable fuera considerablemente mas baja y a la demora de tres a cinco semanas en la sucecion relativa a la que se observó durante 1955-1957 en los 8°45'N, 8°23'W, en el Golfo de Panama. 10. Durante la investigación de marzo, existió una relación inversa bien definida entre la temperatura y la abundancia media de las diatomeas en los diez metros superiores, y entre la transparencia y la abundancia media de las diatomeas. Una relación directa ocurrio entre la salinidad de superficie y la abundancia media de las diatomeas en los diez metros superiores. Estas relaciones se interpretan como indicadoras de que la abundancia de diatomeas refleja primeramente las concentraciones de las sales nutritivas asociadas con una intensidad de afloramiento dada, en lugar de describir relaciones causales. 11. Los resultados de la investigacion indican que la dinamica del fitoplancton observada en los 8°45'N, 79°23'W, desde noviembre de 1954 a mayo de 1957, es generalmente representativa del Golfo de Panama. 12. Durante las investigaciones se observaron las siguientes formas nuevas, las que seran descritas en una publicación posterior: Actinoptychus undulatus f. catenata n.f., Asterionella japonica f. tropicum n.f., Leptocylindrus maximus n. sp., Skeletonema costatum f. tropicum n.f.

Analysis of variance and standardization of longline hook rates of bigeye (Thunnus obesus) and yellowfin (Thunnus albacares) tunas in the Eastern Pacific Ocean during 1975-1987

ENGLISH: Longline hook rates of bigeye and yellowfin tunas in the eastern Pacific Ocean were standardized by maximum depth of fishing, area, and season, using generalized linear models (GLM's). The annual trends of the standardized hook rates differ from the unstandardized, and are more likely to represent the changes in abundance of tunas in the age groups most vulnerable to longliners in the fishing grounds. For both species all of the interactions in the GLM's involving years, depths of fishing, areas, and seasons were significant. This means that the annual trends in hook rates depend on which depths, areas, and seasons are being considered. The overall average hook rates for each were estimated by weighting each 5-degree quadrangle equally and each season by the number of months in it. Since the annual trends in hook rates for each fishing depth category are roughly the same for bigeye, total average annual hook rate estimates are possible with the GLM. For yellowfin, the situation is less clear because of a preponderance of empty cells in the model. The full models explained 55% of the variation in bigeye hook rate and 33% of that of yellowfin. SPANISH: Se estandardizaron las tasas de captura con palangre de atunes patudo y aleta amarilla en el Océano Pacífico oriental por la profunidad máxima de pesca, área, y temporada, usando modelos lineales generalizados (MLG). Las tendencias anuales de las tasas de captura estandardizadas son diferentes a las de las tasas no estandardizadas, y es más que representen los cambios en la abundancia de los atunes en los grupos de edad más vulnerables a los palangreros en las áreas de pesca. Para ambas especies fueron significativas todas las interacciones en los MLG con año, profundidad de pesca, área, y temporada. Esto significa que las tendencias anuales de las tasas de captura dependen de cuál profundidad, área, y temporado se está considerando. Para la estimación de la tasa de captura general media para cada especie se ponderó cada cuadrángulo de 5 grados igualmente y cada temporada por el número de meses que contiene. Ya que las tendencias anuales en las tasas de captura para cada categoría de profundidad de pesca son aproximadamente iguales para el patudo, son posibles estimaciones de la tasa de captura anual media total con el MLG. En el caso del aleta amarilla, la situación es más confusa, debido a una preponderancia de celdas vacías en el modelo. Los modelos completos explican el 55% de la variación de la tasa de captura de patudo y 33% de la del aleta amarilla. (PDF contains 19 pages.)

An evaluation of length-frequency sampling procedures and subsequent data analysis for purse seine-caught yellowfin tuna in the Eastern Pacific Ocean

ENGLISH: A two-stage sampling design is used to estimate the variances of the numbers of yellowfin in different age groups caught in the eastern Pacific Ocean. For purse seiners, the primary sampling unit (n) is a brine well containing fish from a month-area stratum; the number of fish lengths (m) measured from each well are the secondary units. The fish cannot be selected at random from the wells because of practical limitations. The effects of different sampling methods and other factors on the reliability and precision of statistics derived from the length-frequency data were therefore examined. Modifications are recommended where necessary. Lengths of fish measured during the unloading of six test wells revealed two forms of inherent size stratification: 1) short-term disruptions of existing pattern of sizes, and 2) transition zones between long-term trends in sizes. To some degree, all wells exhibited cyclic changes in mean size and variance during unloading. In half of the wells, it was observed that size selection by the unloaders induced a change in mean size. As a result of stratification, the sequence of sizes removed from all wells was non-random, regardless of whether a well contained fish from a single set or from more than one set. The number of modal sizes in a well was not related to the number of sets. In an additional well composed of fish from several sets, an experiment on vertical mixing indicated that a representative sample of the contents may be restricted to the bottom half of the well. The contents of the test wells were used to generate 25 simulated wells and to compare the results of three sampling methods applied to them. The methods were: (1) random sampling (also used as a standard), (2) protracted sampling, in which the selection process was extended over a large portion of a well, and (3) measuring fish consecutively during removal from the well. Repeated sampling by each method and different combinations indicated that, because the principal source of size variation occurred among primary units, increasing n was the most effective way to reduce the variance estimates of both the age-group sizes and the total number of fish in the landings. Protracted sampling largely circumvented the effects of size stratification, and its performance was essentially comparable to that of random sampling. Sampling by this method is recommended. Consecutive-fish sampling produced more biased estimates with greater variances. Analysis of the 1988 length-frequency samples indicated that, for age groups that appear most frequently in the catch, a minimum sampling frequency of one primary unit in six for each month-area stratum would reduce the coefficients of variation (CV) of their size estimates to approximately 10 percent or less. Additional stratification of samples by set type, rather than month-area alone, further reduced the CV's of scarce age groups, such as the recruits, and potentially improved their accuracy. The CV's of recruitment estimates for completely-fished cohorts during the 198184 period were in the vicinity of 3 to 8 percent. Recruitment estimates and their variances were also relatively insensitive to changes in the individual quarterly catches and variances, respectively, of which they were composed. SPANISH: Se usa un diseño de muestreo de dos etapas para estimar las varianzas de los números de aletas amari11as en distintos grupos de edad capturados en el Océano Pacifico oriental. Para barcos cerqueros, la unidad primaria de muestreo (n) es una bodega de salmuera que contenía peces de un estrato de mes-área; el numero de ta11as de peces (m) medidas de cada bodega es la unidad secundaria. Limitaciones de carácter practico impiden la selección aleatoria de peces de las bodegas. Por 10 tanto, fueron examinados los efectos de distintos métodos de muestreo y otros factores sobre la confiabilidad y precisión de las estadísticas derivadas de los datos de frecuencia de ta11a. Se recomiendan modificaciones donde sean necesarias. Las ta11as de peces medidas durante la descarga de seis bodegas de prueba revelaron dos formas de estratificación inherente por ta11a: 1) perturbaciones a corto plazo en la pauta de ta11as existente, y 2) zonas de transición entre las tendencias a largo plazo en las ta11as. En cierto grado, todas las bodegas mostraron cambios cíclicos en ta11a media y varianza durante la descarga. En la mitad de las bodegas, se observo que selección por ta11a por los descargadores indujo un cambio en la ta11a media. Como resultado de la estratificación, la secuencia de ta11as sacadas de todas las bodegas no fue aleatoria, sin considerar si una bodega contenía peces de un solo lance 0 de mas de uno. El numero de ta11as modales en una bodega no estaba relacionado al numero de lances. En una bodega adicional compuesta de peces de varios lances, un experimento de mezcla vertical indico que una muestra representativa del contenido podría estar limitada a la mitad inferior de la bodega. Se uso el contenido de las bodegas de prueba para generar 25 bodegas simuladas y comparar los resultados de tres métodos de muestreo aplicados a estas. Los métodos fueron: (1) muestreo aleatorio (usado también como norma), (2) muestreo extendido, en el cual el proceso de selección fue extendido sobre una porción grande de una bodega, y (3) medición consecutiva de peces durante la descarga de la bodega. EI muestreo repetido con cada método y distintas combinaciones de n y m indico que, puesto que la fuente principal de variación de ta11a ocurría entre las unidades primarias, aumentar n fue la manera mas eficaz de reducir las estimaciones de la varianza de las ta11as de los grupos de edad y el numero total de peces en los desembarcos. El muestreo extendido evito mayormente los efectos de la estratificación por ta11a, y su desempeño fue esencialmente comparable a aquel del muestreo aleatorio. Se recomienda muestrear con este método. El muestreo de peces consecutivos produjo estimaciones mas sesgadas con mayores varianzas. Un análisis de las muestras de frecuencia de ta11a de 1988 indico que, para los grupos de edad que aparecen con mayor frecuencia en la captura, una frecuencia de muestreo minima de una unidad primaria de cada seis para cada estrato de mes-área reduciría los coeficientes de variación (CV) de las estimaciones de ta11a correspondientes a aproximadamente 10% 0 menos. Una estratificación adicional de las muestras por tipo de lance, y no solamente mes-área, redujo aun mas los CV de los grupos de edad escasos, tales como los reclutas, y mejoró potencialmente su precisión. Los CV de las estimaciones del reclutamiento para las cohortes completamente pescadas durante 1981-1984 fueron alrededor de 3-8%. Las estimaciones del reclutamiento y sus varianzas fueron también relativamente insensibles a cambios en las capturas de trimestres individuales y las varianzas, respectivamente, de las cuales fueron derivadas. (PDF contains 70 pages)

Force chains, friction, and flow: Behavior of granular media across length scales

We study the behavior of granular materials at three length scales. At the smallest length scale, the grain-scale, we study inter-particle forces and "force chains". Inter-particle forces are the natural building blocks of constitutive laws for granular materials. Force chains are a key signature of the heterogeneity of granular systems. Despite their fundamental importance for calibrating grain-scale numerical models and elucidating constitutive laws, inter-particle forces have not been fully quantified in natural granular materials. We present a numerical force inference technique for determining inter-particle forces from experimental data and apply the technique to two-dimensional and three-dimensional systems under quasi-static and dynamic load. These experiments validate the technique and provide insight into the quasi-static and dynamic behavior of granular materials.

At a larger length scale, the mesoscale, we study the emergent frictional behavior of a collection of grains. Properties of granular materials at this intermediate scale are crucial inputs for macro-scale continuum models. We derive friction laws for granular materials at the mesoscale by applying averaging techniques to grain-scale quantities. These laws portray the nature of steady-state frictional strength as a competition between steady-state dilation and grain-scale dissipation rates. The laws also directly link the rate of dilation to the non-steady-state frictional strength.

At the macro-scale, we investigate continuum modeling techniques capable of simulating the distinct solid-like, liquid-like, and gas-like behaviors exhibited by granular materials in a single computational domain. We propose a Smoothed Particle Hydrodynamics (SPH) approach for granular materials with a viscoplastic constitutive law. The constitutive law uses a rate-dependent and dilation-dependent friction law. We provide a theoretical basis for a dilation-dependent friction law using similar analysis to that performed at the mesoscale. We provide several qualitative and quantitative validations of the technique and discuss ongoing work aiming to couple the granular flow with gas and fluid flows.

Time-dependent monoenergetic neutron transport in two adjacent semi-infinite media

An exact solution to the monoenergetic Boltzmann equation is obtained for the case of a plane isotropic burst of neutrons introduced at the interface separating two adjacent, dissimilar, semi-infinite media. The method of solution used is to remove the time dependence by a Laplace transformation, solve the transformed equation by the normal mode expansion method, and then invert to recover the time dependence.

The general result is expressed as a sum of definite, multiple integrals, one of which contains the uncollided wave of neutrons originating at the source plane. It is possible to obtain a simplified form for the solution at the interface, and certain numerical calculations are made there.

The interface flux in two adjacent moderators is calculated and plotted as a function of time for several moderator materials. For each case it is found that the flux decay curve has an asymptotic slope given accurately by diffusion theory. Furthermore, the interface current is observed to change directions when the scattering and absorption cross sections of the two moderator materials are related in a certain manner. More specifically, the reflection process in two adjacent moderators appears to depend initially on the scattering properties and for long times on the absorption properties of the media.

This analysis contains both the single infinite and semi-infinite medium problems as special cases. The results in these two special cases provide a check on the accuracy of the general solution since they agree with solutions of these problems obtained by separate analyses.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

Spatiotemporal instabilities in nonlinear Kerr media in the presence of arbitrary higher-order dispersions

Spatiotemporal instabilities in nonlinear Kerr media with arbitrary higher-order dispersions are studied by use of standard linear-stability analysis. A generic expression for instability growth rate that unifies and expands on previous results for temporal, spatial, and spatiotemporal instabilities is obtained. It is shown that all odd-order dispersions contribute nothing to instability, whereas all even-order dispersions not only affect the conventional instability regions but may also lead to the appearance of new instability regions. The role of fourth-order dispersion in spatiotemporal instabilities is studied exemplificatively to demonstrate the generic results. Numerical simulations confirm the obtained analytic results. (C) 2002 Optical Society of America.