Sponge phase transitions from a lattice model

We use Monte Carlo simulations to obtain thermodynamic functions and correlation functions in a lattice model we propose for sponge phases. We demonstrate that the surface-density correlation function dominates the scattering only along the symmetric-sponge (SS) to asymmetric-sponge (AS) phase boundary but not the boundary between the sponge-with-free-edges (SFE) and symmetric-sponge phases. At this second thermodynamic transition the scattering is dominated instead by an edge-density (or seam-density) correlation function. This prediction provides an unambiguous diagnostic for experiments in search of the SS-SFE transition.

Autocalibration of dynamic bacterial growth model for water distribution system using GA

The presence of residual chlorine and organic matter govern the bacterial regrowth within a water distribution system. The bacterial growth model is essential to predict the spatial and temporal variation of all these substances throughout the system. The parameters governing the bacterial growth and biodegradable dissolved organic carbon (BDOC) utilization are difficult to determine by experimentation. In the present study, the estimation of these parameters is addressed by using simulation-optimization procedure. The optimal solution by genetic algorithm (GA) has indicated that the proper combination of parameter values are significant rather than correct individual values. The applicability of the model is illustrated using synthetic data generated by introducing noise in to the error-free measurements. The GA was found to be a potential tool in estimating the parameters controlling the bacterial growth and BDOC utilization. Further, the GA was also used for evaluating the sensitivity issues relating parameter values and objective function. It was observed that mu and k(cl) are more significant and dominating compared to the other parameters. But the magnitude of the parameters is also an important issue in deciding the dominance of a particular parameter. GA is found to be a useful tool in autocalibration of bacterial growth model and a sensitivity study of parameters.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

Mean field lattice model for adsorption isotherms in zeolite NaA

Using a lattice model for adsorption in microporous materials, pure component adsorption isotherms are obtained within a mean field approximation for methane at 300 K and xenon at 300 and 360 K in zeolite NaA. It is argued that the increased repulsive adsorbate-adsorbate interactions at high coverages must play an important role in determining the adsorption behavior. Therefore, this feature is incorporated through a "coverage-dependent interaction'' model, which introduces a free, adjustable parameter. Another important feature, the site volume reduction, has been treated in two ways: a van der Waal model and a 1D hard-rod theory [van Tassel et al., AIChE J. 40, 925 (1994)]; we have also generalized the latter to include all possible adsorbate overlap scenarios. In particular, the 1D hard-rod model, with our coverage-dependent interaction model, is shown to be in best quantitative agreement with the previous grand canonical Monte Carlo isotherms. The expressions for the isosteric heats of adsorption indicate that attractive and repulsive adsorbate-adsorbate interactions increase and decrease the heats of adsorption, respectively. It is concluded that within the mean field approximation, our simple model for repulsive interactions and the 1D hard-rod model for site volume reduction are able to capture most of the important features of adsorption in confined regions. (C) 1999 American Institute of Physics. [S0021-9606(99)70515-5].

Network Error Correction for Unit-Delay, Memory-Free Networks Using Convolutional Codes

A single source network is said to be memory-free if all of the internal nodes (those except the source and the sinks) do not employ memory but merely send linear combinations of the symbols received at their incoming edges on their outgoing edges. In this work, we introduce network-error correction for single source, acyclic, unit-delay, memory-free networks with coherent network coding for multicast. A convolutional code is designed at the source based on the network code in order to correct network- errors that correspond to any of a given set of error patterns, as long as consecutive errors are separated by a certain interval which depends on the convolutional code selected. Bounds on this interval and the field size required for constructing the convolutional code with the required free distance are also obtained. We illustrate the performance of convolutional network error correcting codes (CNECCs) designed for the unit-delay networks using simulations of CNECCs on an example network under a probabilistic error model.

Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation

This work presents a methodology to reconstruct 3D biological organs from image sequences or other scan data using readily available free softwares with the final goal of using the organs (3D solids) for finite element analysis. The methodology deals with issues such as segmentation, conversion to polygonal surface meshes, and finally conversion of these meshes to 3D solids. The user is able to control the detail or the level of complexity of the solid constructed. The methodology is illustrated using 3D reconstruction of a porcine liver as an example. Finally, the reconstructed liver is imported into the commercial software ANSYS, and together with a cyst inside the liver, a nonlinear analysis performed. The results confirm that the methodology can be used for obtaining 3D geometry of biological organs. The results also demonstrate that the geometry obtained by following this methodology can be used for the nonlinear finite element analysis of organs. The methodology (or the procedure) would be of use in surgery planning and surgery simulation since both of these extensively use finite elements for numerical simulations and it is better if these simulations are carried out on patient specific organ geometries. Instead of following the present methodology, it would cost a lot to buy a commercial software which can reconstruct 3D biological organs from scanned image sequences.

