Estimation Of Shear Stress Research Articles

Modeling of Stress-Strain State and Coseismic Effects of Epicentral Zone of Tangshan Earthquake (Southeastern China)

The paper presents the results of numerical modeling and analysis of stress-strain state of the epicentral zone of the strong earthquake in the north-east of China, which occurred on 27.07.1976 with Ms=7.8. Many present-day works continue to discuss the reasons for such a strong earthquake, which occurred in tectonic conditions ‒ far from interplate boundaries, inside the Tangshan tectonic block bounded by tectonic faults. However, published new geodynamic, seismological, geophysical and geodetic data provide confidence in the determining role of fault tectonics in this region. Based on the analysis of the results of modeling of the stress-strain state preceding the Tangshan earthquake with coseismic geophysical and geodetic data, we propose a model of earthquake rupture formation. The results of comparison of independent estimates of shear stresses with the results of modeling in the sources of strong earthquakes suggest that the areas of tectonic stress concentration are localized in the interfault rupture of the Tangshan fault, reaching maximum values at the termination of fault ruptures σi ≈ 50 MPa и τxy ≈ 20 MPa. The hypocenter of the main seismic event (taking into account the error of coordinate determination) is located in the region of stress intensity 35‒50 MPa and the ratio of main stresses σxx/σyy ≈ 8–10. It should be expected that these zones are the starting site of rupture, the extent of which depends on the amount of accumulated elastic potential energy of tectonic stresses in the adjacent region. For the Tangshan earthquake, this area corresponds to a high intensity of stresses exceeding 30 MPa in a band with a length more than 30 km and a width reaching 4.5 km.

Comparative analysis of Zero Pressure Geometry and prestress methods in cardiovascular Fluid-Structure Interaction

Background and Objective:Modelling patient-specific aortic biomechanics with advanced computational techniques, such as Fluid–Structure Interaction (FSI), can be crucial to provide effective decision-making indices to enhance current clinical practices. To effectively simulate Ascending Thoracic Aortic Aneurysms (ATAA), the stress-free configuration must be defined. The Zero Pressure Geometry (ZPG) and the Prestress Tensor (PT) are two of the main approaches to tackle this issue. However, their impact on the numerical results is yet to be analysed. Computed Tomography Angiography (CTA) and Magnetic Resonance Imaging (MRI) data were used to develop patient-specific 2-way FSI frameworks. Methods:Three models were developed considering different tissue prestressing approaches to account for the reference configuration and their numerical results were compared. The selected approaches were: (i) ZPG, (ii) PT and (iii) a combination of the PT approach with a regional mapping of material properties (PTCAL). Results:The pressure fields estimated by all models were equivalent. The estimation of Wall Shear Stress (WSS) based metrics revealed good correspondence between all models except the Relative Residence Time (RRT). Regarding ATAA wall mechanics, the proposed extension to the PT approach presented a closer agreement with the ZPG model than its counterpart. Additionally, the PT and PTCAL approaches required around 60% fewer iterations to achieve cycle-to-cycle convergence than the ZPG algorithm. Conclusion:Using a regional mapping of material properties in combination with the PT method presented a better correspondence with the ZPG approach. The outcomes of this study can pave the way for advancing the accuracy and convergence of ATAA numerical models using the PT methodology.

Coupled CFD-DEM modeling of surface erosion in granular soils: Simulation of erosion function apparatus experiments

This study introduces a numerical modeling approach that couples the computational fluid dynamics (CFD) with the discrete element method (DEM) to simulate grain-scale soil erosion processes induced by water flows. In this modeling framework, CFD simulates fluid flows by solving the volume-averaged Navier-Stokes equations, and uses the k-ω turbulent model for turbulent flows. Simultaneously, DEM computes the displacement of solid particles by incorporating the fluid-particle interactions driven by fluid flows while adhering to Newton's laws of motion. These interactions encompass drag force, buoyancy force, pressure-gradient force, and viscous force exerted by fluid flows and acting on the particles. The coupled CFD-DEM modeling adeptly replicates soil erosion processes, demonstrating good alignment with results obtained from laboratory erosion function apparatus (EFA) tests. In particular, the DEM facilitates the estimation of shear stress acting on the soil surface based on fluid-particle interaction forces, which has been roughly approximated by empirical or semi-empirical models. This study underscores the capability of coupled CFD-DEM in providing valuable insights into the grain-scale behavior of soil particles subjected to fluid flows, with the potential for extension to address soil erosion and fines migration.

