An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Large Scale Numerical Simulations For Multi-Phase Fluid Dynamics With Moving Interfaces

This short communication presents our recent studies to implement numerical simulations for multi-phase flows on top-ranked supercomputer systems with distributed memory architecture. The numerical model is designed so as to make full use of the capacity of the hardware. Satisfactory scalability in terms of both the parallel speed-up rate and the size of the problem has been obtained on two high rank systems with massively parallel processors, the Earth Simulator (Earth simulator research center, Yokohama Kanagawa, Japan) and the TSUBAME (Tokyo Institute of Technology, Tokyo, Japan) supercomputers.

Limiting States Of Internal Wave Rays

Existing models of baroclinic tides are based upon the "traditional approximation'', i. e., neglect of the horizontal component of the Earth's rotation, leading to a well- known conclusion that no freely propagating internal waves can exist beyond the critical latitude and the wave rays are symmetric to the vertical. However, recent studies have contended that the situation may change if both the vertical and horizontal components of the Earth's rotation are taken into account. With the full account of the Coriolis force, characteristics of the internal wavefield generated by tidal flow over uneven topography are investigated. It is found that "nontraditional effects'' profoundly change not only the dynamics of internal waves but also the rate at which the barotropic tidal energy is fed into the internal wavefield. Discarding the traditional approximation, internal waves are proved to be able to generate poleward of the critical latitude, rays of which are no longer symmetric and the limiting values of ray angles become greater or less than 90 degrees, depending on the local latitude and the direction of ray. More importantly, in contrast to the predictions of models based upon the traditional approximation, a substantial conversion occurs in the situations when stratification is so weak that the buoyancy frequency is below the tidal one.

Control dynamics of severe acute respiratory syndrome transmission

Severe acute respiratory syndrome (SARS) is a serious disease with many puzzling features. We present a simple, dynamic model to assess the epidemic potential of SARS and the effectiveness of control measures. With this model, we analysed the SARS epidemic data in Beijing. The data fitting gives the basic case reproduction number of 2.16 leading to the outbreak, and the variation of the effective reproduction number reflecting the control effect. Noticeably, our study shows that the response time and the strength of control measures have significant effects on the scale of the outbreak and the lasting time of the epidemic.

Nonlinear dynamics of atomic force microscopy with intermittent contact

When the atomic force microscopy (AFM) in tapping mode is in intermittent contact with a soft substrate, the contact time can be a significant portion of a cycle, resulting in invalidity of the impact oscillator model, where the contact time is assumed to be infinitely small. Furthermore, we demonstrate that the AFM intermittent contact with soft substrate can induce the motion of higher modes in the AFM dynamic response. Traditional ways of modeling AFM (one degree of freedom (DOF) system or single mode analysis) are shown to have serious mistakes when applied to this kind of problem. A more reasonable displacement criterion on contact is proposed, where the contact time is a function of the mechanical properties of AFM and substrate, driving frequencies/amplitude, initial conditions, etc. Multi-modal analysis is presented and mode coupling is also shown. (c) 2006 Published by Elsevier Ltd.

A dynamic coupling model for hybrid atomistic-continuum computations

A dynamic coupling model is developed for a hybrid atomistic-continuum computation in micro- and nano-fluidics. In the hybrid atomistic-continuum computation, a molecular dynamics (MD) simulation is utilized in one region where the continuum assumption breaks down and the Navier-Stokes (NS) equations are used in another region where the continuum assumption holds. In the overlapping part of these two regions, a constrained particle dynamics is needed to couple the MD simulation and the NS equations. The currently existing coupling models for the constrained particle dynamics have a coupling parameter, which has to be empirically determined. In the present work, a novel dynamic coupling model is introduced where the coupling parameter can be calculated as the computation progresses rather than inputing a priori. The dynamic coupling model is based on the momentum constraint and exhibits a correct relaxation rate. The results from the hybrid simulation on the Couette flow and the Stokes flow are in good agreement with the data from the full MD simulation and the solutions of the NS equations, respectively. (c) 2007 Elsevier Ltd. All rights reserved.

Fiscal imbalances, inflation and sovereign default dynamics

