A combinatorial optimization scheme for parameter structure identification in ground water modeling

This research develops a methodology for parameter structure identification in ground water modeling. For a given set of observations, parameter structure identification seeks to identify the parameter dimension, its corresponding parameter pattern and values. Voronoi tessellation is used to parameterize the unknown distributed parameter into a number of zones. Accordingly, the parameter structure identification problem is equivalent to finding the number and locations as well as the values of the basis points associated with the Voronoi tessellation. A genetic algorithm (GA) is allied with a grid search method and a quasi-Newton algorithm to solve the inverse problem. GA is first used to search for the near-optimal parameter pattern and values. Next, a grid search method and a quasi-Newton algorithm iteratively improve the GA's estimates. Sensitivities of state variables to parameters are calculated by the sensitivity-equation method. MODFLOW and MT3DMS are employed to solve the coupled flow and transport model as well as the derived sensitivity equations. The optimal parameter dimension is determined using criteria based on parameter uncertainty and parameter structure discrimination. Numerical experiments are conducted to demonstrate the proposed methodology, in which the true transmissivity field is characterized by either a continuous distribution or a distribution that can be characterized by zones. We conclude that the optimized transmissivity zones capture the trend and distribution of the true transmissivity field.     

Experimental evaluation of piecewise exact solution for predicting seismic responses of spherical sliding type isolated structures

The use of base isolation for enhancing seismic resistibility has been proven as an efficient method in experimental and theoretical studies. It is usual to insert a flexible device in the horizontal direction to permit the most of relative deformation of a structure at this level. Because the rigidity of the superstructure is much higher than that of the base isolator underneath the structure, the behavior of the superstructure can be idealized as a rigid body during earthquakes. In general, hundreds of degrees of freedom and a step‐by‐step time history analysis are the basic requirements for calculating the seismic response of a base‐isolated structure under earthquakes. In order to develop a simple tool which can be easily adopted for calculating the seismic responses of the spherical sliding type isolated structures, a piecewise exact solution for predicting the seismic responses of base‐isolated structures has been derived in this study. The comparison between the experimental results conducted at the National Center for Research on Earthquake Engineering (NCREE) in Taiwan and the analytical results obtained from the piecewise exact solution show that the formulation derived in this study can predict the seismic responses of the base‐isolated structure with a very high accuracy. Copyright © 2005 John Wiley & Sons, Ltd.     

Analyzing Spatial and Temporal Changes of Aquaculture in Yunlin County,Taiwan*

This article analyzes spatial and temporal changes of aquaculture land use in Yunlin County, Taiwan, and discusses the driving forces for these changes. Digital land‐use maps and satellite imagery for 1982, 1988, 1994, and 2002 were used to extract aquaculture areas. A geographic information system was used to derive changes in aquaculture areas and to compute the spatial clustering statistics over time. A survey of eighty‐seven farmers was conducted to obtain their views and concerns about aquaculture. The results show a rapid expansion in aquaculture during the 1980s, driven mainly by exports to Japan, and a continuous decline since the mid‐1990s because of economic, environmental, and demographic factors.     

A numerical simulation of debris flow and its application

Debris flow is the flow of solid-fluid mixture and was treated as the flow of a continuum in routing in this study. A mathematical model was proposed to describe debris flow including deposition process and then solved numerically with suitable boundary conditions. Laboratory experiments were also conducted for comparison and calibration of the numerical results as well as for investigation of debris flow phenomena. The numerical model was also applied to simulate the debris flow caused by heavy rainfall in Tungman village of Hualien County located in the east of Taiwan on 23 June 1990. The simulated bed topographies in alluvial fan were in good agreement with those obtained from laboratory experiments and field observations.     

Slip, Stress Drop and Ground Motion of Earthquakes: A View from the Perspective of Fractional Brownian Motion

—The characteristics of slip and stress drop distributions accompanying earthquakes are explored from the perspective of fractional Brownian motion (fBm). Slip and stress drop distributions are assumed to be processes of fBm.The Hurst exponent (H), which reveals the roughness of a random process of fBm, is first estimated from ten inferred slip maps for six crustal earthquakes occurring in California. The relationships between the Hurst exponents with respect to static slip (H _u ),stress (H _τ ), static stress drop (H _δσ ) and slip velocity (H _ú ) are then established following Andrews (1980). They are found to be H _δσ = H _τ = H _u -1 = H _ú -0.5. Empirically, H _u is recognized as being about 1 which, according to the theory of fBm, implies that the static slip distribution of an earthquake is just on the margin between being and not being self-similar, depending on the individual case. Cases where H _u is less than 1 (i.e., self-similar) suggest that H _δσ < 0(i.e., the distribution of static stress drop diverges), which is, in light of fBm, invalid. One possible explanation for this paradox is that H _u is less than 1 in crustal earthquake phenomena only over a certain specific bandwidth of wavenumbers, or it could be that the relation H _δσ = H _u -1 is not valid, which implies that static stress drop in the wavenumber domain is not the product of stiffness and slip as described in Andrews (1980). It could be that some different physics apply over this particular bandwidth. In such cases, multi-fractals may be a better way to explore the characteristics of the Hurst exponents of slip. In general, static stress drop and stress distributions are more likely to be self-similar than static slip distribution. H _u ≌ 1 and H _δσ ≌ 0 are good first approxima tions for the slip and stress drop distributions. The spectrum of ground motion displacement falls off as ω ^{-(Hδσ + 2)} with H _δσ ≌ 0, consistent with an ω^-2 model of the earthquake source.     