Mathematical model of COREX melter gasifier: Part I. Steady-state model

The COREX melter gasifier is a countercurrent reactor to produce liquid iron. Directly reduced iron (DRI), noncoking coal, and other additives are charged to the melter gasifier at their respective temperatures, and O-2 is blown through the tuyeres. Functionally, a melter gasifier is divided into three zones: a moving bed, fluidized bed, and free board. A model has been developed for the moving bed, where the tuyere region is two-dimensional (2-D) and the rest is one-dimensional (1-D). It is based on multiphase conservation of mass, momentum, and heat. The fluidized bed has been treated as 1-D. Partial equilibrium is calculated for the free board. The calculated temperature of the hot metal, the top gas, and the chemistry of the top gas agree with the reported plant data. The model has been used to study the effects of bed height, injection of impure O-2, coal chemistry, and reactivity on the process performance.

Two-Stage Extended Kalman Filters with Derivative-Free Local Linearizations

This paper proposes a derivative-free two-stage extended Kalman filter (2-EKF) especially suited for state and parameter identification of mechanical oscillators under Gaussian white noise. Two sources of modeling uncertainties are considered: (1) errors in linearization, and (2) an inadequate system model. The state vector is presently composed of the original dynamical/parameter states plus the so-called bias states accounting for the unmodeled dynamics. An extended Kalman estimation concept is applied within a framework predicated on explicit and derivative-free local linearizations (DLL) of nonlinear drift terms in the governing stochastic differential equations (SDEs). The original and bias states are estimated by two separate filters; the bias filter improves the estimates of the original states. Measurements are artificially generated by corrupting the numerical solutions of the SDEs with noise through an implicit form of a higher-order linearization. Numerical illustrations are provided for a few single- and multidegree-of-freedom nonlinear oscillators, demonstrating the remarkable promise that 2-EKF holds over its more conventional EKF-based counterparts. DOI: 10.1061/(ASCE)EM.1943-7889.0000255. (C) 2011 American Society of Civil Engineers.

A Macroscopic Model for First Return Stroke of Lightning

A lightning return stroke model for a downward flash is proposed. The model includes underlying physical phenomena governing return stroke evolution, namely, electric field due to charge distributed along the leader and cloud, transient enhancement of series channel conductance at the bridging regime, and the nonlinear variation of channel conductance, which supports the return stroke current evolution. Thermal effects of free burning arc at the stroke wave front and its impact on channel conductance are studied. A first-order arc model for determining the dynamic channel conductance along with a field-dependent conductivity for corona sheath is used in the model. The model predicts consistent current propagation along the channel with regard to current amplitude and return stroke velocity. The model is also capable of predicting the remote electromagnetic fields that are consistent with the experimental observations.

Nonlinear Model Predictive Spread Acceleration Guidance with Impact Angle Constraint for Stationary Targets

A new technique named as model predictive spread acceleration guidance (MPSAG) is proposed in this paper. It combines nonlinear model predictive control and spread acceleration guidance philosophies. This technique is then used to design a nonlinear suboptimal guidance law for a constant speed missile against stationary target with impact angle constraint. MPSAG technique can be applied to a class of nonlinear problems, which leads to a closed form solution of the lateral acceleration (latax) history update. Guidance command assumed is the lateral acceleration (latax), applied normal to the velocity vector. The new guidance law is validated by considering the nonlinear kinematics with both lag-free as well as first order autopilot delay. The simulation results show that the proposed technique is quite promising to come up with a nonlinear guidance law that leads to both very small miss distance as well as the desired impact angle.

Nonlinear Model Predictive Spread Acceleration Guidance with Impact Angle Constraint for Slow-moving Targets