Resumen: La cuestión central que este artículo busca responder es como la política monetaria puede afectar el comportamiento de equilibrio de primas por riesgo soberano y cesación de pagos. El artículo se basa en el modelo de “una-tasa-interés”. La deuda pública se hace riesgosa a causa de una política fiscal activa, como en Uribe (2006), reflejando la habilidad limitada de la autoridad fiscal para controlar el superávit primario. El problema de insolvencia es debido a una oleada de mala suerte (shocks negativos que afectan el superávit primario). Pero en contraste a los resultados de Uribe, a medida que aumenta el costo de la deuda soberana (que resulta de un excedente primario débil), la cesación de pagos se anticipa y es reflejada por una creciente prima de riesgo en el país y una probabilidad de cesación de pagos. La cesación de pagos se define como un incumplimiento de un acuerdo contractual y por ende la decisión es tomada por la autoridad fiscal. Mientras tanto, objetivos conflictivos entre la autoridad monetaria y fiscal juegan un rol importante en llevar a la autoridad fiscal a la cesación de pagos sobre sus pasivos. La característica de la política del gobierno necesaria para restaurar el equilibrio después de la cesación de pagos también es analizada.

Vortex dynamics in the studies of looping in tropical cyclone tracks

This paper proposes criteria for predicting the tendency of looping in tropical cyclone tracks using the approach of vortex dynamics. We model the asymmetric structure of a cyclone by a system of vortex patches. The evolution of such system of vortices is simulated by the method of contour dynamics. A new set of exact analytic formulas for contour dynamics calculations is derived, which is shown to be more computationally effective. Based on point-vortex models, we derive analytic formulas for the criteria of looping in a cyclone track. From numerical experiments, the simulated trajectories obtained from the point-vortex system and vortex patch system agree quite well. Hence, the looping criteria obtained from the point-vortex system can be applied by forecasters to stay alert for tendency of looping in a cyclone track. To demonstrate the applicability of the proposed criteria, the trajectory of Typhoon Yancy (9012), whose field data are available from ''TCM-90'', is simulated. The case study shows that the asymmetric structure similar to the pattern of a beta gyre is responsible for its recurvature when Yancy landed Fujian Province, China on 20 August 1990.

Dislocations in the crystallites of nanocrystalline alpha-Fe - A molecular dynamics study

Molecular dynamics simulations are carried out in order to study the atomic structure of crystalline component, of nanocrystalline alpha-Fe when it is consolidated from small grains. A two-dimensional computational block is used to simulate the consolidation process. All the preset dislocations in the original grains glide out of them in the consolidation process, but new dislocations can generate when the grain size is large enough. It shows that dislocations exist in the consolidated material rather than in the original grains. Whether dislocations exist in the crystalline component of the resultant model nana-material depends upon grain size. The critical value of grain size for dislocation generation appears to be about 9 nm. This result agrees with experiments qualitatively.

A unified model for dislocation nucleation, dislocation emission and dislocation free zone

In this paper, a unified model for dislocation nucleation, emission and dislocation free zone is proposed based on the Peierls framework. Three regions are identified ahead of the crack tip. The emitted dislocations, located away from the crack tip in the form of an inverse pileup, define the plastic zone. Between that zone and the cohesive zone immediately ahead of the crack tip, there is a dislocation free zone. With the stress field and the dislocation density field in the cohesive zone and plastic zone being, respectively, expressed in the first and second Chebyshev polynomial series, and the opening and slip displacements in trigonometric series, a set of nonlinear algebraic equations can be obtained and solved with the Newton-Raphson Method. The results of calculations for pure shearing and combined tension and shear loading after dislocation emission are given in detail. An approximate treatment of the dynamic effects of the dislocation emission is also developed in this paper, and the calculation results are in good agreement with those of molecular dynamics simulations.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

The effect of variable currents on internal solitary waves

The effect of variable currents on internal solitary waves is described within the context of a variable coefficient Korteweg-de Vries (KdV) equation, and the approximate slowly varying, solitary-wave solution of this equation. The general theory which leads to the variable coefficient KdV equation is described; a derivation for the special case when the solitary wave and the current are aligned in the same direction is given in the Appendix. Using further simplifications and approximations, a number of analytical expressions are obtained for the variation in the solitary wave amplitude resulting from variable shear in the basic current or from when the basic current is a depth-independent flow which is a simple representation of a geostrophic current, tidal flow or inertial wave.

A model of an isolated magnetic-flux tube in a stratified atmosphere

A two-dimensional model of a magnetic flux tube confined in a gravitational stratified atmosphere is discussed. The magnetic field in the flux tube is assumed to be force-free. By using the approximation of large scale height, the problem of a free boundary with nonlinear conditions may be reduced to one involving a fixed boundary. The two-dimensional features are obtained by applying the perturbation method and adopting the Luest-Schlueter model as the basic state. The results show that the configuration of a flux tube confined in a gravitational stratified atmosphere is divergent, and the more twisted the magnetic field, the more divergent is the flux tube.