A Numerical Approach to Analyzing Shallow Flows over Rough Surfaces

The hydraulic characteristics (such as velocity profiles, near-bed velocity profile, bed shear stress, and resistance coefficients) of shallow flows over rough surfaces were investigated using numerical simulations. A novel method is presented to simulate shallow flows over rough surfaces in a two-dimensional (2D) numerical domain, where the physical numerical domain represents bed topography. Results reveal that the model can accurately predict spatially averaged velocity profiles, turbulence characteristics, shear stresses, and uniform flow depths. The analysis identified two distinct flow regions based on mean and turbulent flow profiles. Results show that the turbulent shear stress profiles provide a more accurate estimation of the bed shear stresses. Resistance coefficients (friction factor or Manning’s roughness coefficient) vary with Froude number and submergence ratio (depth divided by roughness height).

Estimation of friction factor and bed shear stress considering bedform effect in rivers

AbstractThe roughness height plays a crucial role, especially in shallow‐water environments with rough‐bed conditions. Specifically, it is a key parameter for predicting flow velocity and bed shear stress. Therefore, an accurate determination of roughness height is essential for precise predictions of hydraulic phenomena. In this study, we propose a method to estimate form roughness height and friction factor using bathymetry data, with a focus on cross‐sectional data in the transverse direction. To overcome the limitations of using lateral direction data, we derived an empirical formula for estimating bedform length from the bed profiles in the streamwise direction. To assess the validity of bedform analysis and estimated form roughness height based on lateral direction data, we compared the characteristics of bedform and form roughness height analysed using lateral direction data with those analysed using the streamwise direction data. Furthermore, we confirmed the importance of considering the form roughness height in estimating bed shear stress through a comparison with bed shear stress calculations considering only grain roughness height.

A Pseudo-Spectral Method for Wall Shear Stress Estimation from Doppler Ultrasound Imaging in Coronary Arteries.

The Wall Shear Stress (WSS) is the component tangential to the boundary of the normal stress tensor in an incompressible fluid, and it has been recognized as a quantity of primary importance in predicting possible adverse events in cardiovascular diseases, in general, and in coronary diseases, in particular. The quantification of the WSS in patient-specific settings can be achieved by performing a Computational Fluid Dynamics (CFD) analysis based on patient geometry, or it can be retrieved by a numerical approximation based on blood flow velocity data, e.g., ultrasound (US) Doppler measurements. This paper presents a novel method for WSS quantification from 2D vector Doppler measurements. Images were obtained through unfocused plane waves and transverse oscillation to acquire both in-plane velocity components. These velocity components were processed using pseudo-spectral differentiation techniques based on Fourier approximations of the derivatives to compute the WSS. Our Pseudo-Spectral Method (PSM) is tested in two vessel phantoms, straight and stenotic, where a steady flow of 15mL/min is applied. The method is successfully validated against CFD simulations and compared against current techniques based on the assumption of a parabolic velocity profile. The PSM accurately detected Wall Shear Stress (WSS) variations in geometries differing from straight cylinders, and is less sensitive to measurement noise. In particular, when using synthetic data (noise free, e.g., generated by CFD) on cylindrical geometries, the Poiseuille-based methods and PSM have comparable accuracy; on the contrary, when using the data retrieved from US measures, the average error of the WSS obtained with the PSM turned out to be 3 to 9 times smaller than that obtained by state-of-the-art methods. The pseudo-spectral approach allows controlling the approximation errors in the presence of noisy data. This gives a more accurate alternative to the present standard and a less computationally expensive choice compared to CFD, which also requires high-quality data to reconstruct the vessel geometry.

New momentum integral equation applicable to boundary layer flows under arbitrary pressure gradients

