Understanding Micro- and Macro-Mechanics of Soil Liquefaction: A Necessary Step for Field-Scale Assessment

Liquefaction is a devastating instability associated with saturated, loose, and cohesionless soils. It poses a significant risk to distributed infrastructure systems that are vital for the security, economy, safety, health, and welfare of societies. In order to make our cities resilient to the effects of liquefaction, it is important to be able to identify areas that are most susceptible. Some of the prevalent methodologies employed to identify susceptible areas include conventional slope stability analysis and the use of so-called liquefaction charts. However, these methodologies have some limitations, which motivate our research objectives. In this dissertation, we investigate the mechanics of origin of liquefaction in a laboratory test using grain-scale simulations, which helps (i) understand why certain soils liquefy under certain conditions, and (ii) identify a necessary precursor for onset of flow liquefaction. Furthermore, we investigate the mechanics of liquefaction charts using a continuum plasticity model; this can help in modeling the surface hazards of liquefaction following an earthquake. Finally, we also investigate the microscopic definition of soil shear wave velocity, a soil property that is used as an index to quantify liquefaction resistance of soil. We show that anisotropy in fabric, or grain arrangement can be correlated with anisotropy in shear wave velocity. This has the potential to quantify the effects of sample disturbance when a soil specimen is extracted from the field. In conclusion, by developing a more fundamental understanding of soil liquefaction, this dissertation takes necessary steps for a more physical assessment of liquefaction susceptibility at the field-scale.

Simulação de elevação do nível do mar na praia do Leme - RJ.

Os crescentes alertas sobre as mudanças climáticas e suas consequências vêm preocupando a sociedade de modo geral. Geram dúvidas sobre as reais modificações que podem ser ocasionadas, principalmente com relação aos efeitos da elevação do nível médio dos oceanos e seus efeitos nas regiões costeiras. Por vezes, as informações veiculadas contribuem para uma percepção limitada da extensão espaço-temporal dos fenômenos. Em paralelo, a democratização de mapas 2D e 3D está cada vez mais ampla, não sendo mais uma ferramenta direcionada apenas aos profissionais ligados à área da Cartografia e demais Geociências. Este é um dos enormes benefícios proporcionados pela Cartografia Digital. Porém, não teria tanto alcance sem a associação das informações a um poderoso veículo de divulgação como a Web. A Internet é um dos meios de comunicação mais democráticos do mundo e uma ferramenta de grande alcance na sociedade. O uso dessas ferramentas pode ser ainda mais explorado para a construção do conhecimento, permitindo aos cidadãos entenderem mais facilmente um determinado fenômeno por meio da sua visualização/simulação em tempo real. Isto melhora a formação de opinião e o posicionamento da sociedade sobre o caso em questão. Com base nessa argumentação, esta dissertação tem o objetivo de simular cenários de elevações do nível médio do mar e disponibilizá-los em formato de mapas interativos na Web utilizando um Modelo Digital de Elevação suficientemente acurado. Para isto, a área de estudo escolhida foi a Praia do Leme, situada no município do Rio de Janeiro, RJ. Os processos para atingir os objetivos envolveram etapas de coleta de dados de campo, confecção de um banco de dados integrado a dados pré-existentes, interpolações, simulação de cenários da elevação do nível do mar e implementação de uma página Web. Os cenários investigados foram obtidos por simulações de elevações do nível do mar a cada 0,5m até 4,0m acima da cota atual. Os resultados obtidos mostram que a Praia do Leme sofrerá modificações geomorfológicas, com perda de até 18,0m da praia, para uma elevação do nível médio do mar de 2,0m. Entretanto, a orla litorânea sofrerá impactos significativos a partir de uma elevação do nível do mar de 3,0m, em que o avanço planimétrico do mar pode ficar entre 76 e 90m, o que atingiria o calçadão e a região do emissário de águas pluviais do Leme, respectivamente. Isto resultaria em um volume adicional de água na praia de, aproximadamente, 106.000m3 e uma perda de até 41.400m3 de areia. Somente a partir de 3,5m de elevação, o mar atingiria a Avenida Atlântica, sendo que, desta altitude em diante, todo o bairro do Leme seria gradualmente inundado pelo avanço do mar.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.

Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Three-dimensional direct numerical simulations of autoignition in turbulent non-premixed flows with simple and complex chemistry

