Adaptive Edge Detectors for Piecewise Smooth Data Based on the minmod Limiter

We are concerned with the detection of edges—the location and amplitudes of jump discontinuities of piecewise smooth data realized in terms of its discrete grid values. We discuss the interplay between two approaches. One approach, realized in the physical space, is based on local differences and is typically limited to low-order of accuracy. An alternative approach developed in our previous work [Gelb and Tadmor, Appl. Comp. Harmonic Anal., 7, 101–135 (1999)] and realized in the dual Fourier space, is based on concentration factors; with a proper choice of concentration factors one can achieve higher-orders—in fact in [Gelb and Tadmor, SIAM J. Numer. Anal., 38, 1389–1408 (2001)] we constructed exponentially accurate edge detectors. Since the stencil of these highly-accurate detectors is global, an outside threshold parameter is required to avoid oscillations in the immediate neighborhood of discontinuities. In this paper we introduce an adaptive edge detection procedure based on a cross-breading between the local and global detectors. This is achieved by using the minmod limiter to suppress spurious oscillations near discontinuities while retaining high-order accuracy away from the jumps. The resulting method provides a family of robust, parameter-free edge-detectors for piecewise smooth data. We conclude with a series of one- and two-dimensional simulations.To David Gottlieb, on his 60th birthday, with friendship and appreciation.     

Stabilization Methods for Spectral Element Computations of Incompressible Flows

A stabilization method for the spectral element computation of incompressible flow problems is investigated. It is based on a filtering procedure which consists in filtering the velocity field by a spectral vanishing Helmholtz-type operator at each time step. Relationship between this filtering procedure and SVV-stabilization method, introduced recently in [JCP, 2004, 196(2), p680], is established. A number of numerical examples are presented to show the accuracy and stabilization capability of the method.     

Fatigue Behavior and Retrofit Investigation of Distortion-Induced Web Gap Cracking

This paper studies a Kansas Department of Transportation welded plate girder bridge that developed fatigue cracks at small web gaps close to the girder top flange. Repair had been previously performed by softening the connection plate end with a slot retrofit, but cracks were recently found to have reinitiated at some of the repaired details and are again propagating. A comprehensive finite-element method study was performed to investigate the cracking behavior observed in the bridge and to recommend appropriate measures for future bridge retrofit. The analytical results show that stresses developed at the top flange web gaps could exceed yielding under the loading of an HS15 fatigue truck. The current slot repair used in the bridge was found to have introduced higher magnitude fatigue stresses in the web gap. To achieve a permanent repair of the bridge, it is recommended that a welded connection plate to flange attachment be used during future bridge retrofit. The web gap details should be able to withstand unlimited number of load cycles once this additional repair is performed.     

Measurement of the Timoshenko Shear Stiffness. II: Effect of Transverse Compressibility

This is the second of two papers devoted to the issue of measuring the Timoshenko shear stiffness of thin-walled composite beams. In the first paper, the effect of warping on the effective Timoshenko shear stiffness, as measured through bending tests, was studied. The bending test was simulated using finite-element analysis, and the results indicated that the warping effect was minimal. On the other hand, the evidence suggests that transverse flexibility may have a significant influence on the effective Timoshenko shear stiffness, decreasing the effective shear stiffness at shorter test spans. The purpose of the present study is to further investigate this effect and to explore the use of a sandwich theory to predict the measurement error. A higher-order sandwich theory, which captures the transverse strain at concentrated loads and supports, is applied to a commercially available thin-walled composite beam. The results indicate that the sandwich model does capture the decrease in the effective shear stiffness at short spans, and the dependence of the shear stiffness on span-to-depth ratio is similar to that calculated in the first paper, using the finite-element method.     

Self-Irradiation Cascade Simulations in Plutonium Metal: Model Behavior at High Energy