By incorporating the traditionally overlooked advective term in the wall-normal momentum equation, a new momentum integral equation is developed for two-dimensional incompressible turbulent boundary layers under arbitrary pressure gradients. The classical Kármán's integral arises as a special instance of the new momentum integral equation when the pressure gradient is weak. The new momentum integral equation's validity is substantiated by direct numerical simulation data. Unlike the classical Kármán's integral, which is limited to predicting wall shear stress within mild pressure gradients, the new momentum integral equation accurately computes wall shear stress across a broad range of pressure gradients, even in the presence of strong adverse pressure gradients that lead to flow separation. Moreover, a new pressure parameter $\beta _\kappa$ is introduced through examining terms in the new momentum integral equation. This parameter naturally quantifies the pressure gradient's influence on turbulent boundary layers and offers guidance for applying the classical Kármán's integral. Additionally, to facilitate experimental determination of wall shear stress under strong pressure gradients, an approximate integral equation is proposed that relies solely on easily measurable variables. Validation against direct numerical simulation data demonstrates that this simplified equation provides reasonably accurate estimates of wall shear stress in turbulent boundary layers experiencing strong pressure gradients.

The Stress of Measuring Plantar Tissue Stress in People with Diabetes-Related Foot Ulcers: Biomechanical and Feasibility Findings from Two Prospective Cohort Studies.

Reducing high mechanical stress is imperative to heal diabetes-related foot ulcers. We explored the association of cumulative plantar tissue stress (CPTS) and plantar foot ulcer healing, and the feasibility of measuring CPTS, in two prospective cohort studies (Australia (AU) and The Netherlands (NL)). Both studies used multiple sensors to measure factors to determine CPTS: plantar pressures, weight-bearing activities, and adherence to offloading treatments, with thermal stress response also measured to estimate shear stress in the AU-study. The primary outcome was ulcer healing at 12 weeks. Twenty-five participants were recruited: 13 in the AU-study and 12 in the NL-study. CPTS data were complete for five participants (38%) at baseline and one (8%) during follow-up in the AU-study, and one (8%) at baseline and zero (0%) during follow-up in the NL-study. Reasons for low completion at baseline were technical issues (AU-study: 31%, NL-study: 50%), non-adherent participants (15% and 8%) or combinations (15% and 33%); and at follow-up refusal of participants (62% and 25%). These underpowered findings showed that CPTS was non-significantly lower in people who healed compared with non-healed people (457 [117; 727], 679 [312; 1327] MPa·s/day). Current feasibility of CPTS seems low, given technical challenges and non-adherence, which may reflect the burden of treating diabetes-related foot ulcers.

Pro+: Automated protrusion and critical shear stress estimates from 3D point clouds of gravel beds

AbstractThe dimensionless critical shear stress (τ*c) needed for the onset of sediment motion is important for a range of studies from river restoration projects to landscape evolution calculations. Many studies simply assume a τ*c value within the large range of scatter observed in gravel‐bedded rivers because direct field estimates are difficult to obtain. Informed choices of reach‐scale τ*c values could instead be obtained from force balance calculations that include particle‐scale bed structure and flow conditions. Particle‐scale bed structure is also difficult to measure, precluding wide adoption of such force‐balance τ*c values. Recent studies have demonstrated that bed grain size distributions (GSD) can be determined from detailed point clouds (e.g. using G3Point open‐source software). We build on these point cloud methods to introduce Pro+, software that estimates particle‐scale protrusion distributions and τ*c for each grain size and for the entire bed using a force‐balance model. We validated G3Point and Pro+ using two laboratory flume experiments with different grain size distributions and bed topographies. Commonly used definitions of protrusion may not produce representative τ*c distributions, and Pro+ includes new protrusion definitions to better include flow and bed structure influences on particle mobility. The combined G3Point/Pro+ provided accurate grain size, protrusion and τ*c distributions with simple GSD calibration. The largest source of error in protrusion and τ*c distributions were from incorrect grain boundaries and grain locations in G3Point, and calibration of grain software beyond comparing GSD is likely needed. Pro+ can be coupled with grain identifying software and relatively easily obtainable data to provide informed estimates of τ*c. These could replace arbitrary choices of τ*c and potentially improve channel stability and sediment transport estimates.