3D Direct Numerical Simulations (DNS) of autoignition in turbulent non-premixed flows between fuel and hotter air have been carried out using both 1-step and complex chemistry consisting of a 22 species n-heptane mechanism to investigate spontaneous ignition timing and location. The simple chemistry results showed that the previous findings from 2D DNS that ignition occurred at the most reactive mixture fraction (ξMR) and at small values of the conditional scalar dissipation rate (N|ξMR) are valid also for 3D turbulent mixing fields. Performing the same simulation many times with different realizations of the initial velocity field resulted in a very narrow statistical distribution of ignition delay time, consistent with a previous conjecture that the first appearance of ignition is correlated with the low-N content of the conditional probability density function of N. The simulations with complex chemistry for conditions outside the Negative Temperature Coefficient (NTC) regime show behaviour similar to the single-step chemistry simulations. However, in the NTC regime, the most reactive mixture fraction is very rich and ignition seems to occur at high values of scalar dissipation. Copyright © 2006 by ASME.

Towards jet flow LES of conceptual nozzles for acoustic predictions

Progress in simulating chevron nozzle jet flows using ILES/RANS-ILES approaches and using the Ffowcs Williams and Hawkings (FW-H) surface integral method to predict the radiated far field sound is presented in this paper. With the focus on the realistic chevron geometries, SMC001 and SMC006, coarse and fine meshes are generated in the range of 3∼13 million mesh cells. Throughout this work, to minimize numerical dissipation introduced by mesh quality issues, the hexahedral cell type is used. Numerical simulations are then carried out with cell-vertex and cell-centered codes. Despite the modest grids, mean velocities and turbulent statistics are found to be in reasonable accord with measurements. Also, far field sound levels predicted by the FW-H post processor are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

High speed electro-thermal models for inverter simulations

This paper presents the application of advanced compact models of the IGBT and PIN diode to the full electrothermal system simulation of a hybrid electric vehicle converter using a look-up table of device losses. The Fourier-based solution model is used, which takes account of features such as local lifetime control and field-stop technology. Device and circuit parameters are extracted from experimental waveforms and device structural data. Matching of the switching waveforms and the resulting generation of the look-up table is presented. An example of the use of the look-up tables in simulation of inverter device temperatures is also given, for a hypothetical electric vehicle subjected to an urban driving cycle. © 2006 IEEE.

Hot electron field emission via individually transistor-ballasted carbon nanotube arrays.

We present electronically controlled field emission characteristics of arrays of individually ballasted carbon nanotubes synthesized by plasma-enhanced chemical vapor deposition on silicon-on-insulator substrates. By adjusting the source-drain potential we have demonstrated the ability to controllable limit the emission current density by more than 1 order of magnitude. Dynamic control over both the turn-on electric field and field enhancement factor have been noted. A hot electron model is presented. The ballasted nanotubes are populated with hot electrons due to the highly crystalline Si channel and the high local electric field at the nanotube base. This positively shifts the Fermi level and results in a broad energy distribution about this mean, compared to the narrow spread, lower energy thermalized electron population in standard metallic emitters. The proposed vertically aligned carbon nanotube field-emitting electron source offers a viable platform for X-ray emitters and displays applications that require accurate and highly stable control over the emission characteristics.

Response probability distribution of built-up vibro-acoustic systems

The vibro-acoustic response of built-up structures, consisting of stiff components with low modal density and flexible components with high modal density, is sensitive to small imperfections in the flexible components. In this paper, the uncertainty of the response is considered by modeling the low modal density master system as deterministic and the high modal density subsystems in a nonparametric stochastic way, i.e., carrying a diffuse wave field, and by subsequently computing the response probability density function. The master system's mean squared response amplitude follows a singular noncentral complex Wishart distribution conditional on the subsystem energies. For a single degree of freedom, this is equivalent to a chi-square or an exponential distribution, depending on the loading conditions. The subsystem energies follow approximately a chi-square distribution when their relative variance is smaller than unity. The results are validated by application to plate structures, and good agreement with Monte Carlo simulations is found. © 2012 Acoustical Society of America.

