The Pen Branch Fault is a reactivated, high angle, reverse basin border fault that dips to southeast from the basement (Triassic-Paleozoic-Precambrian in age) to near vertical through the Atlantic Coastal Plain sediments (Late Cretaceous to Recent in age) of the Savannah River site in South Carolina. The fault movement has occurred through Late Cretaceous to Tertiary. Faulting might penetrate as shallow as the Dry Branch and the Tobacco Road Sand Formations (Late Eocene). An investigation with high resolution reflection seismic data is undertaken along the Pen Branch Fault. Five of the seismic lines are reprocessed to help in determining the upward depth of penetration of this fault. This is done by utilizing refracted arrivals from the multifold reflection data. The shallowest refractors (4 - 18 m) imaged have an average velocity of 1700 m/s. All of the lines exhibit events that are not flat lying across the data where the fault is believed to be. The lines possess deformation such as offsets, upwarping and channels. Deformation can be related to the Pen Branch Fault as shallow as 4 m from the surface. Displacements along the lines vary from 1 m up to 5 m. Reverse sense of motion is mainly exhibited along the fault zone that is covered by this study. The events resolved portray the Pen Branch Fault in a fault zone of subparallel faults and splays. / Master of Science

Overset, or Chimera, meshes are used to discretize the governing equations within a computational domain using multiple meshes that overlap in an arbitrary manner. The overset meshing technique is most applicable to problems dealing with multiple or moving bodies. Deep water simulations were carried out using both single and overset grid techniques for the evaluation of the overset grid application. These simulations were carried out using the commercial CFD code STAR-CCM+ by CD-adapco. The geometry simulated is that of a SES model (T-Craft) tested at the Naval Surface Warfare Center Carderock Division. The craft is simulated with two degrees of freedom, allowing movement in heave and pitch in response to displacement of the free surface. Agreement between the single and overset grid techniques was deemed reasonable to extend to future shallow water cases. However, due to longer run times of the overset mesh, the traditional or single mesh technique should be employed whenever applicable. In order to extend existing full craft CFD simulations of a surface effect ship (SES) into shallow water and maneuvering cases, an overset mesh is needed. Simulations of the SES were performed and monitored at various depth Froude numbers resulting in subcritical, critical, and supercritical flow regimes. Resistance, pitch response, and free surface response of the SES were compared between the shallow water simulations. The SES produced wider wakes, perpendicular to the craft, at simulations closer to the critical flow regime. Critical flow occurs at a depth Froude number between 0.9 and 0.95. Progression of shallow water effects through the three flow regimes agrees well with shallow water theory. / Master of Science

No / The orientation-sensitive performance of the Schlumberger array, when used to survey narrow, linear features, has long been recognized in geophysical prospecting for geology. However, in spite of frequent use of the array for archaeological survey, particularly in eastern Europe and the Far East, this directional effect is not apparent in the survey of walls and ditches. In order to examine the array's performance some experiments were carried out in a shallow electrolytic tank using insulating and conducting cylinders. Broadside and longitudinal traverses with systematic expansion of the current electrode spacing facilitated the production of pseudosections. The results confirmed the high selectivity of the Schlumberger response to the orientation of the feature. Broadside traverse of the conductor and longitudinal traverse of the insulator produced very large changes: much smaller signals were recorded for the alternative orientations. A subsequent experiment, however, on a simulated ditch in bedrock revealed no signal. The directional effect for a linear insulator was confirmed in field studies of a simple stone-walled structure. Implications for survey of low-contrast linear archaeological features are discussed.

No / Extant literature contains a number of specific case studies on differential decomposition involving adipocere formation or desiccation, but few describe the co-occurrence of these features within a temperate climate. The case of a 65-year-old male, partially buried in a shallow grave for 7 months, is presented in which the soft tissues of the body were outwardly well preserved. The right leg was desiccated, some parts of the body were covered with adipocere (head, neck, right shoulder, upper torso and left leg) and other parts could be classified as in the early stages of decomposition. In this study the taphonomic variables resulting in differential decomposition with desiccation and adipocere formation are discussed.

No

The geometrically nonlinear response of a shallow, circular, cylindrical panel under a midspan line load and induced strain actuation is presented. The panel is a laminate of piezoelectric material perfectly bonded to the convex and concave surfaces of a core of passive material. Since the curved edges are free and the straight edges are pinned a fixed distance apart, the response of the panel is independent of the axial coordinate. Hence, the governing ordinary differential equations are of the same form as for a shallow circular arch. Without induced strain actuation, the panel exhibits snap-through behavior under the midspan load. Induced strain distributions are determined at a constant midspan load to displace the panel to an inverted configuration in a stable manner. This adaptive structure may find application as an electromechanical, nonlinear spring with a digital-like, load-displacement response characteristic. / Master of Science