Surface Depassivation via B-O Dative Bonds Affects the Friction Performance of B-Doped Carbon Coatings.

Boron doping of diamond-like carbon coatings has multiple effects on their tribological properties. While boron typically reduces wear in cutting applications, some B-doped coatings show poor tribological performance compared with undoped films. This is the case of the tribological tests presented in this work in which an alumina ball is placed in frictional contact with different undoped and B-doped amorphous carbon coatings in humid air. With B-doped coatings, a higher friction coefficient at a steady state with respect to their undoped counterparts was observed. Estimates of the average contact shear stress based on experimental friction coefficients, surface topographies, and Persson's contact theory suggest that the increased friction is compatible with the formation of a sparse network of interfacial ether bonds leading to a mild cold-welding friction regime, as documented in the literature. Tight binding and density functional theory simulations were performed to investigate the chemical effect of B-doping on the interfacial properties of the carbon coatings. The results reveal that OH groups that normally passivate carbon surfaces in humid environments can be activated by boron and form B-O dative bonds across the tribological interfaces, leading to a mild cold-welding friction regime. Simulations performed on different tribological pairs suggest that this mechanism could be valid for B-doped carbon surfaces in contact with a variety of materials. In general, this study highlights the impact that subtle modifications in surface and interface chemistry caused by the presence of impurities can have on macroscopic properties, such as friction and wear.

Smoothed particle hydrodynamics based FSI simulation of the native and mechanical heart valves in a patient-specific aortic model

The failure of the aortic heart valve is common, resulting in deterioration of the pumping function of the heart. For the end stage valve failure, bi-leaflet mechanical valve (most popular artificial valve) is implanted. However, due to its non-physiological behaviour, a significant alteration is observed in the normal haemodynamics of the aorta. While in-vivo experimentation of a human heart valve (native and artificial) is a formidable task, in-silico study using computational fluid dynamics (CFD) with fluid structure interaction (FSI) is an effective and economic tool for investigating the haemodynamics of natural and artificial heart valves. In the present work, a haemodynamic model of a natural and mechanical heart valve has been developed using meshless particle-based smoothed particle hydrodynamics (SPH). In order to further enhance its clinical relevance, this study employs a patient-specific vascular geometry and presents a successful validation against traditional finite volume method and 4D magnetic resonance imaging (MRI) data. The results have demonstrated that SPH is ideally suited to simulate the heart valve function due to its Lagrangian description of motion, which is a favourable feature for FSI. In addition, a novel methodology for the estimation of the wall shear stress (WSS) and other related haemodynamic parameters have been proposed from the SPH perspective. Finally, a detailed comparison of the haemodynamic parameters has been carried out for both native and mechanical aortic valve, with a particular emphasis on the clinical risks associated with the mechanical valve.

Subduction plate interface shear stress associated with rapid subduction at deep slow earthquake depths: example from the Sanbagawa belt, southwestern Japan