A new technique named as model predictive spread acceleration guidance (MPSAG) is proposed in this paper. It combines nonlinear model predictive control and spread acceleration guidance philosophies. This technique is then used to design a nonlinear suboptimal guidance law for a constant speed missile against stationary target with impact angle constraint. MPSAG technique can be applied to a class of nonlinear problems, which leads to a closed form solution of the lateral acceleration (latax) history update. Guidance command assumed is the lateral acceleration (latax), applied normal to the velocity vector. The new guidance law is validated by considering the nonlinear kinematics with both lag-free as well as first order autopilot delay. The simulation results show that the proposed technique is quite promising to come up with a nonlinear guidance law that leads to both very small miss distance as well as the desired impact angle.

Docking and free energy simulations to predict conformational domains involved in hCG-LH receptor interactions using recombinant antibodies

Single chain fragment variables (ScFvs) have been extensively employed in studying the protein-protein interactions. ScFvs derived from phage display libraries have an additional advantage of being generated against a native antigen, circumventing loss of information on conformational epitopes. In the present study, an attempt has been made to elucidate human chorionic gonadotropin (hCG)-luteinizing hormone (LH) receptor interactions by using a neutral and two inhibitory ScFvs against hCG. The objective was to dock a computationally derived model of these ScFvs onto the crystal structure of hCG and understand the differential roles of the mapped epitopes in hCG-LH receptor interactions. An anti-hCG ScFv, whose epitope was mapped previously using biochemical tools, served as the positive control for assessing the quality of docking analysis. To evaluate the role of specific side chains at the hCG-ScFv interface, binding free energy as well as residue interaction energies of complexes in solution were calculated using molecular mechanics Poisson-Boltzmann/surface area method after performing the molecular dynamic simulations on the selected hCG-ScFv models and validated using biochemical and SPR analysis. The robustness of these calculations was demonstrated by comparing the theoretically determined binding energies with the experimentally obtained kinetic parameters for hCG-ScFv complexes. Superimposition of hCG-ScFv model onto a model of hCG complexed with the 51-266 residues of LH receptor revealed importance of the residues previously thought to be unimportant for hormone binding and response. This analysis provides an alternate tool for understanding the structure-function analysis of ligand-receptor interactions. Proteins 2011;79:3108-3122. (C) 2011 Wiley-Liss, Inc.

Automatic extraction of free-form surface features (FFSFs)

This paper presents a new algorithm for extracting Free-Form Surface Features (FFSFs) from a surface model. The extraction algorithm is based on a modified taxonomy of FFSFs from that proposed in the literature. A new classification scheme has been proposed for FFSFs to enable their representation and extraction. The paper proposes a separating curve as a signature of FFSFs in a surface model. FFSFs are classified based on the characteristics of the separating curve (number and type) and the influence region (the region enclosed by the separating curve). A method to extract these entities is presented. The algorithm has been implemented and tested for various free-form surface features on different types of free-form surfaces (base surfaces) and is found to correctly identify and represent the features irrespective of the type of underlying surface. The representation and extraction algorithm are both based on topology and geometry. The algorithm is data-driven and does not use any pre-defined templates. The definition presented for a feature is unambiguous and application independent. The proposed classification of FFSFs can be used to develop an ontology to determine semantic equivalences for the feature to be exchanged, mapped and used across PLM applications. (C) 2011 Elsevier Ltd. All rights reserved.

A novel free floating accelerometer force balance system for shock tunnel applications

In order to overcome the interference of the model mounting system with the external aerodynamics of the body during shock tunnel testing, a new free floating internally mountable balance system that ensures unrestrained model motion during testing has been designed, fabricated and tested. Minimal friction ball bearings are used for ensuring the free floating condition of the model during tunnel testing. The drag force acting on a blunt leading edge flat plate at hypersonic Mach number has been measured using the new balance system. Finite element model (FEM) and CFD are exhaustively used in the design as well as for calibrating the new balance system. The experimentally measured drag force on the blunt leading edge flat plate at stagnation enthalpy of 0.7 and 1.2 MJ/kg and nominal Mach number of 5.75 matches well with FEM results. The concept can also be extended for measuring all the three fundamental aerodynamic forces in short duration test facilities like free piston driven shock tunnels.

Geosphere laminations in free groups

We construct for free groups, which are codimension one analogues of geodesic laminations on surfaces. Other analogues that have been constructed by several authors are dimension-one instead of codimension-one. Our main result is that the space of such laminations is compact. This in turn is based on the result that crossing, in the sense of Scott-Swarup, is an open condition. Our construction is based on Hatcher's normal form for spheres in the model manifold.