L’enjeu de la simulation de la dynamique atmosphérique et océanographique a pris ces dernières années une importance accrue avec la question du réchauffement climatique. Le modèle à simuler est complexe. Il combine les équations de la mécanique des fluides avec celles de la thermodynamique. Au 19ème siècle, le mathématicien Adhémar Barré de Saint-Venant formule un système d’équations aux dérivées partielles décrivant les mouvements d’un fluide soumis à la gravité et de faible épaisseur. Il s’agit des équations Shallow Water. L’objectif de cette thèse est de développer et d’analyser un algorithme de résolution des équations Shallow Water sur une sphère en rotation. Dans un premier temps, j’étudie différents aspects mathématiques des opérateurs aux différences finis utilisés par la suite en géométrie sphérique. Les schémas aux différences obtenus sont utilisés pour résoudre l’équation de transport, l’équation des ondes et l’équation de Burgers. Les propriétés de stabilité précision et conservation sont analysées. Dans un second temps, la grille Cubed-Sphere est introduite et analysée. La structure de ce maillage est analogue à celle d’un cube. L'interprétation de la Cubed-Sphere à l’aide de grands cercles permet de construire des opérateurs sphériques discrets gradient, divergence et vorticité d'ordre au moins égal à 3 (en pratique d'ordre 4). La troisième partie de la thèse est dédiée à différents tests pour le système d’équations Shallow Water ainsi que pour l’équation d’advection. Les résultats démontrent une précision proche de celle obtenue par les algorithmes conservatifs d'ordre 4 les plus récents / The problem to obtain accurate simulations of the atmospheric and oceanic equations has become essential in recent years for a proper understanding of the climate change. The full mathematical model to simulate is rather complex. It consists of the coupling of several equations involving fluid dynamics and thermodynamics. In the 19th century, Adhémar Barré de Saint-Venant first formulated the equations describing the dynamic of a fluid subject to gravity and bottom topography. This system is Shallow Water equations. The goal of this thesis is to develop and analyze a numerical scheme to solve the shallow water equation on a rotating sphere. First, a mathematical analsysis of finite difference operators that will be used on the sphere is presented. These schemes are then used to solve various equations in a spehreical setting, in particular the advection equation, the wave equation and the Burgers equation. Stability, accuracy and conservation properties are studied. In a second part, I consider in detail the Cubed-Sphere grid. This particular spherical grid has the mesh topology of a cube. Another interpretation makes use of great circles, this allows to obtain spherical discret operators gradient, divergence and curl of a preved third order. These operators are numercially of 4th order. Numerial results are show in particular for the SW equations an acurracy similar to the one of conservative schemes of 4th order published recently

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第21359号 / 工博第4518号 / 新制||工||1703(附属図書館) / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 岸田 潔, 教授 三村 衛, 教授 木村 亮 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

Flood inundation is a major natural hazard that can have very severe socio-economic consequences. This thesis presents an enhanced numerical model for flood simulation. After setting the context by examining recent large-scale flood events, a literature review is provided on shallow flow numerical models. A new version of the hyperbolic horizontal variable density shallow water equations with source terms in balanced form is used, designed for flows over complicated terrains, suitable for wetting and drying fronts and erodible bed problems. Bed morphodynamics are included in the model by solving a conservation of bed mass equation in conjunction with the variable density shallow water equations. The resulting numerical scheme is based on a Godunov-type finite volume HLLC approximate Riemann solver combined with MUSCL-Hancock time integration and a non-linear slope limiter and is shock-capturing. The model can simulate trans-critical, steep-fronted flows, connecting bodies of water at different elevations. The model is validated for constant density shallow flows using idealised benchmark tests, such as unidirectional and circular dam breaks, damped sloshing in a parabolic tank, dam break flow over a triangular obstacle, and dam break flow over three islands. The simulation results are in excellent agreement with available analytical solutions, alternative numerical predictions, and experimental data. The model is also validated for variable density shallow flows, and a parameter study is undertaken to examine the effects of different density ratios of two adjacent liquids and different hydraulic thrust ratios of species and liquid in mixed flows. The results confirm the ability of the model to simulate shallow water-sediment flows that are of horizontally variable density, while being intensely mixed in the vertical direction. Further validation is undertaken for certain erodible bed cases, including deposition and entrainment of dilute suspended sediment in a flat-bottomed tank with intense mixing, and the results compared against semi-analytical solutions derived by the author. To demonstrate the effectiveness of the model in simulating a complicated variable density shallow flow, the validated numerical model is used to simulate a partial dam-breach flow in an erodible channel. The calibrated model predictions are very similar to experimental data from tests carried out at Tsinghua University. It is believed that the present numerical solver could be useful at describing local horizontal density gradients in sediment laden and debris flows that characterise certain extreme flood events, where sediment deposition is important.

Un des mécanismes proposés pour expliquer l'origine de la variabilité climatique sur des périodes décennales à multidécennales est une oscillation propre de la circulation océanique thermohaline. Son mécanisme s'apparente aux modes de bassin basse fréquence et grande échelle qui résultent de l'interaction entre les ondes lentes planétaires et les ondes rapides de bord au cours du processus d'ajustement du bassin. Toutefois, la plupart des études de ce prototype oscillation décennale ont été menées dans des contextes simplifiés quasi-géostrophiques, à gravité réduite ou à fond plat. On se propose dans ce travail de thèse d'étudier l'effet de la topographie du fond et de la circulation moyenne sur les caractéristiques des modes de bassin baroclines. On utilise un modèle shallow water à deux couches verticales avec surface libre. Différentes bathymétries analytiques type fond plat, dorsale médio-océanique et pentes continentales sont étudiées.L'obtention des vecteurs propres du modèle linéarisé par analyse de stabilité linéaire autour d'un état au repos révèle que (1) la sélection de ces modes à basse résolution s'établit par la dissipation explicite introduite dans le modèle, (2) la période décennale et l'amortissement du mode le moins amorti sont faiblement sensibles à la topographie. Les budgets d'énergie et de vorticité de ces modes sont calculés dans le but de rationaliser le rôle amortisseur de la topographie via la conversion d'énergie qui a lieu entre les modes barotrope et barocline. En effet, une circulation barotrope, absente à fond plat, émerge à travers l'interaction entre le mode barocline à fond plat et la topographie. Toutefois, cette conversion d'énergie sous l'effet JEBAR demeure faible comparée aux processus visqueux.En présence d'une circulation stationnaire forcée par le vent et les flux de chaleur, les intégrations temporelles du modèle nonlinéaire perturbé par des structures baroclines cohérentes type tourbillons gaussiens montrent la forte interaction entre le vortex et la topographie. Cette interaction se manisfeste par : (1) une accélération de la vitesse de phase vers l'ouest par rapport au résultat à fond plat, (2) une circulation barotrope construite par la conversion de l'énergie barocline en barotrope, et (3) un déplacement méridien de l'anomalie dépendant de son signe même en l'absence d'advection nonlinéaire.Par ailleurs, le mode majeur de variabilité barocline, fortement amorti par la topographie et la dissipation dans la configuration non forcée paraît renforcé par l'écoulement stationnaire qui diminue son taux d'amortissement. Sa période d'oscillation développe une dépendance à la migration méridienne de l'advection zonale par l'écoulement moyen: elle est raccourcie (T ̴ 16 ans) pour le forçage par le vent et rallongée (T ̴ 22 ans) pour le forçage par les flux de chaleur. / One of the potential mechanisms at the origin of climatic variability on decadal to multidecadal timescales is the thermohaline oscillation corresponding to large-scale and low frequency basin modes that result from the interaction between long planetary waves and fast inertia-gravity waves during the adjustment process. However, most of the studies dealing with this decadal oscillation were carried out in a simplified flat bottom or reduced-gravity quasigeostrophic context.This dissertation aims to study the effect of bottom topography and mean flow on the characteristics of the gravest baroclinic basin modes in a mid-latitude idealized ocean basin. To that end, we make use of a two-layer shallow water (SW) model. Different bathymetries such as a flat bottom, a mid-ocean ridge and continental slopes are studied. Getting the eigenvectors from the linearized model through linear stability analysis around a state of rest reveals that (1) the selection of these modes is set by the explicit dissipation introduced in the model, (2) the oscillation period and decay rate are weakly sensitive to the form and height of the topography. Vorticity and energy budgets are computed in order to give a rationale for the decaying role of the topography via energy conversion from the baroclinic to the barotropic mode. Indeed, the barotropic flow absent in a flat bottom, results accurately from the interaction of the flat-bottomed baroclinic motion with the topographic height. However, the energy conversion under the JEBAR effect remains weaker with respect to the frictional processes.A stationary circulation is now included through wind or thermal forcing. Temporal integrations of the nonlinear model perturbed by coherent baroclinic structures with a gaussian eddy form show the strong interaction between the vortex and the topography. This interaction implies: (1) a westward acceleration of the zonal phase speed (with respect to the classic flat-bottom result), (2) a barotropic circulation built up by the conversion of the baroclinic energy into a barotropic one, and (3) an eddy sign-dependent meridional migration, even in the absence of nonlinear advection. Moreover, the decadal basin mode strongly damped by the topography and the dissipation shows a decrease of its decay rate by the large scale stationary forcing. Its oscillation period is found to be a strong function of the meridional migration of the eastward advection by the mean flow: it is shortened (T ̴ 16 yrs) in the wind-forced experiment and lengthened (T ̴ 22 yrs) with a thermal forcing.