Abstract. Maximum shear stress along an active deformation zone marking the subduction plate interface is important for understanding earthquake phenomena and is an important input parameter in subduction zone thermomechanical modeling. However, such maximum shear stress is difficult to measure directly at depths more than a few kilometers and is generally estimated by simulation using a range of input parameters with large associated uncertainties. In addition, estimated values generally represent maximum shear stress conditions over short observation timescales, which may not be directly applicable to long-timescale subduction zone modeling. Rocks originally located deep in subduction zones can record information about deformation processes, including maximum shear stress conditions, occurring in regions that cannot be directly accessed. The estimated maximum shear stress is likely to be representative of maximum shear stress experienced over geological timescales and be suitable to use in subduction zone modeling over timescales of millions to tens of millions of years. In this study, we estimated maximum shear stress along a subduction plate interface by using samples from the Sanbagawa metamorphic belt of southwestern (SW) Japan, in which slivers of mantle-wedge-derived serpentinite are widely distributed and in direct contact with metasedimentary rocks derived from the subducted oceanic plate. These areas can be related to the zone of active deformation along the subduction plate interface. To obtain estimates of maximum shear stress at the subduction interface, we focused on the microstructure of quartz-rich metamorphic rocks – quartz is the main component of the rocks we collected and its deformation stress is assumed to be roughly representative of the stress experienced by the surrounding rock and plate interface deformation zone. Maximum shear stress was calculated by applying deformation temperatures estimated by the crystallographic orientation of quartz (the quartz c-axis fabric opening-angle thermometer) and the apparent grain size of dynamically recrystallized quartz in a thin section to an appropriate piezometer. Combined with information on sample deformation depth, estimated from the P–T (pressure–temperature) path and deformation temperatures, it is suggested that there was nearly constant maximum shear stress of 15–41 MPa in the depth range of about 15–30 km, assuming plane stress conditions even when uncertainties related to the measurement direction of thin section and piezometer differences are included. The Sanbagawa belt formed in a warm subduction zone. Deep slow earthquakes are commonly observed in modern-day warm subduction zones such as SW Japan, which has a similar thermal structure to the Sanbagawa belt. In addition, deep slow earthquakes are commonly observed to be concentrated in a domain under the shallow part of the mantle wedge. Samples showed the depth conditions near the mantle wedge, suggesting that these samples were formed in a region with features similar to the deep slow earthquake domain. Estimated maximum shear stress may not only be useful for long-timescale subduction zone modeling but also represent the initial conditions from which slow earthquakes in the same domain nucleated.

Experimental Evaluation of Estimating Local Bed Shear Stress Using Log Law in a Straight River Channel

Bed shear stress is a critical parameter in riverine environments, directly influencing sediment transport, erosion, deposition, and flow dynamics. However, direct measurement of bed shear stress in natural rivers is nearly impossible, and estimating bed shear stress by extrapolating the Reynolds stress profile, though known for accuracy, faces practical limitations due to the challenges of measuring turbulence in natural rivers. As alternatives, indirect methods, such as global bed shear stress estimation based on bed slope and hydraulic radius (assuming uniform flow) and local bed shear stress estimation through regression of velocity profiles using the log law, have been employed. The estimation of local bed shear stress using the log law has primarily been derived from laboratory flow conditions due to the limitations in observing detailed turbulence in field conditions. Consequently, the applicability of the log-law method under field conditions has not been sufficiently evaluated. This study aims to assess local bed shear stress using the log law through a dedicated field experiment in a specific river condition, comparing it with methods utilizing the Reynolds stress profile. The experiment was conducted in a straight river channel with a width of approximately 6.5 m, a depth of 0.61 m, and a mean flow velocity of 0.56 m/s. A total of 197 points of regularly distributed velocity and turbulence data for the given cross-section were acquired using micro-acoustic Doppler velocimeter measurements with a duration of 90 s. We found that the log-law method overestimated local bed shear stress by an overall factor of 2.3 compared to the method using Reynolds stress profiles. Thus, the experimental results conclusively indicate that local bed shear stress estimated from the log law in natural river conditions could be overestimated. However, additional research is needed under more diverse flow and geometric conditions to develop a correction formula for applying the log law in estimating local bed shear stress.

One tune, many tempos: Faults trade off slip in time and space to accommodate relative plate motions

