Integrated parallelization of computations and visualization for large-scale applications

Critical applications like cyclone tracking and earthquake modeling require simultaneous high-performance simulations and online visualization for timely analysis. Faster simulations and simultaneous visualization enable scientists provide real-time guidance to decision makers. In this work, we have developed an integrated user-driven and automated steering framework that simultaneously performs numerical simulations and efficient online remote visualization of critical weather applications in resource-constrained environments. It considers application dynamics like the criticality of the application and resource dynamics like the storage space, network bandwidth and available number of processors to adapt various application and resource parameters like simulation resolution, simulation rate and the frequency of visualization. We formulate the problem of finding an optimal set of simulation parameters as a linear programming problem. This leads to 30% higher simulation rate and 25-50% lesser storage consumption than a naive greedy approach. The framework also provides the user control over various application parameters like region of interest and simulation resolution. We have also devised an adaptive algorithm to reduce the lag between the simulation and visualization times. Using experiments with different network bandwidths, we find that our adaptive algorithm is able to reduce lag as well as visualize the most representative frames.

Modeling the effect of transmission errors on TCP controlled transfers over infrastructure 802.11 wireless LANs

There have been several studies on the performance of TCP controlled transfers over an infrastructure IEEE 802.11 WLAN, assuming perfect channel conditions. In this paper, we develop an analytical model for the throughput of TCP controlled file transfers over the IEEE 802.11 DCF with different packet error probabilities for the stations, accounting for the effect of packet drops on the TCP window. Our analysis proceeds by combining two models: one is an extension of the usual TCP-over-DCF model for an infrastructure WLAN, where the throughput of a station depends on the probability that the head-of-the-line packet at the Access Point belongs to that station; the second is a model for the TCP window process for connections with different drop probabilities. Iterative calculations between these models yields the head-of-the-line probabilities, and then, performance measures such as the throughputs and packet failure probabilities can be derived. We find that, due to MAC layer retransmissions, packet losses are rare even with high channel error probabilities and the stations obtain fair throughputs even when some of them have packet error probabilities as high as 0.1 or 0.2. For some restricted settings we are also able to model tail-drop loss at the AP. Although involving many approximations, the model captures the system behavior quite accurately, as compared with simulations.

Efficient Dendrimer-DNA Complexation and Gene Delivery Vector Properties of Nitrogen-Core Poly(propyl ether imine) Dendrimer in Mammalian Cells

Dendrimers as vectors for gene delivery were established, primarily by utilizing few prominent dendrimer types so far. We report herein studies of DNA complexation efficacies and gene delivery vector properties of a nitrogen-core poly(propyl ether imine) (PETIM) dendrimer, constituted with 22 tertiary amine internal branches and 24 primary amines at the periphery. The interaction of the dendrimer with pEGFPDNA was evaluated through UV-vis, circular dichroism (CD) spectral studies, ethidium bromide fluorescence emission quenching, thermal melting, and gel retardation assays, from which most changes to DNA structure during complexation was found to occur at a weight ratio of dendrimer:DNA similar to 2:1. The zeta potential measurements further confirmed this stoichiometry at electroneutrality. The structure of a DNA oligomer upon dendrimer complexation was simulated through molecular modeling and the simulation showed that the dendrimer enfolded DNA oligomer along both major and minor grooves, without causing DNA deformation, in 1:1 and 2:1 dendrimer-to-DNA complexes. Atomic force microscopy (AFM) studies on dendrimer-pEGFP DNA complex showed an increase in the average z-height as a result of dendrimers decorating the DNA, without causing a distortion of the DNA structure. Cytotoxicity studies involving five different mammalian cell lines, using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide] (MTT) assay, reveal the dendrimer toxicity profile (IC50) values of similar to 400-1000 mu g mL(-1), depending on the cell line tested. Quantitative estimation, using luciferase assay, showed that the gene transfection was at least 100 times higher when compared to poly(ethylene imine) branched polymer, having similar number of cationic sites as the dendrimer. The present study establishes the physicochemical behavior of new nitrogen-core PETIM dendrimer-DNA complexes, their lower toxicities, and efficient gene delivery vector properties.

Simulations and experiments in punching spring-steel devices with sub-millimeter features

This work demonstrates the feasibility of mesoscale (100 μm to mm) punching of multiple holes of intricate shapes in metals. Analytical modeling, finite element (FE)simulation, and experimentations are used in this work. Two dimensional FE simulations in ABAQUS were done with an assumed material modeling and plane-strain condition. A known analytical model was used and compared with the ABAQUS simulation results to understand the effects of clearance between the punch and the die. FE simulation in ABAQUS was done for different clearances and corner radii at punch, die, and holder. A set of punches and dies were used to punch out a miniature spring-steel gripper. Comparison of compliant grippers manufactured by wire-cut electro discharge machining(EDM) and punching shows that realizing sharp interior and re-entrant corners by punching is not easy to achieve. Punching of circular holes with 5 mm and 2.5 mm diameter is achieved. The possibility of realizing meso-scale parts with complicated shapes through punching is demonstrated in this work; and some strategies are suggested for improvement.

Damage assessment of composite plate structures with material and measurement uncertainty

