The A Posteriori Aspects of Estuarine Modeling

This exercise is the application of an analytical method for systematically modeling ecosystems data to observations made on a naturally eutrophic, mesohaline planktonic microcosm. The theory and experimental design are briefly outlined and the particular steps in the acutal modeling process follow. Then there is a discussion as to how the whole endeavor can be refined to culminate in models with predictive capabilities. (PDF has 16 pages.)

Modeling Environmental Stress

The word stress when applied to ecosystems is ambiguous. Stress may be low-level, with accompanying near-linear strain, or it may be of finite magnitude, with nonlinear response and possible disintegration of the system. Since there are practically no widely accepted definitions of ecosystem strain, classification of models of stressed systems is tenuous. Despite appearances, most ecosystem models seem to fall into the low-level linear response category. Although they sometimes simulate systems behavior well, they do not provide necessary and sufficient information about sudden structural changes nor structure after transition. Dynamic models of finiteamplitude response to stress are rare because of analytical difficulties. Some idea as to future transition states can be obtained by regarding the behavior of unperturbed functions under limiting strain conditions. Preliminary work shows that, since community variables do respond in a coherent manner to stress, macroscopic analyses of stressed ecosystems offer possible alternatives to compartmental models.

2D modeling of lower mantle structure with WKM synthetics

The Earth is very heterogeneous, especially in the region close to the surface of the Earth, and in regions close to the core-mantle boundary (CMB). The lowermost mantle (bottom 300km of the mantle) is the place for fast anomaly (3% faster S velocity than PREM, modeled from Scd), for slow anomaly (-3% slower S velocity than PREM, modeled from S,ScS), for extreme anomalous structure (ultra-low velocity zone, 30% lower inS velocity, 10% lower in P velocity). Strong anomaly with larger dimension is also observed beneath Africa and Pacific, originally modeled from travel time of S, SKS and ScS. Given the heterogeneous nature of the earth, more accurate approach (than travel time) has to be applied to study the details of various anomalous structures, and matching waveform with synthetic seismograms has proven effective in constraining the velocity structures. However, it is difficult to make synthetic seismograms in more than 1D cases where no exact analytical solution is possible. Numerical methods like finite difference or finite elements are too time consuming in modeling body waveforms. We developed a 2D synthetic algorithm, which is extended from 1D generalized ray theory (GRT), to make synthetic seismograms efficiently (each seismogram per minutes). This 2D algorithm is related to WKB approximation, but is based on different principles, it is thus named to be WKM, i.e., WKB modified. WKM has been applied to study the variation of fast D" structure beneath the Caribbean sea, to study the plume beneath Africa. WKM is also applied to study PKP precursors which is a very important seismic phase in modeling lower mantle heterogeneity. By matching WKM synthetic seismograms with various data, we discovered and confirmed that (a) The D" beneath Caribbean varies laterally, and the variation is best revealed with Scd+Sab beyond 88 degree where Sed overruns Sab. (b) The low velocity structure beneath Africa is about 1500 km in height, at least 1000km in width, and features 3% reduced S velocity. The low velocity structure is a combination of a relatively thin, low velocity layer (200 km thick or less) beneath the Atlantic, then rising very sharply into mid mantle towards Africa. (c) At the edges of this huge Africa low velocity structures, ULVZs are found by modeling the large separation between S and ScS beyond 100 degree. The ULVZ to the eastern boundary was discovered with SKPdS data, and later is confirmed by PKP precursor data. This is the first time that ULVZ is verified with distinct seismic phase.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

Efficient truncation scheme for modeling deepwater mooring lines