Analysis of incremental slip rates from the four major strike-slip faults of the Marlborough fault system (MFS) of northern South Island, New Zealand, provides a first-ever record at the scale of an entire plate-boundary fault system of how relative plate motions are accommodated in time and space. This record, which spans the past 350–450 m of relative plate motion and ca. 12–14 ky, demonstrates that the fault system as a whole accommodates a steady plate-boundary slip rate, with the MFS faults “keeping up” with the overall rate of relative Pacific-Australia plate motion at relatively short displacement (10 s of meters) and time (102–103 yr) scales. These results affirm the often-assumed but until now unproven assumption that the relative plate-motion rate provides a robust basic constraint on both geodynamical models and analyses of system-level seismic hazard at these scales. In marked contrast, the incremental slip rates of each of the four main Marlborough faults are highly variable through time, marked by coordinated accelerations and decelerations spanning 4–6 earthquakes and several millennia as the faults trade off slip to accommodate a steady relative plate motion rate. These results suggest that (a) the weakest fault in the system will slip faster than average while adjacent mechanically complementary faults slip more slowly, and (b) that these patterns switch back and forth through time, likely reflecting reversible changes in the strength (i.e., resistance to shear) of the individual faults as they collectively accommodate relative plate motion. Interestingly, the periods of fast slip on the MFS faults exhibit ∼20–25 m of displacement, suggesting that these may record periods of fast slip on a weakened fault/ductile shear zone that continues until it uses up all locally stored elastic strain energy, thus potentially approaching local complete stress drop, albeit during a few tens of meters of rapid fault slip during multiple earthquakes, rather than during a single event. This hypothesis is consistent with typical earthquake stress drops of ∼1–10 MPa and estimates of depth-averaged crustal shear stress of a few 10 s of MPa, such as might be released in clusters of 4–6 earthquakes. These results emphasize the need to analyze the collective behavior of the entire fault systems, rather than just individual faults, to understand the mechanics of the system. Moreover, these patterns suggest a potential path forward for more accurate estimation of time-dependent seismic hazard, with the possible incorporation of current position of a fault within a fast- or slow/no-slip period into the probability analysis, as well as a means of potentially estimating crustal shear stress.

Subaqueous silt ripples measured by an echo sounder: Implications for bed roughness, bed shear stress and erosion threshold

Bedforms like ripples are widely distributed in the coastal zone. They influence the bed roughness thus the estimation of bed shear stress and associated sediment transport. An echo sounder was mounted on a bottom-supported tripod intended to measure the erosion and deposition of seabed in the subaqueous Yellow River Delta, China. However, variations in bed elevation are found to be not the net erosion or deposition at the observation site, but the migration of silt ripples which were generated by waves and pushed back and forth or flattened off by the tidal currents. Ripple heights were observed to be within 0.1–0.7 cm and were used for testing the model of Nielsen (1981). The model overestimated the height of silt ripples as it was developed for sands, but the deviation can be well addressed by incorporating a linear modification (R2 = 0.79). Alternatively, a new model specifically for silt ripple height was regressed from the field data with R2 up to 0.78. The existence of silt ripples increases the bed shear stress by 4 times due to the additional bed roughness. A “fluffy layer” overlies the consolidated seabed, therefore, net seabed erosion occurs after the “fluffy layer” is resuspended. A representative critical bed shear stress for net seabed erosion in the study area was found to be 0.8 Pa. The echo sounder can be an alternative tool for observing silt ripples in coastal regions like the Yellow River Delta where the water is too turbid for underwater videos. The proposed model infers silt-ripple features from bed grain size and flow condition and provides a quick estimate for bed roughness improving the understanding on sediment transport.

Experimental study and theoretical analysis of the shear behavior of single-box double-chamber steel-bamboo composite beams

This study reported the findings of experimental and nonlinear analytical investigations of the shear performance of novel single-box double-chamber steel-bamboo composite (SBC) beams. Eight SBC beams with various parameters were prepared, and then the shear tests were performed to evaluate the shear behavior of the beams. The experimental performances, failure modes, and mechanical behavior were observed and analyzed based on the experimental phenomena and data. The influences of different parameters on deformation behavior, failure modes, and shear capacity were revealed, and analytical models were developed for estimating shear stress and shear capacity in the ultimate limit state (ULS). The results demonstrated that (1) the shear behavior of the beams was excellent, and lateral instability and local buckling behaviors were effectively avoided, the steel and bamboo can maintain excellent synergy and function together before the beams failed; (2) two failure modes were identified, primarily determined by the shear span-to-depth ratio; (3) shear resistance was improved by reducing the shear span-to-depth ratio and steel height/thickness ratio, and increasing the flange bamboo width/thickness ratio and height of the web; (4) the analytical models provided accurate predictions for the beams, and the largest discrepancy was controlled within 10%.

Wind, waves, sediment transport, and potential for shoreline development on Earth and on Mars