Composite materials are very useful in structural engineering particularly in weight sensitive applications. Two different test models of the same structure made from composite materials can display very different dynamic behavior due to large uncertainties associated with composite material properties. Also, composite structures can suffer from pre-existing imperfections like delaminations, voids or cracks during fabrication. In this paper, we show that modeling and material uncertainties in composite structures can cause considerable problein in damage assessment. A recently developed C-0 shear deformable locking free refined composite plate element is employed in the numerical simulations to alleviate modeling uncertainty. A qualitative estimate of the impact of modeling uncertainty on the damage detection problem is made. A robust Fuzzy Logic System (FLS) with sliding window defuzzifier is used for delamination damage detection in composite plate type structures. The FLS is designed using variations in modal frequencies due to randomness in material properties. Probabilistic analysis is performed using Monte Carlo Simulation (MCS) on a composite plate finite element model. It is demonstrated that the FLS shows excellent robustness in delamination detection at very high levels of randomness in input data. (C) 2016 Elsevier Ltd. All rights reserved.

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated.We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard–Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Mechanical modeling of grain boundary sliding of polycrystals

Using a variational method, a general three-dimensional solution to the problem of a sliding spherical inclusion embedded in an infinite anisotropic medium is presented in this paper. The inclusion itself is also a general anisotropic elastic medium. The interface is treated as a thin interface layer with interphase anisotropic properties. The displacements in the matrix and the inclusion are expressed as polynomial series of the cartesian coordinate components. Using the virtual work principle, a set of linear algebraic equations about unknown coefficients are obtained. Then the general sliding spherical inclusion problem is accurately solved. Based on this solution, a self-consistent method for sliding polycrystals is proposed. Combining this with a two-dimensional model of an aggregate polycrystal, a systematic analysis of the mechanical behaviour of sliding polycrystals is given in detail. Numerical results are given to show the significant effect of grain boundary sliding on the overall mechanical properties of aggregate polycrystals.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated. We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard-Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Life history and population dynamics of the finetooth shark (Carcharhinus isodon) in the northeastern Gulf of Mexico

The life history and population dynamics of the finetooth shark (Carcharhinus isodon) in the north-eastern Gulf of Mexico were studied by determining age, growth, size-at-maturity, natural mortality, productivity, and elasticity of vital rates of the population. The von Bertalanffy growth model was estimated as Lt=1559 mm TL (1–e–0.24 (t+2.07)) for females and Lt = 1337 mm TL (1–e–0.41 (t+1.39)) for males. For comparison, the Fabens growth equation was also fitted separately to observed size-at-age data, and the fits to the data were found to be similar. The oldest aged specimens were 8.0 and 8.1 yr, and theoretical longevity estimates were 14.4 and 8.5 yr for females and males, respectively. Median length at maturity was 1187 and 1230 mm TL, equivalent to 3.9 and 4.3 yr for males and females, respectively. Two scenarios, based on the results of the two equations used to describe growth, were considered for population modeling and the results were similar. Annual rates of survivorship estimated through five methods ranged from 0.850/yr to 0.607/yr for scenario 1 and from 0.840/yr to 0.590/yr for scenario 2. Productivities were 0.041/yr for scenario 1 and 0.038/yr for scenario 2 when the population level that produces maximum sustain-able yield is assumed to occur at an instantaneous total mortality rate (Z) equaling 1.5 M, and were 0.071/yr and 0.067/yr, when Z=2 M for scenario 1 and 2, respectively. Mean generation time was 6.96 yr and 6.34 yr for scenarios 1 and 2, respectively. Elasticities calculated through simulation of Leslie matrices averaged 12.6% (12.1% for scenario 2) for fertility, 47.7% (46.2% for scenario 2) for juvenile survival, and 39.7% (41.6% for scenario 2) for adult survival. In all, the finetooth shark exhibits life-history and population characteristics intermediate to those of sharks in the small coastal complex and those from some large coastal species, such as the blacktip shark (Carcharhinus limbatus).

Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

Modeling nonlinearities in MEMS oscillators

We present a mathematical model of a microelectromechanical system (MEMS) oscillator that integrates the nonlinearities of the MEMS resonator and the oscillator circuitry in a single numerical modeling environment. This is achieved by transforming the conventional nonlinear mechanical model into the electrical domain while simultaneously considering the prominent nonlinearities of the resonator. The proposed nonlinear electrical model is validated by comparing the simulated amplitude¿frequency response with measurements on an open-loop electrically addressed flexural silicon MEMS resonator driven to large motional amplitudes. Next, the essential nonlinearities in the oscillator circuit are investigated and a mathematical model of a MEMS oscillator is proposed that integrates the nonlinearities of the resonator. The concept is illustrated for MEMS transimpedance-amplifier-based square-wave and sine-wave oscillators. Closed-form expressions of steady-state output power and output frequency are derived for both oscillator models and compared with experimental and simulation results, with a good match in the predicted trends in all three cases. © 1986-2012 IEEE.

Bump-Surface Multicompartment Micelles from a Linear ABC Triblock Copolymer: A Combination Study by Experiment and Computer Simulation

Novel bump-surface multicompartment micelles formed by a linear amphiphilic ABC triblock copolymer via self-assembly in selective solvent were successfully observed both in simulation and experiment. The results revealed that the block A forms the most inner core, and the blocks B and C form the inner and outer layers, respectively, and the bumps were formed by block A and more likely to be born on curving surfaces. Moreover, the micelle shape could be controlled by changing the solvent selectivity of the blocks A and B. Spherical, cylindrical, and discoidal micelles with bumpy surfaces were obtained both in experiment and simulation.