Non-classical damping in dynamic analysis of base-isolated structures with internal equipment

It has been shown that the use of base isolation not only attenuates the response of a primary structural system but also reduces the response of a secondary system mounted on or within the main structure. The isolation system, superstructure and equipment may be made of different materials with significantly different energy dissipation characteristics such that the damping matrix for the combined system is non-classical and can only be approximately expressed by modal damping ratios if the classical mode method is used for analysis. The object of this paper is to evaluate the accuracy of this procedure in approximating the responses of base-isolated structures and internal equipment. The complex mode method can provide exact solutions to problems with non-classical damping and is used here to find the exact response of the isolation-superstructure-equipment system. The entire system is assumed to be linear elastic with viscous damping and the superstructure is assumed to be proportionally damped so that the deformation of the superstructure can be expressed in terms of its classical modes. Recognizing that the ratio of the equipment mass to the structural mass and the ratio of the stiffness of the isolation system to the superstructural stiffness are both small, perturbation methods are used to find the response. This study shows that the response of base-isolated structures can be determined by the classical mode method to some degree of accuracy, but the higher frequency content is distorted. The equipment response derived by the classical mode method is much smaller than the exact solution so that the complex mode method should be applied to find equipment response.     

Seismic response of the superstructure and attached equipment in a base-isolated building

Base isolation can be used both to protect the structure and simultaneously to reduce the response of internal equipment. The seismic response of a base-isolated structure has been studied through the shaking table test or numerical calculation before. The object of this paper is to analyse a base-isolated structure by a different analytical approach—perturbation analysis. Recognizing that the horizontal stiffness of an isolation system is much smaller than that of the superstructure, the mathematical expressions of the modal properties of base-isolated structures are derived by the perturbation method in terms of the modal properties of the superstructure and used to study the dynamic response of superstructure and attached equipment in the base-isolated building. This study shows that the first base-isolated mode not only controls the superstructural response but also dominates the response of high-frequency attachment. The contribution of higher modes to the response of base-isolated structures, which is proportional to the horizontal stiffness of isolation system, is very small.     

FORMOSAT-3/COSMIC observations of the ionospheric auroral oval development

The ionospheric radiance and electron density observed by the tiny ionospheric photometer (TIP) and GPS occultation experiment (GOX) payloads on FORMOSAT-3/COSMIC satellites are applied to determine the boundaries of the auroral oval and its width in the winter nighttime ionosphere for both hemispheres. The TIP collects ionospheric emission at 135.6 nm due to electron impact excitation, while the GOX offers ionospheric electron density profiles with radio occultation (RO) technique. Comparison between them shows similar patterns of the plasma structure in the polar caps. The mean width of the auroral bands ranges between about 2 and 11° latitude in the winter nighttime and it varies with longitudes. The comparison by month suggests that the mean radius of the auroral ovals varies with the intensity of the auroral radiance.     

Design and implementation of an adaptive code discriminator in a DSP/FPGA-based Galileo receiver

In the design of a global navigation satellite system receiver, the tracking performance depends on the code tracking loop and the associated discriminator. An adaptive code discriminator under a multi-correlator architecture has been designed based on a multi-objective principle to achieve a balance among various, sometimes conflicting, design objectives. With the proposed optimization approach and adaptive logic, concerns of large pull-in region and small steady-state error can be addressed. The proposed scheme is implemented in a digital signal processor/field programmable gate array board and an experiment is conducted to process GIOVE-A signals. The test results reveal the advantages of the proposed code tracking architecture and discriminator design.     

Cyclic tests of steel plate shear walls using box‐shape vertical boundary elements with or without infill concrete

The steel plate shear wall (SPSW) is an effective lateral force resisting system in which unstiffened steel infill plates are connected to the horizontal and vertical boundary elements (VBEs) on all sides of the plates. The boundary elements must be designed to resist the tension field force of the infill panels. When the VBEs are made from a steel box section, the flange of each box VBE connected with the infill panels can be pulled out‐of‐plane by significant panel forces, called pull‐out action. This study investigates capacity design methods for box VBEs in SPSWs. Simplified fixed beam and portal frame models aim to estimate the pull‐out responses of the flange of the box sections with and without infill concrete, respectively. In this study, cyclic tests of three full‐scale two‐story SPSWs using box VBEs with or without the infill concrete are conducted. Inelastic pushover analyses of the finite element models are conducted. The tests and analytical results confirm that the proposed design methods, which aimed to prevent the full yield of the flange under the pull‐out action, are applicable. Furthermore, the test and analytical results suggest that the initial yielding of the flange of box VBEs under the collective effects of the pull‐out action on the flange, the gravity load, and the sway action on the SPSW represents a local yielding. A strict prohibition of the initial yielding on the flange under the aforementioned collective effects is not recommended for pursuing a cost‐effective design. Copyright © 2017 John Wiley & Sons, Ltd.