Angular distribution of field emitted electrons from vertically aligned carbon nanotube arrays

Angular field emission (FE) properties of vertically aligned carbon nanotube arrays have been measured on samples grown by plasma enhanced chemical vapor deposition and characterized by scanning electron microscope and I-V measurements. These properties determine the angular divergence of electron beams, a crucial parameter in order to obtain high brilliance FE based cathodes. From angular distributions of the electron beam transmitted through extraction grids of different mesh size and by using ray-tracing simulations, the maximum emission angle from carbon nanotube tips has been determined to be about ± 30 around the tube main axis. © 2012 American Institute of Physics.

Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Ballasted & electrically steerable carbon nanotube field emitters

Here we present our on-going efforts toward the development of stable ballasted carbon nanotube-based field emitters employing hydrothermally synthesized zinc oxide nanowires and thin film silicon-on-insulator substrates. The semiconducting channel in each controllably limits the emission current thereby preventing detrimental burn-out of individual emitters that occurs due to unavoidable statistical variability in emitter characteristics, particularly in their length. Fabrication details and emitter characterization are discussed in addition to their field emission performance. The development of a beam steerable triode electron emitter formed from hexagonal carbon nanotube arrays with central focusing nanotube electrodes, is also described. Numerical ab-initio simulations are presented to account for the empirical emission characteristics. Our engineered ballasted emitters have shown some of the lowest reported lifetime variations (< 0.7%) with on-times of < 1 ms, making them ideally-suited for next-generation displays, environmental lighting and portable x-rays sources. © 2012 SPIE.

Numerical simulation of single-stream jets from a serrated nozzle

Hybrid large-eddy type simulations for cold jet flows from a serrated nozzle are performed at an acoustic Mach number Ma ac = 0.9 and Re = 1.03×10 6. Since the solver being used tends towards having dissipative qualities, the subgrid scale (SGS) model is omitted, giving a numerical type LES (NLES) or implicit LES (ILES) reminiscent procedure. To overcome near wall streak resolution problems a near wall RANS (Reynolds averaged Navier-Stokes) model is smoothly blended to the LES making a hybrid RANS-ILES. The geometric complexity of the serrated nozzle is fully considered without simplification or emulation. An improved but still modest hexahedral multi-block grid with circa 20 million grid points (with respect to 12.5 million in Xia et al.; Int J Heat Fluid Flow 30:1067-1079, 2009) is used. Despite the modest grid size, encouraging and improved results are obtained. Directly resolved mean and second-order fluctuating quantities along the jet centerline and in the jet shear layer compare favorably with measurements. The radiated far-field sound predicted using the Ffowcs Williams and Hawkings (FW-H) surface integral method shows good agreement with the measurements in directivity and sound spectra. © 2011 Springer Science+Business Media B.V.

Large-eddy simulations and measurements of a small-scale high-speed coflowing jet

Measurements and predictions are made of a short-cowl coflowing jet with a bypass ratio of 8:1. The Reynolds number is 300,000, and the inlet Mach numbers are representative of aeroengine conditions. The low Reynolds number of the measurements makes the case well suited to the assessment of large-eddy-simulation-related strategies. The nozzle concentricity is carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both laser Doppler anemometry and particle image velocimetry. The simulations are completed on 6× 106, 12× 106, and 50 × 106 cell meshes. To overcome near-wall modeling problems, a hybrid large-eddy-simulation-Reynolds-averaged-Navier-Stokesrelated method is used. The near-wall Reynolds-averaged-Navier-Stokes layer is helpful in preventing nonphysical separation from the nozzle wall.Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Linear and nonlinear model large-eddy simulations of a plane jet

The nonlinear Kosovic, and mixed Leray and α subgrid scale models are contrasted with linear Smagorinsky and Yoshizawa Large Eddy Simulations for a Re = 4000 plane jet simulation. Comparisons are made with Direct Numerical Simulation data and measurements. Global properties of the jet such as the spreading and centreline velocity decay rates are investigated. The mean-flow and turbulence parameters in the self-similar region are also studied. All models generally give encouraging agreement with the Direct Numerical Simulation data and reliable measurements. Solution differences for the models are relatively minor, none giving clear improvements for all data comparisons.