Modeling the Chesapeake Bay and Tributaries: a synopsis

The last decade has seen the development and application of a spectrum of physical and numerical hydrographic models of the Chesapeake Bay and its tributaries. The success of the James River Hydraulic Model has initiated the construction of an estuarine hydraulic model of the entire Chesapeake System. Numerical analogues for hydrographic behavior and contaminant dispersion in one-, two-, and three dimensional model estuaries exist for various regions of the Bay. From an engineering viewpoint, one dimensional models are sufficiently advanced to be routinely employed in aiding management decisions. Bay investigators are playing leading roles in the development of two- and three-dimensional models of estuarine flows.

Simulation of Natural Gas Production in Hydrate Reservoirs

This paper simulates a one-dimensional physical model of natural gas production from hydrate dissociation in a reservoir by depressurization. According to the principles of solid hydrate decomposition in stratum and flow of natural gas in porous medium, the pressure governing equations for both gas zone and hydrate zone are set up based on the physical production model. Using the approximation reported by N. N. Verigin et al. (1980), the nonlinear governing equations are simplified and the self-similar solutions are obtained. Through calculation, for different reservoir parameters, the distribution characters of pressure are analyzed. The decline trend of natural gas production rate with time is also studied. The simulation results show that production of natural gas from a hydrate reservoir is very sensitive to several reservoir parameters, such as wellbore pressure and stratum porosity and permeability.

A Synoptic View of a Coastal Plain Estuary.

During October, 1972 the Patuxent River Estuary was monitored intensively and synoptically over two tidal cycles to determine the spatial and temporal patterns of various hydrodynamic, chemical and biological features. Forty-one depths at eleven stations along nine transects were sampled simultaneously at hourly intervals for salinity, temperature, dissolved oxygen, chlorohyll a, particulate nitrogen, nitrate, nitrite, total kjeldahl nitrogen, ammonia, particulate carbohydrate, dissolved organic carbon, total hydrolizable phosphorous, dissolved inorganic phosphorous, suspended sediment, particle size distribution, and zooplankton. Tidal velocity was continuously monitored at each depth by recording current meters. Riverine input and meteorological conditions were relatively stable for two weeks preceeding the deployment. This communication describes the calculation of the intrinsic rates of change of the observed variables from their measured distributions in the Estuary. The steady-state, one-dimensional equation of species continuity is employed to separate the advection and tidal dispersion of a hydrodynamically passive substance frbm its intrinsic rate of change at point. A new spatial transform is introduced for the purpose of interpolation and extrapolation of data.The intrinsic rate of change profiles reveal a region of heavy bloom activity in the upper estuary and a secondary bloom near the point in the River that most of the suspended material settles out. The changes in ammonia and nitrates are highly correlated to the productivity patterns. Phosphorous rates are less closely correlated to productivity. The perturbations that the Chalk Point steam electric power plant have on the heat and oxygen balances are easily discernible.

Adiabatic shear localization: occurrence, theories, and applications

Efficient Gate-tunable light-emitting device made of defective boron nitride nanotubes: from ultraviolet to the visible

Boron nitride is a promising material for nanotechnology applications due to its two-dimensional graphene-like, insulating, and highly-resistant structure. Recently it has received a lot of attention as a substrate to grow and isolate graphene as well as for its intrinsic UV lasing response. Similar to carbon, one-dimensional boron nitride nanotubes (BNNTs) have been theoretically predicted and later synthesised. Here we use first principles simulations to unambiguously demonstrate that i) BN nanotubes inherit the highly efficient UV luminescence of hexagonal BN; ii) the application of an external perpendicular field closes the electronic gap keeping the UV lasing with lower yield; iii) defects in BNNTS are responsible for tunable light emission from the UV to the visible controlled by a transverse electric field (TEF). Our present findings pave the road towards optoelectronic applications of BN-nanotube-based devices that are simple to implement because they do not require any special doping or complex growth

Fault-tolerant almost exact state transmission

We show that a category of one-dimensional XY-type models may enable high-fidelity quantum state transmissions, regardless of details of coupling configurations. This observation leads to a fault-tolerant design of a state transmission setup. The setup is fault-tolerant, with specified thresholds, against engineering failures of coupling configurations, fabrication imperfections or defects, and even time-dependent noises. We propose an experimental implementation of the fault-tolerant scheme using hard-core bosons in one-dimensional optical lattices.

Subdiffusion of nonlinear waves in quasiperiodic potentials

19 p.

Pumped double quantum dot with spin-orbit coupling

We study driven by an external electric field quantum orbital and spin dynamics of electron in a one-dimensional double quantum dot with spin-orbit coupling. Two types of external perturbation are considered: a periodic field at the Zeeman frequency and a single half-period pulse. Spin-orbit coupling leads to a nontrivial evolution in the spin and orbital channels and to a strongly spin-dependent probability density distribution. Both the interdot tunneling and the driven motion contribute into the spin evolution. These results can be important for the design of the spin manipulation schemes in semiconductor nanostructures.

On the origin of shear banding instability in metallic glasses

To uncover the physical origin of shear-banding instability in metallic glass (MG), a theoretical description of thermo-mechanical deformation of MG undergoing one-dimensional simple shearing is presented. The coupled thermo-mechanical model takes into account the momentum balance, the energy balance and the dynamics of free volume. The interplay between free-volume production and temperature increase being two potential causes for shear-banding instability is examined on the basis of the homogeneous solution. It is found that the free-volume production facilitates the sudden increase in the temperature before instability and vice versa. A rigorous linear perturbation analysis is used to examine the inhomogeneous deformation, during which the onset criteria and the internal length and time scales for three types of instabilities, namely free-volume softening, thermal softening and coupling softening, are clearly revealed. The shear-banding instability originating from sole free-volume softening takes place easier and faster than that due to sole thermal softening, and dominates in the coupling softening. Furthermore, the coupled thermo-mechanical shear-band analysis does show that an initial slight distribution of local free volume can incur significant strain localization, producing a shear band. During such a localization process, the local free-volume creation occurs indeed prior to the increase in local temperature, indicating that the former is the cause of shear localization, whereas the latter is its consequence. Finally, extension of the above model to include the shear-induced dilatation shows that such dilatation facilitates the shear instability in metallic glasses.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Random response of nonlinear continuous systems

This thesis presents a technique for obtaining the stochastic response of a nonlinear continuous system. First, the general method of nonstationary continuous equivalent linearization is developed. This technique allows replacement of the original nonlinear system with a time-varying linear continuous system. Next, a numerical implementation is described which allows solution of complex problems on a digital computer. In this procedure, the linear replacement system is discretized by the finite element method. Application of this method to systems satisfying the one-dimensional wave equation with two different types of constitutive nonlinearities is described. Results are discussed for nonlinear stress-strain laws of both hardening and softening types.

On the determination of the properties of a medium from its reflection coefficient

This thesis demonstrates how the parameters of a slightly non-homogeneous medium can be derived approximately from the reflection coefficient.

Two types of media are investigated. The first is described by the one-dimensional wave equation, the second by the more complex Timoshenko beam equation. In both cases, the media are assumed to be infinite in extent, with the media parameters becoming homogeneous as the space variable approaches positive or negative infinity.

Much effort is placed in deriving properties of the reflection coefficient for both cases. The wave equation is considered primarily to introduce the techniques used to investigate the more complex Timoshenko equation. In both cases, an approximation is derived for one of the medium parameters involving the reflection coefficient.