The Modified Embedded Atom Method model for Pu metal is revised so that it more accurately captures the behavior of the Ziegler-Biersack-Littmark model of ion-ion interactions. Two revision are tested with somewhat different stiffnesses in the 2-1000 eV range. The revised models show higher damage levels at 20 KeV than an earlier model, suggesting that the behavior of the models above 100 eV is dominating damage production, at least in the earlier stages of the cascade. Work was performed at Los Alamos National Laboratory under the auspices of the US Department of Energy, under contract DE-AC52-06NA25396.     

Reliability Analysis of Single-Degree-of-Freedom Elastoplastic Systems. II: Suboptimal Excitations

A method for constructing an approximation of the critical excitation that drives an elastoplastic system from rest to a target threshold at a specified time instant, referred to as the “suboptimal excitation,” is presented in this paper. It is based on the observations gained from study of the critical excitations in the companion paper. Essentially, for the usual case of interest where the failure time is not small compared to the natural period, the duration of the positive and negative pulses of the critical excitation are roughly equal to half of the natural period. This consideration allows for a simple intuitive approximation of the critical excitation. The amplitudes of the positive and negative pulses are obtained in closed forms using energy balance. Numerical investigations show that the critical excitations are well approximated by the suboptimal excitations.     

Enhancement of Flexural Design of Horizontally Curved Prestressed Bridges

In this paper, flexural behavior of horizontally curved prestressed (posttensioned) box bridges is studied by using three-dimensional and refined finite-element modeling and analysis. Bridge length, section geometry, and material properties are the same in all the models, while angle of curvature varies from 0 to 90°. The results of analysis show that in curved bridges, stress distribution is significantly different in comparison to straight bridges. Also, the level of stresses at some locations of section width is considerably high. It is proposed to vary the distribution of the prestressing tendons across section width in order to optimize the bridge capacity. Results show that by proper redistribution of prestressing in section width, significant reduction in resultant stress is possible.     

Applications of the discrete element method in mechanical engineering

Compared to other fields of engineering, in mechanical engineering, the Discrete Element Method (DEM) is not yet a well known method. Nevertheless, there is a variety of simulation problems where the method has obvious advantages due to its meshless nature. For problems where several free bodies can collide and break after having been largely deformed, the DEM is the method of choice. Neighborhood search and collision detection between bodies as well as the separation of large solids into smaller particles are naturally incorporated in the method. The main DEM algorithm consists of a relatively simple loop that basically contains the three substeps contact detection, force computation and integration. However, there exists a large variety of different algorithms to choose the substeps to compose the optimal method for a given problem. In this contribution, we describe the dynamics of particle systems together with appropriate numerical integration schemes and give an overview over different types of particle interactions that can be composed to adapt the method to fit to a given simulation problem. Surface triangulations are used to model complicated, non-convex bodies in contact with particle systems. The capabilities of the method are finally demonstrated by means of application examples. Commemorative Contribution.     

油气集输系统新增产能建设拓扑优化方法研究

在油气集输新增产能建设中,为充分利用已建管网系统的处理能力,达到降低投资,提高经济效益的目的,进行了油气集输新增产能建设布局优化方法的研究。针对新建集输管网为树状结构的情况,建立了油气集输系统新增产能建设拓扑优化的数学模型,讨论了计算复杂性问题。根据模型的特点,将问题分解为拓扑级优化和几何级优化两个子问题,并通过它们之间的迭代进行求解。实例计算表明,所建立的数学模型正确,求解方法可行,能用于工程实际设计。     

Extension of a combined analytical/numerical initial value problem solver for turbulent mixing with combustion

Here, we describe the development of a reacting flow multi-species/combustion methodology, implemented as an extension to the differential reduced ejector analysis (DREA) computer program [Mathematical and computer modeling, vol. 31, 2000, p. 21; Appl. Math. Model. 25 (2001) 427; Comput. Math. Appl. 43(10–11); NASA Contractor Report, 1998]. Use of the single fluid IVP solver framework that was developed for the original DREA model has been directly coupled into the combustion formulation. With these modifications, the analysis has an elementary single step reaction Fuel+Oxidizer→Product combustion capability. Though approximate in nature, the simplicity and efficiency of the DREA formulation make it suitable for its original niche, namely design and preliminary design environments where more complex and expensive models may be inappropriate.