Artificial retardation of barotropic waves in layered ocean models
Jensen, T.G.
1996-06-01
The effect of reducing the barotropic gravity wave speed in a layered ocean model in order to gain computational speed is explored. In theory the error in the propagation of baroclinic gravity waves typically is less than 3% for a reduction of the external gravity speed by one order of magnitude. This is confirmed in a numerical experiment. For baroclinic Rossby waves, the phase speed error is even less. The barotropic response is limited to the reduced radius of deformation. The method, which we will refer to as gravity wave retardation, is therefore applicable only for oceanic flows where the barotropic mode is of minor importance. It is demonstrated that the method gives very good results for the baroclinic flow of an equatorial jet, spinup of a midlatitude ocean and flow over a midoceanic ridge. The method can be considered as an alternative to multilayer reduced gravity models, and has the advantage that bottom topography can be included. 24 refs., 15 figs.
Mixing by barotropic instability in a nonlinear model
NASA Technical Reports Server (NTRS)
Bowman, Kenneth P.; Chen, Ping
1994-01-01
A global, nonlinear, equivalent barotropic model is used to study the isentropic mixing of passive tracers by barotropic instability. Basic states are analytical zonal-mean jets representative of the zonal-mean flow in the upper stratosphere, where the observed 4-day wave is thought to be a result of barotropic, and possibly baroclinic, instability. As is known from previous studies, the phase speed and growth rate of the unstable waves is fairly sensitive to the shape of the zonal-mean jet; and the dominant wave mode at saturation is not necessarily the fastest growing mode; but the unstable modes share many features of the observed 4-day wave. Lagrangian trajectories computed from model winds are used to characterize the mixing by the flow. For profiles with both midlatitude and polar modes, mixing is stronger in midlatitude than inside the vortex; but there is little exchange of air across the vortex boundary. There is a minimum in the Lyapunov exponents of the flow and the particle dispersion at the jet maximum. For profiles with only polar unstable modes, there is weak mixing inside the vortex, no mixing outside the vortex, and no exchange of air across the vortex boundary. These results support the theoretical arguments that, whether wave disturbances are generated by local instability or propagate from other regions, the mixing properties of the total flow are determined by the locations of the wave critical lines and that strong gradients of potential vorticity are very resistant to mixing.
Bottom friction optimization for a better barotropic tide modelling
NASA Astrophysics Data System (ADS)
Boutet, Martial; Lathuilière, Cyril; Son Hoang, Hong; Baraille, Rémy
2015-04-01
At a regional scale, barotropic tides are the dominant source of variability of currents and water heights. A precise representation of these processes is essential because of their great impacts on human activities (submersion risks, marine renewable energies, ...). Identified sources of error for tide modelling at a regional scale are the followings: bathymetry, boundary forcing and dissipation due to bottom friction. Nevertheless, bathymetric databases are nowadays known with a good accuracy, especially over shelves, and global tide models performances are better than ever. The most promising improvement is thus the bottom friction representation. The method used to estimate bottom friction is the simultaneous perturbation stochastic approximation (SPSA) which consists in the approximation of the gradient based on a fixed number of cost function measurements, regardless of the dimension of the vector to be estimated. Indeed, each cost function measurement is obtained by randomly perturbing every component of the parameter vector. An important feature of SPSA is its relative ease of implementation. In particular, the method does not require the development of tangent linear and adjoint version of the circulation model. Experiments are carried out to estimate bottom friction with the HYbrid Coordinate Ocean Model (HYCOM) in barotropic mode (one isopycnal layer). The study area is the Northeastern Atlantic margin which is characterized by strong currents and an intense dissipation. Bottom friction is parameterized with a quadratic term and friction coefficient is computed with the water height and the bottom roughness. The latter parameter is the one to be estimated. Assimilated data are the available tide gauge observations. First, the bottom roughness is estimated taking into account bottom sediment natures and bathymetric ranges. Then, it is estimated with geographical degrees of freedom. Finally, the impact of the estimation of a mixed quadratic/linear friction
Intraseasonal variability in a barotropic model with seasonal forcing
NASA Technical Reports Server (NTRS)
Strong, Christopher M.; Jin, Fei-Fei; Ghil, Michael
1993-01-01
It has recently been suggested that oscillatory topographic instability could contribute to low-frequency variability over the Northern Hemisphere midlatitudes. A barotropic potential vorticity model, with a hierarchy of forcing and topography configurations on the sphere, is used to investigate the nature of low-frequency oscillations induced by such instabilities. Steady-state solutions of the model include multiple unstable equilibria that sustain oscillatory instabilities with periods of 10 - 15 days, and 150 - 180 days, for a realistic forcing pattern, Time-dependent solutions exhibit chaotic behavior with episodic oscillations, featuring both the intraseasonal (35 - 50 day) and biweekly (10 - 15 day) modes. The former is dominated by standing spatial patterns, the latter by traveling wave patterns. The phases of the intraseasonal oscillation are robust for all cases, exhibiting a clear oscillatory exchange of atmospheric angular momentum with the solid earth via mountain troque. It is demonstrated, through linear stability analysis on the sphere, that the intraseasonal oscillations are induced by topographic instabilities. The role of the seasonal cycle is studied by prescribing an annual cycle in the forcing. In this case, the winter forcing is more favorable than the summer for the occurrence of episodic intraseasonal oscillations. Recent observations are consistent with this model result.
Propagation of signals in basin-scale ocean bottom pressure from a barotropic model
NASA Astrophysics Data System (ADS)
Stepanov, Vladimir N.; Hughes, Chris W.
2006-12-01
The exchange of atmospheric plus oceanic mass between ocean basins is investigated using a global barotropic ocean model. We find two particular cases of exchange between two basins. At periods of 4-6 days, the exchange is between the Atlantic and Pacific basins, and represents a known oscillation forced by atmospheric pressure. This mode represents a failure of the inverse-barometer relationship due to the large scale and high frequency of atmospheric forcing, and the presence of continents. Significant exchange between Atlantic and Pacific also occurs at longer periods. The second case is most prominent at periods longer than 30 days (strongest at periods longer than 100 days), and represents a mass exchange between the Southern Ocean and the Pacific. The Southern Ocean part of this exchange is clearly related to the Southern Mode of fluctuations in Antarctic circumpolar transport, forced by Southern Ocean wind stress. The reason for the exchange being with the Pacific rather than other basins is explored, and is found to be related to the balance of wind stress by form stress in Drake Passage: exchange with the Atlantic and Indian oceans becomes dominant if Drake Passage topography is removed. While recognizing the limitations of a barotropic model, we contend that it is necessary to understand the barotropic adjustment process in order to make sense of longer timescale processes. Accordingly, we end with speculation on the possible importance of the barotropic results for global sea level and tropical dynamics.
Exactly Solvable Wormhole and Cosmological Models with a Barotropic Equation of State
NASA Astrophysics Data System (ADS)
Kuhfittig, P. K. F.
An exact solution of the Einstein field equations given the barotropic equation of state $p=\\omega\\rho$ yields two possible models: (1) if $\\omega <-1$, we obtain the most general possible anisotropic model for wormholes supported by phantom energy and (2) if $\\omega >0$, we obtain a model for galactic rotation curves. Here the equation of state represents a perfect fluid which may include dark matter. These results illustrate the power and usefulness of exact solutions.
An implementation of a barotropic quasigeostrophic model of ocean circulation on the MPP
NASA Technical Reports Server (NTRS)
Grosch, C. E.; Fatoohi, R.
1987-01-01
The implementation on the Massively Parallel Processor (MPP) of a barotropic quasigeostrophic model of ocean circulation is discussed. The mathematical model, including scalings and boundary conditions is discussed. The numerical scheme, which uses compact differencing is also discussed. The implementation of this model on the MPP is then presented. Finally, some performance results are given and compared to results obtained using the VPS-32 and one processor of a CRAY-2.
NASA Technical Reports Server (NTRS)
Bowman, Kenneth P.; Sacks, Jerome; Chang, Yue-Fang
1993-01-01
Methods for the design and analysis of numerical experiments that are especially useful and efficient in multidimensional parameter spaces are presented. The analysis method, which is similar to kriging in the spatial analysis literature, fits a statistical model to the output of the numerical model. The method is applied to a fully nonlinear, global, equivalent-barotropic dynamical model. The statistical model also provides estimates for the uncertainty of predicted numerical model output, which can provide guidance on where in the parameter space to conduct further experiments, if necessary. The method can provide significant improvements in the efficiency with which numerical sensitivity experiments are conducted.
Barotropic-baroclinic time splitting for ocean circulation modeling
Higdon, R.L.; Szoeke, R.A. de
1997-07-15
This paper discusses the following topics on oceanic circulation modeling: time splitting in context of the nonlinear primitive equations; analysis of the stability of this splitting when applied to a linearized flow in two-layer fluid with one horizontal dimensions and a flat lower boundary; analysis of this splitting for horizontal dimensions in a rotating reference frame with constant coriolis parameters; and then numerical tests of the new splitting are discussed.
A balanced diagnostic system compatible with a barotropic prognostic model. [for weather forecasting
NASA Technical Reports Server (NTRS)
Ghil, M.; Shkoller, B.; Yangarber, V.
1977-01-01
A system of diagnostic equations for the velocity field, or wind laws, for a barotropic primitive-equation model of large-scale atmospheric flow is derived. Attention is given to the classical balance equation and its ellipticity condition. Numerical solutions of the diagnostic system are presented, including examples of cases of the mixed elliptic-hyperbolic type and cases with non-zero divergence. Procedures for implementing such a system are outlined, along with a review of factors in using the technique for operational numerical weather prediction.
Brittle damage models in DYNA2D
Faux, D.R.
1997-09-01
DYNA2D is an explicit Lagrangian finite element code used to model dynamic events where stress wave interactions influence the overall response of the system. DYNA2D is often used to model penetration problems involving ductile-to-ductile impacts; however, with the advent of the use of ceramics in the armor-anti-armor community and the need to model damage to laser optics components, good brittle damage models are now needed in DYNA2D. This report will detail the implementation of four brittle damage models in DYNA2D, three scalar damage models and one tensor damage model. These new brittle damage models are then used to predict experimental results from three distinctly different glass damage problems.
NASA Astrophysics Data System (ADS)
Einspigel, D.; Sachl, L.; Martinec, Z.
2014-12-01
We present the DEBOT model, which is a new global barotropic ocean model. The DEBOT model is primarily designed for modelling of ocean flow generated by the tidal attraction of the Moon and the Sun, however it can be used for other ocean applications where the barotropic model is sufficient, for instance, a tsunami wave propagation. The model has been thoroughly tested by several different methods: 1) synthetic example which involves a tsunami-like wave propagation of an initial Gaussian depression and testing of the conservation of integral invariants, 2) a benchmark study with another barotropic model, the LSGbt model, has been performed and 3) results of realistic simulations have been compared with data from tide gauge measurements around the world. The test computations prove the validity of the numerical code and demonstrate the ability of the DEBOT model to simulate the realistic ocean tides. The DEBOT model will be principaly applied in related geophysical disciplines, for instance, in an investigation of an influence of the ocean tides on the geomagnetic field or the Earth's rotation. A module for modelling of the secondary poloidal magnetic field generated by an ocean flow is already implemented in the DEBOT model and preliminary results will be presented. The future aim is to assimilate magnetic data provided by the Swarm satellite mission into the ocean flow model.
Ginsparg, P.
1991-01-01
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
Ginsparg, P.
1991-12-31
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
NASA Astrophysics Data System (ADS)
Hess, G. D.; Garratt, J. R.
Data for the cross-isobaric angle 0, the geostrophic drag coefficient Cg, and the functions A and B of Rossby number similarity theory, obtained from meteorological field experiments, are used to evaluate a range of models of the neutral, barotropic planetary boundary layer. The data give well-defined relationships for 0, Cg, and the integrated dissipation rate over the boundary layer, as a function of the surface Rossby number. Lettau's first-order closure mixing-length model gives an excellent fit to the data; other simple models give reasonable agreement. However more sophisticated models, e.g., higher-order closure, large-eddy simulation, direct numerical simulation and laboratory models, give poor fits to the data. The simplemodels have (at least) one free parameter in their turbulence closure that is matched toatmospheric observations; the more sophisticated models either base their closure onmore general flows or have no free closure parameters. It is suggested that all of theatmospheric experiments that we could locate violate the strict simplifying assumptionsof steady, homogeneous, neutral, barotropic flow required by the sophisticated models.The angle 0 is more sensitive to violations of the assumptions than is Cg.
Experiments on the Assimilation of Patches of Data Using a Barotropic Model
NASA Technical Reports Server (NTRS)
Semazzi, F. H. M.
1985-01-01
Experiments were performed to examine the impact of adopting various procedures in the assimilation of patches of data during the progress of a numerical forecast. A barotropic model described by Takacs and Balgovind, 1983 is used for the experiments. The model is global and its discretization is similar to the GLA Fourth Order GCM. The effect of orography is included in the model. The experiments are of an identical twin kind. The nature run starts from global fields of wind (u and v) and height (h) at 300 mb. They are real data of January, 1979 OOZ which were interpolated from the adjacent sigma levels of the GLA GCM and then balanced by the bounded derivative method. From this initial state a 6 day forecast is performed. At intervals of 3 hours this information is saved over sectors of 60 width in longitude, from the south pole to the north pole, to simulate data coverage and observational frequency of a polar orbiting satellite.
Optimizing internal wave drag in a forward barotropic model with semidiurnal tides
NASA Astrophysics Data System (ADS)
Buijsman, M. C.; Arbic, B. K.; Green, J. A. M.; Helber, R. W.; Richman, J. G.; Shriver, J. F.; Timko, P. G.; Wallcraft, A.
2015-01-01
A global tuning experiment for the semidiurnal tide is performed with a barotropic model. The model is forced with the M2 equilibrium tide and accounts for the self-attraction and loading (SAL) term. In addition to a quadratic drag, various linear internal wave drag terms adjusted by a scale factor of O (1) are applied. The drag terms include the original Nycander (2005) tensor scheme, the Nycander tensor scheme reduced at supercritical slopes, and their scalar sisters, a Nycander scalar scheme computed for additional abyssal hill roughness, and the Jayne and St. Laurent (2001) scalar scheme. The Nycander scheme does not have a tunable parameter, but to obtain the best tidal solutions, it is demonstrated that some tuning is unavoidable. It is shown that the scalar Nycander schemes yield slightly lower root-mean square (RMS) elevation errors vs. the data-assimilative TPXO tide model than the tensor schemes. Although the simulation with the optimally tuned original Nycander scalar yields dissipation rates close to TPXO, the RMS error is among the highest. The RMS error is lowered for the reduced schemes, which place relatively more dissipation in deeper water. The inclusion of abyssal hill roughness improves the regional agreement with TPXO dissipation rates, without changing the RMS errors. It is difficult to have each ocean basin optimally tuned with the application of a constant scale factor. The relatively high RMS error in the Atlantic Ocean is reduced with a spatially varying scale factor with a larger value in the Atlantic. Our best global mean RMS error of 4.4 cm for areas deeper than 1000 m and equatorward of 66° is among the lowest obtained in a forward barotropic tide model.
NASA Astrophysics Data System (ADS)
Ribstein, Bruno; Tissier, Ann-Sophie; Zeitlin, Vladimir
2014-05-01
abstract A detailed linear stability analysis of an easterly barotropic Gaussian jet centered at the equator is performed in the long-wave sector in the framework of one- and two-layer shallow-water models on the equatorial β-plane. It is shown that the dominant instability of the jet is due to phase-locking and resonance between Yanai waves, although the standard barotropic and baroclinic instabilities due to the resonance between Rossby waves are also present. In the one-layer case this dominant instability has non-zero growth rate at zero wavenumber for high enough Rossby and low enough Burger numbers, thus reproducing the classical symmetric inertial instability. Yet its asymmetric counterpart has the highest growth rate [1, 2]. In the two-layer case the dominant instability may be barotropic or baroclinic, the latter being stronger, with the maximum of the growth rate shifting towards smaller downstream wavenumbers as Rossby number increases at fixed Burger number and given thickness and density ratios [1, 2]. At large enough Rossby numbers this instability has a non-zero growth rate limit at zero wavenumber, giving the baroclinic symmetric inertial instability. Again, the maximal growth rate is achieved at small but non-zero wavenumbers, corresponding to the asymmetric inertial instability. At high enough Rossby number and low enough Burger number not only the baroclinic, but also the barotropic symmetric instability appears, as well as higher meridional modes of the baroclinic symmetric instability. Still, all of them are dominated by their asymmetric counterparts. Direct numerical simulations of the saturation of the leading instabilities are performed, showing that the barotropic species of the instability saturates by forming a double vortex street subject to nonlinear oscillations, while the baroclinic, the most vigorous one, saturates by producing strong vertical shears and related dissipation and mixing. Key words: Rotating equatorial β-plane flows
C-Field Cosmological Model for Barotropic Fluid Distribution with Varying Λ in FRW Space Time
NASA Astrophysics Data System (ADS)
Bali, Raj; Saraf, Seema
2013-05-01
A cosmological model for barotropic fluid distribution in creation field cosmology with varying cosmological constant ( Λ) in FRW space-time is investigated. To get the deterministic model satisfying conservation equation, we have assumed \\varLambda = 1/R2 as considered by Chen and Wu (in Phys. Rev. D 41:695, 1990) where R is scale factor. We find that creation field ( C) increases with time, which matches with the result of HN Theory (in Hoyle and Narlikar, Proc. R. Astron. Soc. A 282:178, 1964), \\varLambda˜1/t2, the spatial volume increases with time and deceleration parameter q<0 which shows that universe is accelerating. This result matches with recent observations. The inflationary scenario exists in the models and the results so obtained match with astronomical observations. The various special cases of the model (21) viz. dust filled universe ( γ=0) and radiation dominated era ( ρ=3 p), ρ= p (stiff fluid universe) are also discussed. The models are free from horizon.
Gritsun, A
2013-05-28
The theory of chaotic dynamical systems gives many tools that can be used in climate studies. The widely used ones are the Lyapunov exponents, the Kolmogorov entropy and the attractor dimension characterizing global quantities of a system. Another potentially useful tool from dynamical system theory arises from the fact that the local analysis of a system probability distribution function (PDF) can be accomplished by using a procedure that involves an expansion in terms of unstable periodic orbits (UPOs). The system measure (or its statistical characteristics) is approximated as a weighted sum over the orbits. The weights are inversely proportional to the orbit instability characteristics so that the least unstable orbits make larger contributions to the PDF. Consequently, one can expect some relationship between the least unstable orbits and the local maxima of the system PDF. As a result, the most probable system trajectories (or 'circulation regimes' in some sense) may be explained in terms of orbits. For the special classes of chaotic dynamical systems, there is a strict theory guaranteeing the accuracy of this approach. However, a typical atmospheric model may not qualify for these theorems. In our study, we will try to apply the idea of UPO expansion to the simple atmospheric system based on the barotropic vorticity equation of the sphere. We will check how well orbits approximate the system attractor, its statistical characteristics and PDF. The connection of the most probable states of the system with the least unstable periodic orbits will also be analysed. PMID:23588051
NASA Astrophysics Data System (ADS)
Lamouroux, J.; De Mey, P.; Lyard, F.; Jeansou, E.
2006-07-01
The purpose of this study is to assimilate satellite altimetry and tide-gauge data in the barotropic, free-surface, finite element model MOG2D, covering the Bay of Biscay and nested in a North East Atlantic domain. In a first step, we explore the error sub-space of the model in presence of forcing uncertainties, and especially in presence of high frequency atmospheric forcing errors. This is done by an ensemble modelling approach (Monte-Car lo) in which the atmospheric fields are perturbed in a multivariate way: by generating an a priori ensemble of perturbed atmospheric forcing fields (10-meter wind and surface pressure from ARPEGE the meteorological model), and computing the corresponding a posteriori ensemble of model states, one can approximate the forecast errors of the model by ensemble spread statistics. These statistics are shown to be neither homogeneous over the domain, nor stationary, since they are very dependent on the meteorological forcing. Then, the forecast covariance matrix is modelled through forecast error Ensemble EOFs. These statistics, in form of 3D multivariate EOFs (Sea Level Anomaly, barotropic velocities, surface pressure and wind-stress components), are used in a reduced-order sequential scheme, SEQUOIA, set up in an Optimal Interpolation configuration with the MANTA kernel developed at LEGOS/POC (De Mey, 2005), to constrain the model forecast in the framework of twin experiments. In a reference experiment, the data assimilation system is calibrated and sensitivity tests are conducted. The system provides significant error reduction for all state vector variables, but appears to be sensitive to configuration parameters: particularly, one need to constrain atmospheric forcing fields to achieve an efficient control of the model errors. Finally, the capability of realistic observing networks to reduce the model errors are compared. Frequent and regularly spaced observations, such as tide-gauges (SLA) or HF radars and buoys (velocity
NASA Astrophysics Data System (ADS)
Thompson, Keith R.; Sheng, Jinyu
1997-11-01
Currents measured during the winter of 1985-1986 by four moorings on the inner Scotian Shelf are used to assess the hindcast skill of a three-dimensional circulation model forced by local wind stress and coastal sea level. The model is linear and barotropic. The integration scheme is based on the modification of the Galerkin spectral method proposed recently by Sheng and Thompson [1993]. The skill of the model is measured by the variance of the hindcast errors divided by the variance of the observations (henceforth γ2). The model is most effective within 30 km of shore (0.51<=γ2<=0.66). At the mooring in the Nova Scotia Current, a surface intensified southwestward jet with its center approximately 50 km from shore, the model fails to capture the bulk of the variance at current meters within 70 m of the surface (0.81<=γ2<=0.87). The skill of the model is lowest at the mooring 65 km from shore (0.92<=γ2<=0.94). To put these measures of skill into perspective, the currents are hindcast using a linear statistical model with the same inputs as the circulation model. The statistical model is optimal in the sense no other linear model with these inputs can achieve a lower hindcast error variance. For half of the current meters within 30 km of shore and two in the Nova Scotia Current the skill of the circulation model is not significantly lower than that of the statistical model. The largest discrepancies in the skill of the two models are found at the offshore mooring. We show that the suboptimal performance of the circulation model is due in part to the assumption of a spatially uniform wind field. We speculate that another contributing factor is the assumed form of the sea level profile along the open boundary that is upstream in the sense of coastal trapped wave propagation. Finally, the statistical model is used to estimate the increase in skill of the circulation model that may result from the assimilation of hydrographic data and additional coastal sea levels.
High-resolution barotropic tide modelling in the South China Sea
NASA Astrophysics Data System (ADS)
Luu, Quang-Hung; Tkalich, Pavel
2016-04-01
The South China Sea (SCS) links two of the largest open oceans, the Pacific and the Indian, mainly through the Luzon-Taiwan Straits in the northeast and the Malacca-Karimata Straits in the southwest, respectively. It has a rhino-like shape of 3000-km long, whose belly is contiguous to Vietnam and back leans on the Philippines. The highly irregular topography includes the Gulf of Tonkin in the north, the Gulf Thailand in the southwest, and several small islands in the middle of SCS (i.e., the Spratly and the Paracels) resulting in complicated astronomic tides and tidal dynamics in this region. In this study, we present high-resolution simulation of tides in the SCS using the Semi-Implicit Eulerian-Lagrangian Finite-Element (SELFE) model. We derive the bathymetry from the Shuttle Radar Topography Mission (SMRT) 15-arc second dataset, one of the finest global topography data sources. Our particular interest is to resolve small bathymetry features and islands in the middle of the SCS which we obtained by digitizing very-high resolution satellite images (30-m accuracy). An unstructured triangular mesh comprising of up to 5 million nodes is generated to resolve these features with very high accuracy, while maintaining fairly coarse resolution in rest of the domain. The model is configured to run in barotropic mode by forcing harmonic oscillations from FES2012 global tide predictions along open boundaries, adjusted to account for volume transport at key channels in the SCS. Computed surface elevations and currents agree well with available tide predictions and measurements. Sensitivity study is performed to analyze the role of the small bathymetry features on distorting tides in the SCS.
Modeling the barotropic response of the Mediterranean sea level to atmospheric pressure forcing
NASA Astrophysics Data System (ADS)
Natsiopoulos, Dimitrios A.; Vergos, Georgios S.; Tziavos, Ilias N.
2014-05-01
An important characteristic of the Earth's atmosphere with direct impact on the marine environmental and Earth's gravity field are the variations of atmospheric pressure as it often determines wind and weather patterns across the globe. Variations in atmospheric pressure and especially low atmospheric systems affect the values of radar altimeter sea level anomalies (SLA). This response of sea level is closed to the Inverse Barometer (IB) correction given by the altimeters within their geophysical data records. In this work, altimetric data sets from the satellite remote sensing mission of Jason-2, along with their total IB corrections acquired by the on-board altimeters, have been used for a period of forty days between October and November 2013. This period was characterized by extreme low-pressure fields over the Mediterranean Sea and especially in the area of the Ionian and Adriatic Seas and over the island of Rhodes, Greece. The Jason-2 along-track records of the SLA have been used to study both the sea level response to atmospheric pressure change over short time scales (such as ten days) and examine if the barometer correction (local and global) given by the altimeter is close to the expected response (-1 cm/mbar) of sea level to atmospheric pressure change. For the latter, atmospheric pressure data for the period under study were available from the Live Access Server (LAS) of NOAA, as well, provided at four times per day intervals in a grid format. From the LAS atmospheric pressure data, the IB effect was computed and compared with the one provided by the altimeter for its external evaluation. Finally, a regional multiple regression analysis between sea level anomalies, the LAS atmospheric pressure and wind speed components is carried out to model the barotropic response of the Mediterranean to atmospheric wind and pressure forcing.
NASA Astrophysics Data System (ADS)
He, Ruoying; Wilkin, John L.
2006-08-01
A high-resolution hybrid data assimilative (DA) modeling system is used to study barotropic tides and tidal dynamics on the southeast New England shelf. In situ observations include tidal harmonics of 5 major tidal constituents [M2, S2, N2, O1, and K1] analyzed from coastal sea level and bottom pressure gauges. The DA system consists of both forward and inverse models. The former is the three-dimensional, finite difference, nonlinear Regional Ocean Modeling System (ROMS). The latter is a three-dimensional linearized, frequency domain, finite element model TRUXTON. The DA system assimilates in situ observations via the inversion for the barotropic tidal open boundary conditions (OBCs). Model skill is evaluated by comparing the misfits between the observed and modeled tidal harmonics. The assimilation scheme is found effective and efficient in correcting the tidal OBCs, which in turn improve ROMS tidal solutions. Up to 50% decreases of model/data misfits are achieved after inverse data assimilation. Co-amplitude and co-phase maps and tidal current ellipses for each of 5 tidal constituents are generated, revealing complex tidal variability in this transition region between the tidally amplified Gulf of Maine in the northeast and the tidally much less energetic Middle Atlantic Bight in the southwest. Detailed examinations on the residual circulation, energetics, and momentum balances of the M2 tide reveal the key roles of the unique bottom bathymetry of Nantucket Shoals and the complex coastal geometry in affecting the regional tidal dynamics.
WFR-2D: an analytical model for PWAS-generated 2D ultrasonic guided wave propagation
NASA Astrophysics Data System (ADS)
Shen, Yanfeng; Giurgiutiu, Victor
2014-03-01
This paper presents WaveFormRevealer 2-D (WFR-2D), an analytical predictive tool for the simulation of 2-D ultrasonic guided wave propagation and interaction with damage. The design of structural health monitoring (SHM) systems and self-aware smart structures requires the exploration of a wide range of parameters to achieve best detection and quantification of certain types of damage. Such need for parameter exploration on sensor dimension, location, guided wave characteristics (mode type, frequency, wavelength, etc.) can be best satisfied with analytical models which are fast and efficient. The analytical model was constructed based on the exact 2-D Lamb wave solution using Bessel and Hankel functions. Damage effects were inserted in the model by considering the damage as a secondary wave source with complex-valued directivity scattering coefficients containing both amplitude and phase information from wave-damage interaction. The analytical procedure was coded with MATLAB, and a predictive simulation tool called WaveFormRevealer 2-D was developed. The wave-damage interaction coefficients (WDICs) were extracted from harmonic analysis of local finite element model (FEM) with artificial non-reflective boundaries (NRB). The WFR-2D analytical simulation results were compared and verified with full scale multiphysics finite element models and experiments with scanning laser vibrometer. First, Lamb wave propagation in a pristine aluminum plate was simulated with WFR-2D, compared with finite element results, and verified by experiments. Then, an inhomogeneity was machined into the plate to represent damage. Analytical modeling was carried out, and verified by finite element simulation and experiments. This paper finishes with conclusions and suggestions for future work.
NASA Astrophysics Data System (ADS)
Ngodock, Hans; Wallcraft, Alan; Souopgui, Innocent; Richman, James; Shriver, Jay; Arbic, Brian
2016-04-01
The ocean tidal velocity and elevation can be estimated concurrently with the ocean circulation by adding the astronomical tidal forcing, parameterized topographic internal wave drag, and self-attraction and loading to the general circulation physics. However, the accuracy of these tidal estimates does not yet match accuracies in the best data-assimilative barotropic tidal models. This paper investigates the application of an Augmented State Ensemble Kalman Filter (ASEnKF) to improve the accuracy of M2 barotropic tides embedded in a 1/12.5° three-dimensional ocean general circulation model. The ASEnKF is an alternative to the techniques typically used with linearized tide-only models; such techniques cannot be applied to the embedded tides in a nonlinear eddying circulation. An extra term, meant to correct for errors in the tide model due to imperfectly known topography and damping terms, is introduced into the tidal forcing. Ensembles of the model are created with stochastically generated forcing correction terms. The discrepancies for each ensemble member with TPXO, an existing data-assimilative tide model, are computed. The ASEnKF method yields an optimal estimate of the model forcing correction terms, that minimizes resultant root mean square (RMS) tidal sea surface elevation error with respect to TPXO, as well as an estimate of the tidal elevation. The deep-water, global area-averaged RMS sea surface elevation error of the principal lunar semidiurnal tide M2 is reduced from 4.4 cm in a best-case non-assimilative solution to 2.6 cm. The largest elevation errors in both the non-assimilative and ASEnKF solutions are in the North Atlantic, a highly resonant basin. Possible pathways for achieving further reductions in the RMS error are discussed.
NASA Astrophysics Data System (ADS)
Ngodock, Hans E.; Souopgui, Innocent; Wallcraft, Alan J.; Richman, James G.; Shriver, Jay F.; Arbic, Brian K.
2016-01-01
The ocean tidal velocity and elevation can be estimated concurrently with the ocean circulation by adding the astronomical tidal forcing, parameterized topographic internal wave drag, and self-attraction and loading to the general circulation physics. However, the accuracy of these tidal estimates does not yet match accuracies in the best data-assimilative barotropic tidal models. This paper investigates the application of an augmented state ensemble Kalman Filter (ASEnKF) to improve the accuracy of M2 barotropic tides embedded in a 1/12.5° three-dimensional ocean general circulation model. The ASEnKF is an alternative to the techniques typically used with linearized tide-only models; such techniques cannot be applied to the embedded tides in a nonlinear eddying circulation. An extra term, meant to correct for errors in the tide model due to imperfectly known topography and damping terms, is introduced into the tidal forcing. Ensembles of the model are created with stochastically generated forcing correction terms. The discrepancies for each ensemble member with TPXO, an existing data-assimilative tide model, are computed. The ASEnKF method yields an optimal estimate of the model forcing correction terms, that minimizes resultant root mean square (RMS) tidal sea surface elevation error with respect to TPXO, as well as an estimate of the tidal elevation. The deep-water, global area-averaged RMS sea surface elevation error of the principal lunar semidiurnal tide M2 is reduced from 4.4 cm in a best-case non-assimilative solution to 2.6 cm. The largest elevation errors in both the non-assimilative and ASEnKF solutions are in the North Atlantic, a highly resonant basin. Possible pathways for achieving further reductions in the RMS error are discussed.
Ruge, J.; Li, Y.; McCormick, S.F.
1994-12-31
The formulation and time discretization of problems in meteorology are often tailored to the type of efficient solvers available for use on the discrete problems obtained. A common procedure is to formulate the problem so that a constant (or latitude-dependent) coefficient Poisson-like equation results at each time step, which is then solved using spectral methods. This both limits the scope of problems that can be handled and requires linearization by forward extrapolation of nonlinear terms, which, in turn, requires filtering to control noise. Multigrid methods do not suffer these limitations, and can be applied directly to systems of nonlinear equations with variable coefficients. Here, a global barotropic semi-Lagrangian model, developed by the authors, is presented which results in a system of three coupled nonlinear equations to be solved at each time step. A multigrid method for the solution of these equations is described, and results are presented.
Observational constraints on scalar field models of dark energy with barotropic equation of state
Sergijenko, Olga; Novosyadlyj, Bohdan; Durrer, Ruth E-mail: ruth.durrer@unige.ch
2011-08-01
We constrain the parameters of dynamical dark energy in the form of a classical or tachyonic scalar field with barotropic equation of state jointly with other cosmological parameters using the following datasets: the CMB power spectra from WMAP7, the baryon acoustic oscillations in the space distribution of galaxies from SDSS DR7, the power spectrum of luminous red galaxies from SDSS DR7 and the light curves of SN Ia from 2 different compilations: Union2 (SALT2 light curve fitting) and SDSS (SALT2 and MLCS2k2 light curve fittings). It has been found that the initial value of dark energy equation of state parameter is constrained very weakly by most of the data while the other cosmological parameters are well constrained: their likelihoods and posteriors are similar, their forms are close to Gaussian (or half-Gaussian) and the confidence ranges are narrow. The most reliable determinations of the best-fit value and 1σ confidence range for the initial value of the dark energy equation of state parameter are obtained from the combined datasets including SN Ia data from the full SDSS compilation with MLCS2k2 light curve fitting. In all such cases the best-fit value of this parameter is lower than the value of corresponding parameter for current epoch. Such dark energy loses its repulsive properties and in future the expansion of the Universe changes into contraction. We also perform a forecast for the Planck mock data and show that they narrow significantly the confidence ranges of cosmological parameters values, moreover, their combination with SN SDSS compilation with MLCS2k2 light curve fitting may exclude the fields with initial equation of state parameter > −0.1 at 2σ confidence level.
NASA Astrophysics Data System (ADS)
Kantha, Lakshmi H.; Tierney, Craig; Lopez, Joseph W.; Desai, Shailen D.; Parke, Michael E.; Drexler, Laura
1995-12-01
In this second part, we explore the implications of the tides derived from the high-resolution, data-assimilative, nonlinear barotropic global ocean tidal model described by Kantha (this issue) in altimetric analysis and geophysical applications. It is shown that when applied to the task of removing tidal sea surface height from TOPEX altimetric records, the model performance is comparable to other global tidal models in the open ocean as measured by the reduction in crossover variances. The performance is slightly better than that of the only other high-resolution global tidal model from Grenoble (Le Provost et al., 1994). The results are however mixed in regions shallower than 1000 m and in semienclosed seas such as the Bering Sea, with the model performance slightly worse overall than the Grenoble model. Computations of total power input (and hence total tidal dissipation rate) are shown to be in excellent agreement with recent analyses of TOPEX data and geophysical observations. In addition, distributions of the tidal power input, tidal dissipation, and the power fluxes in the global oceans are shown for the two primary constituents, M2 and K1. Load tides in solid Earth due to ocean tidal loading fluctuations are also computed for the major semidiurnal and diurnal constituents. The load tides are shown to be large in the shallow seas adjacent to the coasts with high tides such as the Patagonian shelf, because of the higher resolution of this global tide model. This has potential implications in geophysical applications.
NASA Astrophysics Data System (ADS)
Ribstein, Bruno; Zeitlin, Vladimir; Tissier, Ann-Sophie
2014-05-01
A detailed linear stability analysis of an easterly barotropic Gaussian jet centered at the equator is performed in the long-wave sector in the framework of one- and two-layer shallow-water models on the equatorial β-plane. It is shown that the dominant instability of the jet is due to phase-locking and resonance between Yanai waves, although the standard barotropic and baroclinic instabilities due to the resonance between Rossby waves are also present. In the one-layer case, this dominant instability has non-zero growth rate at zero wavenumber for high enough Rossby and low enough Burger numbers, thus reproducing the classical symmetric inertial instability. Yet its asymmetric counterpart has the highest growth rate. In the two-layer case, the dominant instability may be barotropic or baroclinic, the latter being stronger, with the maximum of the growth rate shifting towards smaller downstream wavenumbers as Rossby number increases at fixed Burger number, and given thickness and density ratios. At large enough Rossby numbers this instability has a non-zero growth rate limit at zero wavenumber, giving the baroclinic symmetric inertial instability. Again, the maximal growth rate is achieved at small but non-zero wavenumbers, corresponding to the asymmetric inertial instability. At high enough Rossby number and low enough Burger number not only the baroclinic, but also the barotropic symmetric instability appears, as well as higher meridional modes of the baroclinic symmetric instability. Still, all of them are dominated by their asymmetric counterparts. Direct numerical simulations of the saturation of the leading instabilities are performed, showing that the barotropic species of the instability saturates by forming a double vortex street subject to nonlinear oscillations, while the baroclinic, the most vigorous one, saturates by producing strong vertical shears and related dissipation and mixing.
Homogenization models for 2-D grid structures
NASA Technical Reports Server (NTRS)
Banks, H. T.; Cioranescu, D.; Rebnord, D. A.
1992-01-01
In the past several years, we have pursued efforts related to the development of accurate models for the dynamics of flexible structures made of composite materials. Rather than viewing periodicity and sparseness as obstacles to be overcome, we exploit them to our advantage. We consider a variational problem on a domain that has large, periodically distributed holes. Using homogenization techniques we show that the solution to this problem is in some topology 'close' to the solution of a similar problem that holds on a much simpler domain. We study the behavior of the solution of the variational problem as the holes increase in number, but decrease in size in such a way that the total amount of material remains constant. The result is an equation that is in general more complex, but with a domain that is simply connected rather than perforated. We study the limit of the solution as the amount of material goes to zero. This second limit will, in most cases, retrieve much of the simplicity that was lost in the first limit without sacrificing the simplicity of the domain. Finally, we show that these results can be applied to the case of a vibrating Love-Kirchhoff plate with Kelvin-Voigt damping. We rely heavily on earlier results of (Du), (CS) for the static, undamped Love-Kirchhoff equation. Our efforts here result in a modification of those results to include both time dependence and Kelvin-Voigt damping.
NASA Astrophysics Data System (ADS)
Lach, Theodore M.
2003-10-01
The CBM (model) of the nucleus has resulted in the prediction of two new quarks, an "up" quark of mass 237.31 MeV/c2 and a "dn" quark of mass 42.392 MeV/c2. These two new predicted quarks helped to determine that the masses of the quarks and leptons are all related by a geometric progression relationship. The mass of each quark or lepton is just the "geometric mean" of two related elementary particles, either in the same generation or in the same family. This numerology predicts the following masses for the electron family: 0.511000 (electron), 7.74 (predicted), 117.3, 1778.4 (tau), 26950.1 MeV. The geometric ratio of this progression is 15.154 (e to the power e). The mass of the tau in this theory agrees very well with accepted values. This theory suggests that all the "dn like" quarks have a mass of just 10X multiples of 4.24 MeV (the mass of the "d" quark). The first 3 "up like" quark masses are 38, 237.31 and 1500 MeV. This theory also predicts a new heavy generation with a lepton mass of 27 GeV, a "dn like" quark of 42.4 GeV, and an "up like" quark of 65 GeV. Significant evidence already exists for the existence of these new quarks, and lepton. Ref. Masses of the Sub-Nuclear Particles, nucl-th/ 0008026, @ http://xxx.lanl.gov. Infinite Energy, Vol 5, issue 30.
NASA Astrophysics Data System (ADS)
Zhu, D.; Nakamura, N.
2009-12-01
Spontaneous formation of a vorticity staircase and multiple jets is simulated using a one dimensional barotropic model on a beta-plane with parameterized eddy mixing. The model represents nearly inviscid geostrophic turbulence characterized by a uniform forcing of pseudomomentum, nonuniform dissipation due to mixing, and no frictional damping of the mean flow. The dissipation of pseudomomentum (diffusive flux of vorticity) is modeled with the effective diffusivity parameterization proposed recently by Ferrari and Nikurashin(2009). Rossby wave dynamics and upscale energy cascade are not modeled explicitly but implicit in the parameterization. The parameterized effective diffusivity is a decreasing function of squared vorticity gradient, revealing the active role of (potential) vorticity gradient as a barrier to mixing, consistent with the Rossby elasticity idea. Not only does the parameterized diffusivity agree well with the effective diffusivity of a direct numerical simulation, but it allows the 1D model to reproduce other salient features of the direct simulation, most notably the formation of a welldefined vorticity staircase from a uniform vorticity gradient, through inhomogeneous mixing of vorticity. The staircase formation starts as a small-scale, antidiffusive instability in vorticity gradient that develops when the eddy scale is comparable to the Rhines scale. This spawns numerous gaps (barriers) in diffusivity and corresponding small steps in vorticity, but many of them become unstable and disappear later, until a few stable ones remain. The final number of barriers (vorticity steps) is predictable to a good approximation with a few model parameters.
NASA Astrophysics Data System (ADS)
Limbach, P.; Müller, T.; Skoda, R.
2015-12-01
Commonly, for the simulation of cavitation in centrifugal pumps incompressible flow solvers with VOF kind cavitation models are applied. Since the source/sink terms of the void fraction transport equation are based on simplified bubble dynamics, empirical parameters may need to be adjusted to the particular pump operating point. In the present study a barotropic cavitation model, which is based solely on thermodynamic fluid properties and does not include any empirical parameters, is applied on a single flow channel of a pump impeller in combination with a time-explicit viscous compressible flow solver. The suction head curves (head drop) are compared to the results of an incompressible implicit standard industrial CFD tool and are predicted qualitatively correct by the barotropic model.
Technical Review of the UNET2D Hydraulic Model
Perkins, William A.; Richmond, Marshall C.
2009-05-18
The Kansas City District of the US Army Corps of Engineers is engaged in a broad range of river management projects that require knowledge of spatially-varied hydraulic conditions such as velocities and water surface elevations. This information is needed to design new structures, improve existing operations, and assess aquatic habitat. Two-dimensional (2D) depth-averaged numerical hydraulic models are a common tool that can be used to provide velocity and depth information. Kansas City District is currently using a specific 2D model, UNET2D, that has been developed to meet the needs of their river engineering applications. This report documents a tech- nical review of UNET2D.
NASA Astrophysics Data System (ADS)
Lal, A. K.; Pathania, Ankush; Bhalla, Alka; Mohan, C.
2009-12-01
Mohan et al (1992 Astrophys. Space. Sci. 193 69) (1998 Indian J. Pure Appl. Math. 29 199) investigated the problem of equilibrium structures and periods of small adiabatic oscillations of differentially rotating stellar models using a law of differential rotation of the type ω2 = b0 + b1s2 + b2s4 (here ω is a nondimensional measure of the angular velocity of rotation of a fluid element at a distance s from the axis of rotation and b's are suitably chosen constant parameters). This law of differential rotation assumes cylindrical symmetry for the rotating fluid elements. In the present paper, we have extended their study and used a more general law of differential rotation of the type ω2 = b0 + b1s2 + b2s4 + b3z2 + b4z4 + b5z2s2 in which the angular velocity of rotation of a fluid element is assumed to depend both on its distance s from the axis of rotation and on its distance z from the plane through the center of the star perpendicular to the axis of rotation. The main objective of this study has been to investigate whether the dependence of angular velocity of rotation on the parameter z in addition to the parameter s substantially alters the behavior of the eigenfrequencies of small adiabatic barotropic modes of oscillations of differentially rotating stars or not.
An Intercomparison of 2-D Models Within a Common Framework
NASA Technical Reports Server (NTRS)
Weisenstein, Debra K.; Ko, Malcolm K. W.; Scott, Courtney J.; Jackman, Charles H.; Fleming, Eric L.; Considine, David B.; Kinnison, Douglas E.; Connell, Peter S.; Rotman, Douglas A.; Bhartia, P. K. (Technical Monitor)
2002-01-01
A model intercomparison among the Atmospheric and Environmental Research (AER) 2-D model, the Goddard Space Flight Center (GSFC) 2-D model, and the Lawrence Livermore National Laboratory 2-D model allows us to separate differences due to model transport from those due to the model's chemical formulation. This is accomplished by constructing two hybrid models incorporating the transport parameters of the GSFC and LLNL models within the AER model framework. By comparing the results from the native models (AER and e.g. GSFC) with those from the hybrid model (e.g. AER chemistry with GSFC transport), differences due to chemistry and transport can be identified. For the analysis, we examined an inert tracer whose emission pattern is based on emission from a High Speed Civil Transport (HSCT) fleet; distributions of trace species in the 2015 atmosphere; and the response of stratospheric ozone to an HSCT fleet. Differences in NO(y) in the upper stratosphere are found between models with identical transport, implying different model representations of atmospheric chemical processes. The response of O3 concentration to HSCT aircraft emissions differs in the models from both transport-dominated differences in the HSCT-induced perturbations of H2O and NO(y) as well as from differences in the model represent at ions of O3 chemical processes. The model formulations of cold polar processes are found to be the most significant factor in creating large differences in the calculated ozone perturbations
Barotropic Tidal Predictions and Validation in a Relocatable Modeling Environment. Revised
NASA Technical Reports Server (NTRS)
Mehra, Avichal; Passi, Ranjit; Kantha, Lakshmi; Payne, Steven; Brahmachari, Shuvobroto
1998-01-01
Under funding from the Office of Naval Research (ONR), and the Naval Oceanographic Office (NAVOCEANO), the Mississippi State University Center for Air Sea Technology (CAST) has been working on developing a Relocatable Modeling Environment(RME) to provide a uniform and unbiased infrastructure for efficiently configuring numerical models in any geographic/oceanic region. Under Naval Oceanographic Office (NAVO-CEANO) funding, the model was implemented and tested for NAVOCEANO use. With our current emphasis on ocean tidal modeling, CAST has adopted the Colorado University's numerical ocean model, known as CURReNTSS (Colorado University Rapidly Relocatable Nestable Storm Surge) Model, as the model of choice. During the RME development process, CURReNTSS has been relocated to several coastal oceanic regions, providing excellent results that demonstrate its veracity. This report documents the model validation results and provides a brief description of the Graphic user Interface (GUI).
Richardson's Barotropic Forecast: A Reappraisal.
NASA Astrophysics Data System (ADS)
Lynch, Peter
1992-01-01
To elucidate his numerical technique and to examine the effectiveness of geostrophic initial winds, Lewis Fry Richardson carried out an idealized forecast using the linear shallow-water equations and simple analytical pressure and velocity fields. This barotropic forecast has been repeated and extended using a global numerical model, and the results are presented in this paper. Richardson's conclusions regarding the use of geostrophic winds as initial data are reconsidered.An analysis of Richardson's data into normal modes shows that almost 85% of the energy is accounted for by a single eigenmode, the gravest symmetric rotational Hough mode, which travels westward with a period of about five days. This five-day wave has been detected in analyses of stratospheric data. It is striking that the fields chosen by Richardson on considerations of smoothness should so closely resemble a natural oscillation of the atmosphere.The numerical model employed in this study uses an implicit differencing technique, which is stable for large time steps. The numerical instability that would have destroyed Richardson's barotropic forecast, had it been extended, is thereby circumvented. It is sometimes said that computational instability was the cause of the failure of Richardson's baroclinic forecast, for which he obtained a pressure tendency value two orders of magnitude too large. However, the initial tendency is independent of the time step (at least for the explicit scheme used by Richardson). In fact, the spurious tendency resulted from the presence of unrealistically large high-frequency gravity-wave components in the initial fields.High-frequency oscillations are also found in the evolution starting from the idealized data in the barotropic forecast. They are shown to be due to the gravity-wave components of the initial data. These oscillations may be removed by a slight modification of the initial fields. This initialization is effected by means of a simple digital filtering
Aetherizing Lambda: Barotropic fluids as dark energy
Linder, Eric V.; Scherrer, Robert J.
2009-07-15
We examine the class of barotropic fluid models of dark energy, in which the pressure is an explicit function of the density, p=f({rho}). Through general physical considerations we constrain the asymptotic past and future behaviors and show that this class is equivalent to the sum of a cosmological constant and a decelerating perfect fluid, or 'aether', with w{sub AE}{>=}0. Barotropic models give substantially disjoint predictions from quintessence, except in the limit of {lambda}CDM. They are also interesting in that they simultaneously can ameliorate the coincidence problem and yet 'predict' a value of w{approx_equal}-1.
NASA Astrophysics Data System (ADS)
Einšpigel, David; Martinec, Zdeněk
2015-08-01
The purpose of this paper is to present a new global barotropic ocean model-the DEBOT model. The model is based on the shallow water equations which we newly express in geographical coordinates. The derivation includes the boundary conditions and the Reynolds tensor in a form used commonly in oceanography. The numerical model employs finite differences on an Arakawa-C grid in space and a generalized forward-backward scheme in time with a combined third-order Adams-Bashforth and fourth-order Adams-Moulton step. The validity of the model is demonstrated by the tests based on conservation integral invariants. As a practical application, we present ocean circulation simulations generated by the lunisolar tidal force.
Implementation of a barotropic-baroclinic time splitting for isopynic coordinate ocean modeling
Higdon, R.L.
1999-01-20
The motions seen in numerical models of ocean circulation typically include rapidly moving external gravity waves, which are essentially two-dimensional, and a variety of other motions which are much slower and generally three-dimensional. In a number of ocean models, the computational problems resulting from the multiple time scales are addressed by splitting the fast and slow dynamics into separate subproblems that are solved by different techniques. The present paper addresses the practical implementation of such a splitting for the case of isopynic ocean modeling, in which the vertical coordinate is density or some other related quantity. During the implementation of such a splitting for the case of isopynic ocean modeling, in which the vertical coordinate is density or some other related quantity. During the implementation, a major task is to develop formulas suitable for usage with rapidly varying bottom topography. Compared to an earlier splitting, the present splitting has improved stability properties when analyzed in a simple linearized setting. When tested in the same model for which the earlier splitting was developed, the revised splitting substantially reduces a numerically induced sloshing pattern that is seen in the model. This phenomenon appears to be related to a residual term that appears in the momentum equations. In some examples, the residual is one to two orders of magnitude smaller in the case of the revised splitting.
2-D Magnetohydrodynamic Modeling of A Pulsed Plasma Thruster
NASA Technical Reports Server (NTRS)
Thio, Y. C. Francis; Cassibry, J. T.; Wu, S. T.; Rodgers, Stephen L. (Technical Monitor)
2002-01-01
Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) MK-1 pulsed plasma thruster. Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.
2 1/2 -D compressible reconnection model
NASA Astrophysics Data System (ADS)
Skender, M.; Vršnak, B.
The exact solution of the jump conditions on the RD/SMS discontinuity system in a two-and-half-dimensional (2 1/2 -D) symmetrical reconnection model enables one to analyse the outflowing jet characteristics in dependence on the inflow velocity, and to follow changes in transition to the two-dimensional model. Implications arising from the exact solution and its relevance for solar flares are discussed.
NASA Astrophysics Data System (ADS)
Sato, K.; Masuda, A.
2014-12-01
In general, atmospheric motions are categorized into two: One is the quasi-geostrophic (QG) flows including planetary waves (PWs), and the other is gravity waves (GWs). The GWs are largely affected by the QG flows in their generation, propagation, and dissipation, but GWs can modify the QG flows such as the weak wind layer in the upper mesosphere by their ability of momentum transport. In the winter mesosphere, a necessary condition of barotropic and/or baroclinic instability for the QG flow, i.e., negative latitudinal gradient of potential vorticity (PV), is often satisfied. This study examines dynamical mechanism of the formation of such instability condition in boreal winter and discusses the significant role of the GW forcing. We used simulation data from a GW-resolving general circulation model (GCM). As this GCM does not include any GW parameterizations, all waves including GWs are resolved, which allows us to analyze the role of GWs in the momentum budget in the middle atmosphere explicitly. First, a two-dimensional (2-d) analysis using the transformed Eulerian-mean equations was made. It is seen that the negative PV gradient is regarded as an enhanced PV maximum. This maximum is due to the poleward shift of the westerly jet in association with strong EP-flux divergence caused by PWs from the troposphere. Strong GW drag slightly above the westerly jet shifts poleward as well, which can be understood by a selective GW-filtering mechanism. It seems that this GW-drag shift induces strong upwelling in the middle latitudes and adiabatically cools the middle mesosphere. Resultant enhanced static stability is the main cause of the PV maximum in the upper mesosphere. Because of the dominance of PWs during this event, this process may not be zonally uniform. Thus, second, a 3-d analysis was made using a recently derived 3-d transformed Eulerian-mean theory. As expected, the GW drag is distributed depending on the longitude. The zonal structure of PV maximum is
Baroclinic-Barotropic Interaction in Teleconnections
NASA Astrophysics Data System (ADS)
Ji, Xuan
Teleconnections are remote impacts from one location to another, which are transmitted through planetary-scale wave motions, e.g., Rossby waves. Teleconnections can be forced by tropical heat sources, for example, El Nino/Southern Oscillation (ENSO), and result in many societally important impacts in remote ocean basins and continents. It is thus important to understand the complex pathway and wave dynamics in teleconnections. This dissertation discusses the role of baroclinic-barotropic interaction in three different aspects of atmospheric teleconnections. The first major question addressed is what mechanisms control the interhemispheric teleconnections from tropical heat sources. These are investigated using an intermediate complexity model [a Quasi-equilibrium Tropical Circulation Model (QTCM)] and a simple linear two-level model with dry dynamics. Illustrating the interhemispheric teleconnection process with an Atlantic Warm Pool principal case, the heat source directly excites a baroclinic response that spreads across the equator. Three processes involving baroclinic-barotropic interactions---shear advection, surface drag, and vertical advection---then force a cross-equatorial barotropic Rossby wave response. An analysis of these processes in QTCM simulations indicates that: (1) shear advection has a pattern that roughly coincides with the baroclinic signal in the tropics and subtropics; (2) surface drag has large amplitude and spatial extent, and can be very effective in forcing barotropic motions around the globe; (3) vertical advection has a significant contribution locally and remotely where large vertical motions and vertical shear occur. A simple model is modified to perform experiments in which each of these three mechanisms may be included or omitted. By adding surface drag and vertical advection, and comparing each to shear advection, the effects of the three mechanisms on the generation and propagation of the barotropic Rossby waves are shown to be
Modelling RF sources using 2-D PIC codes
Eppley, K.R.
1993-03-01
In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT`S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field (``port approximation``). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.
Modelling RF sources using 2-D PIC codes
Eppley, K.R.
1993-03-01
In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT'S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field ( port approximation''). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.
Unitary matrix models and 2D quantum gravity
Dalley, S. . Joseph Henry Labs.); Johnson, C.V.; Morris, T.R. . Dept. of Physics); Watterstam, A. )
1992-09-21
In this paper the KdV and modified KdV integrable hierarchies are shown to be different descriptions of the same 2D gravitational system - open-closed string theory. Non-perturbative solutions of the multicritical unitary matrix models map to non-singular solutions of the renormalization group equation for the string susceptibility, [P, Q] = Q. The authors also demonstrate that the large-N solutions of unitary matrix integrals in external fields, studied by Gross and Newman, equal the non-singular pure closed-string solutions of [[bar P], Q] = Q.
2D numerical modelling of meandering channel formation
NASA Astrophysics Data System (ADS)
XIAO, Y.; ZHOU, G.; YANG, F. S.
2016-03-01
A 2D depth-averaged model for hydrodynamic sediment transport and river morphological adjustment was established. The sediment transport submodel takes into account the influence of non-uniform sediment with bed surface armoring and considers the impact of secondary flow in the direction of bed-load transport and transverse slope of the river bed. The bank erosion submodel incorporates a simple simulation method for updating bank geometry during either degradational or aggradational bed evolution. Comparison of the results obtained by the extended model with experimental and field data, and numerical predictions validate that the proposed model can simulate grain sorting in river bends and duplicate the characteristics of meandering river and its development. The results illustrate that by using its control factors, the improved numerical model can be applied to simulate channel evolution under different scenarios and improve understanding of patterning processes.
Brane brick models and 2 d (0 , 2) triality
NASA Astrophysics Data System (ADS)
Franco, Sebastián; Lee, Sangmin; Seong, Rak-Kyeong
2016-05-01
We provide a brane realization of 2 d (0 , 2) Gadde-Gukov-Putrov triality in terms of brane brick models. These are Type IIA brane configurations that are T-dual to D1-branes over singular toric Calabi-Yau 4-folds. Triality translates into a local transformation of brane brick models, whose simplest representative is a cube move. We present explicit examples and construct their triality networks. We also argue that the classical mesonic moduli space of brane brick model theories, which corresponds to the probed Calabi-Yau 4-fold, is invariant under triality. Finally, we discuss triality in terms of phase boundaries, which play a central role in connecting Calabi-Yau 4-folds to brane brick models.
Experimental validation of 2D profile photoresist shrinkage model
NASA Astrophysics Data System (ADS)
Bunday, Benjamin; Cordes, Aaron; Self, Andy; Ferry, Lorena; Danilevsky, Alex
2011-03-01
For many years, lithographic resolution has been the main obstacle in allowing the pace of transistor densification to meet Moore's Law. For the 32 nm node and beyond, new lithography techniques will be used, including immersion ArF (iArF) lithography and extreme ultraviolet lithography (EUVL). As in the past, these techniques will use new types of photoresists with the capability to print smaller feature widths and pitches. These smaller feature sizes will also require the use of thinner layers of photoresists, such as under 100 nm. In previous papers, we focused on ArF and iArF photoresist shrinkage. We evaluated the magnitude of shrinkage for both R&D and mature resists as a function of chemical formulation, lithographic sensitivity, scanning electron microscope (SEM) beam condition, and feature size. Shrinkage results were determined by the well accepted methodology described in SEMATECH's CD-SEM Unified Specification. In other associated works, we first developed a 1-D model for resist shrinkage for the bottom linewidth and then a 2-D profile model that accounted for shrinkage of all aspects of a trapezoidal profile along a given linescan. A fundamental understanding of the phenomenology of the shrinkage trends was achieved, including how the shrinkage behaves differently for different sized and shaped features. In the 1-D case, calibration of the parameters to describe the photoresist material and the electron beam was all that was required to fit the models to real shrinkage data, as long as the photoresist was thick enough that the beam could not penetrate the entire layer of resist. The later 2-D model included improvements for solving the CD shrinkage in thin photoresists, which is now of great interest for upcoming realistic lithographic processing to explore the change in resist profile with electron dose and to predict the influence of initial resist profile on shrinkage characteristics. The 2-D model also included shrinkage due to both the primary
A 2D channel-clogging biofilm model.
Winstanley, H F; Chapwanya, M; Fowler, A C; O'Brien, S B G
2015-09-01
We present a model of biofilm growth in a long channel where the biomass is assumed to have the rheology of a viscous polymer solution. We examine the competition between growth and erosion-like surface detachment due to the flow. A particular focus of our investigation is the effect of the biofilm growth on the fluid flow in the pores, and the issue of whether biomass can grow sufficiently to shut off fluid flow through the pores, thus clogging the pore space. Net biofilm growth is coupled along the pore length via flow rate and nutrient transport in the pore flow. Our 2D model extends existing results on stability of 1D steady state biofilm thicknesses to show that, in the case of flows driven by a fixed pressure drop, full clogging of the pore can indeed happen in certain cases dependent on the functional form of the detachment term. PMID:25240390
Mass loss in 2D rotating stellar models
Lovekin, Caterine; Deupree, Bob
2010-10-05
Radiatively driven mass loss is an important factor in the evolution of massive stars . The mass loss rates depend on a number of stellar parameters, including the effective temperature and luminosity. Massive stars are also often rapidly rotating, which affects their structure and evolution. In sufficiently rapidly rotating stars, both the effective temperature and radius vary significantly as a function of latitude, and hence mass loss rates can vary appreciably between the poles and the equator. In this work, we discuss the addition of mass loss to a 2D stellar evolution code (ROTORC) and compare evolution sequences with and without mass loss. Preliminary results indicate that a full 2D calculation of mass loss using the local effective temperature and luminosity can significantly affect the distribution of mass loss in rotating main sequence stars. More mass is lost from the pole than predicted by 1D models, while less mass is lost at the equator. This change in the distribution of mass loss will affect the angular momentum loss, the surface temperature and luminosity, and even the interior structure of the star. After a single mass loss event, these effects are small, but can be expected to accumulate over the course of the main sequence evolution.
2D Quantum Transport Modeling in Nanoscale MOSFETs
NASA Technical Reports Server (NTRS)
Svizhenko, Alexei; Anantram, M. P.; Govindan, T. R.; Biegel, Bryan
2001-01-01
With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density- gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Quantum simulations are focused on MIT 25, 50 and 90 nm "well- tempered" MOSFETs and compared to classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are quantitatively consistent with I D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and sub-threshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.
NASA Astrophysics Data System (ADS)
Salehipour, H.; Peltier, W. R.
2014-12-01
In this paper we will describe the results obtained through integration of a further refined version of the truly global barotropic tidal model of Salehipour et al. (Ocean Modell., 69, 2013) using the most recent reconstruction of ice-age bathymetric conditions as embodied in the recently constructed ICE-6G_C (VM5a) model of Peltier et al. (JGR-Solid Earth, in press, 2014). Our interest is in the spatial and temporal evolution of tidal amplitude, phase and dissipation from the Last Glacial Maximum (LGM) 26,000 years ago until the present. The state-of-the-art higher order nonlinear tidal model of Salehipour et al. (2013) includes a highly parallelized multi-scale framework in which an unstructured tessellation of the global ocean enables extensive local refinement around regions of interest such as the Hawaiian Ridge, the Brazil Basin and the Southern Ocean. At LGM, features such as the Patagonian Shelf were fully exposed land which during the deglaciation process would have been flooded leading to significant changes of tidal range along the evolving coastline. In the further development of this model we have included the fully iterated treatment of the influence of gravitational self-attraction and loading as in, e.g. Egbert et al. (JGR-Oceans, 109, 2004). The treatment of the dissipation of the barotropic tide through dissipation of the internal tide has also been significantly improved. Our paleobathymetry and coastline data sets extend from LGM to present at 500 year intervals and constitute a significant refinement of the widely employed ICE-5G (VM2) model of Peltier (Annu. Rev. Earth Planet. Sci., 32, 2004). Our results will be compared with those recently published by Green & Nycander (JPO, 43, 2013) and Wilmes & Green (JGR-Oceans, 119, 2014) as well as with the earlier results of Griffiths & Peltier (GRL, 35, 2008; J. Clim., 22, 2009).
Cascading rainfall uncertainties into 2D inundation impact models
NASA Astrophysics Data System (ADS)
Souvignet, Maxime; de Almeida, Gustavo; Champion, Adrian; Garcia Pintado, Javier; Neal, Jeff; Freer, Jim; Cloke, Hannah; Odoni, Nick; Coxon, Gemma; Bates, Paul; Mason, David
2013-04-01
Existing precipitation products show differences in their spatial and temporal distribution and several studies have presented how these differences influence the ability to predict hydrological responses. However, an atmospheric-hydrologic-hydraulic uncertainty cascade is seldom explored and how, importantly, input uncertainties propagate through this cascade is still poorly understood. Such a project requires a combination of modelling capabilities, runoff generation predictions based on those rainfall forecasts, and hydraulic flood wave propagation based on the runoff predictions. Accounting for uncertainty in each component is important in decision making for issuing flood warnings, monitoring or planning. We suggest a better understanding of uncertainties in inundation impact modelling must consider these differences in rainfall products. This will improve our understanding of the input uncertainties on our predictive capability. In this paper, we propose to address this issue by i) exploring the effects of errors in rainfall on inundation predictive capacity within an uncertainty framework, i.e. testing inundation uncertainty against different comparable meteorological conditions (i.e. using different rainfall products). Our method cascades rainfall uncertainties into a lumped hydrologic model (FUSE) within the GLUE uncertainty framework. The resultant prediction uncertainties in discharge provide uncertain boundary conditions, which are cascaded into a simplified shallow water 2D hydraulic model (LISFLOOD-FP). Rainfall data captured by three different measurement techniques - rain gauges, gridded data and numerical weather predictions (NWP) models are used to assess the combined input data and model parameter uncertainty. The study is performed in the Severn catchment over the period between June and July 2007, where a series of rainfall events causing record floods in the study area). Changes in flood area extent are compared and the uncertainty envelope is
Duality Between Spin Networks and the 2D Ising Model
NASA Astrophysics Data System (ADS)
Bonzom, Valentin; Costantino, Francesco; Livine, Etera R.
2016-06-01
The goal of this paper is to exhibit a deep relation between the partition function of the Ising model on a planar trivalent graph and the generating series of the spin network evaluations on the same graph. We provide respectively a fermionic and a bosonic Gaussian integral formulation for each of these functions and we show that they are the inverse of each other (up to some explicit constants) by exhibiting a supersymmetry relating the two formulations. We investigate three aspects and applications of this duality. First, we propose higher order supersymmetric theories that couple the geometry of the spin networks to the Ising model and for which supersymmetric localization still holds. Secondly, after interpreting the generating function of spin network evaluations as the projection of a coherent state of loop quantum gravity onto the flat connection state, we find the probability distribution induced by that coherent state on the edge spins and study its stationary phase approximation. It is found that the stationary points correspond to the critical values of the couplings of the 2D Ising model, at least for isoradial graphs. Third, we analyze the mapping of the correlations of the Ising model to spin network observables, and describe the phase transition on those observables on the hexagonal lattice. This opens the door to many new possibilities, especially for the study of the coarse-graining and continuum limit of spin networks in the context of quantum gravity.
Effects of Agent's Repulsion in 2d Flocking Models
NASA Astrophysics Data System (ADS)
Moussa, Najem; Tarras, Iliass; Mazroui, M'hammed; Boughaleb, Yahya
In nature many animal groups, such as fish schools or bird flocks, clearly display structural order and appear to move as a single coherent entity. In order to understand the complex behavior of these systems, many models have been proposed and tested so far. This paper deals with an extension of the Vicsek model, by including a second zone of repulsion, where each agent attempts to maintain a minimum distance from the others. The consideration of this zone in our study seems to play an important role during the travel of agents in the two-dimensional (2D) flocking models. Our numerical investigations show that depending on the basic ingredients such as repulsion radius (R1), effect of density of agents (ρ) and noise (η), our nonequilibrium system can undergo a kinetic phase transition from no transport to finite net transport. For different values of ρ, kinetic phase diagrams in the plane (η ,R1) are found. Implications of these findings are discussed.
2-D Model for Normal and Sickle Cell Blood Microcirculation
NASA Astrophysics Data System (ADS)
Tekleab, Yonatan; Harris, Wesley
2011-11-01
Sickle cell disease (SCD) is a genetic disorder that alters the red blood cell (RBC) structure and function such that hemoglobin (Hb) cannot effectively bind and release oxygen. Previous computational models have been designed to study the microcirculation for insight into blood disorders such as SCD. Our novel 2-D computational model represents a fast, time efficient method developed to analyze flow dynamics, O2 diffusion, and cell deformation in the microcirculation. The model uses a finite difference, Crank-Nicholson scheme to compute the flow and O2 concentration, and the level set computational method to advect the RBC membrane on a staggered grid. Several sets of initial and boundary conditions were tested. Simulation data indicate a few parameters to be significant in the perturbation of the blood flow and O2 concentration profiles. Specifically, the Hill coefficient, arterial O2 partial pressure, O2 partial pressure at 50% Hb saturation, and cell membrane stiffness are significant factors. Results were found to be consistent with those of Le Floch [2010] and Secomb [2006].
Ab initio modeling of 2D layered organohalide lead perovskites.
Fraccarollo, Alberto; Cantatore, Valentina; Boschetto, Gabriele; Marchese, Leonardo; Cossi, Maurizio
2016-04-28
A number of 2D layered perovskites A2PbI4 and BPbI4, with A and B mono- and divalent ammonium and imidazolium cations, have been modeled with different theoretical methods. The periodic structures have been optimized (both in monoclinic and in triclinic systems, corresponding to eclipsed and staggered arrangements of the inorganic layers) at the DFT level, with hybrid functionals, Gaussian-type orbitals and dispersion energy corrections. With the same methods, the various contributions to the solid stabilization energy have been discussed, separating electrostatic and dispersion energies, organic-organic intralayer interactions and H-bonding effects, when applicable. Then the electronic band gaps have been computed with plane waves, at the DFT level with scalar and full relativistic potentials, and including the correlation energy through the GW approximation. Spin orbit coupling and GW effects have been combined in an additive scheme, validated by comparing the computed gap with well known experimental and theoretical results for a model system. Finally, various contributions to the computed band gaps have been discussed on some of the studied systems, by varying some geometrical parameters and by substituting one cation in another's place. PMID:27131557
2-D Inhomogeneous Modeling of the Solar CO Bands
NASA Astrophysics Data System (ADS)
Ayres, T. R.
1996-05-01
The recent discovery of off-limb emissions in the mid-IR ( ~ 5 mu m) vibration-rotation bands of solar carbon monoxide (CO) has sparked new interest in the formation of the molecular lines, and their ability to diagnose thermal conditions at high altitudes. The off-limb extensions of the strong CO lines indicate the penetration of cool material (T ~ 3500 K) several hundred kilometers into the otherwise hot (T ~ 6000 K) chromosphere. The origin of the cool gas, and its role in the thermal energy balance, remain controversial. The interpretation of the CO observations must rely heavily upon numerical modeling, in particular highly-inhomogeneous thermal structures arrayed in a 2-D scheme that can properly treat the geometry of the grazing rays at the solar limb. The radiation transport, itself, is especially simple for the CO off-limb emissions, because the fundamental bands form quite close to LTE (high collision rates; low spontaneous decay rates) and the background continuum is purely thermal as well (f--f transitions in H(-) and H). Thus, the geometrical aspects of the problem can be treated in considerably more detail than would be practical for typical NLTE scattering lines. I describe the recent modeling efforts, and the diagnostic potential of the CO bands for future observational studies of inhomogeneous surface structure on the Sun, and on other stars of late spectral type. This work was supported by NSF grant AST-9218063 to the University of Colorado.
Ab initio modeling of 2D layered organohalide lead perovskites
NASA Astrophysics Data System (ADS)
Fraccarollo, Alberto; Cantatore, Valentina; Boschetto, Gabriele; Marchese, Leonardo; Cossi, Maurizio
2016-04-01
A number of 2D layered perovskites A2PbI4 and BPbI4, with A and B mono- and divalent ammonium and imidazolium cations, have been modeled with different theoretical methods. The periodic structures have been optimized (both in monoclinic and in triclinic systems, corresponding to eclipsed and staggered arrangements of the inorganic layers) at the DFT level, with hybrid functionals, Gaussian-type orbitals and dispersion energy corrections. With the same methods, the various contributions to the solid stabilization energy have been discussed, separating electrostatic and dispersion energies, organic-organic intralayer interactions and H-bonding effects, when applicable. Then the electronic band gaps have been computed with plane waves, at the DFT level with scalar and full relativistic potentials, and including the correlation energy through the GW approximation. Spin orbit coupling and GW effects have been combined in an additive scheme, validated by comparing the computed gap with well known experimental and theoretical results for a model system. Finally, various contributions to the computed band gaps have been discussed on some of the studied systems, by varying some geometrical parameters and by substituting one cation in another's place.
NASA Astrophysics Data System (ADS)
Calafat, Francisco M.; Avgoustoglou, Euripides; Jordà, Gabriel; Flocas, Helena; Tsimplis, Mikis N.
2014-05-01
Storm surges are responsible for considerable damage to coastal property and loss of life every year. Furthermore, the damage from its combined effect with gradual sea level rise constitutes one of the most visible and costly effects of climate change. Understanding the mechanisms of storm surges as well as improving the predictive skill of storm surge models are hence of primary importance. This study investigates the ability of a 2D hydrodynamic model to simulate sea level extremes of meteorological origin in the Mediterranean Sea. For this purpose, the output of the model is compared to hourly sea level observations from a series of coastal tide gauge records covering most regions of the Mediterranean. It is found that the model underestimates the size of the sea level extremes at most stations. In terms of percentage, the underestimation is independent of the value of the extreme and, at some stations like Marseille, can be as large as 65%. Interestingly, in the Adriatic Sea, and more in particular at Trieste, the model can in some cases overestimate the magnitude of the extremes or even generate significant extreme events when they are in fact not present in the tide gauge record. The extent to which higher resolution sea level pressure (SLP) and wind forcing improve the sea level predictions is explored by comparing various simulations with atmospheric forcing of different resolution. Finally, the effect of the rapid deepening during explosive cyclone events on the performance of the ocean model as well as possible departures from the local inverse barometer effect during such events is also investigated. For the study of the events associated with explosive cyclones we use SLP and wind data with a horizontal resolution of about 7 km as simulated by the COSMO.GR model for 23 cases of explosive cyclogenesis spanning the period 2002-2010.
Barotropic instability with divergence - Theory and applications to Venus
NASA Technical Reports Server (NTRS)
Dobrovolskis, Anthony R.; Diner, David J.
1990-01-01
IR images of Venus reveal a curious double-lobed hot spot in the polar region. Elson (1982) has suggested that this dipole represents a barotropic instability associated with a high-latitude jet. Unfortunately, the classical theory of barotropic instability cannot predict temperature variations. This paper generalizes the theory to include horizontal divergence, vertical motions, and temperature variations, and applies it to the stratosphere of Venus. The fastest-growing barotropic instability in the nominal model matches the observed dipole in period and horizontal temperature pattern. The accompanying wind variations are comparable to the speed of the mean jet, indicating strong nonlinear effects. It is concluded that the Venus dipole may represent the self-limited stage of a barotropic instability with divergence.
A 2D electrohydrodynamic model for electrorotation of fluid drops.
Feng, James Q
2002-02-01
A theoretical analysis of spontaneous electrorotation of deformable fluid drops in a DC electric field is presented with a 2D electrohydrodynamic model. The fluids in the system are assumed to be leaky dielectric and Newtonian. If the rotating flow is dominant over the cellular convection type of electrohydrodynamic flow, closed-form solutions for drops of small deformations can be obtained. Because the governing equations are in general nonlinear even when drop deformations are ignored, the general solution for even undeformed drop takes a form of infinite series and can only be evaluated by numerical means. Both closed-form solutions for special cases and numerical solutions for more general cases are obtained here to describe steady-state field variables and first-order drop deformations. In a DC electric field of strength beyond the threshold value, spontaneous electrorotation of a drop is shown to occur when charge relaxation in the surrounding fluid is faster than the fluid inside the drop. With increasing the strength of the applied electric field from the threshold for onset of electrorotation, the axis of drop contraction deviates from from that of the applied electric field in the direction of the rotating flow with an angle increasing with the field strength. PMID:16290391
2D modeling of electromagnetic waves in cold plasmas
Crombé, K.; Van Eester, D.; Koch, R.; Kyrytsya, V.
2014-02-12
The consequences of sheath (rectified) electric fields, resulting from the different mobility of electrons and ions as a response to radio frequency (RF) fields, are a concern for RF antenna design as it can cause damage to antenna parts, limiters and other in-vessel components. As a first step to a more complete description, the usual cold plasma dielectric description has been adopted, and the density profile was assumed to be known as input. Ultimately, the relevant equations describing the wave-particle interaction both on the fast and slow timescale will need to be tackled but prior to doing so was felt as a necessity to get a feeling of the wave dynamics involved. Maxwell's equations are solved for a cold plasma in a 2D antenna box with strongly varying density profiles crossing also lower hybrid and ion-ion hybrid resonance layers. Numerical modelling quickly becomes demanding on computer power, since a fine grid spacing is required to capture the small wavelengths effects of strongly evanescent modes.
Numerical modeling of seismogram envelopes in 2-D random media
NASA Astrophysics Data System (ADS)
Fehler, Michael
2002-11-01
Several portions of seismograms recorded from regional earthquakes cannot be easily explained as resulting from waves propagating along deterministic paths within the Earth. For example, seismic coda, which is the tail portion of the seismogram of an earthquake recorded at distances of less than 100 km, is considered as resulting from waves that are multiply scattered from random heterogeneities in the Earth's lithosphere. At greater distances, observations that the duration of the initial arriving wave packet is much longer than the source-time duration is explained as being due to multiple forward scattering along the path between the source and the receiver. To investigate these phenomena, we use a finite difference method to numerically simulate 2-D scalar-waves that propagate through random media characterized by a von Karman autocorrelation function. Such media are considered to be appropriate models for the random component of the structure of the Earth's lithosphere. We investigate the characteristics of the resulting wavefields and compare them with those of observed seismograms.
VAM2D: Variably saturated analysis model in two dimensions
Huyakorn, P.S.; Kool, J.B.; Wu, Y.S. )
1991-10-01
This report documents a two-dimensional finite element model, VAM2D, developed to simulate water flow and solute transport in variably saturated porous media. Both flow and transport simulation can be handled concurrently or sequentially. The formulation of the governing equations and the numerical procedures used in the code are presented. The flow equation is approximated using the Galerkin finite element method. Nonlinear soil moisture characteristics and atmospheric boundary conditions (e.g., infiltration, evaporation and seepage face), are treated using Picard and Newton-Raphson iterations. Hysteresis effects and anisotropy in the unsaturated hydraulic conductivity can be taken into account if needed. The contaminant transport simulation can account for advection, hydrodynamic dispersion, linear equilibrium sorption, and first-order degradation. Transport of a single component or a multi-component decay chain can be handled. The transport equation is approximated using an upstream weighted residual method. Several test problems are presented to verify the code and demonstrate its utility. These problems range from simple one-dimensional to complex two-dimensional and axisymmetric problems. This document has been produced as a user's manual. It contains detailed information on the code structure along with instructions for input data preparation and sample input and printed output for selected test problems. Also included are instructions for job set up and restarting procedures. 44 refs., 54 figs., 24 tabs.
NASA Astrophysics Data System (ADS)
Calafat, F. M.; Avgoustoglou, E.; Jordà, G.; Flocas, H.; Zodiatis, G.; Tsimplis, M. N.; Kouroutzoglou, J.
2014-11-01
Storm surges are responsible for great damage to coastal property and loss of life every year. Coastal management and adaptation practices are essential to reduce such damage. Numerical models provide a useful tool for informing these practices as they simulate sea level with high spatial resolution. Here we investigate the ability of a barotropic version of the HAMSOM model to simulate sea level extremes of meteorological origin in the Mediterranean Sea, including those caused by explosive cyclones. For this purpose, the output of the model is compared to hourly sea level observations from six tide gauge records (Valencia, Barcelona, Marseille, Civitavecchia, Trieste, and Antalya). It is found that the model underestimates the positive extremes significantly at all stations, in some cases by up to 65%. At Trieste, the model can also sometimes overestimate the extremes significantly. The differences between the model and the residuals are not constant for extremes of a given height, which limits the applicability of the numerical model for storm surge forecasting because calibration is difficult. The 50 and 10 year return levels are reasonably well captured by the model at all stations except Barcelona and Marseille, where they are underestimated by over 30%. The number of exceedances of the 99.9th and 99.95% percentiles over a period of 25 years is severely underestimated by the model at all stations. The skill of the model for predicting the timing and value of the storm surges seems to be higher for the events associated with explosive cyclones at all stations.
Predicting Fracture Using 2D Finite Element Modeling
MacNeil, J.A.M.; Adachi, J.D; Goltzman, D; Josse, R.G; Kovacs, C.S; Prior, J.C; Olszynski, W; Davison, K.S.; Kaiser, S.M
2013-01-01
A decrease in bone density at the hip or spine has been shown to increase the risk of fracture. A limitation of the bone mineral density (BMD) measurement is that it provides only a measure of a bone samples average density when projected onto a 2D surface. Effectively, what determines bone fracture is whether an applied load exceeds ultimate strength, with both bone tissue material properties (can be approximated through bone density), and geometry playing a role. The goal of this project was to use bone geometry and BMD obtained from radiographs and DXA measurements respectively to estimate fracture risk, using a two-dimensional finite element model (FEM) of the sagittal plane of lumbar vertebrae. The Canadian Multicenter Osteoporosis Study (CaMos) data was used for this study. There were 4194 men and women over the age of 50 years, with 786 having fractures. Each subject had BMD testing and radiographs of their lumbar vertebrae. A single two dimensional FEM of the first to fourth lumbar vertebra was automatically generated for each subject. Bone tissue stiffness was assigned based on the BMD of the individual vertebrae, and adjusted for patient age. Axial compression boundary conditions were applied with a force proportional to body mass. The resulting overall strain from the applied force was found. Men and women were analyzed separately. At baseline, the sensitivity of BMD to predict fragility fractures in women and men was 3.77 % and 0.86 %, while the sensitivity of FEM to predict fragility fractures for women and men was 10.8 % and 11.3 %. The FEM ROC curve demonstrated better performance compared to BMD. The relative risk of being considered at high fracture risk using FEM at baseline, was a better predictor of 5 year incident fragility fracture risk compared to BMD. PMID:21959170
A 2-D modeling contribution to river training design
NASA Astrophysics Data System (ADS)
Anselmo, V.; Coccato, M.; Frank, E.; Guiot, E.
2003-04-01
In the last ten years, two major floods (1994 and 2000) occurred in North-western Italy and a few questions arose about the hydraulic behavior of the streams as well about the choice and design of protection works. The River Po Authority is oriented to assign "design flows" in selected cross sections of the main rivers, as a design constraint to land management and river training in the upstream areas. Since the region has been fully developed in the last century and somewhere it is overcrowded, space for spreading flood flows is strongly reduced, while large partially developed areas are prone to flooding and residents ask for being protected. A first question regards the contribution to flood peak reduction of the still existing flood prone undeveloped areas beside the main channels, and a second question is about the best way to improve such a behavior. A 2-D unsteady model (Sobek, originated by Delft Hydraulics) was applied to a 25 km reach of the upper River Po. The effects of major floods was investigated, proving that the reduction of the peak flow is minor mainly because of the rather high slope (0.0015) and of the flood volume (500·106 m3). Aiming to enhance the role of the flooded areas, a few types of river training schemes were checked, with particular attention to the so called "Po system". Depth and extension of compartments are the main variables. Results are interesting, but must be evaluated in front of the cost-benefit analysis. The investigation is being extended to more steep stream reaches (up to 0.01), which are representative of the main upper Po tributaries.
A 2D simulation model for urban flood management
NASA Astrophysics Data System (ADS)
Price, Roland; van der Wielen, Jonathan; Velickov, Slavco; Galvao, Diogo
2014-05-01
The European Floods Directive, which came into force on 26 November 2007, requires member states to assess all their water courses and coast lines for risk of flooding, to map flood extents and assets and humans at risk, and to take adequate and coordinated measures to reduce the flood risk in consultation with the public. Flood Risk Management Plans are to be in place by 2015. There are a number of reasons for the promotion of this Directive, not least because there has been much urban and other infrastructural development in flood plains, which puts many at risk of flooding along with vital societal assets. In addition there is growing awareness that the changing climate appears to be inducing more frequent extremes of rainfall with a consequent increases in the frequency of flooding. Thirdly, the growing urban populations in Europe, and especially in the developing countries, means that more people are being put at risk from a greater frequency of urban flooding in particular. There are urgent needs therefore to assess flood risk accurately and consistently, to reduce this risk where it is important to do so or where the benefit is greater than the damage cost, to improve flood forecasting and warning, to provide where necessary (and possible) flood insurance cover, and to involve all stakeholders in decision making affecting flood protection and flood risk management plans. Key data for assessing risk are water levels achieved or forecasted during a flood. Such levels should of course be monitored, but they also need to be predicted, whether for design or simulation. A 2D simulation model (PriceXD) solving the shallow water wave equations is presented specifically for determining flood risk, assessing flood defense schemes and generating flood forecasts and warnings. The simulation model is required to have a number of important properties: -Solve the full shallow water wave equations using a range of possible solutions; -Automatically adjust the time step and
A 2D simulation model for urban flood management
NASA Astrophysics Data System (ADS)
Price, Roland; van der Wielen, Jonathan; Velickov, Slavco; Galvao, Diogo
2014-05-01
The European Floods Directive, which came into force on 26 November 2007, requires member states to assess all their water courses and coast lines for risk of flooding, to map flood extents and assets and humans at risk, and to take adequate and coordinated measures to reduce the flood risk in consultation with the public. Flood Risk Management Plans are to be in place by 2015. There are a number of reasons for the promotion of this Directive, not least because there has been much urban and other infrastructural development in flood plains, which puts many at risk of flooding along with vital societal assets. In addition there is growing awareness that the changing climate appears to be inducing more frequent extremes of rainfall with a consequent increases in the frequency of flooding. Thirdly, the growing urban populations in Europe, and especially in the developing countries, means that more people are being put at risk from a greater frequency of urban flooding in particular. There are urgent needs therefore to assess flood risk accurately and consistently, to reduce this risk where it is important to do so or where the benefit is greater than the damage cost, to improve flood forecasting and warning, to provide where necessary (and possible) flood insurance cover, and to involve all stakeholders in decision making affecting flood protection and flood risk management plans. Key data for assessing risk are water levels achieved or forecasted during a flood. Such levels should of course be monitored, but they also need to be predicted, whether for design or simulation. A 2D simulation model (PriceXD) solving the shallow water wave equations is presented specifically for determining flood risk, assessing flood defense schemes and generating flood forecasts and warnings. The simulation model is required to have a number of important properties: -Solve the full shallow water wave equations using a range of possible solutions; -Automatically adjust the time step and
The Implementation of C-ID, R2D2 Model on Learning Reading Comprehension
ERIC Educational Resources Information Center
Rayanto, Yudi Hari; Rusmawan, Putu Ngurah
2016-01-01
The purposes of this research are to find out, (1) whether C-ID, R2D2 model is effective to be implemented on learning Reading comprehension, (2) college students' activity during the implementation of C-ID, R2D2 model on learning Reading comprehension, and 3) college students' learning achievement during the implementation of C-ID, R2D2 model on…
NASA Astrophysics Data System (ADS)
Chae, Dongho; Constantin, Peter; Wu, Jiahong
2014-09-01
We give an example of a well posed, finite energy, 2D incompressible active scalar equation with the same scaling as the surface quasi-geostrophic equation and prove that it can produce finite time singularities. In spite of its simplicity, this seems to be the first such example. Further, we construct explicit solutions of the 2D Boussinesq equations whose gradients grow exponentially in time for all time. In addition, we introduce a variant of the 2D Boussinesq equations which is perhaps a more faithful companion of the 3D axisymmetric Euler equations than the usual 2D Boussinesq equations.
NASA Astrophysics Data System (ADS)
Ding, Yang; Bao, Xianwen; Yu, Huaming; Kuang, Liang
2012-04-01
The tides and tidal energetics in the Indonesian seas are simulated using a three-dimensional finite volume coastal ocean model. The high-resolution coastline-fitted model is configured to better resolve the hydrodynamic processes around the numerous barrier islands. A large model domain is adopted to minimize the uncertainty adjacent to open boundaries. The model results with elevation assimilation based on a simple nudge scheme faithfully reproduced the general features of the barotropic tides in the Indonesian Seas. The mean root-mean-square errors between the observed and simulated tidal constants are 2.3, 1.1, 2.4, and 1.5 cm for M2, S2, K1, and O1, respectively. Analysis of the model solutions indicates that the semidiurnal tides in the Indonesian Seas are primarily dominated by the Indian Ocean, whereas the diurnal tides in this region are mainly influenced by the Pacific Ocean, which is consistent with previous studies. Examinations of tidal energy transport reveal that the tidal energy for both of the simulated tidal constituents are transported from the Indian Ocean into the IS mainly through the Lombok Strait and the Timor Sea, whereas only M2 energy enters the Banda Sea and continues northward. The tidal energy dissipates the most in the passages on both sides of Timor Island, with the maximum M2 and K1 tidal energy transport reaching about 750 and 650 kW m-1, respectively. The total energy losses of the four dominant constituents in the IS are nearly 338 GW, with the M2 constituent dissipating 240.8 GW. It is also shown that the bottom dissipation rate for the M2 tide is about 1-2 order of magnitudes larger than that of the other three tidal components in the Indonesian seas.
2-D model of the streamer zone of a leader
NASA Astrophysics Data System (ADS)
Milikh, G. M.; Likhanskii, A. V.; Shneider, M. N.; Raina, A.; George, A.
2016-02-01
Formation of the streamer zone of a leader is an outstanding problem in the physics of electric discharges which is relevant to laboratory leaders, as well as to the leaders formed by lightning. Despite substantial progress in the theoretical understanding of this complicated phenomenon, significant puzzles, such as the low propagation velocity of a leader compared to the fast streamers, remain. The objective of this paper is to present 2-D plasma simulations of the formation and propagation of the streamer zone of a leader. In these simulations we will generate a group of streamers that propagate in a discharge gap while interacting with each other. It is shown that interaction between the streamers significantly reduces their propagation velocity. This explains why the leader, which consists of many streamers, is much slower than a single streamer formed in the same discharge gap. It is shown that the mean velocity suppression of the group of streamers is determined by the inter-streamer distance. The critical value of the packing factor of the streamers at which the interactions between them can be neglected, and thus the discussed process can be treated as caused by a single streamer, is obtained.
2D quantum double models from a 3D perspective
NASA Astrophysics Data System (ADS)
Bernabé Ferreira, Miguel Jorge; Padmanabhan, Pramod; Teotonio-Sobrinho, Paulo
2014-09-01
In this paper we look at three dimensional (3D) lattice models that are generalizations of the state sum model used to define the Kuperberg invariant of 3-manifolds. The partition function is a scalar constructed as a tensor network where the building blocks are tensors given by the structure constants of an involutory Hopf algebra A. These models are very general and are hard to solve in its entire parameter space. One can obtain familiar models, such as ordinary gauge theories, by letting A be the group algebra {C}(G) of a discrete group G and staying on a certain region of the parameter space. We consider the transfer matrix of the model and show that quantum double Hamiltonians are derived from a particular choice of the parameters. Such a construction naturally leads to the star and plaquette operators of the quantum double Hamiltonians, of which the toric code is a special case when A={C}({{{Z}}_{2}}). This formulation is convenient to study ground states of these generalized quantum double models where they can naturally be interpreted as tensor network states. For a surface Σ, the ground state degeneracy is determined by the Kuperberg 3-manifold invariant of \\Sigma \\times {{S}^{1}}. It is also possible to obtain extra models by simply enlarging the allowed parameter space but keeping the solubility of the model. While some of these extra models have appeared before in the literature, our 3D perspective allows for an uniform description of them.
NASA Astrophysics Data System (ADS)
Feng, Jie; Ding, Ruiqiang; Li, Jianping; Liu, Deqiang
2016-09-01
The breeding method has been widely used to generate ensemble perturbations in ensemble forecasting due to its simple concept and low computational cost. This method produces the fastest growing perturbation modes to catch the growing components in analysis errors. However, the bred vectors (BVs) are evolved on the same dynamical flow, which may increase the dependence of perturbations. In contrast, the nonlinear local Lyapunov vector (NLLV) scheme generates flow-dependent perturbations as in the breeding method, but regularly conducts the Gram-Schmidt reorthonormalization processes on the perturbations. The resulting NLLVs span the fast-growing perturbation subspace efficiently, and thus may grasp more components in analysis errors than the BVs. In this paper, the NLLVs are employed to generate initial ensemble perturbations in a barotropic quasi-geostrophic model. The performances of the ensemble forecasts of the NLLV method are systematically compared to those of the random perturbation (RP) technique, and the BV method, as well as its improved version—the ensemble transform Kalman filter (ETKF) method. The results demonstrate that the RP technique has the worst performance in ensemble forecasts, which indicates the importance of a flow-dependent initialization scheme. The ensemble perturbation subspaces of the NLLV and ETKF methods are preliminarily shown to catch similar components of analysis errors, which exceed that of the BVs. However, the NLLV scheme demonstrates slightly higher ensemble forecast skill than the ETKF scheme. In addition, the NLLV scheme involves a significantly simpler algorithm and less computation time than the ETKF method, and both demonstrate better ensemble forecast skill than the BV scheme.
A fully coupled 2D model of equiaxed eutectic solidification
Charbon, Ch.; LeSar, R.
1995-12-31
We propose a model of equiaxed eutectic solidification that couples the macroscopic level of heat diffusion with the microscopic level of nucleation and growth of the eutectic grains. The heat equation with the source term corresponding to the latent heat release due to solidification is calculated numerically by means of an implicit finite difference method. In the time stepping scheme, the evolution of solid fraction is deduced from a stochastic model of nucleation and growth which uses the local temperature (interpolated from the FDM mesh) to determine the local grain density and the local growth rate. The solid-liquid interface of each grain is tracked by using a subdivision of each grain perimeter in a large number of sectors. The state of each sector (i.e. whether it is still in contact with the liquid or already captured by an other grain) and the increase of radius of each grain during one time step allows one to compute the increase of solid fraction. As for deterministic models, the results of the model are the evolution of temperature and of solid fraction at any point of the sample. Moreover the model provides a complete picture of the microstructure, thus not limiting the microstructural information to the average grain density but allowing one to compute any stereological value of interest. We apply the model to the solidification of gray cast iron.
Improvement of a 2D numerical model of lava flows
NASA Astrophysics Data System (ADS)
Ishimine, Y.
2013-12-01
I propose an improved procedure that reduces an improper dependence of lava flow directions on the orientation of Digital Elevation Model (DEM) in two-dimensional simulations based on Ishihara et al. (in Lava Flows and Domes, Fink, JH eds., 1990). The numerical model for lava flow simulations proposed by Ishihara et al. (1990) is based on two-dimensional shallow water model combined with a constitutive equation for a Bingham fluid. It is simple but useful because it properly reproduces distributions of actual lava flows. Thus, it has been regarded as one of pioneer work of numerical simulations of lava flows and it is still now widely used in practical hazard prediction map for civil defense officials in Japan. However, the model include an improper dependence of lava flow directions on the orientation of DEM because the model separately assigns the condition for the lava flow to stop due to yield stress for each of two orthogonal axes of rectangular calculating grid based on DEM. This procedure brings a diamond-shaped distribution as shown in Fig. 1 when calculating a lava flow supplied from a point source on a virtual flat plane although the distribution should be circle-shaped. To improve the drawback, I proposed a modified procedure that uses the absolute value of yield stress derived from both components of two orthogonal directions of the slope steepness to assign the condition for lava flows to stop. This brings a better result as shown in Fig. 2. Fig. 1. (a) Contour plots calculated with the original model of Ishihara et al. (1990). (b) Contour plots calculated with a proposed model.
Stratosphere chemistry in a 2-D model with residual circulation
NASA Technical Reports Server (NTRS)
Guthrie, Paul D.; Jackman, Charles H.
1990-01-01
The objective of this research was to examine the effects of chemical perturbations on the stratosphere using models which can incorporate fully interactive radiative, chemical, and dynamical responses, in the context of a zonally averaged model. Model runs for the unperturbed, chlorine-perturbed and simultaneously chlorine-and CO2-perturbed cases were completed using the JPL-87 chemical kinetics data. The base case was analyzed and submitted for publication. The perturbed cases show substantial sensitivity of the predicted column ozone depletion to the perturbations affecting lower stratosphere temperature, but less to far dynamical perturbations. The column ozone distribution changed substantially when the kinetics data was changed. This implies a greater-than-expected uncertainty in predicted latitude distributions of ozone depletion, due to uncertainty about the accuracy and completeness of the chemical kinetics data set.
Development of CCHE2D embankment break model
Technology Transfer Automated Retrieval System (TEKTRAN)
Earthen embankment breach often results in detrimental impact on downstream residents and infrastructure, especially those located in the flooding zone. Embankment failures are most commonly caused by overtopping or internal erosion. This study is to develop a practical numerical model for simulat...
Phase Structure of the Random Zq Models in 2D
NASA Astrophysics Data System (ADS)
Sasamoto, T.; Nishimori, H.
We discuss the phase diagram of the random Zq models in two dimensions. It is argued that, when q is large enough, there exist three phases in the phase diagram with two axes being the temperature and the strength of randomness. Our conlusions are derived based on the application of the duality arguments for random systems, which have been formulated recently by Maillard et al.
Simulation of subgrid orographic precipitation with an embedded 2-D cloud-resolving model
NASA Astrophysics Data System (ADS)
Jung, Joon-Hee; Arakawa, Akio
2016-03-01
By explicitly resolving cloud-scale processes with embedded two-dimensional (2-D) cloud-resolving models (CRMs), superparameterized global atmospheric models have successfully simulated various atmospheric events over a wide range of time scales. Up to now, however, such models have not included the effects of topography on the CRM grid scale. We have used both 3-D and 2-D CRMs to simulate the effects of topography with prescribed "large-scale" winds. The 3-D CRM is used as a benchmark. The results show that the mean precipitation can be simulated reasonably well by using a 2-D representation of topography as long as the statistics of the topography such as the mean and standard deviation are closely represented. It is also shown that the use of a set of two perpendicular 2-D grids can significantly reduce the error due to a 2-D representation of topography.
Analytic Differentiation of Barlat's 2D Criteria for Inverse Modeling
Endelt, Benny; Nielsen, Karl Brian; Danckert, Joachim
2005-08-05
The demand for alternative identification schemes for identification of constitutive parameters is getting more pronounced as the complexity of the constitutive equations increases, i.e. the number of parameters subject to identification. A general framework for inverse identification of constitutive parameters associated with sheet metal forming is proposed in the article. The inverse problem is solved, through minimization of the least square error between an experimental punch force sampled from a deep drawing and a predicted punch force produced from a coherent finite element model.
NASA Astrophysics Data System (ADS)
Li, Tingting; Li, Xiaofan
2016-06-01
The budget of perturbation kinetic energy is calculated to study barotropic and baroclinic processes associated with convective development in the tropical deep convective regime. The barotropic kinetic-energy conversion is further examined through the derivation and analysis of budget of vertical transport of zonal momentum. The two-dimensional cloud-resolving model simulation data during TOGA COARE are analyzed in this budget study. The analysis of the energy budget shows that while baroclinic conversion from perturbation available potential energy is a major source for the growth of perturbation kinetic energy, barotropic conversion from the mean kinetic energy has significant modification on the variation of perturbation kinetic energy. The general evolution of barotropic conversion is associated with that of the vertical wind shear of the mean circulations, but the significant change in barotropic conversion corresponds to the change in the vertical transport of zonal momentum. Zonal transport of cloud hydrometeor is responsible for the variation of vertical transport of zonal momentum.
Google Earth as a tool in 2-D hydrodynamic modeling
NASA Astrophysics Data System (ADS)
Chien, Nguyen Quang; Keat Tan, Soon
2011-01-01
A method for coupling virtual globes with geophysical hydrodynamic models is presented. Virtual globes such as Google TM Earth can be used as a visualization tool to help users create and enter input data. The authors discuss techniques for representing linear and areal geographical objects with KML (Keyhole Markup Language) files generated using computer codes (scripts). Although virtual globes offer very limited tools for data input, some data of categorical or vector type can be entered by users, and then transformed into inputs for the hydrodynamic program by using appropriate scripts. An application with the AnuGA hydrodynamic model was used as an illustration of the method. Firstly, users draw polygons on the Google Earth screen. These features are then saved in a KML file which is read using a script file written in the Lua programming language. After the hydrodynamic simulation has been performed, another script file is used to convert the resulting output text file to a KML file for visualization, where the depths of inundation are represented by the color of discrete point icons. The visualization of a wind speed vector field was also included as a supplementary example.
Conservation laws and LETKF with 2D Shallow Water Model
NASA Astrophysics Data System (ADS)
Zeng, Yuefei; Janjic, Tijana
2016-04-01
Numerous approaches have been proposed to maintain physical conservation laws in the numerical weather prediction models. However, to achieve a reliable prediction, adequate initial conditions are also necessary, which are produced by a data assimilation algorithm. If an ensemble Kalman filters (EnKF) is used for this purpose, it has been shown that it could yield unphysical analysis ensemble that for example violates principles of mass conservation and positivity preservation (e.g. Janjic et al 2014) . In this presentation, we discuss the selection of conservation criteria for the analysis step, and start with testing the conservation of mass, energy and enstrophy. The simple experiments deal with nonlinear shallow water equations and simulated observations that are assimilated with LETKF (Localized Ensemble Transform Kalman Filter, Hunt et al. 2007). The model is discretized in a specific way to conserve mass, angular momentum, energy and enstrophy. The effects of the data assimilation on the conserved quantities (of mass, energy and enstrophy) depend on observation covarage, localization radius, observed variable and observation operator. Having in mind that Arakawa (1966) and Arakawa and Lamb (1977) showed that the conservation of both kinetic energy and enstrophy by momentum advection schemes in the case of nondivergent flow prevents systematic and unrealistic energy cascade towards high wave numbers, a cause of excessive numerical noise and possible eventual nonlinear instability, we test the effects on prediction depending on the type of errors in the initial condition. The performance with respect to nonlinear energy cascade is assessed as well.
A 2D model to design MHD induction pumps
NASA Astrophysics Data System (ADS)
Stieglitz, R.; Zeininger, J.
2006-09-01
Technical liquid metal systems accompanied by a thermal transfer of energy such as reactor systems, metallurgical processes, metal refinement, casting, etc., require a forced convection of the fluid. The increased temperatures and more often the environmental conditions as, e.g., in a nuclear environment, pumping principles are required, in which rotating parts are absent. Additionally, in many applications a controlled atmosphere is indispensable, in order to ensure the structural integrity of the duct walls. An interesting option to overcome the sealing problem of a mechanical pump towards the surrounding is offered by induction systems. Although their efficiency compared to that of turbo machines is quite low, they have several advantages, which are attractive to the specific requirements in liquid metal applications such as: - low maintenance costs due to the absence of sealings, bearings and moving parts; - low degradation rate of the structural material; - simple replacement of the inductor without cut of the piping system; - fine regulation of flow rate by different inductor connections; - change of pump characteristics without change of the mechanical set-up. Within the article, general design requirements of electromagnetic pumps (EMP) are elaborated. The design of two annular linear induction pumps operating with sodium and lead-bismuth are presented and the calculated pump characteristics and experimentally obtained data are compared. In this context, physical effects leading to deviations between the model and the real data are addressed. Finally, the main results are summarized. Tables 4, Figs 4, Refs 12.
Bond Order Correlations in the 2D Hubbard Model
NASA Astrophysics Data System (ADS)
Moore, Conrad; Abu Asal, Sameer; Yang, Shuxiang; Moreno, Juana; Jarrell, Mark
We use the dynamical cluster approximation to study the bond correlations in the Hubbard model with next nearest neighbor (nnn) hopping to explore the region of the phase diagram where the Fermi liquid phase is separated from the pseudogap phase by the Lifshitz line at zero temperature. We implement the Hirsch-Fye cluster solver that has the advantage of providing direct access to the computation of the bond operators via the decoupling field. In the pseudogap phase, the parallel bond order susceptibility is shown to persist at zero temperature while it vanishes for the Fermi liquid phase which allows the shape of the Lifshitz line to be mapped as a function of filling and nnn hopping. Our cluster solver implements NVIDIA's CUDA language to accelerate the linear algebra of the Quantum Monte Carlo to help alleviate the sign problem by allowing for more Monte Carlo updates to be performed in a reasonable amount of computation time. Work supported by the NSF EPSCoR Cooperative Agreement No. EPS-1003897 with additional support from the Louisiana Board of Regents.
Barotropic and baroclinic instability of Rossby waves on the infinite beta-plane
NASA Astrophysics Data System (ADS)
Anderson, Jeffrey L.
The normal mode instability problem for a stationary Rossby wave of finite amplitude and arbitrary orientation in a two-level model with zonal mean shear in the vertical is examined on the infinite beta-plane using a Floquet technique. The results are compared to the instability of a finite amplitude stationary Rossby wave in a barotropic model. The instability of the barotropic model linearized about the upper level flow of the stationary two-level Rossby wave is also examined; forcing must be introduced to the barotropic equations in this case to balance the time tendency of the basic state. The instabilities of the three models are compared in an attempt to address two questions. First, are the barotropic instabilities related to any of the two-level model instabilities? Second, what are the consequences of introducing forcing in the stability problem? It is found that the barotropic models accurately reproduce those instabilities of the two-level model that have a nearly barotropic structure; these are not necessarily low-frequency instabilities. The introduction of forcing is found to have a significant impact on the barotropic instabilities. Special emphasis is placed on parameter ranges judged to be analogous to the mid-latitude atmosphere. Possible relevance of the results to geophysical flows and the instability problem on the sphere are discussed.
Global barotropic response to a tropical forcing
NASA Technical Reports Server (NTRS)
Zhang, Yuxia; Mcguirk, James P.
1993-01-01
Zonally varying flow has been used to initialize numerical models and has been shown to play an important role in strong localized responses both in extratropics and the tropics. In this study, a climatological 200 mb January mean is used as a steady basic state of a barotropical model which consists of shallow water equations and a mass source centered at 4 deg S/120 deg E to simulate convective heating over Indonesia region. In the experiment, tropical responses appear not only over the western Pacific, where the forcing is located, but also over the eastern Pacific where the response is related to the zonally varying basic state. The westward propagating equatorial Rossby waves excited by the forcing interact with waves out of and into the tropics and the positive and negative phase of the Rossby waves result in blocking circulation over North America and tropical plumes over equatorial eastern Pacific, respectively.
GEO2D - Two-Dimensional Computer Model of a Ground Source Heat Pump System
James Menart
2013-06-07
This file contains a zipped file that contains many files required to run GEO2D. GEO2D is a computer code for simulating ground source heat pump (GSHP) systems in two-dimensions. GEO2D performs a detailed finite difference simulation of the heat transfer occurring within the working fluid, the tube wall, the grout, and the ground. Both horizontal and vertical wells can be simulated with this program, but it should be noted that the vertical wall is modeled as a single tube. This program also models the heat pump in conjunction with the heat transfer occurring. GEO2D simulates the heat pump and ground loop as a system. Many results are produced by GEO2D as a function of time and position, such as heat transfer rates, temperatures and heat pump performance. On top of this information from an economic comparison between the geothermal system simulated and a comparable air heat pump systems or a comparable gas, oil or propane heating systems with a vapor compression air conditioner. The version of GEO2D in the attached file has been coupled to the DOE heating and cooling load software called ENERGYPLUS. This is a great convenience for the user because heating and cooling loads are an input to GEO2D. GEO2D is a user friendly program that uses a graphical user interface for inputs and outputs. These make entering data simple and they produce many plotted results that are easy to understand. In order to run GEO2D access to MATLAB is required. If this program is not available on your computer you can download the program MCRInstaller.exe, the 64 bit version, from the MATLAB website or from this geothermal depository. This is a free download which will enable you to run GEO2D..
Comparison of 1D and 2D modelling with soil erosion model SMODERP
NASA Astrophysics Data System (ADS)
Kavka, Petr; Weyskrabova, Lenka; Zajicek, Jan
2013-04-01
The contribution presents a comparison of a runoff simulated by profile method (1D) and spatially distributed method (2D). Simulation model SMODERP is used for calculation and prediction of soil erosion and surface runoff from agricultural land. SMODERP is physically based model that includes the processes of infiltration (Phillips equation), surface runoff (kinematic wave based equation), surface retention, surface roughness and vegetation impact on runoff. 1D model was developed in past, new 2D model was developed in last two years. The model is being developed at the Department of Irrigation, Drainage and Landscape Engineering, Civil Engineering Faculty, CTU in Prague. 2D model was developed as a tool for widespread GIS software ArcGIS. The physical relations were implemented through Python script. This script uses ArcGIS system tools for raster and vectors treatment of the inputs. Flow direction is calculated by Steepest Descent algorithm in the preliminary version of 2D model. More advanced multiple flow algorithm is planned in the next version. Spatially distributed models enable to estimate not only surface runoff but also flow in the rills. Surface runoff is described in the model by kinematic wave equation. Equation uses Manning roughness coefficient for surface runoff. Parameters for five different soil textures were calibrated on the set of forty measurements performed on the laboratory rainfall simulator. For modelling of the rills a specific sub model was created. This sub model uses Manning formula for flow estimation. Numerical stability of the model is solved by Courant criterion. Spatial scale is fixed. Time step is dynamically changed depending on how flow is generated and developed. SMODERP is meant to be used not only for the research purposes, but mainly for the engineering practice. We also present how the input data can be obtained based on available resources (soil maps and data, land use, terrain models, field research, etc.) and how can
The energy spectrum in a barotropic atmosphere
NASA Astrophysics Data System (ADS)
Kurgansky, M. V.
2008-03-01
In a forced-dissipative barotropic model of the atmosphere on a spherical planet, by following mathematical techniques in (Thompson, P. D.: The equilibrium energy spectrum of randomly forced two-dimensional turbulence, Journal of the Atmospheric Sciences, 30, 1593-1598, 1973) but applying them in a novel context of the discrete spectrum on a rotating sphere, the "minus 2" energy spectrum for wavenumbers much greater than a characteristic wavenumber of the baroclinic forcing has been obtained if the forcing is taken in the simplest and most fundamental form. Some observation-based atmospheric kinetic energy spectra, with their slopes lying between "minus 2" and "minus 3" laws, are discussed from the perspective of the deduced "minus 2" energy spectrum.
Barotropic index w-singularities in cosmology
DaPbrowski, Mariusz P.; Denkiewicz, Tomasz
2009-03-15
We find an explicit cosmological model which allows a special type of cosmological singularity which we call a w-singularity. This singularity has the scale factor finite, the energy density and pressure vanishing, and the only singular behavior appears in a time-dependent barotropic index w(t). It is different from the type IV cosmological singularity in that it does not exhibit the divergence of the higher derivatives of the Hubble parameter and from the big brake since it does not fulfill the anti-Chaplygin gas equation of state. We also find an interesting duality between the w-singularities and the big-bang singularities. Physical examples of w-singularities appear in f(R), scalar field and brane cosmologies.
Barotropic instability of midlatitude zonal jets on Mars, earth and Venus
NASA Technical Reports Server (NTRS)
Michelangeli, Diane V.; Zurek, Richard W.; Elson, Lee S.
1987-01-01
A linearized, nondivergent, barotropic vorticity model on a sphere is used to intercompare the fastest growing, barotropically unstable wave modes computed for zonal jets at high latitudes in the middle atmospheres of Venus, earth, and Mars. The model is briefly described, and the choice of a range of zonal jet parameters - primarily Rossby numbers and the jet width - appropriate to Venus and Mars is discussed. The results are presented and compared with those found by Elson (1982) and Hartmann (1983) for nondivergent, barotropically unstable modes in relatively broad, midlatitude zonal jets in planetary atmospheres. Some specific examples appropriate to Venus and Mars are presented.
Analysis of vegetation effect on waves using a vertical 2-D RANS model
Technology Transfer Automated Retrieval System (TEKTRAN)
A vertical two-dimensional (2-D) model has been applied in the simulation of wave propagation through vegetated water bodies. The model is based on an existing model SOLA-VOF which solves the Reynolds-Averaged Navier-Stokes (RANS) equations with the finite difference method on a staggered rectangula...
Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model
Balakrishnan, Minimol; Chakravarthy, V. Srinivasa; Guhathakurta, Soma
2015-01-01
Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias. PMID:26733873
Fast 2D flood modelling using GPU technology - recent applications and new developments
NASA Astrophysics Data System (ADS)
Crossley, Amanda; Lamb, Rob; Waller, Simon; Dunning, Paul
2010-05-01
In recent years there has been considerable interest amongst scientists and engineers in exploiting the potential of commodity graphics hardware for desktop parallel computing. The Graphics Processing Units (GPUs) that are used in PC graphics cards have now evolved into powerful parallel co-processors that can be used to accelerate the numerical codes used for floodplain inundation modelling. We report in this paper on experience over the past two years in developing and applying two dimensional (2D) flood inundation models using GPUs to achieve significant practical performance benefits. Starting with a solution scheme for the 2D diffusion wave approximation to the 2D Shallow Water Equations (SWEs), we have demonstrated the capability to reduce model run times in ‘real-world' applications using GPU hardware and programming techniques. We then present results from a GPU-based 2D finite volume SWE solver. A series of numerical test cases demonstrate that the model produces outputs that are accurate and consistent with reference results published elsewhere. In comparisons conducted for a real world test case, the GPU-based SWE model was over 100 times faster than the CPU version. We conclude with some discussion of practical experience in using the GPU technology for flood mapping applications, and for research projects investigating use of Monte Carlo simulation methods for the analysis of uncertainty in 2D flood modelling.
Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model.
Balakrishnan, Minimol; Chakravarthy, V Srinivasa; Guhathakurta, Soma
2015-01-01
Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias. PMID:26733873
2D face database diversification based on 3D face modeling
NASA Astrophysics Data System (ADS)
Wang, Qun; Li, Jiang; Asari, Vijayan K.; Karim, Mohammad A.
2011-05-01
Pose and illumination are identified as major problems in 2D face recognition (FR). It has been theoretically proven that the more diversified instances in the training phase, the more accurate and adaptable the FR system appears to be. Based on this common awareness, researchers have developed a large number of photographic face databases to meet the demand for data training purposes. In this paper, we propose a novel scheme for 2D face database diversification based on 3D face modeling and computer graphics techniques, which supplies augmented variances of pose and illumination. Based on the existing samples from identical individuals of the database, a synthesized 3D face model is employed to create composited 2D scenarios with extra light and pose variations. The new model is based on a 3D Morphable Model (3DMM) and genetic type of optimization algorithm. The experimental results show that the complemented instances obviously increase diversification of the existing database.
A multispeed Discrete Boltzmann Model for transcritical 2D shallow water flows
NASA Astrophysics Data System (ADS)
La Rocca, Michele; Montessori, Andrea; Prestininzi, Pietro; Succi, Sauro
2015-03-01
In this work a Discrete Boltzmann Model for the solution of transcritical 2D shallow water flows is presented and validated. In order to provide the model with transcritical capabilities, a particular multispeed velocity set has been employed for the discretization of the Boltzmann equation. It is shown that this particular set naturally yields a simple and closed procedure to determine higher order equilibrium distribution functions needed to simulate transcritical flow. The model is validated through several classical benchmarks and is proven to correctly and accurately simulate both 1D and 2D transitions between the two flow regimes.
A 2-D dynamical model of mesospheric temperature inversions in winter
Hauchecorne, A.; Maillard, A. )
1990-11-01
A 2-D stratospheric and mesospheric dynamical model including drag and diffusion due to gravity wave breaking is used to simulate winter mesospheric temperature inversions similar to those observed by Rayleigh lidar. It is shown that adiabatic heating associated to descending velocities in the mesosphere is the main mechanism involved in the formation of such inversions. Sensitivity tests are performed with the model and confirm this assumption. It is also explained why other previous similar studies with 2-D models did not show mesospheric inversion layers.
Merging of RVR meander with CONCEPTS: Simplified 2D model for long-term meander evolution
Technology Transfer Automated Retrieval System (TEKTRAN)
RVR Meander is a simplified two-dimensional (2D) hydrodynamic and migration model (Abad and Garcia, 2006) while CONCEPTS (CONservational Channel Evolution and Pollutant Transport System) is a one-dimensional (1D) hydrodynamic and morphodynamic model (Langendoen and Alonso, 2008; Langendoen and Simon...
Introducing the R2D2 Model: Online Learning for the Diverse Learners of This World
ERIC Educational Resources Information Center
Bonk, Curtis J.; Zhang, Ke
2006-01-01
The R2D2 method--read, reflect, display, and do--is a new model for designing and delivering distance education, and in particular, online learning. Such a model is especially important to address the diverse preferences of online learners of varied generations and varied Internet familiarity. Four quadrants can be utilized separately or as part…
Evaluation of 2D shallow-water model for spillway flow with a complex geometry
Technology Transfer Automated Retrieval System (TEKTRAN)
Although the two-dimensional (2D) shallow water model is formulated based on several assumptions such as hydrostatic pressure distribution and vertical velocity is negligible, as a simple alternative to the complex 3D model, it has been used to compute water flows in which these assumptions may be ...
Knight shift and spin relaxation in the single band 2D Hubbard model
NASA Astrophysics Data System (ADS)
Leblanc, James; Chen, Xi; Gull, Emanuel
We study in detail the roles of spin and charge fluctuations in the single band 2D Hubbard model. Using dynamical mean field theory and cluster extensions such as the dynamical cluster approximation (DCA), we compute the full two particle susceptibilities in the spin and charge representations. By performing analytic continuations we obtain the temperature and doping dependence of the spin-lattice relaxation (T1- 1) and knight shift in the 2D Hubbard model relevant to NMR results on doped cuprates and connect these to RPA results in weak coupling limits.
NASA Astrophysics Data System (ADS)
Pasternack, Gregory B.; Gilbert, Andrew T.; Wheaton, Joseph M.; Buckland, Evan M.
2006-08-01
SummaryResource managers, scientists, government regulators, and stakeholders are considering sophisticated numerical models for managing complex environmental problems. In this study, observations from a river-rehabilitation experiment involving gravel augmentation and spawning habitat enhancement were used to assess sources and magnitudes of error in depth, velocity, and shear velocity predictions made at the 1-m scale with a commercial two-dimensional (depth-averaged) model. Error in 2D model depth prediction averaged 21%. This error was attributable to topographic survey resolution, which at 1 pt per 1.14 m 2, was inadequate to resolve small humps and depressions influencing point measurements. Error in 2D model velocity prediction averaged 29%. More than half of this error was attributable to depth prediction error. Despite depth and velocity error, 56% of tested 2D model predictions of shear velocity were within the 95% confidence limit of the best field-based estimation method. Ninety percent of the error in shear velocity prediction was explained by velocity prediction error. Multiple field-based estimates of shear velocity differed by up to 160%, so the lower error for the 2D model's predictions suggests such models are at least as accurate as field measurement. 2D models enable detailed, spatially distributed estimates compared to the small number measurable in a field campaign of comparable cost. They also can be used for design evaluation. Although such numerical models are limited to channel types adhering to model assumptions and yield predictions only accurate to ˜20-30%, they can provide a useful tool for river-rehabilitation design and assessment, including spatially diverse habitat heterogeneity as well as for pre- and post-project appraisal.
Validation of DYSTOOL for unsteady aerodynamic modeling of 2D airfoils
NASA Astrophysics Data System (ADS)
González, A.; Gomez-Iradi, S.; Munduate, X.
2014-06-01
From the point of view of wind turbine modeling, an important group of tools is based on blade element momentum (BEM) theory using 2D aerodynamic calculations on the blade elements. Due to the importance of this sectional computation of the blades, the National Renewable Wind Energy Center of Spain (CENER) developed DYSTOOL, an aerodynamic code for 2D airfoil modeling based on the Beddoes-Leishman model. The main focus here is related to the model parameters, whose values depend on the airfoil or the operating conditions. In this work, the values of the parameters are adjusted using available experimental or CFD data. The present document is mainly related to the validation of the results of DYSTOOL for 2D airfoils. The results of the computations have been compared with unsteady experimental data of the S809 and NACA0015 profiles. Some of the cases have also been modeled using the CFD code WMB (Wind Multi Block), within the framework of a collaboration with ACCIONA Windpower. The validation has been performed using pitch oscillations with different reduced frequencies, Reynolds numbers, amplitudes and mean angles of attack. The results have shown a good agreement using the methodology of adjustment for the value of the parameters. DYSTOOL have demonstrated to be a promising tool for 2D airfoil unsteady aerodynamic modeling.
2D-photochemical model for forbidden oxygen line emission for comet 1P/Halley
NASA Astrophysics Data System (ADS)
Cessateur, G.; De Keyser, J.; Maggiolo, R.; Rubin, M.; Gronoff, G.; Gibbons, A.; Jehin, E.; Dhooghe, F.; Gunell, H.; Vaeck, N.; Loreau, J.
2016-08-01
We present here a 2D-model of photochemistry for computing the production and loss mechanisms of the O(1S) and O(1D) states, which are responsible for the emission lines at 577.7 nm, 630 nm, and 636.4 nm, in case of the comet 1P/Halley. The presence of O2 within cometary atmospheres, measured by the in-situ ROSETTA and GIOTTO missions, necessitates a revision of the usual photochemical models. Indeed, the photodissociation of molecular oxygen also leads to a significant production of oxygen in excited electronic states. In order to correctly model the solar UV flux absorption, we consider here a 2D configuration. While the green to red-doublet ratio is not affected by the solar UV flux absorption, estimates of the red-doublet and green lines emissions are, however, overestimated by a factor of two in the 1D model compared to the 2D model. Considering a spherical symmetry, emission maps can be deduced from the 2D model in order to be directly compared to ground and/or in-situ observations.
Oriented Gaussian mixture models for nonrigid 2D/3D coronary artery registration.
Baka, N; Metz, C T; Schultz, C J; van Geuns, R-J; Niessen, W J; van Walsum, T
2014-05-01
2D/3D registration of patient vasculature from preinterventional computed tomography angiography (CTA) to interventional X-ray angiography is of interest to improve guidance in percutaneous coronary interventions. In this paper we present a novel feature based 2D/3D registration framework, that is based on probabilistic point correspondences, and show its usefulness on aligning 3D coronary artery centerlines derived from CTA images with their 2D projection derived from interventional X-ray angiography. The registration framework is an extension of the Gaussian mixture model (GMM) based point-set registration to the 2D/3D setting, with a modified distance metric. We also propose a way to incorporate orientation in the registration, and show its added value for artery registration on patient datasets as well as in simulation experiments. The oriented GMM registration achieved a median accuracy of 1.06 mm, with a convergence rate of 81% for nonrigid vessel centerline registration on 12 patient datasets, using a statistical shape model. The method thereby outperformed the iterative closest point algorithm, the GMM registration without orientation, and two recently published methods on 2D/3D coronary artery registration. PMID:24770908
The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis
Shaheen, Eman De Keyzer, Frederik; Bosmans, Hilde; Ongeval, Chantal Van; Dance, David R.; Young, Kenneth C.
2014-08-15
Purpose: This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images. Methods: Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses. Results: The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly
Impact of high speed civil transports on stratospheric ozone: A 2-D model investigation
Kinnison, D.E.; Connell, P.S.
1996-12-01
This study investigates the effect on stratospheric ozone from a fleet of proposed High Speed Civil Transports (HSCTs). The new LLNL 2-D operator-split chemical-radiative-transport model of the troposphere and stratosphere is used for this HSCT investigation. This model is integrated in a diurnal manner, using an implicit numerical solver. Therefore, rate coefficients are not modified by any sort of diurnal average factor. This model also does not make any assumptions on lumping of chemical species into families. Comparisons to previous model-derived HSCT assessment of ozone change are made, both to the previous LLNL 2-D model and to other models from the international assessment modeling community. The sensitivity to the NO{sub x} emission index and sulfate surface area density is also explored.
Extension and application of the Preissmann slot model to 2D transient mixed flows
NASA Astrophysics Data System (ADS)
Maranzoni, Andrea; Dazzi, Susanna; Aureli, Francesca; Mignosa, Paolo
2015-08-01
This paper presents an extension of the Preissmann slot concept for the modeling of highly transient two-dimensional (2D) mixed flows. The classic conservative formulation of the 2D shallow water equations for free surface flows is adapted by assuming that two fictitious vertical slots, aligned along the two Cartesian plane directions and normally intersecting, are added on the ceiling of each integration element. Accordingly, transitions between free surface and pressurized flow can be handled in a natural and straightforward way by using the same set of governing equations. The opportunity of coupling free surface and pressurized flows is actually useful not only in one-dimensional (1D) problems concerning sewer systems but also for modeling 2D flooding phenomena in which the pressurization of bridges, culverts, or other crossing hydraulic structures can be expected. Numerical simulations are performed by using a shock-capturing MUSCL-Hancock finite volume scheme combined with the FORCE (First-Order Centred) solver for the evaluation of the numerical fluxes. The validation of the mathematical model is accomplished on the basis of both exact solutions of 1D discontinuous initial value problems and reference radial solutions of idealized test cases with cylindrical symmetry. Furthermore, the capability of the model to deal with practical field-scale applications is assessed by simulating the transit of a bore under an arch bridge. Numerical results show that the proposed model is suitable for the prediction of highly transient 2D mixed flows.
Comparison of 3-D finite element model of ashlar masonry with 2-D numerical models of ashlar masonry
NASA Astrophysics Data System (ADS)
Beran, Pavel
2016-06-01
3-D state of stress in heterogeneous ashlar masonry can be also computed by several suitable chosen 2-D numerical models of ashlar masonry. The results obtained from 2-D numerical models well correspond to the results obtained from 3-D numerical model. The character of thermal stress is the same. While using 2-D models the computational time is reduced more than hundredfold and therefore this method could be used for computation of thermal stresses during long time periods with 10 000 of steps.
2D-Raman-THz spectroscopy: A sensitive test of polarizable water models
NASA Astrophysics Data System (ADS)
Hamm, Peter
2014-11-01
In a recent paper, the experimental 2D-Raman-THz response of liquid water at ambient conditions has been presented [J. Savolainen, S. Ahmed, and P. Hamm, Proc. Natl. Acad. Sci. U. S. A. 110, 20402 (2013)]. Here, all-atom molecular dynamics simulations are performed with the goal to reproduce the experimental results. To that end, the molecular response functions are calculated in a first step, and are then convoluted with the laser pulses in order to enable a direct comparison with the experimental results. The molecular dynamics simulation are performed with several different water models: TIP4P/2005, SWM4-NDP, and TL4P. As polarizability is essential to describe the 2D-Raman-THz response, the TIP4P/2005 water molecules are amended with either an isotropic or a anisotropic polarizability a posteriori after the molecular dynamics simulation. In contrast, SWM4-NDP and TL4P are intrinsically polarizable, and hence the 2D-Raman-THz response can be calculated in a self-consistent way, using the same force field as during the molecular dynamics simulation. It is found that the 2D-Raman-THz response depends extremely sensitively on details of the water model, and in particular on details of the description of polarizability. Despite the limited time resolution of the experiment, it could easily distinguish between various water models. Albeit not perfect, the overall best agreement with the experimental data is obtained for the TL4P water model.
Parallelized CCHE2D flow model with CUDA Fortran on Graphics Process Units
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents the CCHE2D implicit flow model parallelized using CUDA Fortran programming technique on Graphics Processing Units (GPUs). A parallelized implicit Alternating Direction Implicit (ADI) solver using Parallel Cyclic Reduction (PCR) algorithm on GPU is developed and tested. This solve...
ERIC Educational Resources Information Center
Park, Elisa L.
2009-01-01
The purpose of this study is to understand the dynamics of Korean students' international mobility to study abroad by using the 2-D Model. The first D, "the driving force factor," explains how and what components of the dissatisfaction with domestic higher education perceived by Korean students drives students' outward mobility to seek foreign…
2D-Raman-THz spectroscopy: A sensitive test of polarizable water models
Hamm, Peter
2014-11-14
In a recent paper, the experimental 2D-Raman-THz response of liquid water at ambient conditions has been presented [J. Savolainen, S. Ahmed, and P. Hamm, Proc. Natl. Acad. Sci. U. S. A. 110, 20402 (2013)]. Here, all-atom molecular dynamics simulations are performed with the goal to reproduce the experimental results. To that end, the molecular response functions are calculated in a first step, and are then convoluted with the laser pulses in order to enable a direct comparison with the experimental results. The molecular dynamics simulation are performed with several different water models: TIP4P/2005, SWM4-NDP, and TL4P. As polarizability is essential to describe the 2D-Raman-THz response, the TIP4P/2005 water molecules are amended with either an isotropic or a anisotropic polarizability a posteriori after the molecular dynamics simulation. In contrast, SWM4-NDP and TL4P are intrinsically polarizable, and hence the 2D-Raman-THz response can be calculated in a self-consistent way, using the same force field as during the molecular dynamics simulation. It is found that the 2D-Raman-THz response depends extremely sensitively on details of the water model, and in particular on details of the description of polarizability. Despite the limited time resolution of the experiment, it could easily distinguish between various water models. Albeit not perfect, the overall best agreement with the experimental data is obtained for the TL4P water model.
Justification for a 2D versus 3D fingertip finite element model during static contact simulations.
Harih, Gregor; Tada, Mitsunori; Dolšak, Bojan
2016-10-01
The biomechanical response of a human hand during contact with various products has not been investigated in details yet. It has been shown that excessive contact pressure on the soft tissue can result in discomfort, pain and also cumulative traumatic disorders. This manuscript explores the benefits and limitations of a simplified two-dimensional vs. an anatomically correct three-dimensional finite element model of a human fingertip. Most authors still use 2D FE fingertip models due to their simplicity and reduced computational costs. However we show that an anatomically correct 3D FE fingertip model can provide additional insight into the biomechanical behaviour. The use of 2D fingertip FE models is justified when observing peak contact pressure values as well as displacement during the contact for the given studied cross-section. On the other hand, an anatomically correct 3D FE fingertip model provides a contact pressure distribution, which reflects the fingertip's anatomy. PMID:26856769
Molecular Dynamics implementation of BN2D or 'Mercedes Benz' water model
NASA Astrophysics Data System (ADS)
Scukins, Arturs; Bardik, Vitaliy; Pavlov, Evgen; Nerukh, Dmitry
2015-05-01
Two-dimensional 'Mercedes Benz' (MB) or BN2D water model (Naim, 1971) is implemented in Molecular Dynamics. It is known that the MB model can capture abnormal properties of real water (high heat capacity, minima of pressure and isothermal compressibility, negative thermal expansion coefficient) (Silverstein et al., 1998). In this work formulas for calculating the thermodynamic, structural and dynamic properties in microcanonical (NVE) and isothermal-isobaric (NPT) ensembles for the model from Molecular Dynamics simulation are derived and verified against known Monte Carlo results. The convergence of the thermodynamic properties and the system's numerical stability are investigated. The results qualitatively reproduce the peculiarities of real water making the model a visually convenient tool that also requires less computational resources, thus allowing simulations of large (hydrodynamic scale) molecular systems. We provide the open source code written in C/C++ for the BN2D water model implementation using Molecular Dynamics.
MODELING THE TRANSVERSE THERMAL CONDUCTIVITY OF 2-D SICF/SIC COMPOSITES MADE WITH WOVEN FABRIC
Youngblood, Gerald E; Senor, David J; Jones, Russell H
2004-06-01
The hierarchical two-layer (H2L) model describes the effective transverse thermal conductivity (Keff) of a 2D-SiCf/SiC composite plate made from stacked and infiltrated woven fabric layers in terms of constituent properties and microstructural and architectural variables. The H2L model includes the effects of fiber-matrix interfacial conductance, high fiber packing fractions within individual tows and the non-uniform nature of 2D fabric/matrix layers that usually include a significant amount of interlayer porosity. Previously, H2L model Keff-predictions were compared to measured values for two versions of 2D Hi-Nicalon/PyC/ICVI-SiC composite, one with a “thin” (0.11m) and the other with a “thick” (1.04m) pyrocarbon (PyC) fiber coating, and for a 2D Tyranno SA/”thin” PyC/FCVI-SIC composite. In this study, H2L model Keff-predictions were compared to measured values for a 2D-SiCf/SiC composite made using the ICVI-process with Hi-Nicalon type S fabric and a “thin” PyC fiber coating. The values of Keff determined for the latter composite were significantly greater than the Keff-values determined for the composites made with either the Hi-Nicalon or the Tyranno SA fabrics. Differences in Keff-values were expected for the different fiber types, but major differences also were due to observed microstructural and architectural variations between the composite systems, and as predicted by the H2L model.
The 2dF Galaxy Redshift Survey: voids and hierarchical scaling models
NASA Astrophysics Data System (ADS)
Croton, Darren J.; Colless, Matthew; Gaztañaga, Enrique; Baugh, Carlton M.; Norberg, Peder; Baldry, I. K.; Bland-Hawthorn, J.; Bridges, T.; Cannon, R.; Cole, S.; Collins, C.; Couch, W.; Dalton, G.; de Propris, R.; Driver, S. P.; Efstathiou, G.; Ellis, R. S.; Frenk, C. S.; Glazebrook, K.; Jackson, C.; Lahav, O.; Lewis, I.; Lumsden, S.; Maddox, S.; Madgwick, D.; Peacock, J. A.; Peterson, B. A.; Sutherland, W.; Taylor, K.
2004-08-01
We measure the redshift-space reduced void probability function (VPF) for 2dFGRS volume-limited galaxy samples covering the absolute magnitude range MbJ-5log10h=-18 to -22. Theoretically, the VPF connects the distribution of voids to the moments of galaxy clustering of all orders, and can be used to discriminate clustering models in the weakly non-linear regime. The reduced VPF measured from the 2dFGRS is in excellent agreement with the paradigm of hierarchical scaling of the galaxy clustering moments. The accuracy of our measurement is such that we can rule out, at a very high significance, popular models for galaxy clustering, including the lognormal distribution. We demonstrate that the negative binomial model gives a very good approximation to the 2dFGRS data over a wide range of scales, out to at least 20 h-1 Mpc. Conversely, the reduced VPF for dark matter in a Λ cold dark matter (ΛCDM) universe does appear to be lognormal on small scales but deviates significantly beyond ~4 h-1 Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy luminosity. Our results hold independently in both the North and South Galactic Pole survey regions.
Simplified 2D Bidomain Model of Whole Heart Electrical Activity and ECG Generation
NASA Astrophysics Data System (ADS)
Sovilj, Siniša; Magjarević, Ratko; Abed, Amr Al; Lovell, Nigel H.; Dokos, Socrates
2014-06-01
The aim of this study was the development of a geometrically simple and highly computationally-efficient two dimensional (2D) biophysical model of whole heart electrical activity, incorporating spontaneous activation of the sinoatrial node (SAN), the specialized conduction system, and realistic surface ECG morphology computed on the torso. The FitzHugh-Nagumo (FHN) equations were incorporated into a bidomain finite element model of cardiac electrical activity, which was comprised of a simplified geometry of the whole heart with the blood cavities, the lungs and the torso as an extracellular volume conductor. To model the ECG, we placed four electrodes on the surface of the torso to simulate three Einthoven leads VI, VII and VIII from the standard 12-lead system. The 2D model was able to reconstruct ECG morphology on the torso from action potentials generated at various regions of the heart, including the sinoatrial node, atria, atrioventricular node, His bundle, bundle branches, Purkinje fibers, and ventricles. Our 2D cardiac model offers a good compromise between computational load and model complexity, and can be used as a first step towards three dimensional (3D) ECG models with more complex, precise and accurate geometry of anatomical structures, to investigate the effect of various cardiac electrophysiological parameters on ECG morphology.
Approximate analytic solutions to 3D unconfined groundwater flow within regional 2D models
NASA Astrophysics Data System (ADS)
Luther, K.; Haitjema, H. M.
2000-04-01
We present methods for finding approximate analytic solutions to three-dimensional (3D) unconfined steady state groundwater flow near partially penetrating and horizontal wells, and for combining those solutions with regional two-dimensional (2D) models. The 3D solutions use distributed singularities (analytic elements) to enforce boundary conditions on the phreatic surface and seepage faces at vertical wells, and to maintain fixed-head boundary conditions, obtained from the 2D model, at the perimeter of the 3D model. The approximate 3D solutions are analytic (continuous and differentiable) everywhere, including on the phreatic surface itself. While continuity of flow is satisfied exactly in the infinite 3D flow domain, water balance errors can occur across the phreatic surface.
2D and 3D numerical models on compositionally buoyant diapirs in the mantle wedge
NASA Astrophysics Data System (ADS)
Hasenclever, Jörg; Morgan, Jason Phipps; Hort, Matthias; Rüpke, Lars H.
2011-11-01
We present 2D and 3D numerical model calculations that focus on the physics of compositionally buoyant diapirs rising within a mantle wedge corner flow. Compositional buoyancy is assumed to arise from slab dehydration during which water-rich volatiles enter the mantle wedge and form a wet, less dense boundary layer on top of the slab. Slab dehydration is prescribed to occur in the 80-180 km deep slab interval, and the water transport is treated as a diffusion-like process. In this study, the mantle's rheology is modeled as being isoviscous for the benefit of easier-to-interpret feedbacks between water migration and buoyant viscous flow of the mantle. We use a simple subduction geometry that does not change during the numerical calculation. In a large set of 2D calculations we have identified that five different flow regimes can form, in which the position, number, and formation time of the diapirs vary as a function of four parameters: subduction angle, subduction rate, water diffusivity (mobility), and mantle viscosity. Using the same numerical method and numerical resolution we also conducted a suite of 3D calculations for 16 selected parameter combinations. Comparing the 2D and 3D results for the same model parameters reveals that the 2D models can only give limited insights into the inherently 3D problem of mantle wedge diapirism. While often correctly predicting the position and onset time of the first diapir(s), the 2D models fail to capture the dynamics of diapir ascent as well as the formation of secondary diapirs that result from boundary layer perturbations caused by previous diapirs. Of greatest importance for physically correct results is the numerical resolution in the region where diapirs nucleate, which must be high enough to accurately capture the growth of the thin wet boundary layer on top of the slab and, subsequently, the formation, morphology, and ascent of diapirs. Here 2D models can be very useful to quantify the required resolution, which we
TRENT2D WG: a smart web infrastructure for debris-flow modelling and hazard assessment
NASA Astrophysics Data System (ADS)
Zorzi, Nadia; Rosatti, Giorgio; Zugliani, Daniel; Rizzi, Alessandro; Piffer, Stefano
2016-04-01
Mountain regions are naturally exposed to geomorphic flows, which involve large amounts of sediments and induce significant morphological modifications. The physical complexity of this class of phenomena represents a challenging issue for modelling, leading to elaborate theoretical frameworks and sophisticated numerical techniques. In general, geomorphic-flows models proved to be valid tools in hazard assessment and management. However, model complexity seems to represent one of the main obstacles to the diffusion of advanced modelling tools between practitioners and stakeholders, although the UE Flood Directive (2007/60/EC) requires risk management and assessment to be based on "best practices and best available technologies". Furthermore, several cutting-edge models are not particularly user-friendly and multiple stand-alone software are needed to pre- and post-process modelling data. For all these reasons, users often resort to quicker and rougher approaches, leading possibly to unreliable results. Therefore, some effort seems to be necessary to overcome these drawbacks, with the purpose of supporting and encouraging a widespread diffusion of the most reliable, although sophisticated, modelling tools. With this aim, this work presents TRENT2D WG, a new smart modelling solution for the state-of-the-art model TRENT2D (Armanini et al., 2009, Rosatti and Begnudelli, 2013), which simulates debris flows and hyperconcentrated flows adopting a two-phase description over a mobile bed. TRENT2D WG is a web infrastructure joining advantages offered by the software-delivering model SaaS (Software as a Service) and by WebGIS technology and hosting a complete and user-friendly working environment for modelling. In order to develop TRENT2D WG, the model TRENT2D was converted into a service and exposed on a cloud server, transferring computational burdens from the user hardware to a high-performing server and reducing computational time. Then, the system was equipped with an
NASA Technical Reports Server (NTRS)
Gao, Shou-Ting; Ping, Fan; Li, Xiao-Fan; Tao, Wei-Kuo
2004-01-01
Although dry/moist potential vorticity is a useful physical quantity for meteorological analysis, it cannot be applied to the analysis of 2D simulations. A convective vorticity vector (CVV) is introduced in this study to analyze 2D cloud-resolving simulation data associated with 2D tropical convection. The cloud model is forced by the vertical velocity, zonal wind, horizontal advection, and sea surface temperature obtained from the TOGA COARE, and is integrated for a selected 10-day period. The CVV has zonal and vertical components in the 2D x-z frame. Analysis of zonally-averaged and mass-integrated quantities shows that the correlation coefficient between the vertical component of the CVV and the sum of the cloud hydrometeor mixing ratios is 0.81, whereas the correlation coefficient between the zonal component and the sum of the mixing ratios is only 0.18. This indicates that the vertical component of the CVV is closely associated with tropical convection. The tendency equation for the vertical component of the CVV is derived and the zonally-averaged and mass-integrated tendency budgets are analyzed. The tendency of the vertical component of the CVV is determined by the interaction between the vorticity and the zonal gradient of cloud heating. The results demonstrate that the vertical component of the CVV is a cloud-linked parameter and can be used to study tropical convection.
Absolute barotropic instability and monsoon depressions
NASA Technical Reports Server (NTRS)
Lindzen, R. S.; Farrell, B.; Rosenthal, A. J.
1983-01-01
Monsoon depressions over the Bay of Bengal develop almost entirely in July and August. After studies conducted by Lindzen et al. (1980) and Stevens and Lindzen (1978), only barotropic instability remains as a mechanism for the development of the wave disturbances associated with monsoon depressions. The present investigation has the objective to show that barotropic instability is able to explain the wave aspects of monsoon depressions, but that normal mode analysis is inadequate. It is found that the local barotropically unstable response to regional perturbations in the Bay of Bengal during July and August will be dominated by the lower troposphere. The analysis clearly identifies the features of the mean flow which lead to monsoon depressions in July. The features include the development of an easterly jet as opposed to semijet structure in the mean flow, and the development of a modest easterly flow at the jet center as opposed to westerly flow.
Parameterising root system growth models using 2D neutron radiography images
NASA Astrophysics Data System (ADS)
Schnepf, Andrea; Felderer, Bernd; Vontobel, Peter; Leitner, Daniel
2013-04-01
Root architecture is a key factor for plant acquisition of water and nutrients from soil. In particular in view of a second green revolution where the below ground parts of agricultural crops are important, it is essential to characterise and quantify root architecture and its effect on plant resource acquisition. Mathematical models can help to understand the processes occurring in the soil-plant system, they can be used to quantify the effect of root and rhizosphere traits on resource acquisition and the response to environmental conditions. In order to do so, root architectural models are coupled with a model of water and solute transport in soil. However, dynamic root architectural models are difficult to parameterise. Novel imaging techniques such as x-ray computed tomography, neutron radiography and magnetic resonance imaging enable the in situ visualisation of plant root systems. Therefore, these images facilitate the parameterisation of dynamic root architecture models. These imaging techniques are capable of producing 3D or 2D images. Moreover, 2D images are also available in the form of hand drawings or from images of standard cameras. While full 3D imaging tools are still limited in resolutions, 2D techniques are a more accurate and less expensive option for observing roots in their environment. However, analysis of 2D images has additional difficulties compared to the 3D case, because of overlapping roots. We present a novel algorithm for the parameterisation of root system growth models based on 2D images of root system. The algorithm analyses dynamic image data. These are a series of 2D images of the root system at different points in time. Image data has already been adjusted for missing links and artefacts and segmentation was performed by applying a matched filter response. From this time series of binary 2D images, we parameterise the dynamic root architecture model in the following way: First, a morphological skeleton is derived from the binary
Validation of a 2-D semi-coupled numerical model for fluid-structure-seabed interaction
NASA Astrophysics Data System (ADS)
Ye, Jianhong; Jeng, Dongsheng; Wang, Ren; Zhu, Changqi
2013-10-01
A 2-D semi-coupled model PORO-WSSI 2D (also be referred as FSSI-CAS 2D) for the Fluid-Structure-Seabed Interaction (FSSI) has been developed by employing RANS equations for wave motion in fluid domain, VARANS equations for porous flow in porous structures; and taking the dynamic Biot's equations (known as "u - p" approximation) for soil as the governing equations. The finite difference two-step projection method and the forward time difference method are adopted to solve the RANS, VARANS equations; and the finite element method is adopted to solve the "u - p" approximation. A data exchange port is developed to couple the RANS, VARANS equations and the dynamic Biot's equations together. The analytical solution proposed by Hsu and Jeng (1994) and some experiments conducted in wave flume or geotechnical centrifuge in which various waves involved are used to validate the developed semi-coupled numerical model. The sandy bed involved in these experiments is poro-elastic or poro-elastoplastic. The inclusion of the interaction between fluid, marine structures and poro-elastoplastic seabed foundation is a special point and highlight in this paper, which is essentially different with other previous coupled models The excellent agreement between the numerical results and the experiment data indicates that the developed coupled model is highly reliablefor the FSSI problem.
2D density model of the Chinese continental lithosphere along a NW-SE transect
NASA Astrophysics Data System (ADS)
Šimonová, Barbora; Bielik, Miroslav; Dérerová, Jana
2015-06-01
This paper presents a 2D density model along a transect from NW to SE China. The model was first constructed by the transformation of seismic velocity to density, revealed by previous deep seismic soundings (DSS) investigations in China. Then, the 2D density model was updated using the GM-SYS software by fitting the computed to the observed gravity data. Based on the density distribution of anomalous layers we divided the Chinese continental crust along the transect into three regions: north-western, central and south-eastern. The first one includes the Junggar Basin, Tianshan and Tarim Basin. The second part consists of the Qilian Orogen, the Qaidam Basin and the Songpan Ganzi Basin. The third region is represented by the Yangtze and the Cathaysia blocks. The low velocity body (vp =5.2 - 6.2 km/s) at the junction of the North-western and Central parts at a depth between 21 - 31 km, which was discovered out by DSS, was also confirmed by our 2D density modelling.
A Neural-FEM tool for the 2-D magnetic hysteresis modeling
NASA Astrophysics Data System (ADS)
Cardelli, E.; Faba, A.; Laudani, A.; Lozito, G. M.; Riganti Fulginei, F.; Salvini, A.
2016-04-01
The aim of this work is to present a new tool for the analysis of magnetic field problems considering 2-D magnetic hysteresis. In particular, this tool makes use of the Finite Element Method to solve the magnetic field problem in real device, and fruitfully exploits a neural network (NN) for the modeling of 2-D magnetic hysteresis of materials. The NS has as input the magnetic inductions components B at the k-th simulation step and returns as output the corresponding values of the magnetic field H corresponding to the input pattern. It is trained by vector measurements performed on the magnetic material to be modeled. This input/output scheme is directly implemented in a FEM code employing the magnetic potential vector A formulation. Validations through measurements on a real device have been performed.
TMRPres2D: high quality visual representation of transmembrane protein models.
Spyropoulos, Ioannis C; Liakopoulos, Theodore D; Bagos, Pantelis G; Hamodrakas, Stavros J
2004-11-22
The 'TransMembrane protein Re-Presentation in 2-Dimensions' (TMRPres2D) tool, automates the creation of uniform, two-dimensional, high analysis graphical images/models of alpha-helical or beta-barrel transmembrane proteins. Protein sequence data and structural information may be acquired from public protein knowledge bases, emanate from prediction algorithms, or even be defined by the user. Several important biological and physical sequence attributes can be embedded in the graphical representation. PMID:15201184
A simple 2-D inundation model for incorporating flood damage in urban drainage planning
NASA Astrophysics Data System (ADS)
Pathirana, A.; Tsegaye, S.; Gersonius, B.; Vairavamoorthy, K.
2008-11-01
In this paper a new inundation model code is developed and coupled with Storm Water Management Model, SWMM, to relate spatial information associated with urban drainage systems as criteria for planning of storm water drainage networks. The prime objective is to achive a model code that is simple and fast enough to be consistently be used in planning stages of urban drainage projects. The formulation for the two-dimensional (2-D) surface flow model algorithms is based on the Navier Stokes equation in two dimensions. An Alternating Direction Implicit (ADI) finite difference numerical scheme is applied to solve the governing equations. This numerical scheme is used to express the partial differential equations with time steps split into two halves. The model algorithm is written using C++ computer programming language. This 2-D surface flow model is then coupled with SWMM for simulation of both pipe flow component and surcharge induced inundation in urban areas. In addition, a damage calculation block is integrated within the inundation model code. The coupled model is shown to be capable of dealing with various flow conditions, as well as being able to simulate wetting and drying processes that will occur as the flood flows over an urban area. It has been applied under idealized and semi-hypothetical cases to determine detailed inundation zones, depths and velocities due to surcharged water on overland surface.
Comparison between 2D and 3D Numerical Modelling of a hot forging simulative test
Croin, M.; Ghiotti, A.; Bruschi, S.
2007-04-07
The paper presents the comparative analysis between 2D and 3D modelling of a simulative experiment, performed in laboratory environment, in which operating conditions approximate hot forging of a turbine aerofoil section. The plane strain deformation was chosen as an ideal case to analyze the process because of the thickness variations in the final section and the consequent distributions of contact pressure and sliding velocity at the interface that are closed to the conditions of the real industrial process. In order to compare the performances of 2D and 3D approaches, two different analyses were performed and compared with the experiments in terms of loads and temperatures peaks at the interface between the dies and the workpiece.
Nested 1D-2D approach for urban surface flood modeling
NASA Astrophysics Data System (ADS)
Murla, Damian; Willems, Patrick
2015-04-01
Floods in urban areas as a consequence of sewer capacity exceedance receive increased attention because of trends in urbanization (increased population density and impermeability of the surface) and climate change. Despite the strong recent developments in numerical modeling of water systems, urban surface flood modeling is still a major challenge. Whereas very advanced and accurate flood modeling systems are in place and operation by many river authorities in support of flood management along rivers, this is not yet the case in urban water management. Reasons include the small scale of the urban inundation processes, the need to have very high resolution topographical information available, and the huge computational demands. Urban drainage related inundation modeling requires a 1D full hydrodynamic model of the sewer network to be coupled with a 2D surface flood model. To reduce the computational times, 0D (flood cones), 1D/quasi-2D surface flood modeling approaches have been developed and applied in some case studies. In this research, a nested 1D/2D hydraulic model has been developed for an urban catchment at the city of Gent (Belgium), linking the underground sewer (minor system) with the overland surface (major system). For the overland surface flood modelling, comparison was made of 0D, 1D/quasi-2D and full 2D approaches. The approaches are advanced by considering nested 1D-2D approaches, including infiltration in the green city areas, and allowing the effects of surface storm water storage to be simulated. An optimal nested combination of three different mesh resolutions was identified; based on a compromise between precision and simulation time for further real-time flood forecasting, warning and control applications. Main streets as mesh zones together with buildings as void regions constitute one of these mesh resolution (3.75m2 - 15m2); they have been included since they channel most of the flood water from the manholes and they improve the accuracy of
Mechanical Modelling of Pultrusion Process: 2D and 3D Numerical Approaches
NASA Astrophysics Data System (ADS)
Baran, Ismet; Hattel, Jesper H.; Akkerman, Remko; Tutum, Cem C.
2015-02-01
The process induced variations such as residual stresses and distortions are a critical issue in pultrusion, since they affect the structural behavior as well as the mechanical properties and geometrical precision of the final product. In order to capture and investigate these variations, a mechanical analysis should be performed. In the present work, the two dimensional (2D) quasi-static plane strain mechanical model for the pultrusion of a thick square profile developed by the authors is further improved using generalized plane strain elements. In addition to that, a more advanced 3D thermo-chemical-mechanical analysis is carried out using 3D quadratic elements which is a novel application for the numerical modelling of the pultrusion process. It is found that the 2D mechanical models give relatively reasonable and accurate stress and displacement evolutions in the transverse direction as compared to the 3D model. Moreover, the generalized plane strain model predicts the longitudinal process induced stresses more similar to the ones calculated in the 3D model as compared with the plane strain model.
Canonical vs. micro-canonical sampling methods in a 2D Ising model
Kepner, J.
1990-12-01
Canonical and micro-canonical Monte Carlo algorithms were implemented on a 2D Ising model. Expressions for the internal energy, U, inverse temperature, Z, and specific heat, C, are given. These quantities were calculated over a range of temperature, lattice sizes, and time steps. Both algorithms accurately simulate the Ising model. To obtain greater than three decimal accuracy from the micro-canonical method requires that the more complicated expression for Z be used. The overall difference between the algorithms is small. The physics of the problem under study should be the deciding factor in determining which algorithm to use. 13 refs., 6 figs., 2 tabs.
Complex zeros of the 2 d Ising model on dynamical random lattices
NASA Astrophysics Data System (ADS)
Ambjørn, J.; Anagnostopoulos, K. N.; Magnea, U.
1998-04-01
We study the zeros in the complex plane of the partition function for the Ising model coupled to 2 d quantum gravity for complex magnetic field and for complex temperature. We compute the zeros by using the exact solution coming from a two matrix model and by Monte Carlo simulations of Ising spins on dynamical triangulations. We present evidence that the zeros form simple one-dimensional patterns in the complex plane, and that the critical behaviour of the system is governed by the scaling of the distribution of singularities near the critical point.
Brane brick models, toric Calabi-Yau 4-folds and 2d (0,2) quivers
NASA Astrophysics Data System (ADS)
Franco, Sebastián; Lee, Sangmin; Seong, Rak-Kyeong
2016-02-01
We introduce brane brick models, a novel type of Type IIA brane configurations consisting of D4-branes ending on an NS5-brane. Brane brick models are T-dual to D1-branes over singular toric Calabi-Yau 4-folds. They fully encode the infinite class of 2 d (generically) {N}=(0,2) gauge theories on the worldvolume of the D1-branes and streamline their connection to the probed geometries. For this purpose, we also introduce new combinatorial procedures for deriving the Calabi-Yau associated to a given gauge theory and vice versa.
2D-3D Registration of CT Vertebra Volume to Fluoroscopy Projection: A Calibration Model Assessment
NASA Astrophysics Data System (ADS)
Bifulco, P.; Cesarelli, M.; Allen, R.; Romano, M.; Fratini, A.; Pasquariello, G.
2009-12-01
This study extends a previous research concerning intervertebral motion registration by means of 2D dynamic fluoroscopy to obtain a more comprehensive 3D description of vertebral kinematics. The problem of estimating the 3D rigid pose of a CT volume of a vertebra from its 2D X-ray fluoroscopy projection is addressed. 2D-3D registration is obtained maximising a measure of similarity between Digitally Reconstructed Radiographs (obtained from the CT volume) and real fluoroscopic projection. X-ray energy correction was performed. To assess the method a calibration model was realised a sheep dry vertebra was rigidly fixed to a frame of reference including metallic markers. Accurate measurement of 3D orientation was obtained via single-camera calibration of the markers and held as true 3D vertebra position; then, vertebra 3D pose was estimated and results compared. Error analysis revealed accuracy of the order of 0.1 degree for the rotation angles of about 1 mm for displacements parallel to the fluoroscopic plane, and of order of 10 mm for the orthogonal displacement.
A 2-D semi-analytical model of double-gate tunnel field-effect transistor
NASA Astrophysics Data System (ADS)
Huifang, Xu; Yuehua, Dai; Ning, Li; Jianbin, Xu
2015-05-01
A 2-D semi-analytical model of double gate (DG) tunneling field-effect transistor (TFET) is proposed. By aid of introducing two rectangular sources located in the gate dielectric layer and the channel, the 2-D Poisson equation is solved by using a semi-analytical method combined with an eigenfunction expansion method. The expression of the surface potential is obtained, which is a special function for the infinite series expressions. The influence of the mobile charges on the potential profile is taken into account in the proposed model. On the basis of the potential profile, the shortest tunneling length and the average electrical field can be derived, and the drain current is then constructed by using Kane's model. In particular, the changes of the tunneling parameters Ak and Bk influenced by the drain—source voltage are also incorporated in the predicted model. The proposed model shows a good agreement with TCAD simulation results under different drain—source voltages, silicon film thicknesses, gate dielectric layer thicknesses, and gate dielectric layer constants. Therefore, it is useful to optimize the DG TFET and this provides a physical insight for circuit level design. Project supported by the National Natural Science Foundation of China (No. 61376106) and the Graduate Innovation Fund of Anhui University.
Heo, Jingu; Savvides, Marios
2012-12-01
In this paper, we propose a novel method for generating a realistic 3D human face from a single 2D face image for the purpose of synthesizing new 2D face images at arbitrary poses using gender and ethnicity specific models. We employ the Generic Elastic Model (GEM) approach, which elastically deforms a generic 3D depth-map based on the sparse observations of an input face image in order to estimate the depth of the face image. Particularly, we show that Gender and Ethnicity specific GEMs (GE-GEMs) can approximate the 3D shape of the input face image more accurately, achieving a better generalization of 3D face modeling and reconstruction compared to the original GEM approach. We qualitatively validate our method using publicly available databases by showing each reconstructed 3D shape generated from a single image and new synthesized poses of the same person at arbitrary angles. For quantitative comparisons, we compare our synthesized results against 3D scanned data and also perform face recognition using synthesized images generated from a single enrollment frontal image. We obtain promising results for handling pose and expression changes based on the proposed method. PMID:22201062
Modeling the Elastic Modulus of 2D Woven CVI SiC Composites
NASA Technical Reports Server (NTRS)
Morscher, Gregory N.
2006-01-01
The use of fiber, interphase, CVI SiC minicomposites as structural elements for 2D-woven SiC fiber reinforced chemically vapor infiltrated (CVI) SiC matrix composites is demonstrated to be a viable approach to model the elastic modulus of these composite systems when tensile loaded in an orthogonal direction. The 0deg (loading direction) and 90deg (perpendicular to loading direction) oriented minicomposites as well as the open porosity and excess SiC associated with CVI SiC composites were all modeled as parallel elements using simple Rule of Mixtures techniques. Excellent agreement for a variety of 2D woven Hi-Nicalon(TradeMark) fiber-reinforced and Sylramic-iBN reinforced CVI SiC matrix composites that differed in numbers of plies, constituent content, thickness, density, and number of woven tows in either direction (i.e, balanced weaves versus unbalanced weaves) was achieved. It was found that elastic modulus was not only dependent on constituent content, but also the degree to which 90deg minicomposites carried load. This depended on the degree of interaction between 90deg and 0deg minicomposites which was quantified to some extent by composite density. The relationships developed here for elastic modulus only necessitated the knowledge of the fractional contents of fiber, interphase and CVI SiC as well as the tow size and shape. It was concluded that such relationships are fairly robust for orthogonally loaded 2D woven CVI SiC composite system and can be implemented by ceramic matrix composite component modelers and designers for modeling the local stiffness in simple or complex parts fabricated with variable constituent contents.
A solidification constitutive model for NIKE2D and NIKE3D
Raboin, P.J.
1994-03-17
This memo updates the current status of a solidification material model development which has been underway for more than a year. Significant modeling goals such as predicting cut-off stresses, thermo-elasto-plasticity, strain rate dependent plasticity and dynamic recovery have been completed. The model is called SOLMAT for solidification material model, and while developed for NIKE2D, it has already been implemented in NIKE3D and NIT03D by B. Maker. This memo details the future development strategy of SOLMAT including liquid and solid constitutive improvements, coupling of deviatoric and dilatational deformation and a plan to switch between constitutive theories. It explains some of the difficulties associated solidification modeling and proposes two experiments to measure properties for using SOLMAT. Due to the sensitive nature of these plans in relation to programmatic and CRADA concerns, this memo should be treated as confidential document.
An Integrative Model of Excitation Driven Fluid Flow in a 2D Uterine Channel
NASA Astrophysics Data System (ADS)
Maggio, Charles; Fauci, Lisa; Chrispell, John
2009-11-01
We present a model of intra-uterine fluid flow in a sagittal cross-section of the uterus by inducing peristalsis in a 2D channel. This is an integrative multiscale computational model that takes as input fluid viscosity, passive tissue properties of the uterine channel and a prescribed wave of membrane depolarization. This voltage pulse is coupled to a model of calcium dynamics inside a uterine smooth muscle cell, which in turn drives a kinetic model of myosin phosphorylation governing contractile muscle forces. Using the immersed boundary method, these muscle forces are communicated to a fluid domain to simulate the contractions which occur in a human uterus. An analysis of the effects of model parameters on the flow properties and emergent geometry of the peristaltic channel will be presented.
Global regularity for the 2D Oldroyd-B model in the corotational case
NASA Astrophysics Data System (ADS)
Ye, Zhuan; Xu, Xiaojing
2016-09-01
This paper is dedicated to the Oldroyd-B model with fractional dissipation $(-\\Delta)^{\\alpha}\\tau$ for any $\\alpha>0$. We establish the global smooth solutions to the Oldroyd-B model in the corotational case with arbitrarily small fractional powers of the Laplacian in two spatial dimensions. The methods described here are quite different from the tedious iterative approach used in recent paper \\cite{XY}. Moreover, in the Appendix we provide some a priori estimates to the Oldroyd-B model in the critical case which may be useful and of interest for future improvement. Finally, the global regularity to to the Oldroyd-B model in the corotational case with $-\\Delta u$ replaced by $(-\\Delta)^{\\gamma}u$ for $\\gamma>1$ are also collected in the Appendix. Therefore our result is more closer to the resolution of the well-known global regularity issue on the critical 2D Oldroyd-B model.
A velocity-dependent anomalous radial transport model for (2-D, 2-V) kinetic transport codes
NASA Astrophysics Data System (ADS)
Bodi, Kowsik; Krasheninnikov, Sergei; Cohen, Ron; Rognlien, Tom
2008-11-01
Plasma turbulence constitutes a significant part of radial plasma transport in magnetically confined plasmas. This turbulent transport is modeled in the form of anomalous convection and diffusion coefficients in fluid transport codes. There is a need to model the same in continuum kinetic edge codes [such as the (2-D, 2-V) transport version of TEMPEST, NEO, and the code being developed by the Edge Simulation Laboratory] with non-Maxwellian distributions. We present an anomalous transport model with velocity-dependent convection and diffusion coefficients leading to a diagonal transport matrix similar to that used in contemporary fluid transport models (e.g., UEDGE). Also presented are results of simulations corresponding to radial transport due to long-wavelength ExB turbulence using a velocity-independent diffusion coefficient. A BGK collision model is used to enable comparison with fluid transport codes.
Evaluation of Hydrus-2D model for solute distribution in subsurface drip
NASA Astrophysics Data System (ADS)
Souza, Claudinei; Bizari, Douglas; Grecco, Katarina
2015-04-01
The competition for water use between agriculture, industry and population has become intense over the years, requiring a rational use of this resource for food production. The subsurface drip irrigation can help producers with the optimization of operating parameters such as frequency and duration of irrigation, flow, spacing and depth of the dripper installation. This information can be obtained by numerical simulations using mathematical models, thus the aim of this study was to evaluate the HYDRUS-2D model from experimental data to predict the size of the wet bulbs generated by emitters of different application rates (1.0 and 1.6 L h-1). The results showed that horizontal displacement (bulb diameter) remained the largest in all the bulbs, observed both in experimental trials and estimated by the model and the correlation between them was high, above 0.90 to below 16% error. We conclude that the HYDRUS-2D model can be used to estimate the dimensions of the wet bulb getting new information on the sizing of the irrigation system.
Exact solution of an anisotropic 2D random walk model with strong memory correlations
NASA Astrophysics Data System (ADS)
Cressoni, J. C.; Viswanathan, G. M.; da Silva, M. A. A.
2013-12-01
Over the last decade, there has been progress in understanding one-dimensional non-Markovian processes via analytic, sometimes exact, solutions. The extension of these ideas and methods to two and higher dimensions is challenging. We report the first exactly solvable two-dimensional (2D) non-Markovian random walk model belonging to the family of the elephant random walk model. In contrast to Lévy walks or fractional Brownian motion, such models incorporate memory effects by keeping an explicit history of the random walk trajectory. We study a memory driven 2D random walk with correlated memory and stops, i.e. pauses in motion. The model has an inherent anisotropy with consequences for its diffusive properties, thereby mixing the dominant regime along one dimension with a subdiffusive walk along a perpendicular dimension. The anomalous diffusion regimes are fully characterized by an exact determination of the Hurst exponent. We discuss the remarkably rich phase diagram, as well as several possible combinations of the independent walks in both directions. The relationship between the exponents of the first and second moments is also unveiled.
Molecular-dynamics of a 2D Model of the Shape Memory Effect
NASA Astrophysics Data System (ADS)
Kastner, Oliver
2006-08-01
This work investigates the thermodynamic properties of a qualitative atomistic model for austenite martensite transitions. The model, still in 2D, employs Lennard-Jones potentials for the determination of the atomic interactions. By use of two atom species it is possible to identify three stable lattice structures in 2D, interpreted as austenite and two variants of martensite. The model is described in the first part of the work [6] in detail. The present work studies the thermodynamic properties of the model concerning a small, 2-dimensional test assembly consisting of 41 atoms. The phase stability is investigated by exploitation of the condition of minimal free energy. The free energy is calculated from the thermal equation of state, which is measured in numerical tensile tests. In the second part of this work a chain of eleven 41-atom assemblies is investigated. The chain is interpreted as an idealized larger body, where the individual crystallites represent crystallographic layers allowing for the creation of micro structure. By use of tensile tests at various temperature conditions we sketch how such chain may exhibit quasi-plasticity, pseudo-elasticity and the shape memory effect.
Momentum Transport: 2D and 3D Cloud Resolving Model Simulations
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo
2001-01-01
The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.
A quasi 2D semianalytical model for the potential profile in hetero and homojunction tunnel FETs
NASA Astrophysics Data System (ADS)
Villani, F.; Gnani, E.; Gnudi, A.; Reggiani, S.; Baccarani, G.
2015-11-01
A quasi 2D semianalytical model for the potential profile in hetero and homojunction tunnel FETs is developed and compared with full-quantum simulation results. It will be shown that the pure analytical solution perfectly matches results at high VDS. However, a coupling with the numerical solution of the 1D Poisson equation in the radial direction is necessary at low VDS, in order to properly account for the charge density in equilibrium with the drain contact. With such an approach we are able to correctly predict the potential profile for both the linear and saturation regimes.
Vector chiral phases in the frustrated 2D XY model and quantum spin chains.
Schenck, H; Pokrovsky, V L; Nattermann, T
2014-04-18
The phase diagram of the frustrated 2D classical and 1D quantum XY models is calculated analytically. Four transitions are found: the vortex unbinding transitions triggered by strong fluctuations occur above and below the chiral transition temperature. Vortex interaction is short range on small and logarithmic on large scales. The chiral transition, though belonging to the Ising universality class by symmetry, has different critical exponents due to nonlocal interaction. In a narrow region close to the Lifshitz point a reentrant phase transition between paramagnetic and quasiferromagnetic phase appears. Applications to antiferromagnetic quantum spin chains and multiferroics are discussed. PMID:24785067
NASA Astrophysics Data System (ADS)
Zou, B.; Li, D. F.; Hu, H. J.; Zhang, H. W.; Lou, L. H.; Chen, M.; Lv, Z. Y.
Based on the verified two dimensional(2D) finite element model for river flow simulation, the effect of estuary training levees on the water flow and sediment movement in the Yellow River estuary is analyzed. For disclosing the effect of setting the two training levees on the flow and sediment motion, the calculation and analysis for the two projects, (one is no levees, the other is setting up two no levees) are given. The results show that when setting up two training levees, water flow is bound by levees and the water flows become more concentrated. As a result, velocity increases in the main channel, sediment carrying capacity of water flow increases correspondingly.
Image restoration using 2D autoregressive texture model and structure curve construction
NASA Astrophysics Data System (ADS)
Voronin, V. V.; Marchuk, V. I.; Petrosov, S. P.; Svirin, I.; Agaian, S.; Egiazarian, K.
2015-05-01
In this paper an image inpainting approach based on the construction of a composite curve for the restoration of the edges of objects in an image using the concepts of parametric and geometric continuity is presented. It is shown that this approach allows to restore the curved edges and provide more flexibility for curve design in damaged image by interpolating the boundaries of objects by cubic splines. After edge restoration stage, a texture restoration using 2D autoregressive texture model is carried out. The image intensity is locally modeled by a first spatial autoregressive model with support in a strongly causal prediction region on the plane. Model parameters are estimated by Yule-Walker method. Several examples considered in this paper show the effectiveness of the proposed approach for large objects removal as well as recovery of small regions on several test images.
Kosaki-Longo index and classification of charges in 2D quantum spin models
NASA Astrophysics Data System (ADS)
Naaijkens, Pieter
2013-08-01
We consider charge superselection sectors of two-dimensional quantum spin models corresponding to cone localisable charges, and prove that the number of equivalence classes of such charges is bounded by the Kosaki-Longo index of an inclusion of certain observable algebras. To demonstrate the power of this result we apply the theory to the toric code on a 2D infinite lattice. For this model we can compute the index of this inclusion, and conclude that there are four distinct irreducible charges in this model, in accordance with the analysis of the toric code model on compact surfaces. We also give a sufficient criterion for the non-degeneracy of the charge sectors, in the sense that Verlinde's matrix S is invertible.
Accelerating numerical modeling of wave propagation through 2-D anisotropic materials using OpenCL.
Molero, Miguel; Iturrarán-Viveros, Ursula
2013-03-01
We present an implementation of the numerical modeling of elastic waves propagation, in 2D anisotropic materials, using the new parallel computing devices (PCDs). Our study is aimed both to model laboratory experiments and explore the capabilities of the emerging PCDs by discussing performance issues. In the experiments a sample plate of an anisotropic material placed inside a water tank is rotated and, for every angle of rotation it is subjected to an ultrasonic wave (produced by a large source transducer) that propagates in the water and through the material producing some reflection and transmission signals that are recording by a "point-like" receiver. This experiment is numerically modeled by running a finite difference code covering a set of angles θ∈[-50°, 50°], and recorded the signals for the transmission and reflection results. Transversely anisotropic and weakly orthorhombic materials are considered. We accelerated the computation using an open-source toolkit called PyOpenCL, which lets one to easily access the OpenCL parallel computation API's from the high-level programming environment of Python. A speedup factor over 19 using the GPU is obtained when compared with the execution of the same program in parallel using a CPU multi-core (in this case we use the 4-cores that has the CPU). The performance for different graphic cards and operating systems is included together with the full 2-D finite difference code with PyOpenCL. PMID:23290584
Self-Organization in 2D Traffic Flow Model with Jam-Avoiding Drive
NASA Astrophysics Data System (ADS)
Nagatani, Takashi
1995-04-01
A stochastic cellular automaton (CA) model is presented to investigate the traffic jam by self-organization in the two-dimensional (2D) traffic flow. The CA model is the extended version of the 2D asymmetric exclusion model to take into account jam-avoiding drive. Each site contains either a car moving to the up, a car moving to the right, or is empty. A up car can shift right with probability p ja if it is blocked ahead by other cars. It is shown that the three phases (the low-density phase, the intermediate-density phase and the high-density phase) appear in the traffic flow. The intermediate-density phase is characterized by the right moving of up cars. The jamming transition to the high-density jamming phase occurs with higher density of cars than that without jam-avoiding drive. The jamming transition point p 2c increases with the shifting probability p ja. In the deterministic limit of p ja=1, it is found that a new jamming transition occurs from the low-density synchronized-shifting phase to the high-density moving phase with increasing density of cars. In the synchronized-shifting phase, all up cars do not move to the up but shift to the right by synchronizing with the move of right cars. We show that the jam-avoiding drive has an important effect on the dynamical jamming transition.
Stochastic dynamics of phase singularities under ventricular fibrillation in 2D Beeler-Reuter model
NASA Astrophysics Data System (ADS)
Suzuki, Akio; Konno, Hidetoshi
2011-09-01
The dynamics of ventricular fibrillation (VF) has been studied extensively, and the initiation mechanism of VF has been elucidated to some extent. However, the stochastic dynamical nature of sustained VF remains unclear so far due to the complexity of high dimensional chaos in a heterogeneous system. In this paper, various statistical mechanical properties of sustained VF are studied numerically in 2D Beeler-Reuter-Drouhard-Roberge (BRDR) model with normal and modified ionic current conductance. The nature of sustained VF is analyzed by measuring various fluctuations of spatial phase singularity (PS) such as velocity, lifetime, the rates of birth and death. It is found that the probability density function (pdf) for lifetime of PSs is independent of system size. It is also found that the hyper-Gamma distribution serves as a universal pdf for the counting number of PSs for various system sizes and various parameters of our model tissue under VF. Further, it is demonstrated that the nonlinear Langevin equation associated with a hyper-Gamma process can mimic the pdf and temporal variation of the number of PSs in the 2D BRDR model.
NASA Astrophysics Data System (ADS)
Yan, Bo; Li, Yuguo; Liu, Ying
2016-07-01
In this paper, we present an adaptive finite element (FE) algorithm for direct current (DC) resistivity modeling in 2-D generally anisotropic conductivity structures. Our algorithm is implemented on an unstructured triangular mesh that readily accommodates complex structures such as topography and dipping layers and so on. We implement a self-adaptive, goal-oriented grid refinement algorithm in which the finite element analysis is performed on a sequence of refined grids. The grid refinement process is guided by an a posteriori error estimator. The problem is formulated in terms of total potentials where mixed boundary conditions are incorporated. This type of boundary condition is superior to the Dirichlet type of conditions and improves numerical accuracy considerably according to model calculations. We have verified the adaptive finite element algorithm using a two-layered earth with azimuthal anisotropy. The FE algorithm with incorporation of mixed boundary conditions achieves high accuracy. The relative error between the numerical and analytical solutions is less than 1% except in the vicinity of the current source location, where the relative error is up to 2.4%. A 2-D anisotropic model is used to demonstrate the effects of anisotropy upon the apparent resistivity in DC soundings.
Qu, Qiang; Qu, Jian; Han, Lu; Zhan, Min; Wu, Lan-xiang; Zhang, Yi-wen; Zhang, Wei; Zhou, Hong-hao
2014-01-01
Aim: Herbal products have been widely used, and the safety of herb-drug interactions has aroused intensive concerns. This study aimed to investigate the effects of phytochemicals on the catalytic activities of human CYP2D6*1 and CYP2D6*10 in vitro. Methods: HepG2 cells were stably transfected with CYP2D6*1 and CYP2D6*10 expression vectors. The metabolic kinetics of the enzymes was studied using HPLC and fluorimetry. Results: HepG2-CYP2D6*1 and HepG2-CYP2D6*10 cell lines were successfully constructed. Among the 63 phytochemicals screened, 6 compounds, including coptisine sulfate, bilobalide, schizandrin B, luteolin, schizandrin A and puerarin, at 100 μmol/L inhibited CYP2D6*1- and CYP2D6*10-mediated O-demethylation of a coumarin compound AMMC by more than 50%. Furthermore, the inhibition by these compounds was dose-dependent. Eadie-Hofstee plots demonstrated that these compounds competitively inhibited CYP2D6*1 and CYP2D6*10. However, their Ki values for CYP2D6*1 and CYP2D6*10 were very close, suggesting that genotype-dependent herb-drug inhibition was similar between the two variants. Conclusion: Six phytochemicals inhibit CYP2D6*1 and CYP2D6*10-mediated catalytic activities in a dose-dependent manner in vitro. Thus herbal products containing these phytochemicals may inhibit the in vivo metabolism of co-administered drugs whose primary route of elimination is CYP2D6. PMID:24786236
Modeling and Control of 2-D Grasping of an Object with Arbitrary Shape under Rolling Contact
NASA Astrophysics Data System (ADS)
Arimoto, Suguru; Yoshida, Morio; Sekimoto, Masahiro; Tahara, Kenji
Modeling, control, and stabilization of dynamics of two-dimensional object grasping by using a pair of multi-joint robot fingers are investigated under rolling contact constraints and an arbitrary geometry of the object and fingertips. First, a fundamental testbed problem of modeling and control of rolling motion between 2-D rigid bodies with an arbitrary shape is treated under the assumption that the two contour curves coincide at the contact point and share the same tangent. The rolling constraint induces the Euler equation of motion that is parameterized by a common arclength parameter and constrained onto the kernel space orthogonally complemented to the image space spanned from the constraint gradient. By extending the analysis to the problem of stable grasp of a 2-D object with an arbitrary shape by a pair of robot fingers, the Euler-Lagrange equation of motion of the overall fingers/object system parametrized by arclength parameters is derived, together with a couple of first-order differential equations that express evolutions of contact points in terms of the second fundamental form. It is shown that 2-D rolling constraints are integrable in the sense of Frobonius even if their Pfaffian forms are characterized by arclength parameters. A control signal called “blind grasping” is introduced and shown to be effective in stabilization of grasping without using the details of the object shape and parameters or external sensing. An extension of the Dirichlet-Lagrange stability theorem to a class of systems with DOF-redundancy under constraints is suggested by using a Morse-Bott-Lyapunov function.
Uncertainties in modelling Mt. Pinatubo eruption with 2-D AER model and CCM SOCOL
NASA Astrophysics Data System (ADS)
Kenzelmann, P.; Weisenstein, D.; Peter, T.; Luo, B. P.; Rozanov, E.; Fueglistaler, S.; Thomason, L. W.
2009-04-01
Large volcanic eruptions may introduce a strong forcing on climate. They challenge the skills of climate models. In addition to the short time attenuation of solar light by ashes the formation of stratospheric sulphate aerosols, due to volcanic sulphur dioxide injection into the lower stratosphere, may lead to a significant enhancement of the global albedo. The sulphate aerosols have a residence time of about 2 years. As a consequence of the enhanced sulphate aerosol concentration both the stratospheric chemistry and dynamics are strongly affected. Due to absorption of longwave and near infrared radiation the temperature in the lower stratosphere increases. So far chemistry climate models overestimate this warming [Eyring et al. 2006]. We present an extensive validation of extinction measurements and model runs of the eruption of Mt. Pinatubo in 1991. Even if Mt. Pinatubo eruption has been the best quantified volcanic eruption of this magnitude, the measurements show considerable uncertainties. For instance the total amount of sulphur emitted to the stratosphere ranges from 5-12 Mt sulphur [e.g. Guo et al. 2004, McCormick, 1992]. The largest uncertainties are in the specification of the main aerosol cloud. SAGE II, for instance, could not measure the peak of the aerosol extinction for about 1.5 years, because optical termination was reached. The gap-filling of the SAGE II [Thomason and Peter, 2006] using lidar measurements underestimates the total extinctions in the tropics for the first half year after the eruption by 30% compared to AVHRR [Rusell et. al 1992]. The same applies to the optical dataset described by Stenchikov et al. [1998]. We compare these extinction data derived from measurements with extinctions derived from AER 2D aerosol model calculations [Weisenstein et al., 2007]. Full microphysical calculations with injections of 14, 17, 20 and 26 Mt SO2 in the lower stratosphere were performed. The optical aerosol properties derived from SAGE II
Longtime Well-posedness for the 2D Groma-Balogh Model
NASA Astrophysics Data System (ADS)
Wan, Renhui; Chen, Jiecheng
2016-07-01
In this paper, we consider the cauchy problem for the 2D Groma-Balogh model (Acta Mater 47:3647-3654, 1999). From the works Cannone et al. (Arch Ration Mech Anal 196:71-96, 2010) and El Hajj (Ann Inst Henri Poincaré Anal Nonlinéaire 27:21-35, 2010), one can see global well-posedness for this model is an open question. However, we can prove longtime well-posedness. In particular, we show that this model admits a unique solution with the lifespan T^star satisfying T^star log ^2(1+T^star )≳ ɛ ^{-2} if the initial data is of size ɛ . To achieve this, we first establish some new decay estimates concerning the operator e^{-{R}_{12}^2t} . Then, we prove the longtime well-posedness by utilizing the weak dissipation to deal with the nonlinear terms.
Well-posedness and generalized plane waves simulations of a 2D mode conversion model
NASA Astrophysics Data System (ADS)
Imbert-Gérard, Lise-Marie
2015-12-01
Certain types of electro-magnetic waves propagating in a plasma can undergo a mode conversion process. In magnetic confinement fusion, this phenomenon is very useful to heat the plasma, since it permits to transfer the heat at or near the plasma center. This work focuses on a mathematical model of wave propagation around the mode conversion region, from both theoretical and numerical points of view. It aims at developing, for a well-posed equation, specific basis functions to study a wave mode conversion process. These basis functions, called generalized plane waves, are intrinsically based on variable coefficients. As such, they are particularly adapted to the mode conversion problem. The design of generalized plane waves for the proposed model is described in detail. Their implementation within a discontinuous Galerkin method then provides numerical simulations of the process. These first 2D simulations for this model agree with qualitative aspects studied in previous works.
The concept models and implementations of multiport neural net associative memory for 2D patterns
NASA Astrophysics Data System (ADS)
Krasilenko, Vladimir G.; Nikolskyy, Aleksandr I.; Yatskovskaya, Rimma A.; Yatskovsky, Victor I.
2011-04-01
The paper considers neural net models and training and recognizing algorithms with base neurobiologic operations: p-step autoequivalence and non-equivalenc The Modified equivalently models (MEMs) of multiport neural net associative memory (MNNAM) are offered with double adaptive - equivalently weighing (DAEW) for recognition of 2D-patterns (images). It is shown, the computing process in MNNAM under using the proposed MEMs, is reduced to two-step and multi-step algorithms and step-by-step matrix-matrix (tensor-tensor) procedures. The given results of computer simulations confirmed the perspective of such models. Besides the result was received when MNNAM capacity on base of MEMs exceeded the amount of neurons.
Modeling floods in a dense urban area using 2D shallow water equations
NASA Astrophysics Data System (ADS)
Mignot, E.; Paquier, A.; Haider, S.
2006-07-01
SummaryA code solving the 2D shallow water equations by an explicit second-order scheme is used to simulate the severe October 1988 flood in the Richelieu urban locality of the French city of Nîmes. A reference calculation using a detailed description of the street network and of the cross-sections of the streets, considering impervious residence blocks and neglecting the flow interaction with the sewer network provides a mean peak water elevation 0.13 m lower than the measured flood marks with a standard deviation between the measured and computed water depths of 0.53 m. Sensitivity analysis of various topographical and numerical parameters shows that globally, the results keep the same level of accuracy, which reflects both the stability of the calculation method and the smoothening of results. However, the local flow modifications due to change of parameter values can drastically modify the local water depths, especially when the local flow regime is modified. Furthermore, the flow distribution to the downstream parts of the city can be altered depending on the set of parameters used. Finally, a second event, the 2002 flood, was simulated with the calibrated model providing results similar to 1988 flood calculation. Thus, the article shows that, after calibration, a 2D model can be used to help planning mitigation measures in a dense urban area.
Modeling and 2-D discrete simulation of dislocation dynamics for plastic deformation of metal
NASA Astrophysics Data System (ADS)
Liu, Juan; Cui, Zhenshan; Ou, Hengan; Ruan, Liqun
2013-05-01
Two methods are employed in this paper to investigate the dislocation evolution during plastic deformation of metal. One method is dislocation dynamic simulation of two-dimensional discrete dislocation dynamics (2D-DDD), and the other is dislocation dynamics modeling by means of nonlinear analysis. As screw dislocation is prone to disappear by cross-slip, only edge dislocation is taken into account in simulation. First, an approach of 2D-DDD is used to graphically simulate and exhibit the collective motion of a large number of discrete dislocations. In the beginning, initial grains are generated in the simulation cells according to the mechanism of grain growth and the initial dislocation is randomly distributed in grains and relaxed under the internal stress. During the simulation process, the externally imposed stress, the long range stress contribution of all dislocations and the short range stress caused by the grain boundaries are calculated. Under the action of these forces, dislocations begin to glide, climb, multiply, annihilate and react with each other. Besides, thermal activation process is included. Through the simulation, the distribution of dislocation and the stress-strain curves can be obtained. On the other hand, based on the classic dislocation theory, the variation of the dislocation density with time is described by nonlinear differential equations. Finite difference method (FDM) is used to solve the built differential equations. The dislocation evolution at a constant strain rate is taken as an example to verify the rationality of the model.
2D and 3D shape based segmentation using deformable models.
El-Baz, Ayman; Yuksel, Seniha E; Shi, Hongjian; Farag, Aly A; El-Ghar, Mohamed A; Eldiasty, Tarek; Ghoneim, Mohamed A
2005-01-01
A novel shape based segmentation approach is proposed by modifying the external energy component of a deformable model. The proposed external energy component depends not only on the gray level of the images but also on the shape information which is obtained from the signed distance maps of objects in a given data set. The gray level distribution and the signed distance map of the points inside and outside the object of interest are accurately estimated by modelling the empirical density function with a linear combination of discrete Gaussians (LCDG) with positive and negative components. Experimental results on the segmentation of the kidneys from low-contrast DCE-MRI and on the segmentation of the ventricles from brain MRI's show how the approach is accurate in segmenting 2-D and 3-D data sets. The 2D results for the kidney segmentation have been validated by a radiologist and the 3D results of the ventricle segmentation have been validated with a geometrical phantom. PMID:16686036
Interpretation of gravity data using 2-D continuous wavelet transformation and 3-D inverse modeling
NASA Astrophysics Data System (ADS)
Roshandel Kahoo, Amin; Nejati Kalateh, Ali; Salajegheh, Farshad
2015-10-01
Recently the continuous wavelet transform has been proposed for interpretation of potential field anomalies. In this paper, we introduced a 2D wavelet based method that uses a new mother wavelet for determination of the location and the depth to the top and base of gravity anomaly. The new wavelet is the first horizontal derivatives of gravity anomaly of a buried cube with unit dimensions. The effectiveness of the proposed method is compared with Li and Oldenburg inversion algorithm and is demonstrated with synthetics and real gravity data. The real gravity data is taken over the Mobrun massive sulfide ore body in Noranda, Quebec, Canada. The obtained results of the 2D wavelet based algorithm and Li and Oldenburg inversion on the Mobrun ore body had desired similarities to the drill-hole depth information. In all of the inversion algorithms the model non-uniqueness is the challenging problem. Proposed method is based on a simple theory and there is no model non-uniqueness on it.
A 2D Electromechanical Model of Human Atrial Tissue Using the Discrete Element Method
Brocklehurst, Paul; Adeniran, Ismail; Yang, Dongmin; Sheng, Yong; Zhang, Henggui; Ye, Jianqiao
2015-01-01
Cardiac tissue is a syncytium of coupled cells with pronounced intrinsic discrete nature. Previous models of cardiac electromechanics often ignore such discrete properties and treat cardiac tissue as a continuous medium, which has fundamental limitations. In the present study, we introduce a 2D electromechanical model for human atrial tissue based on the discrete element method (DEM). In the model, single-cell dynamics are governed by strongly coupling the electrophysiological model of Courtemanche et al. to the myofilament model of Rice et al. with two-way feedbacks. Each cell is treated as a viscoelastic body, which is physically represented by a clump of nine particles. Cell aggregations are arranged so that the anisotropic nature of cardiac tissue due to fibre orientations can be modelled. Each cell is electrically coupled to neighbouring cells, allowing excitation waves to propagate through the tissue. Cell-to-cell mechanical interactions are modelled using a linear contact bond model in DEM. By coupling cardiac electrophysiology with mechanics via the intracellular Ca2+ concentration, the DEM model successfully simulates the conduction of cardiac electrical waves and the tissue's corresponding mechanical contractions. The developed DEM model is numerically stable and provides a powerful method for studying the electromechanical coupling problem in the heart. PMID:26583141
Locally adaptive 2D-3D registration using vascular structure model for liver catheterization.
Kim, Jihye; Lee, Jeongjin; Chung, Jin Wook; Shin, Yeong-Gil
2016-03-01
Two-dimensional-three-dimensional (2D-3D) registration between intra-operative 2D digital subtraction angiography (DSA) and pre-operative 3D computed tomography angiography (CTA) can be used for roadmapping purposes. However, through the projection of 3D vessels, incorrect intersections and overlaps between vessels are produced because of the complex vascular structure, which makes it difficult to obtain the correct solution of 2D-3D registration. To overcome these problems, we propose a registration method that selects a suitable part of a 3D vascular structure for a given DSA image and finds the optimized solution to the partial 3D structure. The proposed algorithm can reduce the registration errors because it restricts the range of the 3D vascular structure for the registration by using only the relevant 3D vessels with the given DSA. To search for the appropriate 3D partial structure, we first construct a tree model of the 3D vascular structure and divide it into several subtrees in accordance with the connectivity. Then, the best matched subtree with the given DSA image is selected using the results from the coarse registration between each subtree and the vessels in the DSA image. Finally, a fine registration is conducted to minimize the difference between the selected subtree and the vessels of the DSA image. In experimental results obtained using 10 clinical datasets, the average distance errors in the case of the proposed method were 2.34±1.94mm. The proposed algorithm converges faster and produces more correct results than the conventional method in evaluations on patient datasets. PMID:26824922
2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach
NASA Technical Reports Server (NTRS)
Svizhenko, Alexei; Anantram, M. P.; Govindan, T. R.; Biegel, Bryan
2001-01-01
We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.
Be2D: A model to understand the distribution of meteoric 10Be in soilscapes
NASA Astrophysics Data System (ADS)
Campforts, Benjamin; Vanacker, Veerle; Vanderborght, Jan; Govers, Gerard
2016-04-01
Cosmogenic nuclides have revolutionised our understanding of earth surface process rates. They have become one of the standard tools to quantify soil production by weathering, soil redistribution and erosion. Especially Beryllium-10 has gained much attention due to its long half-live and propensity to be relatively conservative in the landscape. The latter makes 10Be an excellent tool to assess denudation rates over the last 1000 to 100 × 103 years, bridging the anthropogenic and geological time scale. Nevertheless, the mobility of meteoric 10Be in soil systems makes translation of meteoric 10Be inventories into erosion and deposition rates difficult. Here we present a coupled soil hillslope model, Be2D, that is applied to synthetic and real topography to address the following three research questions. (i) What is the influence of vertical meteoric Be10 mobility, caused by chemical mobility, clay translocation and bioturbation, on its lateral redistribution over the soilscape, (ii) How does vertical mobility influence erosion rates and soil residence times inferred from meteoric 10Be inventories and (iii) To what extent can a tracer with a half-life of 1.36 Myr be used to distinguish between natural and human-disturbed soil redistribution rates? The model architecture of Be2D is designed to answer these research questions. Be2D is a dynamic model including physical processes such as soil formation, physical weathering, clay migration, bioturbation, creep, overland flow and tillage erosion. Pathways of meteoric 10Be mobility are simulated using a two step approach which is updated each timestep. First, advective and diffusive mobility of meteoric 10Be is simulated within the soil profile and second, lateral redistribution because of lateral soil fluxes is calculated. The performance and functionality of the model is demonstrated through a number of synthetic and real model runs using existing datasets of meteoric 10Be from case-studies in southeastern US. Brute
Estimating nitrogen losses in furrow irrigated soil amended by compost using HYDRUS-2D model
NASA Astrophysics Data System (ADS)
Iqbal, Shahid; Guber, Andrey; Zaman Khan, Haroon; ullah, Ehsan
2014-05-01
Furrow irrigation commonly results in high nitrogen (N) losses from soil profile via deep infiltration. Estimation of such losses and their reduction is not a trivial task because furrow irrigation creates highly nonuniform distribution of soil water that leads to preferential water and N fluxes in soil profile. Direct measurements of such fluxes are impractical. The objective of this study was to assess applicability of HYDRUS-2D model for estimating nitrogen balance in manure amended soil under furrow irrigation. Field experiments were conducted in a sandy loam soil amended by poultry manure compost (PMC) and pressmud compost (PrMC) fertilizers. The PMC and PrMC contained 2.5% and 0.9% N and were applied at 5 rates: 2, 4, 6, 8 and 10 ton/ha. Plots were irrigated starting from 26th day from planting using furrows with 1x1 ridge to furrow aspect ratio. Irrigation depths were 7.5 cm and time interval between irrigations varied from 8 to 15 days. Results of the field experiments showed that approximately the same corn yield was obtained with considerably higher N application rates using PMC than using PrMC as a fertilizer. HYDRUS-2D model was implemented to evaluate N fluxes in soil amended by PMC and PrMC fertilizers. Nitrogen exchange between two pools of organic N (compost and soil) and two pools of mineral N (soil NH4-N and soil NO3-N) was modeled using mineralization and nitrification reactions. Sources of mineral N losses from soil profile included denitrification, root N uptake and leaching with deep infiltration of water. HYDRUS-2D simulations showed that the observed increases in N root water uptake and corn yields associated with compost application could not be explained by the amount of N added to soil profile with the compost. Predicted N uptake by roots significantly underestimated the field data. Good agreement between simulated and field-estimated values of N root uptake was achieved when the rate of organic N mineralization was increased
The dynamics of barotropic vortex merging
NASA Astrophysics Data System (ADS)
Kieu, Chanh
2016-08-01
The merging of multiple vortices is a fundamental process of the dynamics of Earth's atmosphere and oceans. In this study, the interaction of like-signed vortices is analytically and numerically examined in a framework of two-dimensional inviscid barotropic flows. It is shown that barotropic vortex interaction turns out to be more intricate than simple merging scenarios often assumed in previous studies. Some particular configurations exist in which the vortex merging process is never complete despite strong interaction of like-signed vortices, regardless of the strengths or distances between the vortices. While the conditions for a complete vortex merging process introduced in this study appear to be too strict for most practical applications, this study suggests that careful criteria for vortex mergers should be properly defined when simulating the interaction of vortices, because the merging may not always result in a final enhanced circulation at the end of the interaction, as usually assumed in the literature.
NASA Astrophysics Data System (ADS)
Bezzeccheri, E.; Colasanti, S.; Falco, A.; Liguori, R.; Rubino, A.; Lugli, P.
2016-05-01
Vertical Organic Transistors and Phototransistors have been proven to be promising technologies due to the advantages of reduced channel length and larger sensitive area with respect to planar devices. Nevertheless, a real improvement of their performance is subordinate to the quantitative description of their operation mechanisms. In this work, we present a comparative study on the modeling of vertical and planar Organic Phototransistor (OPT) structures. Computer-based simulations of the devices have been carried out with Synopsys Sentaurus TCAD in a 2D Drift-Diffusion framework. The photoactive semiconductor material has been modeled using the virtual semiconductor approach as the archetypal P3HT:PC61BM bulk heterojunction. It has been found that both simulated devices have comparable electrical and optical characteristics, accordingly to recent experimental reports on the subject.
An investigation of DTNS2D for use as an incompressible turbulence modelling test-bed
NASA Technical Reports Server (NTRS)
Steffen, Christopher J., Jr.
1992-01-01
This paper documents an investigation of a two dimensional, incompressible Navier-Stokes solver for use as a test-bed for turbulence modelling. DTNS2D is the code under consideration for use at the Center for Modelling of Turbulence and Transition (CMOTT). This code was created by Gorski at the David Taylor Research Center and incorporates the pseudo compressibility method. Two laminar benchmark flows are used to measure the performance and implementation of the method. The classical solution of the Blasius boundary layer is used for validating the flat plate flow, while experimental data is incorporated in the validation of backward facing step flow. Velocity profiles, convergence histories, and reattachment lengths are used to quantify these calculations. The organization and adaptability of the code are also examined in light of the role as a numerical test-bed.
Singularities of the Partition Function for the Ising Model Coupled to 2D Quantum Gravity
NASA Astrophysics Data System (ADS)
Ambjørn, J.; Anagnostopoulos, K. N.; Magnea, U.
We study the zeros in the complex plane of the partition function for the Ising model coupled to 2D quantum gravity for complex magnetic field and real temperature, and for complex temperature and real magnetic field, respectively. We compute the zeros by using the exact solution coming from a two-matrix model and by Monte-Carlo simulations of Ising spins on dynamical triangulations. We present evidence that the zeros form simple one-dimensional curves in the complex plane, and that the critical behaviour of the system is governed by the scaling of the distribution of the singularities near the critical point. Despite the small size of the systems studied, we can obtain a reasonable estimate of the (known) critical exponents.
A guide to using material model No. 11 in NIKE2D: An internal variable, viscoplasticity model
Flower, E.C.; Nikkel, D.J. Jr.
1990-10-30
The need to accurately model the superplastic forming process which is highly rate and temperature dependent motivated the evaluation of Bammann's internal variable, viscoplasticity material model. The model is based upon the concepts of unified creep plasticity, but employs a yield surface for efficient implementation into large-scale numerical computer codes. It has proven elsewhere to be quite successful in describing large strain, thermal-mechanical behavior of crystalline materials. Features of the model enable it to simulate the apparent strain-rate behavior exhibited by many metals above one half the melt temperature. It is the efficient incorporation of features that make the model attractive for use in finite element modeling of metal deformation processes. Although this model was implemented into the Lawrence Livermore National Laboratory's NIKE2D finite element program in 1986, there have been no known reports of successful use by NIKE2D users. The purpose of this report is to provide the user the proper format to input model parameters, a procedure for determining appropriate values for material constants from experimental data, and supplemental information on the model relevant to the implementation in the NIKE2D finite element program. Detailed accounts of the theoretical aspects of the model can be found in the cited references. 4 refs., 8 figs.
Neutrino-electron Scattering in 2-D Models of Supernova Convection
NASA Astrophysics Data System (ADS)
DeNisco, K. R.; Swesty, F. D.; Myra, E. S.
2005-12-01
We present results from 2-D supernova simulations which include the effects of neutrino-electron scattering. The importance of neutrino-electron scattering in stellar collapse has been known for two decades. Yet it has often been neglected in multidimensional simulations due to the difficulty of implementing it consistently. The inclusion of this process is numerically challenging because of the extremely short scattering timescales involved. The stiffness resulting from this short timescale precludes an explicit numerical treatment of this phenomenon, such as those that have recently been utilized in some 2-D models. We describe our fully-implicit treatment of this process and present our initial results. This work was performed at the State University of New York at Stony Brook as part of the TeraScale Supernova Initiative, and is funded by SciDAC grant DE-FC02-01ER41185 from the U.S. Department of Energy, Office of Science High-Energy, Nuclear, and Advanced Scientific Computing Research Programs. We gratefully acknowledge support of the National Energy Research Scientific Computing Center (NERSC) for computational and consulting support.
Implications of lack-of-ergodicity in 2D Potts model
NASA Astrophysics Data System (ADS)
Ota, Smita
2015-03-01
Microcanonical Monte Carlo simulation is used to study two dimensional (2D) q state Potts model. We consider a 2D square lattice having NxN spins with periodic boundary condition and simulated the system with N =15 and q =10. The demon energy distribution is found to be exponential for high system energy and large system size. For smaller system size and above the first order transition the demon energy distribution is found to deviate from exp(- βED) and has the form exp(- βED + γ ED2). Here β = 1/kBT and kB is the Boltzmann constant. It is found that γ is finite at higher temperatures. As the system energy is reduced γ becomes zero near the first order transition. It is found that during cooling γ changes sign from negative to positive and then to negative again near the 1st order transition. Therefore the demon energy distribution becomes exp(- βED) (or ergodic) at two values of system energy near the 1st order transition. Further cooling or at still lower temperatures the system shows lack of ergodicity. However, difference in heating cooling curves are apparent in E vs γ. The system energies for which γ is zero during cooling can represent the 'ergodic' states. This can be related to the two-level systems observed in glasses at low temperatures.
Nutter, C.
1980-11-01
GRAV2D is an interactive computer program used for modeling 2-1/2 dimensional gravity data. A forward algorithm is used to give the theoretical attraction of gravity intensity at a station due to a perturbing body given by the initial model. The resultant model can then be adjusted for a better fit by a combination of manual adjustment, one-dimensional automatic search, and Marquardt inversion. GRAV2D has an interactive data management system for data manipulation and display built around subroutines to do a forward problem, a one-dimensional direct search and an inversion. This is a user's guide and documentation for GRAV2D.
Instability of magnetic equilibria in barotropic stars
NASA Astrophysics Data System (ADS)
Mitchell, J. P.; Braithwaite, J.; Reisenegger, A.; Spruit, H.; Valdivia, J. A.; Langer, N.
2015-02-01
In stably stratified stars, numerical magnetohydrodynamics simulations have shown that arbitrary initial magnetic fields evolve into stable equilibrium configurations, usually containing nearly axisymmetric, linked poloidal and toroidal fields that stabilize each other. In this work, we test the hypothesis that stable stratification is a requirement for the existence of such stable equilibria. For this purpose, we follow numerically the evolution of magnetic fields in barotropic (and thus neutrally stable) stars, starting from two different types of initial conditions, namely random disordered magnetic fields, as well as linked poloidal-toroidal configurations resembling the previously found equilibria. With many trials, we always find a decay of the magnetic field over a few Alfvén times, never a stable equilibrium. This strongly suggests that there are no stable equilibria in barotropic stars, thus clearly invalidating the assumption of barotropic equations of state often imposed on the search of magnetic equilibria. It also supports the hypothesis that, as dissipative processes erode the stable stratification, they might destabilize previously stable magnetic field configurations, leading to their decay.
2D time-domain finite-difference modeling for viscoelastic seismic wave propagation
NASA Astrophysics Data System (ADS)
Fan, Na; Zhao, Lian-Feng; Xie, Xiao-Bi; Ge, Zengxi; Yao, Zhen-Xing
2016-07-01
Real Earth media are not perfectly elastic. Instead, they attenuate propagating mechanical waves. This anelastic phenomenon in wave propagation can be modeled by a viscoelastic mechanical model consisting of several standard linear solids. Using this viscoelastic model, we approximate a constant Q over a frequency band of interest. We use a four-element viscoelastic model with a tradeoff between accuracy and computational costs to incorporate Q into 2D time-domain first-order velocity-stress wave equations. To improve the computational efficiency, we limit the Q in the model to a list of discrete values between 2 and 1000. The related stress and strain relaxation times that characterize the viscoelastic model are pre-calculated and stored in a database for use by the finite-difference calculation. A viscoelastic finite-difference scheme that is second-order in time and fourth-order in space is developed based on the MacCormack algorithm. The new method is validated by comparing the numerical result with analytical solutions that are calculated using the generalized reflection/transmission coefficient method. The synthetic seismograms exhibit greater than 95 per cent consistency in a two-layer viscoelastic model. The dispersion generated from the simulation is consistent with the Kolsky-Futterman dispersion relationship.
Transforming 2d Cadastral Data Into a Dynamic Smart 3d Model
NASA Astrophysics Data System (ADS)
Tsiliakou, E.; Labropoulos, T.; Dimopoulou, E.
2013-08-01
3D property registration has become an imperative need in order to optimally reflect all complex cases of the multilayer reality of property rights and restrictions, revealing their vertical component. This paper refers to the potentials and multiple applications of 3D cadastral systems and explores the current state-of-the art, especially the available software with which 3D visualization can be achieved. Within this context, the Hellenic Cadastre's current state is investigated, in particular its data modeling frame. Presenting the methodologies and specifications addressing the registration of 3D properties, the operating cadastral system's shortcomings and merits are pointed out. Nonetheless, current technological advances as well as the availability of sophisticated software packages (proprietary or open source) call for 3D modeling. In order to register and visualize the complex reality in 3D, Esri's CityEngine modeling software has been used, which is specialized in the generation of 3D urban environments, transforming 2D GIS Data into Smart 3D City Models. The application of the 3D model concerns the Campus of the National Technical University of Athens, in which a complex ownership status is established along with approved special zoning regulations. The 3D model was built using different parameters based on input data, derived from cadastral and urban planning datasets, as well as legal documents and architectural plans. The process resulted in a final 3D model, optimally describing the cadastral situation and built environment and proved to be a good practice example of 3D visualization.
NASA Astrophysics Data System (ADS)
Mo, Yike; Greenhalgh, Stewart A.; Robertsson, Johan O. A.; Karaman, Hakki
2015-05-01
Lateral velocity variations and low velocity near-surface layers can produce strong scattered and guided waves which interfere with reflections and lead to severe imaging problems in seismic exploration. In order to investigate these specific problems by laboratory seismic modelling, a simple 2D ultrasonic model facility has been recently assembled within the Wave Propagation Lab at ETH Zurich. The simulated geological structures are constructed from 2 mm thick metal and plastic sheets, cut and bonded together. The experiments entail the use of a piezoelectric source driven by a pulse amplifier at ultrasonic frequencies to generate Lamb waves in the plate, which are detected by piezoelectric receivers and recorded digitally on a National Instruments recording system, under LabVIEW software control. The 2D models employed were constructed in-house in full recognition of the similitude relations. The first heterogeneous model features a flat uniform low velocity near-surface layer and deeper dipping and flat interfaces separating different materials. The second model is comparable but also incorporates two rectangular shaped inserts, one of low velocity, the other of high velocity. The third model is identical to the second other than it has an irregular low velocity surface layer of variable thickness. Reflection as well as transmission experiments (crosshole & vertical seismic profiling) were performed on each model. The two dominant Lamb waves recorded are the fundamental symmetric mode (non-dispersive) and the fundamental antisymmetric (flexural) dispersive mode, the latter normally being absent when the source transducer is located on a model edge but dominant when it is on the flat planar surface of the plate. Experimental group and phase velocity dispersion curves were determined and plotted for both modes in a uniform aluminium plate. For the reflection seismic data, various processing techniques were applied, as far as pre-stack Kirchhoff migration. The
Turbulence modeling for subsonic separated flows over 2-D airfoils and 3-D wings
NASA Astrophysics Data System (ADS)
Rosen, Aaron M.
Accurate predictions of turbulent boundary layers and flow separation through computational fluid dynamics (CFD) are becoming more and more essential for the prediction of loads in the design of aerodynamic flight components. Standard eddy viscosity models used in many commercial codes today do not capture the nonequilibrium effects seen in a separated flow and thus do not generally make accurate separation predictions. Part of the reason for this is that under nonequilibrium conditions such as a strong adverse pressure gradient, the history effects of the flow play an important role in the growth and decay of turbulence. More recent turbulence models such as Olsen and Coakley's Lag model and Lillard's lagRST model seek to simulate these effects by lagging the turbulent variables when nonequilibrium effects become important. The purpose of the current research is to assess how these nonequilibrium turbulence models capture the separated regions on various 2-D airfoils and 3-D wings. Nonequilibrium models including the Lag model and the lagRST model are evaluated in comparison with three baseline models (Spalart-Allmaras, Wilcox's k-omega, and Menter's SST) using a modified version of the OVERFLOW code. Tuning the model coefficients of the Lag and lagRST models is also explored. Results show that the various lagRST formulations display an improvement in velocity profile predictions over the standard RANS models, but have trouble capturing the edge of the boundary layer. Experimental separation location measurements were not available, but several trends are noted which may be useful to tuning the model coefficients in the future.
2D-photochemical modeling of Saturn’s stratosphere: hydrocarbon and water distributions
NASA Astrophysics Data System (ADS)
Hue, Vincent; Cavalié, Thibault; Hersant, Franck; Dobrijevic, Michel; Greathouse, Thomas; Lellouch, Emmanuel; Hartogh, Paul; Cassidy, Timothy; Spiga, Aymeric; Guerlet, Sandrine; Sylvestre, Melody
2014-11-01
Saturn’s axial tilt of 27° produces seasons in a similar way as on Earth. The seasonal forcing over Saturn’s 30 years period influences the production/loss of the major atmospheric absorbers and coolants through photochemistry, and influences therefore Saturn’s stratospheric temperatures. We have developed a 2D time-dependent photochemical model of Saturn’s atmosphere [Hue et al., in prep.], coupled to a radiative-climate model [Greathouse et al., 2008] to study seasonal effects on its atmospheric composition. Cassini spacecraft has revealed that the distribution of hydrocarbons in Saturn’s stratosphere [Guerlet et al., 2009] differs from pure photochemical predictions, i.e. without meridional transport [Moses et al., 2005]. Differences between the observed distribution of hydrocarbons and 2D-photochemical predictions are likely to be an indicator of dynamical forcing.Disentangling the origin of water in the stratosphere of this planet has been a long-term issue. Due to Saturn’s cold tropopause trap, which acts as a transport barrier, the water vapor observed by the Infrared Space Observatory (ISO) [Feuchtgruber et al., 1997] has an external origin. Three external sources have been identified: (i) permanent flux from interplanetary dust particles, (ii) local sources form planetary environments (rings, satellites), (iii) large cometary impacts, similar to Shoemaker-Levy 9 on Jupiter. Previous observations of Saturn with Herschel’s Hsso program [Hartogh et al., 2009] led to the detection of a water torus around Saturn [Hartogh et al., 2011], fed by Enceladus’ geysers. A substantial fraction of this torus is predicted to be a local source of water for Saturn’s and its satellites, as it will spread in this system [Cassidy et al., 2010]. Using the new 2D-photochemical model, we test here the validity of Enceladus’ torus as the source of Saturn’s stratospheric water.References : Hue et al., in prep. Greathouse et al., 2008. AGU Fall Meeting
On the assimilation of flood extension images into 2D shallow-water models
NASA Astrophysics Data System (ADS)
Monnier, J.; Couderc, F.; Dartus, D.; Madec, R.; Vila, J.
2012-12-01
In river hydraulics, assimilation of water level measurements at gauging stations is well controlled, while assimilation of images (e.g. from satellite) is still delicate. In the present talk, we address the richness of satellite information to constraint a 2D shallow-water model, and present also related difficulties. A preliminary study done on Mosel river is presented in [LaMo] [HoLaMoPu]. On selected parts of the image, an 0th order model flow allows to obtain some reliable water levels with quantified uncertainties (C. Puech et al.). Next, variationnal sensitivities (based on a gradient computation and adjoint equations) reveal some difficulties that a model designer have to tackle (e.g. roughness parameters at open boundaries), and allow to better understand both the model and the flow. Next, a variational data assimilation algorithm (4D-var) shows that such data lead to a better calibration of the model (e.g. roughness coefficients) and potentially allows to identify the incoming and/or outgoing flow at open boundaries, [LaMo] [HoLaMoPu]. On the other side, the flood dynamic extension is difficult to represent accurately using a 2D SW model since the wet-dry front dynamics is difficult to compute. We compare some 2nd order finite volume solvers and obtain an accurate and stable scheme at wet-dry front. Then, we present some basic rules of compatibility between data and mesh resolution in order to be reliable enough to constraint the model with flood extension data, [CoMaMoViDa]. All the algorithms are implemented into DassFlow software (Fortran, MPI, adjoint) [Da]. [CoMaMoViDa] F. Couderc, R. Madec, J. Monnier, J.-P. Vila, D. Dartus. "Sensitivity analysis and variational data assimilation for geophysical shallow water flows". Submitted. [Da] DassFlow - Data Assimilation for Free Surface Flows. Open-source computational software http://www-gmm.insa-toulouse.fr/~monnier/DassFlow/ [HoLaMoPu] R. Hostache, X. Lai, J. Monnier, C. Puech. "Assimilation of spatial
The combined effect of attraction and orientation zones in 2D flocking models
NASA Astrophysics Data System (ADS)
Iliass, Tarras; Cambui, Dorilson
2016-01-01
In nature, many animal groups, such as fish schools or bird flocks, clearly display structural order and appear to move as a single coherent entity. In order to understand the complex motion of these systems, we study the Vicsek model of self-propelled particles (SPP) which is an important tool to investigate the behavior of collective motion of live organisms. This model reproduces the biological behavior patterns in the two-dimensional (2D) space. Within the framework of this model, the particles move with the same absolute velocity and interact locally in the zone of orientation by trying to align their direction with that of the neighbors. In this paper, we model the collective movement of SPP using an agent-based model which follows biologically motivated behavioral rules, by adding a second region called the attraction zone, where each particles move towards each other avoiding being isolated. Our main goal is to present a detailed numerical study on the effect of the zone of attraction on the kinetic phase transition of our system. In our study, the consideration of this zone seems to play an important role in the cohesion. Consequently, in the directional orientation, the zone that we added forms the compact particle group. In our simulation, we show clearly that the model proposed here can produce two collective behavior patterns: torus and dynamic parallel group. Implications of these findings are discussed.
Coronary arteries motion modeling on 2D x-ray images
NASA Astrophysics Data System (ADS)
Gao, Yang; Sundar, Hari
2012-02-01
During interventional procedures, 3D imaging modalities like CT and MRI are not commonly used due to interference with the surgery and radiation exposure concerns. Therefore, real-time information is usually limited and building models of cardiac motion are difficult. In such case, vessel motion modeling based on 2-D angiography images become indispensable. Due to issues with existing vessel segmentation algorithms and the lack of contrast in occluded vessels, manual segmentation of certain branches is usually necessary. In addition, such occluded branches are the most important vessels during coronary interventions and obtaining motion models for these can greatly help in reducing the procedure time and radiation exposure. Segmenting different cardiac phases independently does not guarantee temporal consistency and is not efficient for occluded branches required manual segmentation. In this paper, we propose a coronary motion modeling system which extracts the coronary tree for every cardiac phase, maintaining the segmentation by tracking the coronary tree during the cardiac cycle. It is able to map every frame to the specific cardiac phase, thereby inferring the shape information of the coronary arteries using the model corresponding to its phase. Our experiments show that our motion modeling system can achieve promising results with real-time performance.
2-D modeling of laterally acoustically coupled thin film bulk acoustic wave resonator filters.
Pensala, Tuomas; Meltaus, Johanna; Kokkonen, Kimmo; Ylilammi, Markku
2010-11-01
A 2-D model is developed for calculating lateral acoustical coupling between adjacent thin film BAW resonators forming an electrical N-port. The model is based on solution and superposition of lateral eigenmodes and eigenfrequencies in a structure consisting of adjacent regions with known plate wave dispersion properties. Mechanical and electrical response of the device are calculated as a superposition of eigenmodes according to voltage drive at one electrical port at a time while extracting current induced in the other ports, leading to a full Y-parameter description of the device. Exemplary cases are simulated to show the usefulness of the model in the study of the basic design rules of laterally coupled thin film BAW resonator filters. Model predictions are compared to an experimental 1.9-GHz band-pass filter based on aluminum nitride thin film technology and lateral acoustical coupling. Good agreement is obtained in prediction of passband behavior. The eigenmode-based model forms a useful tool for fast simulation of laterally coupled acoustic devices. It allows one to gain insight into basic device physics in a very intuitive fashion compared with more detailed but heavier finite element method. Shortcomings of this model and possible improvements are discussed. PMID:21041141
Models Ion Trajectories in 2D and 3D Electrostatic and Magnetic Fields
Energy Science and Technology Software Center (ESTSC)
2000-02-21
SIMION3D7.0REV is a C based ion optics simulation program that can model complex problems using Laplace equation solutions for potential fields. The program uses an ion optics workbench that can hold up to 200 2D and/or 3D electrostatic/magnetic potential arrays. Arrays can have up to 50,000,000 points. SIMION3D7.0''s 32 bit virtual Graphics User Interface provides a highly interactive advanced user environment. All potential arrays are visualized as 3D objects that the user can cut awaymore » to inspect ion trajectories and potential energy surfaces. User programs allow the user to customize the program for specific simulations. A geometry file option supports the definition of highly complex array geometry. Algorithm modifications have improved this version''s computational speed and accuracy.« less
Surface delta interaction in the g7/2 - d5/2 model space
NASA Astrophysics Data System (ADS)
Yu, Xiaofei; Zamick, Larry
2016-05-01
Using an attractive surface delta interaction we obtain wave functions for 2 neutrons (or neutron holes) in the g7/2 -d5/2 model space. If we take the single particle energies to be degenerate we find that the g factors for I = 2 , 4 and 6 are all the same G (J) =gl, the orbital g factor of the nucleon. For a free neutron gl = 0, so in this case all 2 particles or 2 holes' g factors are equal to zero. Only the orbital part of the g-factors contributes - the spin part cancels out. We then consider the effects of introducing a single energy splitting between the 2 orbits. We make a linear approximation for all other n values.
3D prostate boundary segmentation from ultrasound images using 2D active shape models.
Hodge, Adam C; Ladak, Hanif M
2006-01-01
Boundary outlining, or segmentation, of the prostate is an important task in diagnosis and treatment planning for prostate cancer. This paper describes an algorithm for semi-automatic, three-dimensional (3D) segmentation of the prostate boundary from ultrasound images based on two-dimensional (2D) active shape models (ASM) and rotation-based slicing. Evaluation of the algorithm used distance- and volume-based error metrics to compare algorithm generated boundary outlines to gold standard (manually generated) boundary outlines. The mean absolute distance between the algorithm and gold standard boundaries was 1.09+/-0.49 mm, the average percent absolute volume difference was 3.28+/-3.16%, and a 5x speed increase as compared manual planimetry was achieved. PMID:17946106
Optical fiber poling by induction: analysis by 2D numerical modeling.
De Lucia, F; Huang, D; Corbari, C; Healy, N; Sazio, P J A
2016-04-15
Since their first demonstration some 25 years ago, thermally poled silica fibers have been used to realize device functions such as electro-optic modulation, switching, polarization-entangled photons, and optical frequency conversion with a number of advantages over bulk free-space components. We have recently developed an innovative induction poling technique that could allow for the development of complex microstructured fiber geometries for highly efficient χ^{(2)}-based device applications. To systematically implement these more advanced poled fiber designs, we report here the development of comprehensive numerical models of the induction poling mechanism itself via two-dimensional (2D) simulations of ion migration and space-charge region formation using finite element analysis. PMID:27082323
Models Ion Trajectories in 2D and 3D Electrostatic and Magnetic Fields
Dahl, David
2000-02-21
SIMION3D7.0REV is a C based ion optics simulation program that can model complex problems using Laplace equation solutions for potential fields. The program uses an ion optics workbench that can hold up to 200 2D and/or 3D electrostatic/magnetic potential arrays. Arrays can have up to 50,000,000 points. SIMION3D7.0''s 32 bit virtual Graphics User Interface provides a highly interactive advanced user environment. All potential arrays are visualized as 3D objects that the user can cut away to inspect ion trajectories and potential energy surfaces. User programs allow the user to customize the program for specific simulations. A geometry file option supports the definition of highly complex array geometry. Algorithm modifications have improved this version''s computational speed and accuracy.
Robust autonomous model learning from 2D and 3D data sets.
Langs, Georg; Donner, René; Peloschek, Philipp; Bischof, Horst
2007-01-01
In this paper we propose a weakly supervised learning algorithm for appearance models based on the minimum description length (MDL) principle. From a set of training images or volumes depicting examples of an anatomical structure, correspondences for a set of landmarks are established by group-wise registration. The approach does not require any annotation. In contrast to existing methods no assumptions about the topology of the data are made, and the topology can change throughout the data set. Instead of a continuous representation of the volumes or images, only sparse finite sets of interest points are used to represent the examples during optimization. This enables the algorithm to efficiently use distinctive points, and to handle texture variations robustly. In contrast to standard elasticity based deformation constraints the MDL criterion accounts for systematic deformations typical for training sets stemming from medical image data. Experimental results are reported for five different 2D and 3D data sets. PMID:18051152
Calibration Of 2D Hydraulic Inundation Models In The Floodplain Region Of The Lower Tagus River
NASA Astrophysics Data System (ADS)
Pestanana, R.; Matias, M.; Canelas, R.; Araujo, A.; Roque, D.; Van Zeller, E.; Trigo-Teixeira, A.; Ferreira, R.; Oliveira, R.; Heleno, S.
2013-12-01
In terms of inundated area, the largest floods in Portugal occur in the Lower Tagus River. On average, the river overflows every 2.5 years, at times blocking roads and causing important agricultural damages. This paper focus on the calibration of 2D-horizontal flood simulation models for the floods of 2001 and 2006 on a 70-km stretch of the Lower Tagus River. Flood extent maps, derived from ERS SAR and ENVISAT ASAR imagery were compared with the flood extent maps obtained for each simulation, to calibrate roughness coefficients. The combination of the calibration results from the 2001 and 2006 floods provided a preliminary Manning coefficient map of the study area.
Universality Class of the Nishimori Point in the 2D +/-J Random-Bond Ising Model
NASA Astrophysics Data System (ADS)
Honecker, A.; Picco, M.; Pujol, P.
2001-07-01
We study the universality class of the Nishimori point in the 2D +/-J random-bond Ising model by means of the numerical transfer-matrix method. Using the domain-wall free energy, we locate the position of the fixed point along the Nishimori line at the critical concentration value pc = 0.1094+/-0.0002 and estimate ν = 1.33+/-0.03. Then, we obtain the exponents for the moments of the spin-spin correlation functions as well as the value for the central charge c = 0.464+/-0.004. The main qualitative result is the fact that percolation is now excluded as a candidate for describing the universality class of this fixed point.
Universality class of the Nishimori point in the 2D +/- J random-bond Ising model.
Honecker, A; Picco, M; Pujol, P
2001-07-23
We study the universality class of the Nishimori point in the 2D +/- J random-bond Ising model by means of the numerical transfer-matrix method. Using the domain-wall free energy, we locate the position of the fixed point along the Nishimori line at the critical concentration value p(c) = 0.1094 +/- 0.0002 and estimate nu = 1.33 +/- 0.03. Then, we obtain the exponents for the moments of the spin-spin correlation functions as well as the value for the central charge c = 0.464 +/- 0.004. The main qualitative result is the fact that percolation is now excluded as a candidate for describing the universality class of this fixed point. PMID:11461639
Structure-approximating inverse protein folding problem in the 2D HP model.
Gupta, Arvind; Manuch, Ján; Stacho, Ladislav
2005-12-01
The inverse protein folding problem is that of designing an amino acid sequence which has a particular native protein fold. This problem arises in drug design where a particular structure is necessary to ensure proper protein-protein interactions. In this paper, we show that in the 2D HP model of Dill it is possible to solve this problem for a broad class of structures. These structures can be used to closely approximate any given structure. One of the most important properties of a good protein (in drug design) is its stability--the aptitude not to fold simultaneously into other structures. We show that for a number of basic structures, our sequences have a unique fold. PMID:16379538
An application of the distributed hydrologic model CASC2D to a tropical montane watershed
NASA Astrophysics Data System (ADS)
Marsik, Matt; Waylen, Peter
2006-11-01
SummaryIncreased stormflow in the Quebrada Estero watershed (2.5 km 2), in the northwestern Central Valley tectonic depression of Costa Rica, reportedly has caused flooding of the city of San Ramón in recent decades. Although scientifically untested, urban expansion was deemed the cause and remedial measures were recommended by the Programa de Investigación en Desarrollo Humano Sostenible (ProDUS). CASC2D, a physically-based, spatially explicit hydrologic model, was constructed and calibrated to a June 10th 2002 storm that delivered 110.5 mm of precipitation in 4.5 h visibly exceeded the bankfull stage (0.9 m) of the Quebrada flooding portions of San Ramón. The calibrated hydrograph showed a peak discharge 16.68% (2.5 m 3 s -1) higher, an above flood stage duration 20% shorter, and time to peak discharge 11 min later than the same observed discharge hydrograph characteristics. Simulations of changing land cover conditions from 1979 to 1999 showed an increase also in the peak discharge, above flood stage duration, and time to peak discharge. Analysis using a modified location quotient identified increased urbanization in lower portions of the watershed over the time period studied. These results suggest that increased urbanization in the Quebrada Estero watershed have increased flooding peaks, and durations above threshold, confirming the ProDUS report. These results and the CASC2D model offer an easy-to-use, pragmatic planning tool for policymakers in San Ramón to assess future development scenarios and their potential flooding impacts to San Ramón.
2D Distinct Element Method (DEM) models of the initiation, propagation and saturation of rock joints
NASA Astrophysics Data System (ADS)
Arslan, A.; Schöpfer, M. P.; Walsh, J. J.; Childs, C.
2009-12-01
In layered sequences, rock joints usually best develop within the more brittle layers and commonly display a regular spacing that scales with layer thickness. A variety of conceptual and mechanical models have been developed for these observations. A limitation of previous approaches, however, is that fracture initiation and associated interface slip are not explicitly simulated; instead, fractures were predefined and interfaces were welded. To surmount this problem, we have modelled the formation and growth of joints in layered sequences by using the two-dimensional Distinct Element Method (DEM) as implemented in the Particle Flow Code (PFC-2D). In PFC-2D, rock is represented by an assemblage of circular particles that are bonded at particle-particle contacts. Failure occurs if either the tensile or shear strength of a bond is exceeded. The models comprise a central brittle layer with high Young’s modulus, which is embedded in a low Young’s modulus matrix. The interfaces between the layers are defined by ‘smooth joint’ contacts, a modelling feature that eliminates interparticle bumpiness and associated interlocking friction. Consequently, this feature allows the user to assign macroscopic properties such as friction and cohesion along layer interfaces in a controlled manner. Layer parallel extension is applied by assigning a velocity to particles at the lateral boundaries of the model while maintaining a constant vertical confining pressure. Models were extended until joint saturation in the central layer was reached. We thereby explored the impact of confining pressure and interface properties (friction, cohesion) on joint spacing. A number of important conclusions can be drawn from our models: (i) The distributions of average horizontal normal stress within the layer and of shear stress at the interface are consistent with analytical solutions (stress-transfer theory). (ii) At low interfacial shear strength, new joints form preferentially midway between
Assessment of the Impacts of Compensation Flow Changes Upon Instream Habitat Using 2D Modelling
NASA Astrophysics Data System (ADS)
Mould, D. C.; Lane, S. N.; Christmas, M.
2004-05-01
Many millstone-grit rivers in northern England are impounded. In such cases the water company in the area has to release compensation flows from the reservoirs, traditionally to meet industrial needs: these flows are rarely set with ecology in mind; and have commonly involved constant flow. Dam overtopping may create spates, but spawning in many fish species is prompted by a spate flow in the early autumn when dams are rarely full enough to overtop. Such flows are important for fine sediment flushing and controlling the wetted useable area for spawning. Classical physical habitat modelling for instream habitat has been largely reliant upon 1D approaches, such as the Instream Flow Incremental Methodology (IFIM). Here we use a 2D finite element model (FESWMS), to simulate changes in instream habitat with variations in the compensation flow regimes. The spatial resolution of 2D models can be adapted to the scale of fish habitats so providing better representation of the reach-scale flow processes (such as slack water in the margins, wetting and drying) than the 1D case. The model is applied to the Rivers Rivelin and Loxley in Sheffield, Northern England. At the confluence of the two rivers, the compensation flow level is set at 30.6 Thousand Cubic Metres per Day (TCMD). Due to historical reasons, the compensation is not divided equally, as the Loxley receives 28 TCMD whilst the Rivelin receives only 2.6 TCMD. The model is used to simulate a transfer of 6 TCMD from the Loxley to the Rivelin. After validation, model predictions are combined with available habitat requirement data (e.g. velocity and depth needs) to develop an index of change in habitat suitability in terms of first order variables (e.g. velocity, depth and wetted useable area). This suggests that the change in compensation may significantly improve instream ecology in relation to macroinvertebrates, brown trout (Salmo trutta) and bullhead (Cottus gobio) in the Rivelin without causing detrimental impacts
Spin Circuit Model for Spin Orbit Torques in 2D Channels
NASA Astrophysics Data System (ADS)
Hong, Seokmin
2015-03-01
Recently, the unique coupling between charge and spin in topological insulators has been explored through various types of electrical measurements, which could have interesting applications. In this talk, we present a spin circuit model for spin orbit torques in topological insulator surface states and other 2D channels. We show with a simple example that results from the circuit model agree well with those obtained from nonequilibrium Green's function (NEGF) based quantum transport simulation. Some predictions of our model have already received experimental support and we hope this model can provide a unifying framework that can be used to critically evaluate experimental results, to explore new types of devices as well as to answer fundamental questions regarding these materials. The model for spin-orbit torques described here can be incorporated into a broader spin-circuit approach which, we believe, provides a natural platform for multi-physics, multi-component spintronic devices. This work was supported by FAME, a Semiconductor Research Corporation program sponsored by MARCO and DARPA.
Spin Circuit Model for 2D Channels with Spin-Orbit Coupling
Hong, Seokmin; Sayed, Shehrin; Datta, Supriyo
2016-01-01
In this paper we present a general theory for an arbitrary 2D channel with “spin momentum locking” due to spin-orbit coupling. It is based on a semiclassical model that classifies all the channel electronic states into four groups based on the sign of the z-component of the spin (up (U), down (D)) and the sign of the x-component of the velocity (+, −). This could be viewed as an extension of the standard spin diffusion model which uses two separate electrochemical potentials for U and D states. Our model uses four: U+, D+, U−, and D−. We use this formulation to develop an equivalent spin circuit that is also benchmarked against a full non-equilibrium Green’s function (NEGF) model. The circuit representation can be used to interpret experiments and estimate important quantities of interest like the charge to spin conversion ratio or the maximum spin current that can be extracted. The model should be applicable to topological insulator surface states with parallel channels as well as to other layered structures with interfacial spin-orbit coupling. PMID:26932563
Thermochemical Nonequilibrium 2D Modeling of Nitrogen Inductively Coupled Plasma Flow
NASA Astrophysics Data System (ADS)
Yu, Minghao; Yusuke, Takahashi; Hisashi, Kihara; Ken-ichi, Abe; Kazuhiko, Yamada; Takashi, Abe; Satoshi, Miyatani
2015-09-01
Two-dimensional (2D) numerical simulations of thermochemical nonequilibrium inductively coupled plasma (ICP) flows inside a 10-kW inductively coupled plasma wind tunnel (ICPWT) were carried out with nitrogen as the working gas. Compressible axisymmetric Navier-Stokes (N-S) equations coupled with magnetic vector potential equations were solved. A four-temperature model including an improved electron-vibration relaxation time was used to model the internal energy exchange between electron and heavy particles. The third-order accuracy electron transport properties (3rd AETP) were applied to the simulations. A hybrid chemical kinetic model was adopted to model the chemical nonequilibrium process. The flow characteristics such as thermal nonequilibrium, inductive discharge, effects of Lorentz force were made clear through the present study. It was clarified that the thermal nonequilibrium model played an important role in properly predicting the temperature field. The prediction accuracy can be improved by applying the 3rd AETP to the simulation for this ICPWT. supported by Grant-in-Aid for Scientific Research (No. 23560954), sponsored by the Japan Society for the Promotion of Science
Incorporating a Turbulence Transport Model into 2-D Hybrid Hall Thruster Simulations
NASA Astrophysics Data System (ADS)
Cha, Eunsun; Cappelli, Mark A.; Fernandez, Eduardo
2014-10-01
2-D hybrid simulations of Hall plasma thrusters that do not resolve cross-field transport-generating fluctuations require a model to capture how electrons migrate across the magnetic field. We describe the results of integrating a turbulent electron transport model into simulations of plasma behavior in a plane spanned by the E and B field vectors. The simulations treat the electrons as a fluid and the heavy species (ions/neutrals) as discrete particles. The transport model assumes that the turbulent eddy cascade in the electron fluid to smaller scales is the primary means of electron energy dissipation. Using this model, we compare simulations to experimental measurements made on a laboratory Hall discharge over a range of discharge voltage. Both the current-voltage trends as well as the plasma properties such as plasma temperature, electron density, and ion velocities seem agree favorably with experiments, where a simple Bohm transport model tends to perform poorly in capturing much of the discharge behavior.
Spin Circuit Model for 2D Channels with Spin-Orbit Coupling.
Hong, Seokmin; Sayed, Shehrin; Datta, Supriyo
2016-01-01
In this paper we present a general theory for an arbitrary 2D channel with "spin momentum locking" due to spin-orbit coupling. It is based on a semiclassical model that classifies all the channel electronic states into four groups based on the sign of the z-component of the spin (up (U), down (D)) and the sign of the x-component of the velocity (+, -). This could be viewed as an extension of the standard spin diffusion model which uses two separate electrochemical potentials for U and D states. Our model uses four: U+, D+, U-, and D-. We use this formulation to develop an equivalent spin circuit that is also benchmarked against a full non-equilibrium Green's function (NEGF) model. The circuit representation can be used to interpret experiments and estimate important quantities of interest like the charge to spin conversion ratio or the maximum spin current that can be extracted. The model should be applicable to topological insulator surface states with parallel channels as well as to other layered structures with interfacial spin-orbit coupling. PMID:26932563
Spin Circuit Model for 2D Channels with Spin-Orbit Coupling
NASA Astrophysics Data System (ADS)
Hong, Seokmin; Sayed, Shehrin; Datta, Supriyo
2016-03-01
In this paper we present a general theory for an arbitrary 2D channel with “spin momentum locking” due to spin-orbit coupling. It is based on a semiclassical model that classifies all the channel electronic states into four groups based on the sign of the z-component of the spin (up (U), down (D)) and the sign of the x-component of the velocity (+, -). This could be viewed as an extension of the standard spin diffusion model which uses two separate electrochemical potentials for U and D states. Our model uses four: U+, D+, U-, and D-. We use this formulation to develop an equivalent spin circuit that is also benchmarked against a full non-equilibrium Green’s function (NEGF) model. The circuit representation can be used to interpret experiments and estimate important quantities of interest like the charge to spin conversion ratio or the maximum spin current that can be extracted. The model should be applicable to topological insulator surface states with parallel channels as well as to other layered structures with interfacial spin-orbit coupling.
Field Evaluation of a Novel 2D Preferential Flow Snowpack Hydrology Model
NASA Astrophysics Data System (ADS)
Leroux, N.; Pomeroy, J. W.; Kinar, N. J.
2015-12-01
Accurate estimation of snowmelt flux is of primary importance for runoff hydrograph prediction, which is used for water management and flood forecasting. Lateral flows and preferential flow pathways in porous media flow have proven critical for improving soil and groundwater flow models, but though many physically-based layered snowmelt models have been developed, only 1D matrix flow is accounted for in these models. Therefore, there is a need for snowmelt models that include these processes so as to examine the potential to improve snowmelt hydrological modelling. A 2D model is proposed that enables an improved understanding of energy and water flows within deep heterogeneous snowpacks, including those on slopes. A dual pathway theory is presented that simulates the formation of preferential flow paths, vertical and lateral water flows through the snow matrix and flow fingers, internal energy fluxes, melt, wet snow metamorphism, and internal refreezing. The dual pathway model utilizes an explicit finite volume method to solve for the energy and water flux equations over a non-orthogonal grid. It was run and evaluated using in-situ data collected from snowpit - accessed gravimetric, thermometric, photographic, and dielectric observations and novel non-invasive acoustic observations of layering, temperature, flowpath geometry, density and wetness at the Fortress Mountain Snow Laboratory, Alberta, Canada. The melt of a natural snowpack was artificially generated after detailed observation of snowpack initial conditions such as snow layer properties, temperature, and liquid water content. Snowpack ablation and liquid water content distribution over time were then measured and used for model parameterization and validation. Energy available at the snow surface and soil slope angle were set as mondel inputs. Model verification was based on snowpack property evolution. The heterogeneous flow model can be an important tool to help understand snowmelt flow processes, how
NASA Astrophysics Data System (ADS)
Arbic, B. K.; Wetzel, A.; Cerovecki, I.; Hendershott, M.; Karsten, R. H.; Molinari, J.; Mueller, M.
2012-12-01
Recently there has been a great deal of interest in ocean mixing, primarily because ocean mixing is thought to control the strength of the oceanic meridional overturning circulation. An important source of ocean mixing is breaking internal waves, and an important energy source for oceanic internal waves is the astronomical tide. The ocean climate process team discussed in this session aims to develop parameterizations of internal-wave driven mixing which are suitable for the altered climates of our anticipated future as well as for the present climate. For this reason it is of interest to ask how climatic changes in the oceanic stratification have altered and will alter the ocean tides. By solving the Sturm-Liouville vertical mode problem with a present-day profile of oceanic temperature and salinity as well as the same profile altered by the observed temperature changes occuring over the last 135 years (Roemmich et al. 2012), we find that the square of the phase speed of the first mode baroclinic mode internal tide has changed by about 5-6% over that time frame. We alter the g' and layer depths of a two-layer numerical model of the global barotropic and baroclinic tides such that the phase speed changes by this amount. We will show the changes taking place in the baroclinic tide as well as in the barotropic tide. The barotropic tide changes because the phase speed of the barotropic tide depends on stratification (Gill 1982); this in turn affects the resonance of barotropic tides. The results of the numerical simulations are compared to a two-layer analytical model of tides in a forced-damped closed basin, which also predicts a sensitivity of both barotropic and baroclinic tides to stratification and to climatic perturbations to the stratification. Finally, the numerical results are compared to tide gauge observations of secular changes in tides occurring over the last century or so.
NASA Astrophysics Data System (ADS)
Shen, D.; Wang, J.; Cheng, X.; Rui, Y.; Ye, S.
2015-02-01
The rapid progress of Light Detection And Ranging (LiDAR) technology has made acquirement and application of high-resolution digital elevation model (DEM) data increasingly popular, especially with regards to the study of floodplain flow modeling. High-resolution DEM data include many redundant interpolation points, needs a high amount of calculation, and does not match the size of computational mesh. These disadvantages are a common problem for floodplain flow modeling studies. Two-dimensional (2-D) hydraulic modeling, a popular method of analyzing floodplain flow, offers high precision of elevation parameterization for computational mesh while ignoring much micro-topographic information of the DEM data itself. We offer a flood simulation method that integrates 2-D hydraulic model results and high-resolution DEM data, enabling the calculation of flood water levels in DEM grid cells through local inverse distance weighted interpolation. To get rid of the false inundation areas during interpolation, it employs the run-length encoding method to mark the inundated DEM grid cells and determine the real inundation areas through the run-length boundary tracing technique, which solves the complicated problem of the connectivity between DEM grid cells. We constructed a 2-D hydraulic model for the Gongshuangcha polder, a flood storage area of Dongting Lake, using our integrated method to simulate the floodplain flow. The results demonstrate that this method can solve DEM associated problems efficiently and simulate flooding processes with greater accuracy than DEM only simulations.
A new model for two-dimensional numerical simulation of pseudo-2D gas-solids fluidized beds
Li, Tingwen; Zhang, Yongmin
2013-10-11
Pseudo-two dimensional (pseudo-2D) fluidized beds, for which the thickness of the system is much smaller than the other two dimensions, is widely used to perform fundamental studies on bubble behavior, solids mixing, or clustering phenomenon in different gas-solids fluidization systems. The abundant data from such experimental systems are very useful for numerical model development and validation. However, it has been reported that two-dimensional (2D) computational fluid dynamic (CFD) simulations of pseudo-2D gas-solids fluidized beds usually predict poor quantitative agreement with the experimental data, especially for the solids velocity field. In this paper, a new model is proposed to improve the 2D numerical simulations of pseudo-2D gas-solids fluidized beds by properly accounting for the frictional effect of the front and back walls. Two previously reported pseudo-2D experimental systems were simulated with this model. Compared to the traditional 2D simulations, significant improvements in the numerical predictions have been observed and the predicted results are in better agreement with the available experimental data.
NASA Astrophysics Data System (ADS)
Kaspi, Y.; Davighi, J. E.; Galanti, E.; Hubbard, W. B.
2016-09-01
The upcoming Juno and Cassini gravity measurements of Jupiter and Saturn, respectively, will allow probing the internal dynamics of these planets through accurate analysis of their gravity spectra. To date, two general approaches have been suggested for relating the flow velocities and gravity fields. In the first, barotropic potential surface models, which naturally take into account the oblateness of the planet, are used to calculate the gravity field. However, barotropicity restricts the flows to be constant along cylinders parallel to the rotation axis. The second approach, calculated in the reference frame of the rotating planet, assumes that due to the large scale and rapid rotation of these planets, the winds are to leading order in geostrophic balance. Therefore, thermal wind balance relates the wind shear to the density gradients. While this approach can take into account any internal flow structure, it is limited to only calculating the dynamical gravity contributions, and has traditionally assumed spherical symmetry. This study comes to relate the two approaches both from a theoretical perspective, showing that they are analytically identical in the barotropic limit, and numerically, through systematically comparing the different model solutions for the gravity harmonics. For the barotropic potential surface models we employ two independent solution methods - the potential-theory and Maclaurin spheroid methods. We find that despite the sphericity assumption, in the barotropic limit the thermal wind solutions match well the barotropic oblate potential-surface solutions.
A case study of fluid flow in fractured rock mass based on 2-D DFN modeling
NASA Astrophysics Data System (ADS)
Han, Jisu; Noh, Young-Hwan; Um, Jeong-Gi; Choi, Yosoon
2014-05-01
A two dimensional steady-state fluid flow through fractured rock mass of an abandoned copper mine in Korea is addressed based on discrete fracture network modeling. An injection well and three observation wells were installed at the field site to monitor the variations of total heads induced by injection of fresh water. A series of packer tests were performed to estimate the rock mass permeability. First, the two dimensional stochastic fracture network model was built and validated for a granitic rock mass using the geometrical and statistical data obtained from surface exposures and borehole logs. This validated fracture network model was combined with the fracture data observed on boreholes to generate a stochastic-deterministic fracture network system. Estimated apertures for each of the fracture sets using permeability data obtained from borehole packer tests were discussed next. Finally, a systematic procedure for fluid flow modeling in fractured rock mass in two dimensional domain was presented to estimate the conductance, flow quantity and nodal head in 2-D conceptual linear pipe channel network. The results obtained in this study clearly show that fracture geometry parameters (orientation, density and size) play an important role in the hydraulic behavior of fractured rock masses.
Distributed and coupled 2D electro-thermal model of power semiconductor devices
NASA Astrophysics Data System (ADS)
Belkacem, Ghania; Lefebvre, Stéphane; Joubert, Pierre-Yves; Bouarroudj-Berkani, Mounira; Labrousse, Denis; Rostaing, Gilles
2014-05-01
The development of power electronics in the field of transportations (automotive, aeronautics) requires the use of power semiconductor devices providing protection and diagnostic functions. In the case of series protections power semiconductor devices which provide protection may operate in shortcircuit and act as a current limiting device. This mode of operations is very constraining due to the large dissipation of power. In these particular conditions of operation, electro-thermal models of power semiconductor devices are of key importance in order to optimize their thermal design and increase their reliability. The development of such an electro-thermal model for power MOSFET transistors based on the coupling between two computation softwares (Matlab and Cast3M) is described in this paper. The 2D electro-thermal model is able to predict (i) the temperature distribution on chip surface well as in the volume under short-circuit operations, (ii) the effect of the temperature on the distribution of the current flowing within the die and (iii) the effects of the ageing of the metallization layer on the current density and the temperature. In this paper, the electrical and thermal models are described as well as the implemented coupling scheme.
Studies of Arctic stratospheric ozone in a 2-D model including some effects of zonal asymmetries
Isaksen, I.S.A.; Rognerud, B.; Stordal, F. ); Coffey, M.T.; Mankin, W.G. )
1990-03-01
A two-dimensional (2-D) zonally averaged chemistry-transport model of the stratosphere has been extended to include some zonally asymmetric effects to study the chemically disturbed conditions in the Arctic winter during the occurrence of polar stratospheric clouds (PSCs). The model allows air parcels that have been in PSCs in the polar night to be exposed to sunlight during the passage south through a wave trough. Large enhancements of ClO are estimated as well as significant ozone reductions, most pronounced around the 20 km height level. The ozone depletions maximize in late March, about one month after the cease in PSC activity in the model, and amount to 5-8% in column ozone at 70{degree}N. In agreement with column measurements made from the DC-8, the model estimates an increase in the columns of HNO{sub 3} and ClONO{sub 2}, and a decrease in the HCl column within the polar vortex.
LBQ2D, Extending the Line Broadened Quasilinear Model to TAE-EP Interaction
NASA Astrophysics Data System (ADS)
Ghantous, Katy; Gorelenkov, Nikolai; Berk, Herbert
2012-10-01
The line broadened quasilinear model was proposed and tested on the one dimensional electrostatic case of the bump on tailfootnotetextH.L Berk, B. Breizman and J. Fitzpatrick, Nucl. Fusion, 35:1661, 1995 to study the wave particle interaction. In conventional quasilinear theory, the sea of overlapping modes evolve with time as the particle distribution function self consistently undergo diffusion in phase space. The line broadened quasilinear model is an extension to the conventional theory in a way that allows treatment of isolated modes as well as overlapping modes by broadening the resonant line in phase space. This makes it possible to treat the evolution of modes self consistently from onset to saturation in either case. We describe here the model denoted by LBQ2D which is an extension of the proposed one dimensional line broadened quasilinear model to the case of TAEs interacting with energetic particles in two dimensional phase space, energy as well as canonical angular momentum. We study the saturation of isolated modes in various regimes and present the analytical derivation and numerical results. Finally, we present, using ITER parameters, the case where multiple modes overlap and describe the techniques used for the numerical treatment.
Obtaining instantaneous water levels relative to a geoid with a 2D storm surge model
NASA Astrophysics Data System (ADS)
Slobbe, D. C.; Verlaan, M.; Klees, R.; Gerritsen, H.
2013-01-01
Current and new applications of 2D storm surge models such as the Dutch Continental Shelf Model (DCSM) require that the models provide proper estimates of the instantaneous water levels expressed relative to a particular geoid, rather than only the tide and surge components expressed relative to the ill-defined model's zero height surface. For DCSM, this is realized by adding the depth-averaged horizontal baroclinic pressure gradients to the model equations, which are derived from 4D salinity and temperature fields provided by the Proudman Oceanographic Laboratory hydrodynamic model (POL's hindcast). The vertical datum of the extended model is fixed to that of the European Gravimetric Geoid model 2008 (EGG08). This is done by an adjustment of the model parameters that depend on the choice of the reference surface (e.g., bathymetry) and by referring the water levels along the open boundaries to this reference surface. Using different numerical experiments we investigate the effects on the water levels of several approximations we have made during the implementation. The ability of the model to reproduce both the mean sea level (MSL) and instantaneous water levels is assessed by a comparison with the MSL derived from POL's hindcast as well as with instantaneous water levels acquired by various radar altimeter satellites. From this comparison we conclude that our modeled MSL is in good agreement with the MSL derived from POL's hindcast; the standard deviation of the differences is below 2 cm. However, larger differences in MSL are observed when comparing the model output with the MSL derived from radar altimeter data. They are attributed to either geoid errors or errors in the used salinity and temperature fields. The root mean squared (rms) differences between observed and modeled instantaneous water levels over the entire model domain varies from 9 cm for data acquired by the TOPEX satellite to 11 cm for data acquired by the GFO-1 satellite. These numbers improve to
A New 2D-Advection-Diffusion Model Simulating Trace Gas Distributions in the Lowermost Stratosphere
NASA Astrophysics Data System (ADS)
Hegglin, M. I.; Brunner, D.; Peter, T.; Wirth, V.; Fischer, H.; Hoor, P.
2004-12-01
Tracer distributions in the lowermost stratosphere are affected by both, transport (advective and non-advective) and in situ sources and sinks. They influence ozone photochemistry, radiative forcing, and heating budgets. In-situ measurements of long-lived species during eight measurement campaigns revealed relatively simple behavior of the tracers in the lowermost stratosphere when represented in an equivalent-latitude versus potential temperature framework. We here present a new 2D-advection-diffusion model that simulates the main transport pathways influencing the tracer distributions in the lowermost stratosphere. The model includes slow diabatic descent of aged stratospheric air and vertical and/or horizontal diffusion across the tropopause and within the lowermost stratosphere. The diffusion coefficients used in the model represent the combined effects of different processes with the potential of mixing tropospheric air into the lowermost stratosphere such as breaking Rossby and gravity waves, deep convection penetrating the tropopause, turbulent diffusion, radiatively driven upwelling etc. They were specified by matching model simulations to observed distributions of long-lived trace gases such as CO and N2O obtained during the project SPURT. The seasonally conducted campaigns allow us to study the seasonal dependency of the diffusion coefficients. Despite its simplicity the model yields a surprisingly good description of the small scale features of the measurements and in particular of the observed tracer gradients at the tropopause. The correlation coefficients between modeled and measured trace gas distributions were up to 0.95. Moreover, mixing across isentropes appears to be more important than mixing across surfaces of constant equivalent latitude (or PV). With the aid of the model, the distribution of the fraction of tropospheric air in the lowermost stratosphere can be determined.
Transectional heat transfer in thermoregulating bigeye tuna (Thunnus obesus) - a 2D heat flux model.
Boye, Jess; Musyl, Michael; Brill, Richard; Malte, Hans
2009-11-01
We developed a 2D heat flux model to elucidate routes and rates of heat transfer within bigeye tuna Thunnus obesus Lowe 1839 in both steady-state and time-dependent settings. In modeling the former situation, we adjusted the efficiencies of heat conservation in the red and the white muscle so as to make the output of the model agree as closely as possible with observed cross-sectional isotherms. In modeling the latter situation, we applied the heat exchanger efficiencies from the steady-state model to predict the distribution of temperature and heat fluxes in bigeye tuna during their extensive daily vertical excursions. The simulations yielded a close match to the data recorded in free-swimming fish and strongly point to the importance of the heat-producing and heat-conserving properties of the white muscle. The best correspondence between model output and observed data was obtained when the countercurrent heat exchangers in the blood flow pathways to the red and white muscle retained 99% and 96% (respectively) of the heat produced in these tissues. Our model confirms that the ability of bigeye tuna to maintain elevated muscle temperatures during their extensive daily vertical movements depends on their ability to rapidly modulate heating and cooling rates. This study shows that the differential cooling and heating rates could be fully accounted for by a mechanism where blood flow to the swimming muscles is either exclusively through the heat exchangers or completely shunted around them, depending on the ambient temperature relative to the body temperature. Our results therefore strongly suggest that such a mechanism is involved in the extensive physiological thermoregulatory abilities of endothermic bigeye tuna. PMID:19880733
Lobe dynamics applied to barotropic Rossby-wave breaking
NASA Astrophysics Data System (ADS)
Koh, Tieh-Yong; Plumb, R. Alan
2000-06-01
We applied the methods of lobe dynamics to the problem of transport across the edge of a barotropic vortex patch. The model used captures the essential dynamics of filament shedding in the wintertime stratospheric polar vortex. Two approaches were adopted for the problem: (1) the dominant periodic component of the vortical flow was identified and conventional lobe dynamics methods for periodic dynamical systems were applied to it; (2) the full aperiodic, dynamically consistent flow was retained and a modified brand of lobe dynamics was used to quantify the transport. Our results show that in the periodic case, much reversible transport occurs across the lobe dynamical boundary due to overlapping intruding and extruding lobes. In the aperiodic case, a small amount of intrusion was noted, contrary to the well-established fact that potential vorticity shedding in barotropic vortices is uniquely outwards. In our discussion, we argue that while lobe dynamics provides a rigorous framework for quantifying transport across the lobe dynamical boundary, this boundary may not be appropriate for quantifying transport across internal transport barriers, such as the stratospheric polar vortex edge.
Comparing a 2D fluid model of the DC planar magnetron cathode to experiments
Garcia, M.
1996-05-01
Planar magnetron cathodes have arching magnetic field lines which concentrate plasma density near the electrode surface. This enhances the ion bombardment of the surface and the yield of sputtered atoms. Magnetron cathodes are used in the Plasma Electrode Pockels Cell (PEPC) devices of the Laser Program because they provide for significantly higher conduction than do glow discharges. An essential feature of magnetron cathodes is that the vector product of the perpendicular electric field, E[sub y], with the parallel component of the magnetic field, B[sub x], forms a closed track with a circulating current along the cathode surface. An analytical, 2D, two component, quasi-neutral, continuum model yields formulas for the plasma density, the total and component current densities, the electric field, and the positive electrical potential, between the cathode surface and a distant, uniform plasma. For a specific gas, the free parameters are electron temperature, gas number density, and total current. The model is applied to the interpretation of experimental data from the PEPC device, as well as a small vacuum facility for testing magnetron cathodes. Finally, the model has been applied to generate cross sectional views of a PEPC magnetron cathode track.
Doubled CO2 Effects on NO(y) in a Coupled 2D Model
NASA Technical Reports Server (NTRS)
Rosenfield, J. E.; Douglass, A. R.
1998-01-01
Changes in temperature and ozone have been the main focus of studies of the stratospheric impact of doubled CO2. Increased CO2 is expected to cool the stratosphere, which will result in increases in stratospheric ozone through temperature dependent loss rates. Less attention has been paid to changes in minor constituents which affect the O3 balance and which may provide additional feedbacks. Stratospheric NO(y) fields calculated using the GSFC 2D interactive chemistry-radiation-dynamics model show significant sensitivity to the model CO2. Modeled upper stratospheric NO(y) decreases by about 15% in response to CO2 doubling, mainly due to the temperature decrease calculated to result from increased cooling. The abundance of atomic nitrogen, N, increases because the rate of the strongly temperature dependent reaction N + O2 yields NO + O decreases at lower temperatures. Increased N leads to an increase in the loss of NO(y) which is controlled by the reaction N + NO yields N2 + O. The NO(y) reduction is shown to be sensitive to the NO photolysis rate. The decrease in the O3 loss rate due to the NO(y) changes is significant when compared to the decrease in the O3 loss rate due to the temperature changes.
Simulation of abrasive flow machining process for 2D and 3D mixture models
NASA Astrophysics Data System (ADS)
Dash, Rupalika; Maity, Kalipada
2015-12-01
Improvement of surface finish and material removal has been quite a challenge in a finishing operation such as abrasive flow machining (AFM). Factors that affect the surface finish and material removal are media viscosity, extrusion pressure, piston velocity, and particle size in abrasive flow machining process. Performing experiments for all the parameters and accurately obtaining an optimized parameter in a short time are difficult to accomplish because the operation requires a precise finish. Computational fluid dynamics (CFD) simulation was employed to accurately determine optimum parameters. In the current work, a 2D model was designed, and the flow analysis, force calculation, and material removal prediction were performed and compared with the available experimental data. Another 3D model for a swaging die finishing using AFM was simulated at different viscosities of the media to study the effects on the controlling parameters. A CFD simulation was performed by using commercially available ANSYS FLUENT. Two phases were considered for the flow analysis, and multiphase mixture model was taken into account. The fluid was considered to be a
Complex 2D matrix model and geometrical map on the complex-Nc plane
NASA Astrophysics Data System (ADS)
Nawa, Kanabu; Ozaki, Sho; Nagahiro, Hideko; Jido, Daisuke; Hosaka, Atsushi
2013-08-01
We study the parameter dependence of the internal structure of resonance states by formulating a complex two-dimensional (2D) matrix model, where the two dimensions represent two levels of resonances. We calculate a critical value of the parameter at which a "nature transition" with character exchange occurs between two resonance states, from the viewpoint of geometry on complex-parameter space. Such a critical value is useful for identifying the internal structure of resonance states with variation of the parameter in the system. We apply the model to analyze the internal structure of hadrons with variation of the color number N_c from infty to a realistic value 3. By regarding 1/N_c as the variable parameter in our model, we calculate a critical color number of the nature transition between hadronic states in terms of a quark-antiquark pair and a mesonic molecule as exotics from the geometry on the complex-N_c plane. For large-N_c effective theory, we employ the chiral Lagrangian induced by holographic QCD with a D4/D8/overline {D8} multi-D brane system in type IIA superstring theory.
Field-induced magnetization jumps and quantum criticality in the 2D J-Q model
NASA Astrophysics Data System (ADS)
Iaizzi, Adam; Sandvik, Anders
The J-Q model is a `designer hamiltonian' formed by adding a four spin `Q' term to the standard antiferromagnetic S = 1 / 2 Heisenberg model. The Q term drives a quantum phase transition to a valence-bond solid (VBS) state: a non-magnetic state with a pattern of local singlets which breaks lattice symmetries. The elementary excitations of the VBS are triplons, i.e. gapped S=1 quasiparticles. There is considerable interest in the quantum phase transition between the Néel and VBS states as an example of deconfined quantum criticality. Near the phase boundary, triplons deconfine into pairs of bosonic spin-1/2 excitations known as spinons. Using exact diagonalization and the stochastic series expansion quantum monte carlo method, we study the 2D J-Q model in the presence of an external magnetic field. We use the field to force a nonzero density of magnetic excitations at T=0 and look for signatures of Bose-Einstein condensation of spinons. At higher magnetic fields, there is a jump in the induced magnetization caused by the onset of an effective attractive interaction between magnons on a ferromagnetic background. We characterize the first order quantum phase transition and determine the minimum value of the coupling ratio q ≡ Q / J required to produce this jump. Funded by NSF DMR-1410126.
2d Affine XY-Spin Model/4d Gauge Theory Duality and Deconfinement
Anber, Mohamed M.; Poppitz, Erich; Unsal, Mithat; /SLAC /Stanford U., Phys. Dept. /San Francisco State U.
2012-08-16
We introduce a duality between two-dimensional XY-spin models with symmetry-breaking perturbations and certain four-dimensional SU(2) and SU(2) = Z{sub 2} gauge theories, compactified on a small spatial circle R{sup 1,2} x S{sup 1}, and considered at temperatures near the deconfinement transition. In a Euclidean set up, the theory is defined on R{sup 2} x T{sup 2}. Similarly, thermal gauge theories of higher rank are dual to new families of 'affine' XY-spin models with perturbations. For rank two, these are related to models used to describe the melting of a 2d crystal with a triangular lattice. The connection is made through a multi-component electric-magnetic Coulomb gas representation for both systems. Perturbations in the spin system map to topological defects in the gauge theory, such as monopole-instantons or magnetic bions, and the vortices in the spin system map to the electrically charged W-bosons in field theory (or vice versa, depending on the duality frame). The duality permits one to use the two-dimensional technology of spin systems to study the thermal deconfinement and discrete chiral transitions in four-dimensional SU(N{sub c}) gauge theories with n{sub f} {ge} 1 adjoint Weyl fermions.
Prominence fine-structure dynamics as inferred from 2D non-LTE models
NASA Astrophysics Data System (ADS)
Gunar, Stanislav; Schmieder, Brigitte; Mein, Pierre; Heinzel, Petr
2012-07-01
2D multi-thread prominence fine structure models are able to produce synthetic Lyman spectra in very good agreement with spectral observations by SOHO/SUMER including the spectral line asymmetries. The synthetic differential emission measure curves derived from these models are also in a good agreement with observations. Now we show that these models are also able to produce synthetic H-alpha line profiles in very good agreement with observations which allows us to analyze not only the physical parameters of the prominence fine-structure plasma but also some aspects of its dynamical behaviour. We compare the synthetic H-alpha spectra with the observed spectra of the April 26, 2007 prominence using three statistical parameters: the line integrated intensity, the line full-width at the half-maximum (FWHM), and the Doppler velocity derived from shifts of the line profiles. This statistical analysis allows us to conclude that the overall statistical distribution of the LOS velocities in the April 26, 2007 prominence at the time of the observations was below +/-15 km/s and in the prominence core was close to +/-10 km/s. In combination with the analysis of the Lyman spectra we determine several physical parameters of the observed prominence fine structures which show that the April 26, 2007 prominence was relatively less massive. We are also able to put some constrains on the prominence core temperature that might be relatively low, reaching values below 6000 K.
Automatic 3D high-fidelity traffic interchange modeling using 2D road GIS data
NASA Astrophysics Data System (ADS)
Wang, Jie; Shen, Yuzhong
2011-03-01
3D road models are widely used in many computer applications such as racing games and driving simulations. However, almost all high-fidelity 3D road models were generated manually by professional artists at the expense of intensive labor. There are very few existing methods for automatically generating 3D high-fidelity road networks, especially for those existing in the real world. Real road network contains various elements such as road segments, road intersections and traffic interchanges. Among them, traffic interchanges present the most challenges to model due to their complexity and the lack of height information (vertical position) of traffic interchanges in existing road GIS data. This paper proposes a novel approach that can automatically produce 3D high-fidelity road network models, including traffic interchange models, from real 2D road GIS data that mainly contain road centerline information. The proposed method consists of several steps. The raw road GIS data are first preprocessed to extract road network topology, merge redundant links, and classify road types. Then overlapped points in the interchanges are detected and their elevations are determined based on a set of level estimation rules. Parametric representations of the road centerlines are then generated through link segmentation and fitting, and they have the advantages of arbitrary levels of detail with reduced memory usage. Finally a set of civil engineering rules for road design (e.g., cross slope, superelevation) are selected and used to generate realistic road surfaces. In addition to traffic interchange modeling, the proposed method also applies to other more general road elements. Preliminary results show that the proposed method is highly effective and useful in many applications.
Govind Rajan, Ananth; Warner, Jamie H; Blankschtein, Daniel; Strano, Michael S
2016-04-26
Transition metal dichalcogenides (TMDs) like molybdenum disulfide (MoS2) and tungsten disulfide (WS2) are layered materials capable of growth to one monolayer thickness via chemical vapor deposition (CVD). Such CVD methods, while powerful, are notoriously difficult to extend across different reactor types and conditions, with subtle variations often confounding reproducibility, particularly for 2D TMD growth. In this work, we formulate the first generalized TMD synthetic theory by constructing a thermodynamic and kinetic growth mechanism linked to CVD reactor parameters that is predictive of specific geometric shape, size, and aspect ratio from triangular to hexagonal growth, depending on specific CVD reactor conditions. We validate our model using experimental data from Wang et al. (Chem. Mater. 2014, 26, 6371-6379) that demonstrate the systemic evolution of MoS2 morphology down the length of a flow CVD reactor where variations in gas phase concentrations can be accurately estimated using a transport model (CSulfur = 9-965 μmol/m(3); CMoO3 = 15-16 mmol/m(3)) under otherwise isothermal conditions (700 °C). A stochastic model which utilizes a site-dependent activation energy barrier based on the intrinsic TMD bond energies and a series of Evans-Polanyi relations leads to remarkable, quantitative agreement with both shape and size evolution along the reactor. The model is shown to extend to the growth of WS2 at 800 °C and MoS2 under varied process conditions. Finally, a simplified theory is developed to translate the model into a "kinetic phase diagram" of the growth process. The predictive capability of this model and its extension to other TMD systems promise to significantly increase the controlled synthesis of such materials. PMID:26937889
2D dynamical magma propagation modeling: application to the 2001 Mount Etna eruption
NASA Astrophysics Data System (ADS)
Pinel, Virginie; Carrara, Alexandre; Maccaferri, Francesco; Rivalta, Eleonora; Corbi, Fabio
2016-04-01
Numerical and analog studies of dike propagation in a stress field induced by volcanic edifice construction have shown that surface loading tends both to attract the magma and to reduce its velocity. Available numerical models can either calculate the trajectory or the velocity of the ascending dikes, but not both of them simultaneously. We developed a hybrid model of dyke propagation in two dimensions solving both for the magma trajectory and velocity as a function of the source overpressure, the magma physical properties (density and viscosity) as well as the crustal density and stress field. We first calculate a dyke trajectory in 2D and secondly run a 1D dynamical model of dyke propagation along this trajectory taken into account the influence of the stress field seen by the magma along this path. This model is used to characterize the influence of surface load on magma migration towards the surface and compared to previous results obtained by analog modeling.We find that the amplitude of dyke deflection and magma velocity variation depend on the ratio between the dyke driving pressure (source overpressure as well buoyancy) and the stress field perturbation. Our model is then applied to the July 2001 eruption of Etna, where the final dyke deflection had been previously interpreted as due to the topographic load by Bonaccorso et al. [2010]. We show that the velocity decrease observed during the last stage of the propagation can also be attributed to the local stress field. We use the dyke propagation duration to estimate the magma overpressure at the dyke bottom to be less than 4 MPa.
Wu, Chao; Jiang, Xi-Ling; Shen, Hong-Wu; Yu, Ai-Ming
2009-01-01
Harmaline is a β-carboline alkaloid showing neuroprotective and neurotoxic properties. Our recent studies have revealed an important role for cytochrome P450 2D6 (CYP2D6) in harmaline O-demethylation. This study, therefore, aimed to delineate the effects of CYP2D6 phenotype/genotype on harmaline metabolism, pharmacokinetics (PK) and pharmacodynamics (PD), and to develop a pharmacogenetics mechanism-based compartmental PK model. In vitro kinetic studies on metabolite formation in human CYP2D6 extensive metabolizer (EM) and poor metabolizer (PM) hepatocytes indicated that harmaline O-demethylase activity (Vmax/Km) was about 9-fold higher in EM hepatocytes. Substrate depletion showed mono-exponential decay trait, and estimated in vitro harmaline clearance (CLint, μL/min/106 cells) was significantly lower in PM hepatocytes (28.5) than EM hepatocytes (71.1). In vivo studies in CYP2D6-humanized and wild-type mouse models showed that wild-type mice were subjected to higher and longer exposure to harmaline (5 and 15 mg/kg; i.v. and i.p.), and more severe hypothermic responses. The PK/PD data were nicely described by our pharmacogenetics-based PK model involving the clearance of drug by CYP2D6 (CLCYP2D6) and other mechanisms (CLother), and an indirect response PD model, respectively. Wild-type mice were also more sensitive to harmaline in marble-burying tests, as manifested by significantly lower ED50 and steeper Hill slope. These findings suggest that distinct CYP2D6 status may cause considerable variations in harmaline metabolism, PK and PD. In addition, the pharmacogenetics-based PK model may be extended to define PK difference caused by other polymorphic drug-metabolizing enzyme in different populations. PMID:19445902
On the assimilation of SWOT type data into 2D shallow-water models
NASA Astrophysics Data System (ADS)
Frédéric, Couderc; Denis, Dartus; Pierre-André, Garambois; Ronan, Madec; Jérôme, Monnier; Jean-Paul, Villa
2013-04-01
In river hydraulics, assimilation of water level measurements at gauging stations is well controlled, while assimilation of images is still delicate. In the present talk, we address the richness of satellite mapped information to constrain a 2D shallow-water model, but also related difficulties. 2D shallow models may be necessary for small scale modelling in particular for low-water and flood plain flows. Since in both cases, the dynamics of the wet-dry front is essential, one has to elaborate robust and accurate solvers. In this contribution we introduce robust second order, stable finite volume scheme [CoMaMoViDaLa]. Comparisons of real like tests cases with more classical solvers highlight the importance of an accurate flood plain modelling. A preliminary inverse study is presented in a flood plain flow case, [LaMo] [HoLaMoPu]. As a first step, a 0th order data processing model improves observation operator and produces more reliable water level derived from rough measurements [PuRa]. Then, both model and flow behaviours can be better understood thanks to variational sensitivities based on a gradient computation and adjoint equations. It can reveal several difficulties that a model designer has to tackle. Next, a 4D-Var data assimilation algorithm used with spatialized data leads to improved model calibration and potentially leads to identify river discharges. All the algorithms are implemented into DassFlow software (Fortran, MPI, adjoint) [Da]. All these results and experiments (accurate wet-dry front dynamics, sensitivities analysis, identification of discharges and calibration of model) are currently performed in view to use data from the future SWOT mission. [CoMaMoViDaLa] F. Couderc, R. Madec, J. Monnier, J.-P. Vila, D. Dartus, K. Larnier. "Sensitivity analysis and variational data assimilation for geophysical shallow water flows". Submitted. [Da] DassFlow - Data Assimilation for Free Surface Flows. Computational software http
MAST-2D diffusive model for flood prediction on domains with triangular Delaunay unstructured meshes
NASA Astrophysics Data System (ADS)
Aricò, C.; Sinagra, M.; Begnudelli, L.; Tucciarelli, T.
2011-11-01
A new methodology for the solution of the 2D diffusive shallow water equations over Delaunay unstructured triangular meshes is presented. Before developing the new algorithm, the following question is addressed: it is worth developing and using a simplified shallow water model, when well established algorithms for the solution of the complete one do exist? The governing Partial Differential Equations are discretized using a procedure similar to the linear conforming Finite Element Galerkin scheme, with a different flux formulation and a special flux treatment that requires Delaunay triangulation but entire solution monotonicity. A simple mesh adjustment is suggested, that attains the Delaunay condition for all the triangle sides without changing the original nodes location and also maintains the internal boundaries. The original governing system is solved applying a fractional time step procedure, that solves consecutively a convective prediction system and a diffusive correction system. The non linear components of the problem are concentrated in the prediction step, while the correction step leads to the solution of a linear system of the order of the number of computational cells. A semi-analytical procedure is applied for the solution of the prediction step. The discretized formulation of the governing equations allows to handle also wetting and drying processes without any additional specific treatment. Local energy dissipations, mainly the effect of vertical walls and hydraulic jumps, can be easily included in the model. Several numerical experiments have been carried out in order to test (1) the stability of the proposed model with regard to the size of the Courant number and to the mesh irregularity, (2) its computational performance, (3) the convergence order by means of mesh refinement. The model results are also compared with the results obtained by a fully dynamic model. Finally, the application to a real field case with a Venturi channel is presented.
Basic Brackets of a 2D Model for the Hodge Theory Without its Canonical Conjugate Momenta
NASA Astrophysics Data System (ADS)
Kumar, R.; Gupta, S.; Malik, R. P.
2016-06-01
We deduce the canonical brackets for a two (1+1)-dimensional (2D) free Abelian 1-form gauge theory by exploiting the beauty and strength of the continuous symmetries of a Becchi-Rouet-Stora-Tyutin (BRST) invariant Lagrangian density that respects, in totality, six continuous symmetries. These symmetries entail upon this model to become a field theoretic example of Hodge theory. Taken together, these symmetries enforce the existence of exactly the same canonical brackets amongst the creation and annihilation operators that are found to exist within the standard canonical quantization scheme. These creation and annihilation operators appear in the normal mode expansion of the basic fields of this theory. In other words, we provide an alternative to the canonical method of quantization for our present model of Hodge theory where the continuous internal symmetries play a decisive role. We conjecture that our method of quantization is valid for a class of field theories that are tractable physical examples for the Hodge theory. This statement is true in any arbitrary dimension of spacetime.
Learning structured models for segmentation of 2-D and 3-D imagery.
Lucchi, Aurelien; Marquez-Neila, Pablo; Becker, Carlos; Li, Yunpeng; Smith, Kevin; Knott, Graham; Fua, Pascal
2015-05-01
Efficient and accurate segmentation of cellular structures in microscopic data is an essential task in medical imaging. Many state-of-the-art approaches to image segmentation use structured models whose parameters must be carefully chosen for optimal performance. A popular choice is to learn them using a large-margin framework and more specifically structured support vector machines (SSVM). Although SSVMs are appealing, they suffer from certain limitations. First, they are restricted in practice to linear kernels because the more powerful nonlinear kernels cause the learning to become prohibitively expensive. Second, they require iteratively finding the most violated constraints, which is often intractable for the loopy graphical models used in image segmentation. This requires approximation that can lead to reduced quality of learning. In this paper, we propose three novel techniques to overcome these limitations. We first introduce a method to "kernelize" the features so that a linear SSVM framework can leverage the power of nonlinear kernels without incurring much additional computational cost. Moreover, we employ a working set of constraints to increase the reliability of approximate subgradient methods and introduce a new way to select a suitable step size at each iteration. We demonstrate the strength of our approach on both 2-D and 3-D electron microscopic (EM) image data and show consistent performance improvement over state-of-the-art approaches. PMID:25438309
T-duality without isometry via extended gauge symmetries of 2D sigma models
NASA Astrophysics Data System (ADS)
Chatzistavrakidis, Athanasios; Deser, Andreas; Jonke, Larisa
2016-01-01
Target space duality is one of the most profound properties of string theory. However it customarily requires that the background fields satisfy certain invariance conditions in order to perform it consistently; for instance the vector fields along the directions that T-duality is performed have to generate isometries. In the present paper we examine in detail the possibility to perform T-duality along non-isometric directions. In particular, based on a recent work of Kotov and Strobl, we study gauged 2D sigma models where gauge invariance for an extended set of gauge transformations imposes weaker constraints than in the standard case, notably the corresponding vector fields are not Killing. This formulation enables us to follow a procedure analogous to the derivation of the Buscher rules and obtain two dual models, by integrating out once the Lagrange multipliers and once the gauge fields. We show that this construction indeed works in non-trivial cases by examining an explicit class of examples based on step 2 nilmanifolds.
A computational model of the short-cut rule for 2D shape decomposition.
Luo, Lei; Shen, Chunhua; Liu, Xinwang; Zhang, Chunyuan
2015-01-01
We propose a new 2D shape decomposition method based on the short-cut rule. The short-cut rule originates from cognition research, and states that the human visual system prefers to partition an object into parts using the shortest possible cuts. We propose and implement a computational model for the short-cut rule and apply it to the problem of shape decomposition. The model we proposed generates a set of cut hypotheses passing through the points on the silhouette, which represent the negative minima of curvature. We then show that most part-cut hypotheses can be eliminated by analysis of local properties of each. Finally, the remaining hypotheses are evaluated in ascending length order, which guarantees that of any pair of conflicting cuts only the shortest will be accepted. We demonstrate that, compared with state-of-the-art shape decomposition methods, the proposed approach achieves decomposition results, which better correspond to human intuition as revealed in psychological experiments. PMID:25438318
Optimal implicit 2-D finite differences to model wave propagation in poroelastic media
NASA Astrophysics Data System (ADS)
Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.
2016-08-01
Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite differences (FD) to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10 kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P and slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High-order explicit FD can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here, we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.