A Computational Fluid Dynamic (CFD) model study of wave and structure interactions on an elevated residential building under various air gap and surge/wave conditions was performed using the olaFlow, an open-source program using the OpenFOAM (Open-source Fields Operation And Manipulation) platform. The numerical model results, including free surface elevation, wave velocity, and vertical pressures on the underside of the elevated structure, showed a good agreement with the measured time-series data from the 1:6 scale hydraulic experiment (Duncan et al., 2021). The numerical simulations were used to extend the physical model tests by computing the vertical distribution of the pressure and resulting wave-induced horizontal forces/pressures, which were not measured in the physical model studies. The simulated results indicate that the pattern of pressure distributions at the frontal face of the elevated structure was controlled by water depth and wavebreaking types (nonbreaking, breaking, and broken waves). The wave induced-vertical force on the elevated structure strongly depends on wave height and the air gap, which is a net elevation from the still water level to the bottom of the structure, but the horizontal force shows complicated patterns due to the varied surge levels (flow depth), wave heights and air gaps. The new dimensionless parameter, α′/h, comprised of the air gap, incident wave height, and flow depth, is introduced and utilized to predict the horizontal forces on the elevated structure.
more »
« less
Smoothed particle hydrodynamics modeling of elevated structures impacted by tsunami-like waves
Accurate characterization of the response of coastal structures when subjected to tsunami-like waves is important for structural engineering assessment and design. The weakly-compressive Smoothed Particle Hydrodynamic (SPH) model can theoretically investigate such phenomena in both horizontal and vertical directions. Yet, the convergence of the solutions is sensitive to physical and numerical parameters used in the modeling. In this paper, multiple three- and two-dimensional SPH models are used to study the numerical convergence of free-surface elevation solutions for various initial inter-particle distances, domain locations along the flume and vicinity of the structure, and unbroken and broken wave flow conditions. The results are used to infer on the trade-offs between the accuracy of the SPH solutions and computational costs of the simulations, including computing time and data storage requirements. Two-dimensional models and an approximate ratio of ten particles per wave height can reasonably predict the nonturbulent unbroken wave case. The broken wave case requires three-dimensional models and four times the ratio of particles per wave height. A correlation between experimental and numerical results is then performed, showing adequacy of the free-surface elevation converged SPH models to capture global force responses. The distribution of horizontal and vertical pressures exerted on the elevated structure are characterized and compared with an analytical equation derived from the experimental dataset, highlighting the symbiotic relationship between experimental data, for calibration of the models, and numerical insights, for physical setup design. For example, additional instruments should be placed at strategic locations in future experimental programs to further validate numerical local responses, such as pressures near the edges and corners of structures. Such insights are important to support future work and development of updated US and European guidelines for the design of overland built infrastructures.
more »
« less
- Award ID(s):
- 2103713
- PAR ID:
- 10579323
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Engineering Structures
- Volume:
- 270
- Issue:
- C
- ISSN:
- 0141-0296
- Page Range / eLocation ID:
- 114851
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Seismic design of water retaining structures relies heavily on the response of the retained water to shaking. The water dynamic response has been evaluated by means of analytical, numerical, and experimental approaches. In practice, it is common to use simplified code‐based methods to evaluate the added demands imposed by water sloshing. Yet, such methods were developed with an inherent set of assumptions that might limit their application. Alternatively, numerical modeling methods offer a more accurate way of quantifying the water response and have been commonly validated using 1 g shake table experiments. In this study, a unique series of five centrifuge tests was conducted with the goal of investigating the hydrodynamic behavior of water by varying its height and length. Moreover, sine wave and earthquake motions were applied to examine the water response at different types and levels of excitation. Arbitrary Lagrangian‐Eulerian finite element models were then developed to reproduce 1 g shake table experiments available in the literature in addition to the centrifuge tests conducted in this study. The results of the numerical simulations as well as the simplified and analytical methods were compared to the experimental measurements, in terms of free surface elevation and hydrodynamic pressures, to evaluate their applicability and limitations. The comparison showed that the numerical models were able to reasonably capture the water response of all configurations both under earthquake and sine wave motions. The analytical solutions performed well except for cases with resonance under harmonic motions. As for the simplified methods, they provided acceptable results for the peak responses under earthquake motions. However, under sine wave motions, where convective sloshing is significant, they underpredict the response. Also, beyond peak ground accelerations of 0.5 g., a mild nonlinear increase in peak dynamic pressures was measured which deviates from assumed linear response in the simplified methods. The study confirmed the reliability of numerical models in capturing water dynamic responses, demonstrating their broad applicability for use in complex problems of fluid‐structure‐soil interaction.more » « less
-
Particulate matter in the environment, such as sediment, marine debris and plankton, is transported by surface waves. The transport of these inertial particles is different from that of fluid parcels described by Stokes drift. In this study, we consider the transport of negatively buoyant particles that settle in flow induced by surface waves as described by linear wave theory in arbitrary depth. We consider particles that fall under both a linear drag regime in the low Reynolds number limit and in a nonlinear drag regime in the transitional Reynolds number range. Based on an analysis of typical applications, we find that the nonlinear regime is the most widely applicable. From an expansion in the particle Stokes number, we find kinematic expressions for inertial particle motion in waves, and from a multiscale expansion in the dimensionless wave amplitude, we find expressions for the wave-averaged drift velocities. These drift velocities are analogous to Stokes drift and can be used in large-scale models that do not resolve surface waves. We find that the horizontal drift velocity is reduced relative to the Stokes drift of fluid parcels and that the vertical drift velocity is enhanced relative to the particle terminal settling velocity. We also demonstrate that a cloud of settling particles released simultaneously will disperse in the horizontal direction. Finally, we discuss the accuracy of our expressions by comparing against numerical simulations, which show excellent agreement, and against experimental data, which show the same trends.more » « less
-
Abstract In shallow coastal oceans, turbulent flows driven by surface winds and waves and constrained by a solid bottom disperse particles. This work examines the mechanisms driving horizontal and vertical dispersion of buoyant and sinking particles for times much greater than turbulent integral time scales. Turbulent fields are modeled using a wind‐stress driven large eddy simulation (LES), incorporating wave‐driven Langmuir turbulence, surface breaking wave turbulent kinetic energy inputs, and a solid bottom boundary. A Lagrangian stochastic model is paired to the LES to incorporate Lagrangian particle tracking. Within a subset of intermediate buoyant rise velocities, particles experience synergistic vertical mixing in which breaking waves (BW) inject particles into Langmuir downwelling velocities sufficient to drive deep mixing. Along‐wind dispersion is controlled by vertical shear in mean along‐wind velocities. Wind and bottom friction‐driven vertical shear enhances dispersion of buoyant and sinking particles, while energetic turbulent mixing, such as from BW, dampens shear dispersion. Strongly rising and sinking particles trapped at the ocean surface and bottom, respectively, experience no vertical shear, resulting in low rates of along‐wind dispersion. Crosswind dispersion is shaped by particle advection in wind‐aligned fields of counter‐rotating Langmuir and Couette roll cells. Langmuir cells enhance crosswind dispersion in neutrally to intermediately buoyant particles through enhanced cell hopping. Surface trapping restricts particles to Langmuir convergence regions, strongly inhibiting crosswind dispersion. In shallow coastal systems, particle dispersion depends heavily on particle buoyancy and wave‐dependent turbulent effects.more » « less
-
Abstract In many electrochemical processes, the transport of charged species is governed by the Nernst–Planck equation, which includes terms for both diffusion and electrochemical migration. In this work, a multi-physics, multi-species model based on the smoothed particle hydrodynamics (SPH) method is presented to model the Nernst–Planck equation in systems with electrodeposition. Electrodeposition occurs when ions are deposited onto an electrode. These deposits create complex boundary geometries, which can be challenging for numerical methods to resolve. SPH is a particularly effective numerical method for systems with moving and deforming boundaries due to its particle nature. This paper discusses the SPH implementation of the Nernst–Planck equations with electrodeposition and verifies the model with an analytical solution and a numerical integrator. A convergence study of migration and precipitation is presented to illustrate the model’s accuracy, along with comparisons of the deposition growth front to experimental results.more » « less
An official website of the United States government