This paper is aimed at enabling the confident use of existing model test facilities for ultra deepwater application without having to compromise on the widely accepted range of scales currently used by the floating production industry. Passive line truncation has traditionally been the preferred method of creating an equivalent numerical model at reduced depth; however, these techniques tend to suffer in capturing accurately line dynamic response and so reproducing peak tensions. In an attempt to improve credibility of model test data the proposed truncation procedure sets up the truncated model, based on line dynamic response rather than quasi-static system stiffness. The upper sections of each line are modeled in detail, capturing the wave action zone and all coupling effects with the vessel. These terminate to an approximate analytical model that aims to simulate the remainder of the line. Stages 1 & 2 are used to derive a water depth truncation ratio. Here vibration decay of transverse elastic waves is assessed and it is found that below a certain length criterion, the transverse vibrational characteristics for each line are inertia driven, hence with respect to these motions the truncated model can assume a linear damper whose coefficient depends on the local line properties and vibration frequency. Stage 3 endeavors to match the individual line stiffness between the full depth and truncated models. In deepwater it is likely that taut polyester moorings will be used which are predominantly straight and have high axial stiffness that provides the principal restoring force to static and low frequency vessel motions. Consequently, it means that the natural frequencies of axial vibrations are above the typical wave frequency range allowing for a quasi-static solution. In cases of exceptionally large wave frequency vessel motions, localized curvature at the chain seabed segment and tangential skin drag on the polyester rope can increase dynamic peak tensions considerably. The focus of this paper is to develop an efficient scheme based on analytic formulation, for replicating these forces at the truncation. The paper will close with an example case study of a single mooring under extreme conditions that replicates exactly the static and dynamic characteristics of the full depth line. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

3D modeling of high-T c superconductors by finite element software

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E-J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.

FPGA, physics-based modeling of IGBT and pin diode for hardware co-simulation of complex power electronic converters and systems

This paper presents the steps and the challenges for implementing analytical, physics-based models for the insulated gate bipolar transistor (IGBT) and the PIN diode in hardware and more specifically in field programmable gate arrays (FPGAs). The models can be utilised in hardware co-simulation of complex power electronic converters and entire power systems in order to reduce the simulation time without compromising the accuracy of results. Such a co-simulation allows reliable prediction of the system's performance as well as accurate investigation of the power devices' behaviour during operation. Ultimately, this will allow application-specific optimisation of the devices' structure, circuit topologies as well as enhancement of the control and/or protection schemes.

Design of brushless doubly-fed machines based: On magnetic circuit modeling

This paper proposes a magnetic circuit model (MCM) for the design of a brushless doubly-fed machine (BDFM). The BDFM possesses advantages in terms of high reliability and reduced gearbox stages, and it requires a fractionally-rated power converter. This makes it suitable for utilization in offshore wind turbines. It is difficult for conventional design methods to calculate the flux in the stator because the two sets of stator windings, which have different pole number, form a complex flux pattern which is not easily determined using common analytical approaches. However, it is advantageous to predict the flux density in the teeth and air-gap at the initial design stage for sizing purposes without recourse finite element analysis. Therefore a magnetic circuit model is developed in this paper to calculate the flux density. A BDFM is used as a case study with FEA validation. © 1965-2012 IEEE.

Analytical field-effect method for extraction of subgap states in thin-film transistors

We present an analytical field-effect method to extract the density of subgap states (subgap DOS) in amorphous semiconductor thin-film transistors (TFTs), using a closed-form relationship between surface potential and gate voltage. By accounting the interface states in the subthreshold characteristics, the subgap DOS is retrieved, leading to a reasonably accurate description of field-effect mobility and its gate voltage dependence. The method proposed here is very useful not only in extracting device performance but also in physically based compact TFT modeling for circuit simulation. © 2012 IEEE.

Corrugated all-composite sandwich structures. Part 1: Modeling

An analytical model for the compressive and shear response of monolithic and hierarchical corrugated composite cores has been developed. The stiffness model considers the contribution in stiffness from the bending- and the shear deformations of the core members in addition to the stretching deformation. The strength model is based on the normal stress and shear stress distribution over each core member when subjected to a shear or compressive load condition. The strength model also accounts for initial imperfections. In part 1 of this series, the analytical model is described and the results are compared to finite element predictions. In part 2, the analytical model is compared to experimental results and the behaviour of the corrugated structures is investigated more thoroughly using failure mechanism maps. © 2008 Elsevier Ltd. All rights reserved.