The potential for shoreline formation surrounding ancient Martian lake basins is evaluated by reviewing the fundamental processes of how waves form, propagate, break, and run up a beach, while assessing their potential to initiate sediment transport upon reaching a shore, first on Earth and then on Mars. The modeling approach uses a series of well-established semi-empirical equations to simulate the transfer of wind energy to a water surface, grow and propagate wind waves across that surface according to fetch length, and upon reaching the shore, convert the velocity of the breaking waves to the shear stress required to move particles of different sizes on beaches with a range of slopes. This approach is tested by using instrumental wind records to generate wave fields on the Great Salt Lake, Utah and compare the shear stress generated upon wave breaking to the particle size distributions of nine different beaches that were formed during 1986–1987. The critical threshold shear stress required to move the coarse particles (7–28 cm) of these beaches matches reasonably well to the estimated shear stresses derived from the modeled wave climate. This procedure is then applied to ancient Martian beaches, after adjusting for the lower strength of gravity, to estimate the probable size of waves and their ability to transport coarse mixtures of sediment and form shorelines. The results indicate that, given reasonable fetch lengths (25–200 km) and wind velocities (5–25 m/s) during ice-free conditions, large waves (1–10 m) would have likely traversed ancient Martian lakes. Upon reaching a shore these waves would have generated breaking wave velocities of about 1–2 m/s, and been capable of transporting cobbles and small boulders (5–50 cm). The results further imply that shorelines would likely have formed in Martian lake basins.

Comparison of CFD-Derived Estimates of Endothelial Shear Stress to Invasive Coronary Assessment

Comparison of CFD-Derived Estimates of Endothelial Shear Stress to Invasive Coronary Assessment

Semi-analytical method to estimate boundary shear stress in smooth rectangular and trapezoidal open channels

The boundary shear stress is obtained from the continuity and momentum equations in smooth prismatic open channels. Shear stress is a function of gravity, secondary flow, and gradient velocity. In this research, semi-analytical equations were obtained to estimate the average boundary shear stress in smooth rectangular and trapezoidal open channels by 0.5, 1, 1.5 and 2 side slopes using conformal mapping techniques. After dividing the channel cross section into the bed and wall sub-sections, the present study used the method by Guo and Julien (2005) to estimate the average shear stress of the bed and the wall in open channels. For the studied case, by the first shear stress estimation, secondary current and constant eddy viscosity were neglected. After recommending two secondary current and eddy viscosity correction factors, a second approximation well matched the experimental measurements of the average shear stress for all aspect ratios.

Estimation of Mouse Carotid Arterial Wall Shear Stress Using High-Frequency Ultrasound Imaging.

Wall shear stress (WSS) is a crucial hemodynamic factor that promotes atherosclerosis (plaque) development in arteries; although the relationship between WSS and arterial atherosclerosis has been explored in many animal studies, it is not fully understood. No suitable tool, however, exists for rapidly estimating dynamic WSS in small-animal studies. This study proposes a 40-MHz high-frequency ultrasound (HFUS) imaging system for dynamic WSS estimation based on mouse carotid artery blood flow velocity gradient measurements by vector Doppler imaging (VDI). Aliasing reduces the accuracy of Doppler measurements, which can be prevented by increasing the imaging frame rate. Conventionally, imaging is performed at two tilted angles by alternating between the angles; in the proposed method, the frame rate was doubled by imaging at each tilted angle sequentially and by then temporally aligning the sequences based on pulsatile flow characteristics. Velocity estimation using this method had low errors for both a steady-flow straight-tube and pulsatile flow 60%-stenosis phantom. The method was tested for wild-type (WT) C57BL/6 mice at 16 weeks old and apolipoprotein E knockout (ApoE KO) mice at 16 and 24 weeks old; differences in time-averaged and oscillatory WSS were observed, and histology confirmed that the 24-week ApoE KO mice with the highest oscillatory WSS had the greatest plaque formation. The proposed HFUS WSS imaging method can predict the location and extent of plaque development; thus, this method is useful for small-animal studies investigating the WSS effect on atherosclerotic plaque development.