Design and modeling of an integrated device for acoustic resonance spectroscopy

This paper investigates the design and modelling of an integrated device for acoustic resonance spectroscopy (ARS). Miniaturisation of such platforms can be achieved using MEMS technology thereby enabling scaling of device dimensions to investigate smaller specimens while simultaneously operating at higher frequencies. We propose an integrated device where the transducers are mounted in close proximity with the specimen to be analysed (e.g. by integrating ultrasound transducers within a microfluidic channel). A finite element (FE) model and a simplified analytical model have been constructed to predict the acoustic response of a sample embedded in such a device configuration. A FE simulation is performed in COMSOL by embedding the piezoelectric transducers in representative fluid media. Resonant frequencies associated with the measurement can be extracted from this data. The response of various media modelled through FEA matches with analytical predictions for a range of biological media. A variety of biological media may be identified by using the measured resonant frequencies as a signature of relevant physical characteristics. The paper establishes the modelling basis of an integrated acoustic resonant spectrometer that is then applied to examine the impact of geometrical scaling on system resolution. © 2013 IEEE.

Advanced modeling of periodic structures in SEA

The SEA properties of a periodic structure are computed from the FE analysis of a single periodic cell. The periodic theory is used in conjunction with FE so that any geometry can be considered. Some efficient algorithms have been implemented to get the subsystems intrinsic properties (modal density, damping, and equivalent mass), as well as the coupling properties of the subsystem with acoustic subsystems (radiation and transmission). Comparisons with analytical results validate the method. © (2006) by the Katholieke Universiteit Leuven Department of Mechanical Engineering All rights reserved.

Generalized Viscoplastic Modeling of Debris Flow

Various concepts have been proposed or used in the development of rheological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold’s concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non-Newtonian fluid mechanics. Although the generality of Bagnold’s model is still at issue, debris-flow modelers in Japan have generally accepted Takahashi’s formulas derived from Bagnold’s model. Some efforts have recently been made by theoreticians in non-Newtonian fluid mechanics to modify or improve Bagnold’s concept or model. A viable rheological model should consist both of a rate-independent part and a rate-dependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris-flow modeling. In fact, Bagnold’s model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold’s simplified assumption of constant grain concentration.

Formation and modeling of disinfection by-products in drinking water of six cities in China

Funding and support for this project was provided by NSFC (Grant No. 40771015), and Key International Science and Technology Cooperation Projects (Grant No. 22007DFC20180). The authors also gratefully acknowledge the support of Key Projects in the National Science & Technology Pillar Program in the Eleventh Five-year Plan Period (Grant No. 2006BAD01B06-02). The authors thank the CDCs of Daqing, Beijing, Tianjin, Zhengzhou, Changsha and Shenzhen cities for field and laboratory technical support.

Modeling of protein adsorption in membrane affinity chromatography

For the design of affinity membranes, protein adsorption in membrane affinity chromatography (MAC) was studied by frontal analysis. According to fast mass transfer, small thickness of affinity membranes and high affinity between the protein and the ligand, an ideal adsorption (IA) model was proposed for MAC and was used together with equilibrium-dispersive (E-D) model to describe the adsorption of bovine serum albumin (BSA) onto cellulose diacetate/polyethyleneimine (CA/PEI) blend membranes with and without Cu2+ chelating. E-D model was found to better describe the initial region of experimental breakthrough curves. The influence of axial dispersion was revealed and it showed the importance of design of the module to homogenously distribute feed solution. IA model was found to be better for the whole experimental breakthrough curve. According to it, the capacity of affinity membranes and the specificity of the interaction are of equal importance for the design of affinity membranes. An optimum feed concentration was also found in the operation of MAC. The discrepancy between experimental optimum feed concentrations and predicted ones from IA model may be due to the ignorance of some experimental effects such as axial dispersion.