Combined convective and diffusive modeling of the ring current and radiation belt electron dynamics using the VERB-4D code

Science.gov (United States)

Aseev, N.; Shprits, Y.; Drozdov, A.; Kellerman, A. C.; Wang, D.

2017-12-01

Ring current and radiation belts are key elements in the global dynamics of the Earth's magnetosphere. Comprehensive mathematical models are useful tools that allow us to understand the multiscale dynamics of these charged particle populations. In this work, we present results of simulations of combined ring current - radiation belt electron dynamics using the four-dimensional Versatile Electron Radiation Belt (VERB-4D) code. The VERB-4D code solves the modified Fokker-Planck equation including convective terms and models simultaneously ring current (1 - 100 keV) and radiation belt (100 keV - several MeV) electron dynamics. We apply the code to the number of geomagnetic storms that occurred in the past, compare the results with different satellite observations, and show how low-energy particles can affect the high-energy populations. Particularly, we use data from Polar Operational Environmental Satellite (POES) mission that provides a very good MLT coverage with 1.5-hour time resolution. The POES data allow us to validate the approach of the VERB-4D code for modeling MLT-dependent processes such as electron drift, wave-particle interactions, and magnetopause shadowing. We also show how different simulation parameters and empirical models can affect the results, making a particular emphasis on the electric and magnetic field models. This work will help us reveal advantages and disadvantages of the approach behind the code and determine its prediction efficiency.

3D instantaneous dynamics modeling of present-day Aegean subduction

Science.gov (United States)

Glerum, Anne; Spakman, Wim; van Hinsbergen, Douwe; Pranger, Casper

2017-04-01

To study the sensitivity of surface observables to subduction and mantle flow, i.e. the coupling of crustal tectonics and the underlying mantle dynamics, we have developed 3D numerical models of the instantaneous crust-mantle dynamics of the eastern Mediterranean. These models comprise both a realistic crust-lithosphere system and the underlying mantle. The focus for this presentation lies on the regional crustal flow response to the present-day Aegean subduction system. Our curved model domain measures 40°x40°x2900km with the Aegean subduction system taken as the geographic center. Model set-ups are based on geological and geophysical data of the eastern Mediterranean. We first create a 3D synthetic geometry of the crust-lithosphere system in a stand-alone program, including the present-day configuration of the plates in the region and crust and lithosphere thickness variations abstracted from Moho and LAB maps (Faccenna et al., 2014, Carafa et al., 2015). In addition we construct the geometry of the Aegean slab from a seismic tomography model (UU-P07; Amaru, 2007) and earthquake hypocenters (NCEDC, 2014). Geometries are then imported into the finite element code ASPECT (Kronbichler et al., 2012) using specially designed plugins. The mantle initial temperature conditions can include deviations from an adiabatic profile obtained from conversion of the UU-P07 seismic velocity anomalies to temperature anomalies using a depth-dependent scaling (Karato, 2008). We model compressible mantle flow for which material properties are obtained from thermodynamics P-T lookup-tables (Perple_X, Connolly, 2009) in combination with nonlinear viscoplastic rheology laws. Sublithospheric flow through the lateral model boundaries is left free via open boundary conditions (Chertova et al., 2012), while plate motion is prescribed at the model sides in terms of relative as well as absolute plate motion velocities (e.g. Doubrovine et al., 2012). So far, we used a free-slip surface, but

The northern tidal dynamic of Aceh waters: A 3D numerical model

Science.gov (United States)

Irham, M.; Miswar, E.; Ilhamsyah, Y.; Setiawan, I.

2018-05-01

The northern tidal dynamic of Aceh waters studied by employing three-dimensional (3D) numerical hydrodynamic model. The purpose of this study is to understand the phenomena and the characteristic of the northern tidal dynamic of Aceh waters. The research used the explicit-splitting scheme numerical model of Navier-Stokes formulation. The result displays that the vertical rotation of flow movement (vertical eddy) at a depth of 15 to 25 meter eastern part of the study area. Hence, the result also informs that the current circulation identically to the upwelling in the western region of Aceh during the wet season and vice versa. However, during the transitional season, the flow circulation depends on how the tidal dynamic occurs in the area.

Dynamical localization in the 3-D kicked Rydberg atom

International Nuclear Information System (INIS)

Persson, E.; Yoshida, S.; Tong, X.-M.; Reinhold, C.; Burgdoerfer, J.

2001-01-01

Full text: The dynamical localization for the 3D periodically kicked Rydberg atom is analyzed. For the 1D kicked atom, earlier work shows dynamical localization as the quantum suppression of classically fast ionization associated with unbounded chaotic trajectories. The corresponding wave functions localize around unstable periodic orbits. For the experimental observation, the crucial question is the dependence of the dynamical localization on the dimension. As the first step, we simulate the full 3D evolution of an extreme parabolic initial state elongated in the direction of the unidirectional kicks. We compare this simulation with the 1D model and find signatures of localization also in 3D. We further examine the dependence of quantum localization on the parabolic quantum number of the initial state. In the limit of high kick frequencies, the origin of the localization can be understood in terms of Stark states in the average field. We discuss conditions for where an experimental observation of the localization is most likely. (author)

Nested 1D-2D approach for urban surface flood modeling

Science.gov (United States)

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

Simulation of Missing Pellet Surface thermal behavior with 3D dynamic gap element

International Nuclear Information System (INIS)

Kim, Hyo Chan; Yang, Yong Sik; Koo, Yang Hyun; Kang, Chang Hak; Lee Sung Uk; Yang, Dong Yol

2014-01-01

Most of the fuel performance codes that are able to simulate a multidimensional analysis are used to calculate the radial temperature distribution and perform a multidimensional mechanical analysis based on a one-dimensional (1D) temperature result. The FRAPCON-FRAPTRAN code system incorporates a 1D thermal module and two-dimensional (2D) mechanical module when FEM option is activated. In this method, the multidimensional gap conductance model is not required because one-dimensional thermal analysis is carried out. On the other hand, a gap conductance model for a multi-dimension should be developed in the code to perform a multidimensional thermal analysis. ALCYONE developed by CEA introduces an equivalent heat convection coefficient that represents the multidimensional gap conductance. However, the code does not employ dynamic gap conductance which is a function of gap thickness and gap characteristics in direct. The BISON code, which has been developed by INL (Idaho National Laboratory), employed a thermo-mechanical contact method that is specifically designed for tightly-coupled implicit solutions that employ Jacobian-free solution methods. Owing to tightly-coupled implicit solutions, the BISON code solves gap conductance and gap thickness simultaneously with given boundary conditions. In this paper, 3D dynamic gap element has been proposed to resolve convergence issue and nonlinear characteristic of multidimensional gap conductance. To evaluate 3D dynamic gap element module, 3D thermomechanical module using FORTRAN77 has been implemented incorporating 3D dynamic gap element. To demonstrate effect of 3D dynamic gap element, thermal behavior of missing pellet surface (MPS) has been simulated by the developed module. LWR fuel performance codes should incorporate thermo-mechanical loop to solve gap conductance problem, iteratively. However, gap conductance in multidimensional model is difficult issue owing to its nonlinearity and convergence characteristics. In

Axial turbomachine modelling with a 1D axisymmetric approach

International Nuclear Information System (INIS)

Tauveron, Nicolas; Saez, Manuel; Ferrand, Pascal; Leboeuf, Francis

2007-01-01

This work concerns the design and safety analysis of direct cycle gas cooled reactor. The estimation of compressor and turbine performances in transient operations is of high importance for the designer. The first goal of this study is to provide a description of compressor behaviour in unstable conditions with a better understanding than the models based on performance maps ('traditional' 0D approach). A supplementary objective is to provide a coherent description of the turbine behaviour. The turbomachine modelling approach consists in the solution of 1D axisymmetric Navier-Stokes equations on an axial grid inside the turbomachine: mass, axial momentum, circumferential momentum and total-enthalpy balances are written. Blade forces are taken into account by using compressor or turbine blade cascade steady correlations. A particular effort has been developed to generate or test correlations in low mass flow and negative mass flow regimes, based on experimental data. The model is tested on open literature cases of the gas turbine aircraft community. For compressor and turbine, steady situations are fairly described, especially for medium and high mass flow rate. The dynamic behaviour of compressor is also quite well described, even in unstable operation (surge): qualitative tendencies (role of plenum volume and role of throttle) and some quantitative characteristics (frequency) are in a good agreement with experimental data. The application to transient simulations of gas cooled nuclear reactors is concentrated on the hypothetical 10 in. break accident. The results point out the importance of the location of the pipe rupture in a hypothetical break event. In some detailed cases, compressor surge and back flow through the circuit can occur. In order to be used in a design phase, a simplified model of surge has also been developed. This simplified model is applied to the gas fast reactor (GFR) and compared quite favourably with 1D axisymmetric simulation results

1D + 3D two-phase flow numerical model of a proton exchange membrane fuel cell

International Nuclear Information System (INIS)

Ferreira, Rui B.; Falcão, D.S.; Oliveira, V.B.; Pinto, A.M.F.R.

2017-01-01

Highlights: •A 1D + 3D model of a PEM fuel cell is described and experimentally validated. •VOF method tracks the two-phase flow and electrochemical reactions are considered. •Water dynamics inside a serpentine channel is analyzed for different voltages. •Water content in different regions of channel is quantified. •Important issues on coupling of the VOF model with electrochemical reactions are addressed. -- Abstract: In this work, a numerical model of a proton exchange membrane (PEM) fuel cell is presented. The volume of fluid (VOF) method is employed to simulate the air-water two-phase flow in the cathode gas channel, at the same time that the cell electrochemical performance is predicted. The model is validated against an experimental polarization curve and through the visualization of water distribution inside a transparent fuel cell. The water dynamics inside a serpentine gas channel is numerically analyzed under different operating voltages. Moreover, water content in different regions of the channel is quantified. Current density and water generation rate spatial distributions are also displayed and it is shown how they affect the process of water emergence into the gas channel. Important issues on the simulation of the PEM fuel cells two-phase flow are addressed, especially concerning the coupling of the VOF technique with electrochemical reactions. Both the model and the numerical results aim to contribute to a better understanding of the two-phase flow phenomenon that occurs in these devices.

3D Dynamic Modeling of the Head-Neck Complex for Fast Eye and Head Orientation Movements Research

Directory of Open Access Journals (Sweden)

Daniel A. Sierra

2011-01-01

Full Text Available A 3D dynamic computer model for the movement of the head-neck complex is presented. It incorporates anatomically correct information about the diverse elements forming the system. The skeleton is considered as a set of interconnected rigid 3D bodies following the Newton-Euler laws of movement. The muscles are modeled using Enderle's linear model, which shows equivalent dynamic characteristics to Loeb's virtual muscle model. The soft tissues, namely, the ligaments, intervertebral disks, and facet joints, are modeled considering their physiological roles and dynamics. In contrast with other head and neck models developed for safety research, the model is aimed to study the neural control of the complex during fast eye and head movements, such as saccades and gaze shifts. In particular, the time-optimal hypothesis and the feedback control ones are discussed.

A 3-D model of tumor progression based on complex automata driven by particle dynamics.

Science.gov (United States)

Wcisło, Rafał; Dzwinel, Witold; Yuen, David A; Dudek, Arkadiusz Z

2009-12-01

The dynamics of a growing tumor involving mechanical remodeling of healthy tissue and vasculature is neglected in most of the existing tumor models. This is due to the lack of efficient computational framework allowing for simulation of mechanical interactions. Meanwhile, just these interactions trigger critical changes in tumor growth dynamics and are responsible for its volumetric and directional progression. We describe here a novel 3-D model of tumor growth, which combines particle dynamics with cellular automata concept. The particles represent both tissue cells and fragments of the vascular network. They interact with their closest neighbors via semi-harmonic central forces simulating mechanical resistance of the cell walls. The particle dynamics is governed by both the Newtonian laws of motion and the cellular automata rules. These rules can represent cell life-cycle and other biological interactions involving smaller spatio-temporal scales. We show that our complex automata, particle based model can reproduce realistic 3-D dynamics of the entire system consisting of the tumor, normal tissue cells, blood vessels and blood flow. It can explain phenomena such as the inward cell motion in avascular tumor, stabilization of tumor growth by the external pressure, tumor vascularization due to the process of angiogenesis, trapping of healthy cells by invading tumor, and influence of external (boundary) conditions on the direction of tumor progression. We conclude that the particle model can serve as a general framework for designing advanced multiscale models of tumor dynamics and it is very competitive to the modeling approaches presented before.

In silico modelling and molecular dynamics simulation studies of thiazolidine based PTP1B inhibitors.

Science.gov (United States)

Mahapatra, Manoj Kumar; Bera, Krishnendu; Singh, Durg Vijay; Kumar, Rajnish; Kumar, Manoj

2018-04-01

Protein tyrosine phosphatase 1B (PTP1B) has been identified as a negative regulator of insulin and leptin signalling pathway; hence, it can be considered as a new therapeutic target of intervention for the treatment of type 2 diabetes. Inhibition of this molecular target takes care of both diabetes and obesity, i.e. diabestiy. In order to get more information on identification and optimization of lead, pharmacophore modelling, atom-based 3D QSAR, docking and molecular dynamics studies were carried out on a set of ligands containing thiazolidine scaffold. A six-point pharmacophore model consisting of three hydrogen bond acceptor (A), one negative ionic (N) and two aromatic rings (R) with discrete geometries as pharmacophoric features were developed for a predictive 3D QSAR model. The probable binding conformation of the ligands within the active site was studied through molecular docking. The molecular interactions and the structural features responsible for PTP1B inhibition and selectivity were further supplemented by molecular dynamics simulation study for a time scale of 30 ns. The present investigation has identified some of the indispensible structural features of thiazolidine analogues which can further be explored to optimize PTP1B inhibitors.

Multi-coupling dynamic model and 3d simulation program for in-situ leaching of uranium mining

International Nuclear Information System (INIS)

Tan Kaixuan; Zeng Sheng; Sang Xiao; Sun Bing

2010-01-01

The in-situ leaching of uranium mining is a very complicated non-linear dynamic system, which involves couplings and positive/negative feedback among many factors and processes. A comprehensive, coupled multi-factors and processes dynamic model and simulation method was established to study the in-situ leaching of uranium mining. The model accounts for most coupling among various processes as following: (1) rock texture mechanics and its evolution, (2)the incremental stress rheology of rock deformation, (3) 3-D viscoelastic/ plastic multi-deformation processes, (4) hydrofracturing, (5) tensorial (anisotropic) fracture and rock permeability, (6) water-rock interactions and mass-transport (both advective and diffusive), (7) dissolution-induced chemical compaction, (8) multi-phase fluid flow. A 3-D simulation program was compiled based on Fortran and C++. An example illustrating the application of this model to simulating acidification, production and terminal stage of in situ leaching of uranium mining is presented for the some mine in Xinjiang, China. This model and program can be used for theoretical study, mine design, production management, the study of contaminant transport and restoration in groundwater of in-situ leaching of uranium mining. (authors)

Developing a Dynamic Pharmacophore Model for HIV-1 Integrase

International Nuclear Information System (INIS)

Carlson, Heather A.; Masukawa, Keven M.; Rubins, Kathleen; Bushman, Frederic; Jorgensen, William L.; Lins, Roberto; Briggs, James; Mccammon, Andy

2000-01-01

We present the first receptor-based pharmacophore model for HIV-1 integrase. The development of ''dynamic'' pharmacophore models is a new method that accounts for the inherent flexibility of the active site and aims to reduce the entropic penalties associated with binding a ligand. Furthermore, this new drug discovery method overcomes the limitation of an incomplete crystal structure of the target protein. A molecular dynamics (MD) simulation describes the flexibility of the uncomplexed protein. Many conformational models of the protein are saved from the MD simulations and used in a series of multi-unit search for interacting conformers (MUSIC) simulations. MUSIC is a multiple-copy minimization method, available in the BOSS program; it is used to determine binding regions for probe molecules containing functional groups that complement the active site. All protein conformations from the MD are overlaid, and conserved binding regions for the probe molecules are identified. Those conserved binding regions define the dynamic pharmacophore model. Here, the dynamic model is compared to known inhibitors of the integrase as well as a three-point, ligand-based pharmacophore model from the literature. Also, a ''static'' pharmacophore model was determined in the standard fashion, using a single crystal structure. Inhibitors thought to bind in the active site of HIV-1 integrase fit the dynamic model but not the static model. Finally, we have identified a set of compounds from the Available Chemicals Directory that fit the dynamic pharmacophore model, and experimental testing of the compounds has confirmed several new inhibitors

Comparison between a coupled 1D-2D model and a fully 2D model for supercritical flow simulation in crossroads

KAUST Repository

Ghostine, Rabih

2014-12-01

In open channel networks, flow is usually approximated by the one-dimensional (1D) Saint-Venant equations coupled with an empirical junction model. In this work, a comparison in terms of accuracy and computational cost between a coupled 1D-2D shallow water model and a fully two-dimensional (2D) model is presented. The paper explores the ability of a coupled model to simulate the flow processes during supercritical flows in crossroads. This combination leads to a significant reduction in the computational time, as a 1D approach is used in branches and a 2D approach is employed in selected areas only where detailed flow information is essential. Overall, the numerical results suggest that the coupled model is able to accurately simulate the main flow processes. In particular, hydraulic jumps, recirculation zones, and discharge distribution are reasonably well reproduced and clearly identified. Overall, the proposed model leads to a 30% reduction in run times. © 2014 International Association for Hydro-Environment Engineering and Research.

Boundary Control of Linear Uncertain 1-D Parabolic PDE Using Approximate Dynamic Programming.

Science.gov (United States)

Talaei, Behzad; Jagannathan, Sarangapani; Singler, John

2018-04-01

This paper develops a near optimal boundary control method for distributed parameter systems governed by uncertain linear 1-D parabolic partial differential equations (PDE) by using approximate dynamic programming. A quadratic surface integral is proposed to express the optimal cost functional for the infinite-dimensional state space. Accordingly, the Hamilton-Jacobi-Bellman (HJB) equation is formulated in the infinite-dimensional domain without using any model reduction. Subsequently, a neural network identifier is developed to estimate the unknown spatially varying coefficient in PDE dynamics. Novel tuning law is proposed to guarantee the boundedness of identifier approximation error in the PDE domain. A radial basis network (RBN) is subsequently proposed to generate an approximate solution for the optimal surface kernel function online. The tuning law for near optimal RBN weights is created, such that the HJB equation error is minimized while the dynamics are identified and closed-loop system remains stable. Ultimate boundedness (UB) of the closed-loop system is verified by using the Lyapunov theory. The performance of the proposed controller is successfully confirmed by simulation on an unstable diffusion-reaction process.

Regulation of dopamine D1 receptor dynamics within the postsynaptic density of hippocampal glutamate synapses.

Directory of Open Access Journals (Sweden)

Laurent Ladepeche

Full Text Available Dopamine receptor potently modulates glutamate signalling, synaptic plasticity and neuronal network adaptations in various pathophysiological processes. Although key intracellular signalling cascades have been identified, the cellular mechanism by which dopamine and glutamate receptor-mediated signalling interplay at glutamate synapse remain poorly understood. Among the cellular mechanisms proposed to aggregate D1R in glutamate synapses, the direct interaction between D1R and the scaffold protein PSD95 or the direct interaction with the glutamate NMDA receptor (NMDAR have been proposed. To tackle this question we here used high-resolution single nanoparticle imaging since it provides a powerful way to investigate at the sub-micron resolution the dynamic interaction between these partners in live synapses. We demonstrate in hippocampal neuronal networks that dopamine D1 receptors (D1R laterally diffuse within glutamate synapses, in which their diffusion is reduced. Disrupting the interaction between D1R and PSD95, through genetical manipulation and competing peptide, did not affect D1R dynamics in glutamatergic synapses. However, preventing the physical interaction between D1R and the GluN1 subunit of NMDAR abolished the synaptic stabilization of diffusing D1R. Together, these data provide direct evidence that the interaction between D1R and NMDAR in synapses participate in the building of the dopamine-receptor-mediated signalling, and most likely to the glutamate-dopamine cross-talk.

3D dynamic pituitary MR imaging with CAIPIRINHA: Initial experience and comparison with 2D dynamic MR imaging

Energy Technology Data Exchange (ETDEWEB)

Fushimi, Yasutaka, E-mail: yfushimi@kuhp.kyoto-u.ac.jp [Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan); Okada, Tomohisa; Kanagaki, Mitsunori; Yamamoto, Akira; Kanda, Yumiko; Sakamoto, Ryo [Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan); Hojo, Masato; Takahashi, Jun C.; Miyamoto, Susumu [Department of Neurosurgery, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan); Togashi, Kaori [Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan)

2014-10-15

Objectives: To evaluate the validity of 3D dynamic pituitary MR imaging with controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA), with special emphasis on demarcation of pituitary posterior lobe and stalk. Methods: Participants comprised 32 patients who underwent dynamic pituitary MR imaging due to pituitary or parasellar lesions. 3D dynamic MR with CAIPIRINHA was performed at 3 T with 20-s-interval, precontrast, 1st to 5th dynamic images. Normalized values and enhanced ratios (dynamic postcontrast image values divided by precontrast ones) were compared between 3D and 2D dynamic MR imaging for patients with visual identification of posterior lobe and stalk. Results: In 3D, stalk was identified in 29 patients and unidentified in 3, and posterior lobe was identified in 28 and unidentified in 4. In 2D, stalk was identified in 26 patients and unidentified in 6 patients, and posterior lobe was identified in 15 and unidentified in 17. Normalized values of pituitary posterior lobe and stalk were higher in 3D than 2D (P < 0.001). No significant difference in enhancement ratio was seen between 3D and 2D. Conclusions: 3D dynamic pituitary MR provided better identification and higher normalized values of pituitary posterior lobe and stalk than 2D.

Development of a 3D consistent 1D neutronics model for reactor core simulation

International Nuclear Information System (INIS)

Lee, Ki Bog; Joo, Han Gyu; Cho, Byung Oh; Zee, Sung Quun

2001-02-01

In this report a 3D consistent 1D model based on nonlinear analytic nodal method is developed to reproduce the 3D results. During the derivation, the current conservation factor (CCF) is introduced which guarantees the same axial neutron currents obtained from the 1D equation as the 3D reference values. Furthermore in order to properly use 1D group constants, a new 1D group constants representation scheme employing tables for the fuel temperature, moderator density and boron concentration is developed and functionalized for the control rod tip position. To test the 1D kinetics model with CCF, several steady state and transient calculations were performed and compared with 3D reference values. The errors of K-eff values were reduced about one tenth when using CCF without significant computational overhead. And the errors of power distribution were decreased to the range of one fifth or tenth at steady state calculation. The 1D kinetics model with CCF and the 1D group constant functionalization employing tables as a function of control rod tip position can provide preciser results at the steady state and transient calculation. Thus it is expected that the 1D kinetics model derived in this report can be used in the safety analysis, reactor real time simulation coupled with system analysis code, operator support system etc.

System dynamics model developement for R and D budget allocation strategy

International Nuclear Information System (INIS)

Kwack, S. M.; Kim, D. H.; Lee, Y. S.; Jung, M. T.

2003-01-01

A computer simulation model was developed for R and D budget allocation problems that usually contain very complex and non-linear social issues. The System Dynamics approach was employed, which is proper to complex and non-linear social problem modeling. An issue of budget allocation to each step in five research areas was analyzed for an application example. The base scenario, that assumes to allocate a large portion of budget to demonstration step, was found to have a weakness in long-term sense. To overcome this weakness, some other better alternatives were recommended through the analysis. In addition, this paper suggests the ways to utilize the updated model in the future

Analytic observations for the d=1+ 1 bridge site (or single-step) deposition model

International Nuclear Information System (INIS)

Evans, J.W.; Kang, H.C.

1991-01-01

Some exact results for a reversible version of the d=1+1 bridge site (or single-step) deposition model are presented. Exact steady-state properties are determined directly for finite systems with various mean slopes. These show explicitly how the asymptotic growth velocity and fluctuations are quenched as the slope approaches its maximum allowed value. Next, exact hierarchial equations for the dynamics are presented. For the special case of ''equilibrium growth,'' these are analyzed exactly at the pair-correlation level directly for an infinite system. This provided further insight into asymptotic scaling behavior. Finally, the above hierarchy is compared with one generated from a discrete form of the Kardar--Parisi--Zhang equations. Some differences are described

Does one-dimensional (1D) adatom and cluster diffusion of Pt on the Pt(110)-(1 x 2) surface lead to 1D ripening?

International Nuclear Information System (INIS)

Linderoth, T R; Horch, S; Petersen, L; Laegsgaard, E; Stensgaard, I; Besenbacher, F

2005-01-01

The technique of scanning tunnelling microscopy (STM) uniquely allows dynamic processes on surfaces to be followed directly in real space and at atomic resolution. Results for the 551225 surface diffusion of Pt adatoms and clusters on the anisotropic, missing row reconstructed Pt(110)-(1 x 2) surface are briefly reviewed. Mass transport in this system is entirely one-dimensional (1D) since, at low adatom coverage, atoms and clusters are confined to the missing row troughs. In this paper, we therefore address the question if Pt/Pt(110)-(1 x 2) is a 1D model system to study late stage growth phenomena such as island ripening? From STM measurements, we quantify the morphology changes resulting from annealing a surface configuration with small 1D Pt islands in the missing row troughs to temperatures in the interval 369-395 K. Interestingly, the resulting increase in island sizes (ripening) cannot be accounted for by the known island and adatom mobilities within a 1D model. An explanation is provided from dynamic, time-resolved 'STM-movies' that directly reveal two novel island-mediated mechanisms for inter-trough mass transport which cause the Pt/Pt(110)-(1 x 2) system not to be purely 1D at the higher surface coverage used in the annealing experiments

Comparison between 3D and 1D simulations of a regenerative blower for fuel cell applications

International Nuclear Information System (INIS)

Badami, M.; Mura, M.

2012-01-01

Highlights: ► A hydrogen recirculation blower for automotive fuel cells applications is studied. ► A 3D CFD analysis has been carried out to better understand the internal flows of the machine. ► The CFD results are compared to a 1D model set up by the authors in previous works. ► The main hypotheses put forward for the theoretical 1D model are compatible with the 3D analysis. - Abstract: A 3D Computational Fluid Dynamics (CFD) analysis has been carried out to better understand the internal fluid dynamics of a regenerative blower used for hydrogen recirculation in a Proton Exchange Membrane (PEM), Fuel Cell (FC) utilized for automotive applications. The obtained results are used to highlight the motion of the fluid in the vanes and in the side channel of the machine and to verify the main hypotheses put forward concerning the theoretical 1D model set up by the authors in previous works on the basis of the momentum exchange theory. Finally, the CFD analysis has been used to point out the effect of the slope of the vanes on the performance of the regenerative blower, and the results have been compared with those obtained using of the 1D model.

NanoSail - D Orbital and Attitude Dynamics

Science.gov (United States)

Heaton, Andrew F.; Faller, Brent F.; Katan, Chelsea K.

2013-01-01

NanoSail-D unfurled January 20th, 2011 and successfully demonstrated the deployment and deorbit capability of a solar sail in low Earth orbit. The orbit was strongly perturbed by solar radiation pressure, aerodynamic drag, and oblate gravity which were modeled using STK HPOP. A comparison of the ballistic coefficient history to the orbit parameters exhibits a strong relationship between orbital lighting, the decay rate of the mean semi-major axis and mean eccentricity. A similar comparison of mean solar area using the STK HPOP solar radiation pressure model exhibits a strong correlation of solar radiation pressure to mean eccentricity and mean argument of perigee. NanoSail-D was not actively controlled and had no capability on-board for attitude or orbit determination. To estimate attitude dynamics we created a 3-DOF attitude dynamics simulation that incorporated highly realistic estimates of perturbing forces into NanoSail-D torque models. By comparing the results of this simulation to the orbital behavior and ground observations of NanoSail-D, we conclude that there is a coupling between the orbit and attitude dynamics as well as establish approximate limits on the location of the NanoSail-D solar center of pressure. Both of these observations contribute valuable data for future solar sail designs and missions.

In Silico Exploration of 1,7-Diazacarbazole Analogs as Checkpoint Kinase 1 Inhibitors by Using 3D QSAR, Molecular Docking Study, and Molecular Dynamics Simulations

Directory of Open Access Journals (Sweden)

Xiaodong Gao

2016-05-01

Full Text Available Checkpoint kinase 1 (Chk1 is an important serine/threonine kinase with a self-protection function. The combination of Chk1 inhibitors and anti-cancer drugs can enhance the selectivity of tumor therapy. In this work, a set of 1,7-diazacarbazole analogs were identified as potent Chk1 inhibitors through a series of computer-aided drug design processes, including three-dimensional quantitative structure–activity relationship (3D-QSAR modeling, molecular docking, and molecular dynamics simulations. The optimal QSAR models showed significant cross-validated correlation q2 values (0.531, 0.726, fitted correlation r2 coefficients (higher than 0.90, and standard error of prediction (less than 0.250. These results suggested that the developed models possess good predictive ability. Moreover, molecular docking and molecular dynamics simulations were applied to highlight the important interactions between the ligand and the Chk1 receptor protein. This study shows that hydrogen bonding and electrostatic forces are key interactions that confer bioactivity.

Dynamic origins of fermionic D -terms

Science.gov (United States)

Hudson, Jonathan; Schweitzer, Peter

2018-03-01

The D -term is defined through matrix elements of the energy-momentum tensor, similarly to mass and spin, yet this important particle property is experimentally not known any fermion. In this work we show that the D -term of a spin 1/2 fermion is of dynamical origin: it vanishes for a free fermion. This is in pronounced contrast to the bosonic case where already a free spin-0 boson has a non-zero intrinsic D -term. We illustrate in two simple models how interactions generate the D -term of a fermion with an internal structure, the nucleon. All known matter is composed of elementary fermions. This indicates the importance to study this interesting particle property in more detail, which will provide novel insights especially on the structure of the nucleon.

IMPAC-2, Dynamic Impact Analysis for 1-D Nonlinear Spring Shipping Container Model

International Nuclear Information System (INIS)

Payne, J. B.; Counts, J.

1980-01-01

1 - Description of problem or function: IMPAC2 solves the equations of motion for a one-dimensional, lumped-mass, nonlinear spring mathematical container model. The program was designed to analyze the dynamic response of metallic shipping containers impacting an unyielding surface. The container may consist of several hollow concentric cylinders, each of a different material and length. 2 - Method of solution: An iteration technique (Newmark-Beta) is used to solve the equations of motion for displacement, velocity, and acceleration of the masses at each time-step. 3 - Restrictions on the complexity of the problem: The maximum number of masses is 100. The program is written for a maximum of 10 different materials. The first seven have been given the following designations: 1 steel, 2 clad lead, 3 uranium, 4 pine, 5 lead, 6 balsa, 7 marine plywood. Stress-strain characteristics for materials 1, 3, 4, and 6 are approximated by a bilinear curve fit which requires as input the modulus of elasticity, the stress at the intersection of the two lines and the slope of the second line

NASTRAN analysis of the 1/8-scale space shuttle dynamic model

Science.gov (United States)

Bernstein, M.; Mason, P. W.; Zalesak, J.; Gregory, D. J.; Levy, A.

1973-01-01

The space shuttle configuration has more complex structural dynamic characteristics than previous launch vehicles primarily because of the high model density at low frequencies and the high degree of coupling between the lateral and longitudinal motions. An accurate analytical representation of these characteristics is a primary means for treating structural dynamics problems during the design phase of the shuttle program. The 1/8-scale model program was developed to explore the adequacy of available analytical modeling technology and to provide the means for investigating problems which are more readily treated experimentally. The basic objectives of the 1/8-scale model program are: (1) to provide early verification of analytical modeling procedures on a shuttle-like structure, (2) to demonstrate important vehicle dynamic characteristics of a typical shuttle design, (3) to disclose any previously unanticipated structural dynamic characteristics, and (4) to provide for development and demonstration of cost effective prototype testing procedures.

Research Summary 3-D Computational Fluid Dynamics (CFD) Model Of The Human Respiratory System

Science.gov (United States)

The U.S. EPA’s Office of Research and Development (ORD) has developed a 3-D computational fluid dynamics (CFD) model of the human respiratory system that allows for the simulation of particulate based contaminant deposition and clearance, while being adaptable for age, ethnicity,...

Development of 1D Liner Compression Code for IDL

Science.gov (United States)

Shimazu, Akihisa; Slough, John; Pancotti, Anthony

2015-11-01

A 1D liner compression code is developed to model liner implosion dynamics in the Inductively Driven Liner Experiment (IDL) where FRC plasmoid is compressed via inductively-driven metal liners. The driver circuit, magnetic field, joule heating, and liner dynamics calculations are performed at each time step in sequence to couple these effects in the code. To obtain more realistic magnetic field results for a given drive coil geometry, 2D and 3D effects are incorporated into the 1D field calculation through use of correction factor table lookup approach. Commercial low-frequency electromagnetic fields solver, ANSYS Maxwell 3D, is used to solve the magnetic field profile for static liner condition at various liner radius in order to derive correction factors for the 1D field calculation in the code. The liner dynamics results from the code is verified to be in good agreement with the results from commercial explicit dynamics solver, ANSYS Explicit Dynamics, and previous liner experiment. The developed code is used to optimize the capacitor bank and driver coil design for better energy transfer and coupling. FRC gain calculations are also performed using the liner compression data from the code for the conceptual design of the reactor sized system for fusion energy gains.

Dynamically accumulated dose and 4D accumulated dose for moving tumors

International Nuclear Information System (INIS)

Li Heng; Li Yupeng; Zhang Xiaodong; Li Xiaoqiang; Liu Wei; Gillin, Michael T.; Zhu, X. Ronald

2012-01-01

Purpose: The purpose of this work was to investigate the relationship between dynamically accumulated dose (dynamic dose) and 4D accumulated dose (4D dose) for irradiation of moving tumors, and to quantify the dose uncertainty induced by tumor motion. Methods: The authors established that regardless of treatment modality and delivery properties, the dynamic dose will converge to the 4D dose, instead of the 3D static dose, after multiple deliveries. The bounds of dynamic dose, or the maximum estimation error using 4D or static dose, were established for the 4D and static doses, respectively. Numerical simulations were performed (1) to prove the principle that for each phase, after multiple deliveries, the average number of deliveries for any given time converges to the total number of fractions (K) over the number of phases (N); (2) to investigate the dose difference between the 4D and dynamic doses as a function of the number of deliveries for deliveries of a “pulsed beam”; and (3) to investigate the dose difference between 4D dose and dynamic doses as a function of delivery time for deliveries of a “continuous beam.” A Poisson model was developed to estimate the mean dose error as a function of number of deliveries or delivered time for both pulsed beam and continuous beam. Results: The numerical simulations confirmed that the number of deliveries for each phase converges to K/N, assuming a random starting phase. Simulations for the pulsed beam and continuous beam also suggested that the dose error is a strong function of the number of deliveries and/or total deliver time and could be a function of the breathing cycle, depending on the mode of delivery. The Poisson model agrees well with the simulation. Conclusions: Dynamically accumulated dose will converge to the 4D accumulated dose after multiple deliveries, regardless of treatment modality. Bounds of the dynamic dose could be determined using quantities derived from 4D doses, and the mean dose

Comparison between a coupled 1D-2D model and a fully 2D model for supercritical flow simulation in crossroads

KAUST Repository

Ghostine, Rabih; Hoteit, Ibrahim; Vazquez, Jose; Terfous, Abdelali; Ghenaim, Abdellah; Mose, Robert

2014-01-01

In open channel networks, flow is usually approximated by the one-dimensional (1D) Saint-Venant equations coupled with an empirical junction model. In this work, a comparison in terms of accuracy and computational cost between a coupled 1D-2D

Modeling tree crown dynamics with 3D partial differential equations.

Science.gov (United States)

Beyer, Robert; Letort, Véronique; Cournède, Paul-Henry

2014-01-01

We characterize a tree's spatial foliage distribution by the local leaf area density. Considering this spatially continuous variable allows to describe the spatiotemporal evolution of the tree crown by means of 3D partial differential equations. These offer a framework to rigorously take locally and adaptively acting effects into account, notably the growth toward light. Biomass production through photosynthesis and the allocation to foliage and wood are readily included in this model framework. The system of equations stands out due to its inherent dynamic property of self-organization and spontaneous adaptation, generating complex behavior from even only a few parameters. The density-based approach yields spatially structured tree crowns without relying on detailed geometry. We present the methodological fundamentals of such a modeling approach and discuss further prospects and applications.

Electro-mechanical dynamics of spiral waves in a discrete 2D model of human atrial tissue.

Directory of Open Access Journals (Sweden)

Paul Brocklehurst

Full Text Available We investigate the effect of mechano-electrical feedback and atrial fibrillation induced electrical remodelling (AFER of cellular ion channel properties on the dynamics of spiral waves in a discrete 2D model of human atrial tissue. The tissue electro-mechanics are modelled using the discrete element method (DEM. Millions of bonded DEM particles form a network of coupled atrial cells representing 2D cardiac tissue, allowing simulations of the dynamic behaviour of electrical excitation waves and mechanical contraction in the tissue. In the tissue model, each cell is modelled by nine particles, accounting for the features of individual cellular geometry; and discrete inter-cellular spatial arrangement of cells is also considered. The electro-mechanical model of a human atrial single-cell was constructed by strongly coupling the electrophysiological model of Colman et al. to the mechanical myofilament model of Rice et al., with parameters modified based on experimental data. A stretch-activated channel was incorporated into the model to simulate the mechano-electrical feedback. In order to investigate the effect of mechano-electrical feedback on the dynamics of spiral waves, simulations of spiral waves were conducted in both the electromechanical model and the electrical-only model in normal and AFER conditions, to allow direct comparison of the results between the models. Dynamics of spiral waves were characterized by tracing their tip trajectories, stability, excitation frequencies and meandering range of tip trajectories. It was shown that the developed DEM method provides a stable and efficient model of human atrial tissue with considerations of the intrinsically discrete and anisotropic properties of the atrial tissue, which are challenges to handle in traditional continuum mechanics models. This study provides mechanistic insights into the complex behaviours of spiral waves and the genesis of atrial fibrillation by showing an important role of

Electro-mechanical dynamics of spiral waves in a discrete 2D model of human atrial tissue.

Science.gov (United States)

Brocklehurst, Paul; Ni, Haibo; Zhang, Henggui; Ye, Jianqiao

2017-01-01

We investigate the effect of mechano-electrical feedback and atrial fibrillation induced electrical remodelling (AFER) of cellular ion channel properties on the dynamics of spiral waves in a discrete 2D model of human atrial tissue. The tissue electro-mechanics are modelled using the discrete element method (DEM). Millions of bonded DEM particles form a network of coupled atrial cells representing 2D cardiac tissue, allowing simulations of the dynamic behaviour of electrical excitation waves and mechanical contraction in the tissue. In the tissue model, each cell is modelled by nine particles, accounting for the features of individual cellular geometry; and discrete inter-cellular spatial arrangement of cells is also considered. The electro-mechanical model of a human atrial single-cell was constructed by strongly coupling the electrophysiological model of Colman et al. to the mechanical myofilament model of Rice et al., with parameters modified based on experimental data. A stretch-activated channel was incorporated into the model to simulate the mechano-electrical feedback. In order to investigate the effect of mechano-electrical feedback on the dynamics of spiral waves, simulations of spiral waves were conducted in both the electromechanical model and the electrical-only model in normal and AFER conditions, to allow direct comparison of the results between the models. Dynamics of spiral waves were characterized by tracing their tip trajectories, stability, excitation frequencies and meandering range of tip trajectories. It was shown that the developed DEM method provides a stable and efficient model of human atrial tissue with considerations of the intrinsically discrete and anisotropic properties of the atrial tissue, which are challenges to handle in traditional continuum mechanics models. This study provides mechanistic insights into the complex behaviours of spiral waves and the genesis of atrial fibrillation by showing an important role of the mechano

Structure and dynamics of a [1:1] drug-DNA complex: Analysis of 2D NMR data using molecular mechanics and molecular dynamics calculations

International Nuclear Information System (INIS)

Sarma, R.H.; Sarma, M.H.; Umemoto, K.

1990-01-01

1D/2D NMR studies are reported for a [1:1] complex of d(GA 4 T 4 C) 2 and Dst2 (an analogue of distamycin A). Full- Matrix NOESY Simulations, Molecular Mechanics and Molecular Dynamics Calculations are performed to analyze the NMR data. Results show that drug-DNA complex formation is driven by static features like H-bonding and steric interactions in the minor-groove of DNA. As a consequence of drug binding, a non-linear oscillatory mode is activated. In this mode the molecule samples equilibrium structural states of difference degrees of bending. It is noted that these structures belong to three distinctly different energy wells that satisfy the same NMR data. 14 refs., 4 figs., 2 tabs

FLOWPLOT2, 2-D, 3-D Fluid Dynamic Plots

International Nuclear Information System (INIS)

Cobb, C.K.; Tunstall, J.N.

1989-01-01

1 - Description of program or function: FLOWPLOT2 is a plotting program used with numerical or analytical fluid dynamics codes to create velocity vector plots, contour plots of up to three fluid parameters (e.g. pressure, density, and temperature), two-dimensional profile plots, three-dimensional curve plots, and/or three-dimensional surface plots for either the u or v velocity components. If the fluid dynamics code computes a transient or simulated time related solution, FLOWPLOT2 can also be used to generate these plots for any specified time interval. Multiple cases generating different plots for different time intervals may be run in one execution of the program. In addition, plots can be created for selected two- dimensional planes of three-dimensional steady-state problems. The user has the option of producing plots on CalComp or Versatec plotters or microfiche and of creating a compressed dataset before plotting. 2 - Method of solution: FLOWPLOT2 reads a dataset written by the fluid dynamics code. This dataset must be written in a specified format and must contain parametric data at the nodal points of a uniform or non-uniform rectangular grid formed by the intersection of the grid lines of the model. 3 - Restrictions on the complexity of the problem - Maxima of: 2500 nodes, 40 y-values for 2-D profile plots and 3-D curve plots, 20 contour values, 3 fluid parameters

Dynamic 3D MR-defecography

Energy Technology Data Exchange (ETDEWEB)

Ratz, V.; Wech, T.; Schindele, A.; Dierks, A.; Sauer, A.; Reibetanz, J.; Borzi, A.; Bley, T.; Koestler, H.

2016-09-15

Epidemiological studies have estimated the incidence of chronic constipation to be up to 27% of the general population. The gold standard to evaluate affected patients is the dynamic entero-colpo-cysto-defecography. In the clinical routine 2 D MR-defecography is also performed, but only one to three 2 D slices at a temporal footprint of about one second are acquired. To improve the detection of lateral localized pathologies, we developed and implemented dynamic 3 D MR-defecography. Each 3 D block consisted of seven slices with an in-plane spatial resolution of 1.3 x 1.3 mm{sup 2} to 2.3 x 2.3 mm{sup 2} and an image update rate between 0.8 s and 1.3 s. We used a fast bSSFP sequence with a modified stack-of-stars sampling scheme for data acquisition and a modified FISTA compressed sensing algorithm to reconstruct the undersampled datasets. We performed a study including 6 patients to optimize the acquisition parameters with respect to image quality.

Dynamic 3D MR-defecography

International Nuclear Information System (INIS)

Ratz, V.; Wech, T.; Schindele, A.; Dierks, A.; Sauer, A.; Reibetanz, J.; Borzi, A.; Bley, T.; Koestler, H.

2016-01-01

Epidemiological studies have estimated the incidence of chronic constipation to be up to 27% of the general population. The gold standard to evaluate affected patients is the dynamic entero-colpo-cysto-defecography. In the clinical routine 2 D MR-defecography is also performed, but only one to three 2 D slices at a temporal footprint of about one second are acquired. To improve the detection of lateral localized pathologies, we developed and implemented dynamic 3 D MR-defecography. Each 3 D block consisted of seven slices with an in-plane spatial resolution of 1.3 x 1.3 mm 2 to 2.3 x 2.3 mm 2 and an image update rate between 0.8 s and 1.3 s. We used a fast bSSFP sequence with a modified stack-of-stars sampling scheme for data acquisition and a modified FISTA compressed sensing algorithm to reconstruct the undersampled datasets. We performed a study including 6 patients to optimize the acquisition parameters with respect to image quality.

Comparison and analysis of 1D/2D/3D neutronics modeling for a fusion reactor

International Nuclear Information System (INIS)

Li, J.; Zeng, Q.; Chen, M.; Jiang, J.; Wu, Y.

2007-01-01

During the course of analyzing the characteristics for fusion reactors, the refined calculations are needed to confirm that the nuclear design requirements are met. Since the long computational time is consumed, the refined three-dimensional (3D) representation has been used primarily for establishing the baseline reference values, analyzing problems which cannot be reduced by symmetry considerations to lower dimensions, or where a high level of accuracy is desired locally. The two-dimensional (2D) or one-dimensional (1D) description leads itself readily to resolve many problems, such as the studies for the material fraction optimization, or for the blanket size optimization. The purpose of this paper is to find out the differences among different geometric descriptions, which can guide the way to approximate and simplify the computational model. The fusion power reactor named FDS-II was designed as an advanced fusion power reactor to demonstrate and validate the commercialization of fusion power by Institute of Plasma Physics, Chinese Academy of Science. In this contribution, the dual-cooled lithium lead (DLL) blanket of FDS-II was used as a reference for neutronics comparisons and analyses. The geometric descriptions include 1D concentric sphere model, 1D, 2D and 3D cylinder models. The home-developed multi-functional neutronics analysis code system VisualBUS, the Monte Carlo transport code MCNP and nuclear data library HENDL have been used for these analyses. The neutron wall loading distribution, tritium breeding ratio (TBR) and nuclear heat were calculated to evaluate the nuclear performance. The 3D calculation has been used as a comparison reference because it has the least errors in the treatment of geometry. It is suggested that the value of TBR calculated by the 1D approach should be greater than 1.3 to satisfy the practical need of tritium self-sufficiency. The distribution of nuclear heat based on the 2D and 3D models were similar since they all consider

Stochastic dynamics of phase singularities under ventricular fibrillation in 2D Beeler-Reuter model

Directory of Open Access Journals (Sweden)

Akio Suzuki

2011-09-01

Full Text Available 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.

Modeling of impurity spectroscopy in the divertor and SOL of DIII-D using the 1D multifluid model NEWT1D

International Nuclear Information System (INIS)

West, W.P.; Evans, T.E.; Brooks, N.H.

1996-10-01

NEWT1D, a one dimensional multifluid model of the scrape-off layer and divertor plasma, has been used to model the plasma including the distribution of carbon ionization states in the SOL and divertor of ELMing H-mode at two injected power levels in DIII-D. Comparison of the code predictions to the measured divertor and scrape-off layer (SOL) plasma density and temperature shows good agreement. Comparison of the predicted line emissions to the spectroscopic data suggests that physically sputtered carbon from the strike point is not transported up the flux tube; a distributed source of carbon a few centimeters up the flux tube is required to achieve reasonable agreement

Application of an unstructured 3D finite volume numerical model to flows and salinity dynamics in the San Francisco Bay-Delta

Science.gov (United States)

Martyr-Koller, R.C.; Kernkamp, H.W.J.; Van Dam, Anne A.; Mick van der Wegen,; Lucas, Lisa; Knowles, N.; Jaffe, B.; Fregoso, T.A.

2017-01-01

A linked modeling approach has been undertaken to understand the impacts of climate and infrastructure on aquatic ecology and water quality in the San Francisco Bay-Delta region. The Delft3D Flexible Mesh modeling suite is used in this effort for its 3D hydrodynamics, salinity, temperature and sediment dynamics, phytoplankton and water-quality coupling infrastructure, and linkage to a habitat suitability model. The hydrodynamic model component of the suite is D-Flow FM, a new 3D unstructured finite-volume model based on the Delft3D model. In this paper, D-Flow FM is applied to the San Francisco Bay-Delta to investigate tidal, seasonal and annual dynamics of water levels, river flows and salinity under historical environmental and infrastructural conditions. The model is driven by historical winds, tides, ocean salinity, and river flows, and includes federal, state, and local freshwater withdrawals, and regional gate and barrier operations. The model is calibrated over a 9-month period, and subsequently validated for water levels, flows, and 3D salinity dynamics over a 2 year period.Model performance was quantified using several model assessment metrics and visualized through target diagrams. These metrics indicate that the model accurately estimated water levels, flows, and salinity over wide-ranging tidal and fluvial conditions, and the model can be used to investigate detailed circulation and salinity patterns throughout the Bay-Delta. The hydrodynamics produced through this effort will be used to drive affiliated sediment, phytoplankton, and contaminant hindcast efforts and habitat suitability assessments for fish and bivalves. The modeling framework applied here will serve as a baseline to ultimately shed light on potential ecosystem change over the current century.

Application of an unstructured 3D finite volume numerical model to flows and salinity dynamics in the San Francisco Bay-Delta

Science.gov (United States)

Martyr-Koller, R. C.; Kernkamp, H. W. J.; van Dam, A.; van der Wegen, M.; Lucas, L. V.; Knowles, N.; Jaffe, B.; Fregoso, T. A.

2017-06-01

A linked modeling approach has been undertaken to understand the impacts of climate and infrastructure on aquatic ecology and water quality in the San Francisco Bay-Delta region. The Delft3D Flexible Mesh modeling suite is used in this effort for its 3D hydrodynamics, salinity, temperature and sediment dynamics, phytoplankton and water-quality coupling infrastructure, and linkage to a habitat suitability model. The hydrodynamic model component of the suite is D-Flow FM, a new 3D unstructured finite-volume model based on the Delft3D model. In this paper, D-Flow FM is applied to the San Francisco Bay-Delta to investigate tidal, seasonal and annual dynamics of water levels, river flows and salinity under historical environmental and infrastructural conditions. The model is driven by historical winds, tides, ocean salinity, and river flows, and includes federal, state, and local freshwater withdrawals, and regional gate and barrier operations. The model is calibrated over a 9-month period, and subsequently validated for water levels, flows, and 3D salinity dynamics over a 2 year period. Model performance was quantified using several model assessment metrics and visualized through target diagrams. These metrics indicate that the model accurately estimated water levels, flows, and salinity over wide-ranging tidal and fluvial conditions, and the model can be used to investigate detailed circulation and salinity patterns throughout the Bay-Delta. The hydrodynamics produced through this effort will be used to drive affiliated sediment, phytoplankton, and contaminant hindcast efforts and habitat suitability assessments for fish and bivalves. The modeling framework applied here will serve as a baseline to ultimately shed light on potential ecosystem change over the current century.

Screening masses in quenched (2+1)d Yang-Mills theory: Universality from dynamics?

International Nuclear Information System (INIS)

Frigori, Rafael B.

2010-01-01

We compute the spectrum of gluonic screening-masses in the 0 ++ channel of quenched 3d Yang-Mills theory near the phase-transition. Our finite-temperature lattice simulations are performed at scaling region, using state-of-art techniques for thermalization and spectroscopy, which allows for thorough data extrapolations to thermodynamic limit. Ratios among mass-excitations with the same quantum numbers on the gauge theory, 2d Ising and λφ 4 models are compared, resulting in a nice agreement with predictions from universality. In addition, a gauge-to-scalar mapping, previously employed to fit QCD Green's functions at deep IR, is verified to dynamically describe these universal spectroscopic patterns.

Modelling atmospheric transport of persistent organic pollutants in the Northern Hemisphere with a 3-D dynamical model: DEHM-POP

OpenAIRE

Hansen , K. M.; Christensen , J. H.; Brandt , J.; Frohn , L. M.; Geels , C.

2004-01-01

International audience; The Danish Eulerian Hemispheric Model (DEHM) is a 3-D dynamical atmospheric transport model originally developed to describe the atmospheric transport of sulphur into the Arctic. A new version of the model, DEHM-POP, developed to study the atmospheric transport and environmental fate of persistent organic pollutants (POPs) is presented. During environmental cycling, POPs can be deposited and re-emitted several times before reaching a final destination. A description of...

Floodplain simulation for Musi River using integrated 1D/2D hydrodynamic model

Directory of Open Access Journals (Sweden)

Al Amin Muhammad B.

2017-01-01

Full Text Available This paper presents the simulation of floodplain at Musi River using integrated 1D and 2D hydrodynamic model. The 1D flow simulation was applied for the river channel with flow hydrograph as upstream boundary condition. The result of 1D flow simulation was integrated into 2D flow simulation in order to know the area and characteristics of flood inundation. The input data of digital terrain model which was used in this research had grid resolution of 10m×10m, but for 2D simulation the resolution was with grid resolution 50 m × 50 m so as to limit simulation time since the model size was big enough. The result of the simulation showed that the inundated area surrounding Musi River is about 107.44 km2 with maximum flood depth is 3.24 m, water surface velocity ranges from 0.00 to 0.83 m/s. Most of floodplain areas varied from middle to high flood hazard level, and only few areas had very high level of flood hazard especially on river side. The structural flood control measurement to be recommended to Palembang is to construct flood dike and flood gate. The non structural measurement one is to improve watershed management and socialization of flood awareness.

A 2-D process-based model for suspended sediment dynamics: a first step towards ecological modeling

Science.gov (United States)

Achete, F. M.; van der Wegen, M.; Roelvink, D.; Jaffe, B.

2015-06-01

In estuaries suspended sediment concentration (SSC) is one of the most important contributors to turbidity, which influences habitat conditions and ecological functions of the system. Sediment dynamics differs depending on sediment supply and hydrodynamic forcing conditions that vary over space and over time. A robust sediment transport model is a first step in developing a chain of models enabling simulations of contaminants, phytoplankton and habitat conditions. This works aims to determine turbidity levels in the complex-geometry delta of the San Francisco estuary using a process-based approach (Delft3D Flexible Mesh software). Our approach includes a detailed calibration against measured SSC levels, a sensitivity analysis on model parameters and the determination of a yearly sediment budget as well as an assessment of model results in terms of turbidity levels for a single year, water year (WY) 2011. Model results show that our process-based approach is a valuable tool in assessing sediment dynamics and their related ecological parameters over a range of spatial and temporal scales. The model may act as the base model for a chain of ecological models assessing the impact of climate change and management scenarios. Here we present a modeling approach that, with limited data, produces reliable predictions and can be useful for estuaries without a large amount of processes data.

A 2-D process-based model for suspended sediment dynamics: A first step towards ecological modeling

Science.gov (United States)

Achete, F. M.; van der Wegen, M.; Roelvink, D.; Jaffe, B.

2015-01-01

In estuaries suspended sediment concentration (SSC) is one of the most important contributors to turbidity, which influences habitat conditions and ecological functions of the system. Sediment dynamics differs depending on sediment supply and hydrodynamic forcing conditions that vary over space and over time. A robust sediment transport model is a first step in developing a chain of models enabling simulations of contaminants, phytoplankton and habitat conditions. This works aims to determine turbidity levels in the complex-geometry delta of the San Francisco estuary using a process-based approach (Delft3D Flexible Mesh software). Our approach includes a detailed calibration against measured SSC levels, a sensitivity analysis on model parameters and the determination of a yearly sediment budget as well as an assessment of model results in terms of turbidity levels for a single year, water year (WY) 2011. Model results show that our process-based approach is a valuable tool in assessing sediment dynamics and their related ecological parameters over a range of spatial and temporal scales. The model may act as the base model for a chain of ecological models assessing the impact of climate change and management scenarios. Here we present a modeling approach that, with limited data, produces reliable predictions and can be useful for estuaries without a large amount of processes data.

[Dynamic study of the female levator ani muscle using MRI 3D vectorial modeling].

Science.gov (United States)

Delmas, Vincent; Ami, Olivier; Iba-Zizen, Marie-Thérèse

2010-06-01

The levator ani muscle has a major role in the female pelvic floor, and is involved in the pathophysiology of pelvic prolapse and stress urinary incontinence. We conducted an anatomical and morphological study of this muscle using dynamic 3D vectorial reconstruction MRI, in order to analyze the contraction of two major components of the levator ani: the iliococcygeus and pubococcygeus. Three volunteer healthy continent nulliparous women aged from 19 to 22 underwent dynamic pelvic MRI. Coronal T2-weighted pelvic images were obtained in the supine position, at rest, holding back, and during Valsalva stress effort. 3D vectorial models were reconstructed by manual segmentation of the source images, and were set up on bony anatomic marks. Iliococcygeus and pubococcygeus volumes were measured in the three positions. Volumetrics, displacement and dynamic morphing changes were analyzed with 3D vectorial animation software. The urogenital hiatus extended more holding back (mean +4.31 mm) than on effort (mean +2.78 mm). The iliococcygeus lowered (mean -3.95 mm) and deviated outward (mean +3.01 mm). The basic tone of the iliococcygeus muscle gives it a dome shape, and its reflex contraction against abdominal strain ensures anal and urinary continence The levator ani is more than a pelvic diaphragm: it is a truly dynamic pelvic floor. Its points of support on the stiff osseous frame allow it to retain the pelvic organs. The levator ani muscle seems to prevent anal prolapse during stress strain.

Energy transfer and reaction dynamics of matrix-isolated 1,2-difluoroethane-d4

Science.gov (United States)

Raff, Lionel M.

1990-09-01

The molecular dynamics of vibrationally excited 1,2-difluoroethane-d4 isolated in Ar, Kr, and Xe matrices at 12 K are investigated using trajectory methods. The matrix model is an fcc crystal containing 125 unit cells with 666 atoms in a cubic (5×5×5) arrangement. It is assumed that 1,2-difluoroethane-d4 is held interstitially within the volume bounded by the innermost unit cell of the crystal. The transport effects of the bulk are simulated using the velocity reset method introduced by Riley, Coltrin, and Diestler [J. Chem. Phys. 88, 5934 (1988)]. The system potential is written as the separable sum of a lattice potential, a lattice-molecule interaction and a gas-phase potential for 1,2-difluoroethane. The first two of these are assumed to have pairwise form while the molecular potential is a modified form of the global potential previously developed for 1,2-difluoroethane [J. Phys. Chem. 91, 3266 (1987)]. Calculated sublimation energies for the pure crystals are in good accord with the experimental data. The distribution of metastable-state energies for matrix-isolated 1,2-difluoroethane-d4 is Gaussian in form. In krypton, the full width at half maximum for the distribution is 0.37 eV. For a total excitation energy of 6.314 eV, the observed dynamic processes are vibrational relaxation, orientational exchange, and four-center DF elimination reactions. The first of these processes is characterized by a near linear, first-order decay curve with rate coefficients in the range 1.30-1.48×1011 s-1. The average rates in krypton and xenon are nearly equal. The process is slightly slower in argon. The decay curves exhibit characteristic high-frequency oscillations that are generally seen in energy transfer studies. It is demonstrated that these oscillations are associated with the frequencies for intramolecular energy transfer so that the entire frequency spectrum for such transfer processes can be obtained from the Fourier transform of the decay curve. Orientational

Coupled 0D-1D CFD Modeling of Right Heart and Pulmonary Artery Morphometry Tree

Science.gov (United States)

Dong, Melody; Yang, Weiguang; Feinstein, Jeffrey A.; Marsden, Alison

2017-11-01

Pulmonary arterial hypertension (PAH) is characterized by elevated pulmonary artery (PA) pressure and remodeling of the distal PAs resulting in right ventricular (RV) dysfunction and failure. It is hypothesized that patients with untreated ventricular septal defects (VSD) may develop PAH due to elevated flows and pressures in the PAs. Wall shear stress (WSS), due to elevated flows, and circumferential stress, due to elevated pressures, are known to play a role in vascular mechanobiology. Thus, simulating VSD hemodynamics and wall mechanics may facilitate our understanding of mechanical stimuli leading to PAH initiation and progression. Although 3D CFD models can capture detailed hemodynamics in the proximal PAs, they cannot easily model hemodynamics and wave propagation in the distal PAs, where remodeling occurs. To improve current PA models, we will present a new method that couples distal PA hemodynamics with RV function. Our model couples a 0D lumped parameter model of the RV to a 1D model of the PA tree, based on human PA morphometry data, to characterize RV performance and WSS changes in the PA tree. We will compare a VSD 0D-1D model and a 0D-3D model coupled to a mathematical morphometry tree model to quantify WSS in the entire PA vascular tree.

Model-based dynamic resistive wall mode identification and feedback control in the DIII-D tokamak

International Nuclear Information System (INIS)

In, Y.; Kim, J.S.; Edgell, D.H.; Strait, E.J.; Humphreys, D.A.; Walker, M.L.; Jackson, G.L.; Chu, M.S.; Johnson, R.; La Haye, R.J.; Okabayashi, M.; Garofalo, A.M.; Reimerdes, H.

2006-01-01

A new model-based dynamic resistive wall mode (RWM) identification and feedback control algorithm has been developed. While the overall RWM structure can be detected by a model-based matched filter in a similar manner to a conventional sensor-based scheme, it is significantly influenced by edge-localized-modes (ELMs). A recent study suggested that such ELM noise might cause the RWM control system to respond in an undesirable way. Thus, an advanced algorithm to discriminate ELMs from RWM has been incorporated into this model-based control scheme, dynamic Kalman filter. Specifically, the DIII-D [J. L. Luxon, Nucl. Fusion 42, 614 (2002)] resistive vessel wall was modeled in two ways: picture frame model or eigenmode treatment. Based on the picture frame model, the first real-time, closed-loop test results of the Kalman filter algorithms during DIII-D experimental operation are presented. The Kalman filtering scheme was experimentally confirmed to be effective in discriminating ELMs from RWM. As a result, the actuator coils (I-coils) were rarely excited during ELMs, while retaining the sensitivity to RWM. However, finding an optimized set of operating parameters for the control algorithm requires further analysis and design. Meanwhile, a more advanced Kalman filter based on a more accurate eigenmode model has been developed. According to this eigenmode approach, significant improvement in terms of control performance has been predicted, while maintaining good ELM discrimination

Comparison of 2D and 3D modeled tumor motion estimation/prediction for dynamic tumor tracking during arc radiotherapy

Science.gov (United States)

Liu, Wu; Ma, Xiangyu; Yan, Huagang; Chen, Zhe; Nath, Ravinder; Li, Haiyun

2017-05-01

Many real-time imaging techniques have been developed to localize a target in 3D space or in a 2D beam’s eye view (BEV) plane for intrafraction motion tracking in radiation therapy. With tracking system latency, the 3D-modeled method is expected to be more accurate even in terms of 2D BEV tracking error. No quantitative analysis, however, has been reported. In this study, we simulated co-planar arc deliveries using respiratory motion data acquired from 42 patients to quantitatively compare the accuracy between 2D BEV and 3D-modeled tracking in arc therapy and to determine whether 3D information is needed for motion tracking. We used our previously developed low kV dose adaptive MV-kV imaging and motion compensation framework as a representative of 3D-modeled methods. It optimizes the balance between additional kV imaging dose and 3D tracking accuracy and solves the MLC blockage issue. With simulated Gaussian marker detection errors (zero mean and 0.39 mm standard deviation) and ~155/310/460 ms tracking system latencies, the mean percentage of time that the target moved  >2 mm from the predicted 2D BEV position are 1.1%/4.0%/7.8% and 1.3%/5.8%/11.6% for the 3D-modeled and 2D-only tracking, respectively. The corresponding average BEV RMS errors are 0.67/0.90/1.13 mm and 0.79/1.10/1.37 mm. Compared to the 2D method, the 3D method reduced the average RMS unresolved motion along the beam direction from ~3 mm to ~1 mm, resulting in on average only  <1% dosimetric advantage in the depth direction. Only for a small fraction of the patients, when tracking latency is long, the 3D-modeled method showed significant improvement of BEV tracking accuracy, indicating potential dosimetric advantage. However, if the tracking latency is short (~150 ms or less), those improvements are limited. Therefore, 2D BEV tracking has sufficient targeting accuracy for most clinical cases. The 3D technique is, however, still important in solving the MLC blockage problem

Simplification of one-dimensional hydraulic networks by automated processes evaluated on 1D/2D deterministic flood models

DEFF Research Database (Denmark)

Davidsen, Steffen; Löwe, Roland; Thrysøe, Cecilie

2017-01-01

Evaluation of pluvial flood risk is often based on computations using 1D/2D urban flood models. However, guidelines on choice of model complexity are missing, especially for one-dimensional (1D) network models. This study presents a new automatic approach for simplification of 1D hydraulic networ...

REVIEW One-Dimensional Dynamical Modeling of Earthquakes: A Review

Directory of Open Access Journals (Sweden)

Jeen-Hwa Wang

2008-01-01

Full Text Available Studies of the power-law relations of seismicity and earthquake source parameters based on the one-dimensional (1-D Burridge-Knopoff¡¦s (BK dynamical lattice model, especially those studies conducted by Taiwan¡¦s scientists, are reviewed in this article. In general, velocity- and/or state-dependent friction is considered to control faulting. A uniform distribution of breaking strengths (i.e., the static friction strength is taken into account in some studies, and inhomogeneous distributions in others. The scaling relations in these studies include: Omori¡¦s law, the magnitude-frequency or energy-frequency relation, the relation between source duration time and seismic moment, the relation between rupture length and seismic moment, the frequency-length relation, and the source power spectra. The main parameters of the one-dimensional (1-D Burridge-Knopoff¡¦s (BK dynamical lattice model include: the decreasing rate (r of dynamic friction strength with sliding velocity; the type and degree of heterogeneous distribution of the breaking strengths, the stiffness ratio (i.e., the ratio between the stiffness of the coil spring connecting two mass elements and that of the leaf spring linking a mass element and the moving plate; the frictional drop ratio of the minimum dynamic friction strength to the breaking strength; and the maximum breaking strength. For some authors, the distribution of the breaking strengths was considered to be a fractal function. Hence, the fractal dimension of such a distribution is also a significant parameter. Comparison between observed scaling laws and simulation results shows that the 1-D BK dynamical lattice model acceptably approaches fault dynamics.

Activated sludge model No. 2d, ASM2d

DEFF Research Database (Denmark)

Henze, M.

1999-01-01

The Activated Sludge Model No. 2d (ASM2d) presents a model for biological phosphorus removal with simultaneous nitrification-denitrification in activated sludge systems. ASM2d is based on ASM2 and is expanded to include the denitrifying activity of the phosphorus accumulating organisms (PAOs......). This extension of ASM2 allows for improved modeling of the processes, especially with respect to the dynamics of nitrate and phosphate. (C) 1999 IAWQ Published by Elsevier Science Ltd. All rights reserved....

Improved Large-Scale Inundation Modelling by 1D-2D Coupling and Consideration of Hydrologic and Hydrodynamic Processes - a Case Study in the Amazon

Science.gov (United States)

Hoch, J. M.; Bierkens, M. F.; Van Beek, R.; Winsemius, H.; Haag, A.

2015-12-01

Understanding the dynamics of fluvial floods is paramount to accurate flood hazard and risk modeling. Currently, economic losses due to flooding constitute about one third of all damage resulting from natural hazards. Given future projections of climate change, the anticipated increase in the World's population and the associated implications, sound knowledge of flood hazard and related risk is crucial. Fluvial floods are cross-border phenomena that need to be addressed accordingly. Yet, only few studies model floods at the large-scale which is preferable to tiling the output of small-scale models. Most models cannot realistically model flood wave propagation due to a lack of either detailed channel and floodplain geometry or the absence of hydrologic processes. This study aims to develop a large-scale modeling tool that accounts for both hydrologic and hydrodynamic processes, to find and understand possible sources of errors and improvements and to assess how the added hydrodynamics affect flood wave propagation. Flood wave propagation is simulated by DELFT3D-FM (FM), a hydrodynamic model using a flexible mesh to schematize the study area. It is coupled to PCR-GLOBWB (PCR), a macro-scale hydrological model, that has its own simpler 1D routing scheme (DynRout) which has already been used for global inundation modeling and flood risk assessments (GLOFRIS; Winsemius et al., 2013). A number of model set-ups are compared and benchmarked for the simulation period 1986-1996: (0) PCR with DynRout; (1) using a FM 2D flexible mesh forced with PCR output and (2) as in (1) but discriminating between 1D channels and 2D floodplains, and, for comparison, (3) and (4) the same set-ups as (1) and (2) but forced with observed GRDC discharge values. Outputs are subsequently validated against observed GRDC data at Óbidos and flood extent maps from the Dartmouth Flood Observatory. The present research constitutes a first step into a globally applicable approach to fully couple

Coupling of the computational fluid dynamics code ANSYS CFX with the 3D neutron kinetic core model DYN3D

International Nuclear Information System (INIS)

Kliem, S.; Grahn, A.; Rohde, U.; Schuetze, J.; Frank, Th.

2010-01-01

The computational fluid dynamics code ANSYS CFX has been coupled with the neutron-kinetic core model DYN3D. ANSYS CFX calculates the fluid dynamics and related transport phenomena in the reactors coolant and provides the corresponding data to DYN3D. In the fluid flow simulation of the coolant, the core itself is modeled within the porous body approach. DYN3D calculates the neutron kinetics and the fuel behavior including the heat transfer to the coolant. The physical data interface between the codes is the volumetric heat release rate into the coolant. In the prototype that is currently available, the coupling is restricted to single-phase flow problems. In the time domain an explicit coupling of the codes has been implemented so far. Steady-state and transient verification calculations for two small-size test problems confirm the correctness of the implementation of the prototype coupling. The first test problem was a mini-core consisting of nine real-size fuel assemblies with quadratic cross section. Comparison was performed with the DYN3D stand-alone code. In the steady state, the effective multiplication factor obtained by the DYN3D/ANSYS CFX codes hows a deviation of 9.8 pcm from the DYN3D stand-alone solution. This difference can be attributed to the use of different water property packages in the two codes. The transient test case simulated the withdrawal of the control rod from the central fuel assembly at hot zero power in the same mini-core. Power increase during the introduction of positive reactivity and power reduction due to fuel temperature increase are calculated in the same manner by the coupled and the stand-alone codes. The maximum values reached during the power rise differ by about 1 MW at a power level of 50 MW. Beside the different water property packages, these differences are caused by the use of different flow solvers. The same calculations were carried for a mini-core with seven real-size fuel assemblies with hexagonal cross section in

Flood Inundation Modelling in the Kuantan River Basin using 1D-2D Flood Modeller coupled with ASTER-GDEM

Science.gov (United States)

Ng, Z. F.; Gisen, J. I.; Akbari, A.

2018-03-01

Topography dataset is an important input in performing flood inundation modelling. However, it is always difficult to obtain high resolution topography that provide accurate elevation information. Fortunately, there are some open source topography datasets available with reasonable resolution such as SRTM and ASTER-GDEM. In Malaysia particularly in Kuantan, the modelling research on the floodplain area is still lacking. This research aims to: a) to investigate the suitability of ASTER-GDEM to be applied in the 1D-2D flood inundation modelling for the Kuantan River Basin; b) to generate flood inundation map for Kuantan river basin. The topography dataset used in this study is ASTER-GDEM to generate physical characteristics of watershed in the basin. It is used to perform rainfall runoff modelling for hydrological studies and to delineate flood inundation area in the Flood Modeller. The results obtained have shown that a 30m resolution ASTER-GDEM is applicable as an input for the 1D-2D flood modelling. The simulated water level in 2013 has NSE of 0.644 and RSME of 1.259. As a conclusion, ASTER-GDEM can be used as one alternative topography datasets for flood inundation modelling. However, the flood level obtained from the hydraulic modelling shows low accuracy at flat urban areas.

Two-phase 1D+1D model of a DMFC: development and validation on extensive operating conditions range

Energy Technology Data Exchange (ETDEWEB)

Casalegno, A.; Marchesi, R.; Parenti, D. [Dipartimento di Energetica, Politecnico di Milano (Italy)

2008-02-15

A two-phase 1D+1D model of a direct methanol fuel cell (DMFC) is developed, considering overall mass balance, methanol transport in gas phase through anode diffusion layer, methanol and water crossover. The model is quantitatively validated on an extensive range of operating conditions, 24 polarisation curves. The model accurately reproduces DMFC performance in the validation range and, outside this, it is able to predict values under feasible operating conditions. Finally, the estimations of methanol crossover flux are qualitatively and quantitatively similar to experimental measures and the main local quantities' trends are coherent with results obtained with more complex models. (Abstract Copyright [2008], Wiley Periodicals, Inc.)

A 3D dynamical model of the colliding winds in binary systems

Science.gov (United States)

Parkin, E. R.; Pittard, J. M.

2008-08-01

We present a three-dimensional (3D) dynamical model of the orbital-induced curvature of the wind-wind collision region in binary star systems. Momentum balance equations are used to determine the position and shape of the contact discontinuity between the stars, while further downstream the gas is assumed to behave ballistically. An Archimedean spiral structure is formed by the motion of the stars, with clear resemblance to high-resolution images of the so-called `pinwheel nebulae'. A key advantage of this approach over grid or smoothed particle hydrodynamic models is its significantly reduced computational cost, while it also allows the study of the structure obtained in an eccentric orbit. The model is relevant to symbiotic systems and γ-ray binaries, as well as systems with O-type and Wolf-Rayet stars. As an example application, we simulate the X-ray emission from hypothetical O+O and WR+O star binaries, and describe a method of ray tracing through the 3D spiral structure to account for absorption by the circumstellar material in the system. Such calculations may be easily adapted to study observations at wavelengths ranging from the radio to γ-ray.

Livrable D1.2 of the PERSEE project : Perceptual-Modelling-Definition-of-the-Models

OpenAIRE

Wang , Junle; Bosc , Emilie; Li , Jing; Ricordel , Vincent

2011-01-01

Livrable D1.2 du projet ANR PERSEE; Ce rapport a été réalisé dans le cadre du projet ANR PERSEE (n° ANR-09-BLAN-0170). Exactement il correspond au livrable D1.2 du projet. Son titre : Perceptual-Modelling-Definition-of-the-Models

How Spatial Abilities and Dynamic Visualizations Interplay When Learning Functional Anatomy with 3D Anatomical Models

Science.gov (United States)

Berney, Sandra; Bétrancourt, Mireille; Molinari, Gaëlle; Hoyek, Nady

2015-01-01

The emergence of dynamic visualizations of three-dimensional (3D) models in anatomy curricula may be an adequate solution for spatial difficulties encountered with traditional static learning, as they provide direct visualization of change throughout the viewpoints. However, little research has explored the interplay between learning material…

Black hole dynamics at large D

CERN Multimedia

CERN. Geneva

2016-01-01

We demonstrate that the classical dynamics of black holes can be reformulated as a dynamical problem of a codimension one membrane moving in flat space. This membrane - roughly the black hole event horizon - carries a conserved charge current and stress tensor which source radiation. This `membrane paradigm' may be viewed as a simplification of the equations of general relativity at large D, and suggests the possibility of using 1/D as a useful expansion parameter in the analysis of complicated four dimensional solutions of general relativity, for instance the collision between two black holes.

Motor function is associated with 1,25(OH)(2)D and indices of insulin-glucose dynamics in non-diabetic older adults.

Science.gov (United States)

Justice, Jamie N; Pierpoint, Lauren A; Mani, Diba; Schwartz, Robert S; Enoka, Roger M

2014-06-01

Advancing age is accompanied by changes in metabolic characteristics, such as reduced insulin sensitivity and low levels of vitamin D, which may exacerbate age-related declines in physical function. The aim of the present study was to determine the associations between insulin-glucose dynamics, vitamin D metabolites, and performance on a battery of motor tasks in healthy, non-diabetic older adults. Sixty-nine community-dwelling men and women (65-90 years) were recruited. Insulin-glucose dynamics were determined by an intravenous glucose tolerance test, and vitamin D metabolites were measured. Motor function was characterized by the time to walk 500 m, chair-rise time, lower body strength, dorsiflexor steadiness and endurance time, and muscle coactivation. Significant unadjusted correlations were found between insulin-glucose dynamics and 1,25-dihydroxyvitamin D [1,25(OH)2D] with walk time, strength, steadiness, endurance time, and muscle activation (p glucose (R (2) = 0.55, p glucose dynamics and the bio-active vitamin D metabolite 1,25(OH)2D. Walking endurance and strength were explained by 1,25(OH)2D and fasting blood glucose and insulin, even after adjusting for age, sex, and body fat. Motor function in a relatively small sample of non-diabetic older men and women was associated with metabolic factors that increase in prevalence with aging.

1D and 2D urban dam-break flood modelling in Istanbul, Turkey

Science.gov (United States)

Ozdemir, Hasan; Neal, Jeffrey; Bates, Paul; Döker, Fatih

2014-05-01

Urban flood events are increasing in frequency and severity as a consequence of several factors such as reduced infiltration capacities due to continued watershed development, increased construction in flood prone areas due to population growth, the possible amplification of rainfall intensity due to climate change, sea level rise which threatens coastal development, and poorly engineered flood control infrastructure (Gallegos et al., 2009). These factors will contribute to increased urban flood risk in the future, and as a result improved modelling of urban flooding according to different causative factor has been identified as a research priority (Gallegos et al., 2009; Ozdemir et al. 2013). The flooding disaster caused by dam failures is always a threat against lives and properties especially in urban environments. Therefore, the prediction of dynamics of dam-break flows plays a vital role in the forecast and evaluation of flooding disasters, and is of long-standing interest for researchers. Flooding occurred on the Ayamama River (Istanbul-Turkey) due to high intensity rainfall and dam-breaching of Ata Pond in 9th September 2009. The settlements, industrial areas and transportation system on the floodplain of the Ayamama River were inundated. Therefore, 32 people were dead and millions of Euros economic loses were occurred. The aim of this study is 1 and 2-Dimensional flood modelling of the Ata Pond breaching using HEC-RAS and LISFLOOD-Roe models and comparison of the model results using the real flood extent. The HEC-RAS model solves the full 1-D Saint Venant equations for unsteady open channel flow whereas LISFLOOD-Roe is the 2-D shallow water model which calculates the flow according to the complete Saint Venant formulation (Villanueva and Wright, 2006; Neal et al., 2011). The model consists a shock capturing Godunov-type scheme based on the Roe Riemann solver (Roe, 1981). 3 m high resolution Digital Surface Model (DSM), natural characteristics of the pond

The dynamical core of the Aeolus 1.0 statistical-dynamical atmosphere model: validation and parameter optimization

Science.gov (United States)

Totz, Sonja; Eliseev, Alexey V.; Petri, Stefan; Flechsig, Michael; Caesar, Levke; Petoukhov, Vladimir; Coumou, Dim

2018-02-01

We present and validate a set of equations for representing the atmosphere's large-scale general circulation in an Earth system model of intermediate complexity (EMIC). These dynamical equations have been implemented in Aeolus 1.0, which is a statistical-dynamical atmosphere model (SDAM) and includes radiative transfer and cloud modules (Coumou et al., 2011; Eliseev et al., 2013). The statistical dynamical approach is computationally efficient and thus enables us to perform climate simulations at multimillennia timescales, which is a prime aim of our model development. Further, this computational efficiency enables us to scan large and high-dimensional parameter space to tune the model parameters, e.g., for sensitivity studies.Here, we present novel equations for the large-scale zonal-mean wind as well as those for planetary waves. Together with synoptic parameterization (as presented by Coumou et al., 2011), these form the mathematical description of the dynamical core of Aeolus 1.0.We optimize the dynamical core parameter values by tuning all relevant dynamical fields to ERA-Interim reanalysis data (1983-2009) forcing the dynamical core with prescribed surface temperature, surface humidity and cumulus cloud fraction. We test the model's performance in reproducing the seasonal cycle and the influence of the El Niño-Southern Oscillation (ENSO). We use a simulated annealing optimization algorithm, which approximates the global minimum of a high-dimensional function.With non-tuned parameter values, the model performs reasonably in terms of its representation of zonal-mean circulation, planetary waves and storm tracks. The simulated annealing optimization improves in particular the model's representation of the Northern Hemisphere jet stream and storm tracks as well as the Hadley circulation.The regions of high azonal wind velocities (planetary waves) are accurately captured for all validation experiments. The zonal-mean zonal wind and the integrated lower

Self-Consistent 3D Modeling of Electron Cloud Dynamics and Beam Response

International Nuclear Information System (INIS)

Furman, Miguel; Furman, M.A.; Celata, C.M.; Kireeff-Covo, M.; Sonnad, K.G.; Vay, J.-L.; Venturini, M.; Cohen, R.; Friedman, A.; Grote, D.; Molvik, A.; Stoltz, P.

2007-01-01

We present recent advances in the modeling of beam electron-cloud dynamics, including surface effects such as secondary electron emission, gas desorption, etc, and volumetric effects such as ionization of residual gas and charge-exchange reactions. Simulations for the HCX facility with the code WARP/POSINST will be described and their validity demonstrated by benchmarks against measurements. The code models a wide range of physical processes and uses a number of novel techniques, including a large-timestep electron mover that smoothly interpolates between direct orbit calculation and guiding-center drift equations, and a new computational technique, based on a Lorentz transformation to a moving frame, that allows the cost of a fully 3D simulation to be reduced to that of a quasi-static approximation

3D Protein Dynamics in the Cell Nucleus.

Science.gov (United States)

Singh, Anand P; Galland, Rémi; Finch-Edmondson, Megan L; Grenci, Gianluca; Sibarita, Jean-Baptiste; Studer, Vincent; Viasnoff, Virgile; Saunders, Timothy E

2017-01-10

The three-dimensional (3D) architecture of the cell nucleus plays an important role in protein dynamics and in regulating gene expression. However, protein dynamics within the 3D nucleus are poorly understood. Here, we present, to our knowledge, a novel combination of 1) single-objective based light-sheet microscopy, 2) photoconvertible proteins, and 3) fluorescence correlation microscopy, to quantitatively measure 3D protein dynamics in the nucleus. We are able to acquire >3400 autocorrelation functions at multiple spatial positions within a nucleus, without significant photobleaching, allowing us to make reliable estimates of diffusion dynamics. Using this tool, we demonstrate spatial heterogeneity in Polymerase II dynamics in live U2OS cells. Further, we provide detailed measurements of human-Yes-associated protein diffusion dynamics in a human gastric cancer epithelial cell line. Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.

Ckmeans.1d.dp: Optimal k-means Clustering in One Dimension by Dynamic Programming.

Science.gov (United States)

Wang, Haizhou; Song, Mingzhou

2011-12-01

The heuristic k -means algorithm, widely used for cluster analysis, does not guarantee optimality. We developed a dynamic programming algorithm for optimal one-dimensional clustering. The algorithm is implemented as an R package called Ckmeans.1d.dp . We demonstrate its advantage in optimality and runtime over the standard iterative k -means algorithm.

A classical approach in simple nuclear fusion reaction 1H2+1H3 using two-dimension granular molecular dynamics model

International Nuclear Information System (INIS)

Viridi, S.; Kurniadi, R.; Waris, A.; Perkasa, Y. S.

2012-01-01

Molecular dynamics in 2-D accompanied by granular model provides an opportunity to investigate binding between nuclei particles and its properties that arises during collision in a fusion reaction. A fully classical approach is used to observe the influence of initial angle of nucleus orientation to the product yielded by the reaction. As an example, a simplest fusion reaction between 1 H 2 and 1 H 3 is observed. Several products of the fusion reaction have been obtained, even the unreported ones, including temporary 2 He 4 nucleus.

Linking 1D coastal ocean modelling to environmental management: an ensemble approach

Science.gov (United States)

Mussap, Giulia; Zavatarelli, Marco; Pinardi, Nadia

2017-12-01

The use of a one-dimensional interdisciplinary numerical model of the coastal ocean as a tool contributing to the formulation of ecosystem-based management (EBM) is explored. The focus is on the definition of an experimental design based on ensemble simulations, integrating variability linked to scenarios (characterised by changes in the system forcing) and to the concurrent variation of selected, and poorly constrained, model parameters. The modelling system used was previously specifically designed for the use in "data-rich" areas, so that horizontal dynamics can be resolved by a diagnostic approach and external inputs can be parameterised by nudging schemes properly calibrated. Ensembles determined by changes in the simulated environmental (physical and biogeochemical) dynamics, under joint forcing and parameterisation variations, highlight the uncertainties associated to the application of specific scenarios that are relevant to EBM, providing an assessment of the reliability of the predicted changes. The work has been carried out by implementing the coupled modelling system BFM-POM1D in an area of Gulf of Trieste (northern Adriatic Sea), considered homogeneous from the point of view of hydrological properties, and forcing it by changing climatic (warming) and anthropogenic (reduction of the land-based nutrient input) pressure. Model parameters affected by considerable uncertainties (due to the lack of relevant observations) were varied jointly with the scenarios of change. The resulting large set of ensemble simulations provided a general estimation of the model uncertainties related to the joint variation of pressures and model parameters. The information of the model result variability aimed at conveying efficiently and comprehensibly the information on the uncertainties/reliability of the model results to non-technical EBM planners and stakeholders, in order to have the model-based information effectively contributing to EBM.

Detection of phasic dopamine by D1 and D2 striatal medium spiny neurons.

Science.gov (United States)

Yapo, Cedric; Nair, Anu G; Clement, Lorna; Castro, Liliana R; Hellgren Kotaleski, Jeanette; Vincent, Pierre

2017-12-15

Brief dopamine events are critical actors of reward-mediated learning in the striatum; the intracellular cAMP-protein kinase A (PKA) response of striatal medium spiny neurons to such events was studied dynamically using a combination of biosensor imaging in mouse brain slices and in silico simulations. Both D1 and D2 medium spiny neurons can sense brief dopamine transients in the sub-micromolar range. While dopamine transients profoundly change cAMP levels in both types of medium spiny neurons, the PKA-dependent phosphorylation level remains unaffected in D2 neurons. At the level of PKA-dependent phosphorylation, D2 unresponsiveness depends on protein phosphatase-1 (PP1) inhibition by DARPP-32. Simulations suggest that D2 medium spiny neurons could detect transient dips in dopamine level. The phasic release of dopamine in the striatum determines various aspects of reward and action selection, but the dynamics of the dopamine effect on intracellular signalling remains poorly understood. We used genetically encoded FRET biosensors in striatal brain slices to quantify the effect of transient dopamine on cAMP or PKA-dependent phosphorylation levels, and computational modelling to further explore the dynamics of this signalling pathway. Medium-sized spiny neurons (MSNs), which express either D 1 or D 2 dopamine receptors, responded to dopamine by an increase or a decrease in cAMP, respectively. Transient dopamine showed similar sub-micromolar efficacies on cAMP in both D1 and D2 MSNs, thus challenging the commonly accepted notion that dopamine efficacy is much higher on D 2 than on D 1 receptors. However, in D2 MSNs, the large decrease in cAMP level triggered by transient dopamine did not translate to a decrease in PKA-dependent phosphorylation level, owing to the efficient inhibition of protein phosphatase 1 by DARPP-32. Simulations further suggested that D2 MSNs can also operate in a 'tone-sensing' mode, allowing them to detect transient dips in basal dopamine

Dynamic magnetic hysteresis behavior and dynamic phase transition in the spin-1 Blume-Capel model

Energy Technology Data Exchange (ETDEWEB)

Deviren, Bayram, E-mail: bayram.deviren@nevsehir.edu.tr [Department of Physics, Nevsehir University, 50300 Nevsehir (Turkey); Keskin, Mustafa [Department of Physics, Erciyes University, 38039 Kayseri (Turkey)

2012-03-15

The nature (time variation) of response magnetization m(wt) of the spin-1 Blume-Capel model in the presence of a periodically varying external magnetic field h(wt) is studied by employing the effective-field theory (EFT) with correlations as well as the Glauber-type stochastic dynamics. We determine the time variations of m(wt) and h(wt) for various temperatures, and investigate the dynamic magnetic hysteresis behavior. We also investigate the temperature dependence of the dynamic magnetization, hysteresis loop area and correlation near the transition point in order to characterize the nature (first- or second-order) of the dynamic transitions as well as obtain the dynamic phase transition temperatures. The hysteresis loops are obtained for different reduced temperatures and we find that the areas of the loops are decreasing with the increasing of the reduced temperatures. We also present the dynamic phase diagrams and compare the results of the EFT with the results of the dynamic mean-field approximation. The phase diagrams exhibit many dynamic critical points, such as tricritical ( Bullet ), zero-temperature critical (Z), triple (TP) and multicritical (A) points. According to values of Hamiltonian parameters, besides the paramagnetic (P), ferromagnetic (F) fundamental phases, one coexistence or mixed phase region, (F+P) and the reentrant behavior exist in the system. The results are in good agreement with some experimental and theoretical results. - Highlights: Black-Right-Pointing-Pointer Kinetic spin-1 Blume-Capel model is studied using the effective-field theory. Black-Right-Pointing-Pointer We investigated the dynamic magnetic hysteresis behavior. Black-Right-Pointing-Pointer Dynamic magnetization, hysteresis loop area and correlation are investigated. Black-Right-Pointing-Pointer System exhibits tricritical, zero-temperature, triple and multicritical points. Black-Right-Pointing-Pointer We present the dynamic phase diagrams and compare the results of the EFT

Mineral vein dynamics modeling (FRACS). Phase 1

Energy Technology Data Exchange (ETDEWEB)

Urai, J.; Virgo, S.; Arndt, M. [RWTH Aachen (Germany). Geologie-Endogene Dynamik] [and others

2013-07-15

The Mineral Vein Dynamics Modeling group ''FRACS'' is a team of 7 research groups from the Universities of Mainz, Aachen, Tuebingen, Karlsruhe, Bayreuth, ETH Zuerich and Glasgow working on an understanding of the dynamic development of fracturing, fluid flow and fracture sealing. World-class field laboratories, especially carbonate sequences from the Oman Mountains are studied and classified. State of the art numerical programs are written, expanded and used to simulate the dynamic interaction of fracturing, flow and resealing and the results are compared with the natural examples. Newest analytical technologies including laser scanning, high resolution X-ray microtomography, fluid inclusion and isotope analysis are performed to understand and compare the results of simulations with natural examples. A new statistical program was developed to classify the natural fracture and vein systems and compare them with dynamic numerical simulations and analytical models. The results of the first project phase are extremely promising. Most of the numerical models have been developed up to the stage where they can be used to simulate the natural examples. The models allow a definition of the first proxies for high fluid pressure and tectonic stresses. It was found out that the Oman Mountains are a complex and very dynamic system that constantly fractures and reseals from the scale of small veins up to the scale of large normal and strike slip faults. The numerical simulations also indicate that the permeability of such systems is not a constant but that the system adjusts to the driving force, for ex-ample high fluid pressure. When the system reseals fast a fluctuating behavior can be observed in the models where the system constantly fractures and reseals, which is in accordance with the observation of the natural laboratory.

Comparison of 1D and 2D CSR Models with Application to the FERMI(at)ELETTRA Bunch Compressors

International Nuclear Information System (INIS)

Bassi, G.; Ellison, J.A.; Heinemann, K.

2011-01-01

We compare our 2D mean field (Vlasov-Maxwell) treatment of coherent synchrotron radiation (CSR) effects with 1D approximations of the CSR force which are commonly implemented in CSR codes. In our model we track particles in 4D phase space and calculate 2D forces (1). The major cost in our calculation is the computation of the 2D force. To speed up the computation and improve 1D models we also investigate approximations to our exact 2D force. As an application, we present numerical results for the Fermi(at)Elettra first bunch compressor with the configuration described in (1).

HELIOS-CR - A 1-D radiation-magnetohydrodynamics code with inline atomic kinetics modeling

International Nuclear Information System (INIS)

MacFarlane, J.J.; Golovkin, I.E.; Woodruff, P.R.

2006-01-01

HELIOS-CR is a user-oriented 1D radiation-magnetohydrodynamics code to simulate the dynamic evolution of laser-produced plasmas and z-pinch plasmas. It includes an in-line collisional-radiative (CR) model for computing non-LTE atomic level populations at each time step of the hydrodynamics simulation. HELIOS-CR has been designed for ease of use, and is well-suited for experimentalists, as well as graduate and undergraduate student researchers. The energy equations employed include models for laser energy deposition, radiation from external sources, and high-current discharges. Radiative transport can be calculated using either a multi-frequency flux-limited diffusion model, or a multi-frequency, multi-angle short characteristics model. HELIOS-CR supports the use of SESAME equation of state (EOS) tables, PROPACEOS EOS/multi-group opacity data tables, and non-LTE plasma properties computed using the inline CR modeling. Time-, space-, and frequency-dependent results from HELIOS-CR calculations are readily displayed with the HydroPLOT graphics tool. In addition, the results of HELIOS simulations can be post-processed using the SPECT3D Imaging and Spectral Analysis Suite to generate images and spectra that can be directly compared with experimental measurements. The HELIOS-CR package runs on Windows, Linux, and Mac OSX platforms, and includes online documentation. We will discuss the major features of HELIOS-CR, and present example results from simulations

Critical bifurcation surfaces of 3D discrete dynamics

Directory of Open Access Journals (Sweden)

Michael Sonis

2000-01-01

Full Text Available This paper deals with the analytical representation of bifurcations of each 3D discrete dynamics depending on the set of bifurcation parameters. The procedure of bifurcation analysis proposed in this paper represents the 3D elaboration and specification of the general algorithm of the n-dimensional linear bifurcation analysis proposed by the author earlier. It is proven that 3D domain of asymptotic stability (attraction of the fixed point for a given 3D discrete dynamics is bounded by three critical bifurcation surfaces: the divergence, flip and flutter surfaces. The analytical construction of these surfaces is achieved with the help of classical Routh–Hurvitz conditions of asymptotic stability. As an application the adjustment process proposed by T. Puu for the Cournot oligopoly model is considered in detail.

Towards high resolution mapping of 3-D mesoscale dynamics from observations

Directory of Open Access Journals (Sweden)

B. Buongiorno Nardelli

2012-10-01

Full Text Available The MyOcean R&D project MESCLA (MEsoSCaLe dynamical Analysis through combined model, satellite and in situ data was devoted to the high resolution 3-D retrieval of tracer and velocity fields in the oceans, based on the combination of in situ and satellite observations and quasi-geostrophic dynamical models. The retrieval techniques were also tested and compared with the output of a primitive equation model, with particular attention to the accuracy of the vertical velocity field as estimated through the Q vector formulation of the omega equation. The project focused on a test case, covering the region where the Gulf Stream separates from the US East Coast. This work demonstrated that innovative methods for the high resolution mapping of 3-D mesoscale dynamics from observations can be used to build the next generations of operational observation-based products.

D2D-Enabled Small Cell Network Control Scheme Based on the Dynamic Stackelberg Game

Directory of Open Access Journals (Sweden)

Sungwook Kim

2017-01-01

Full Text Available For current and future cellular networks, small cell structure with licensed and unlicensed bandwidth, caching content provisioning, and device-to-device (D2D communications is seen as a necessary architecture. Recently, a series of control methods have been developed to address a myriad of challenges in next-generation small cell networks. In this study, we focus on the design of novel D2D-enabled small cell network control scheme by allowing caching and unlicensed D2D communications. Motivated by game theory and learning algorithm, the proposed scheme adaptively selects caching contents and splits the available bandwidth for licensed and unlicensed communications. Under dynamically changing network environments, we capture the dynamics of the network system and design a new dynamic Stackelberg game model. Based on a hierarchical and feedback based control manner, small base stations and users can be leaders or followers dynamically while improving 5G network performance. Simulations and performance analysis verify the efficiency of the proposed scheme, showing that our approach can outperform existing schemes by about 5%~15% in terms of bandwidth utilization, cache hit ratio, and system throughput.

Evaluation of the scale dependent dynamic SGS model in the open source code caffa3d.MBRi in wall-bounded flows

Science.gov (United States)

Draper, Martin; Usera, Gabriel

2015-04-01

The Scale Dependent Dynamic Model (SDDM) has been widely validated in large-eddy simulations using pseudo-spectral codes [1][2][3]. The scale dependency, particularly the potential law, has been proved also in a priori studies [4][5]. To the authors' knowledge there have been only few attempts to use the SDDM in finite difference (FD) and finite volume (FV) codes [6][7], finding some improvements with the dynamic procedures (scale independent or scale dependent approach), but not showing the behavior of the scale-dependence parameter when using the SDDM. The aim of the present paper is to evaluate the SDDM in the open source code caffa3d.MBRi, an updated version of the code presented in [8]. caffa3d.MBRi is a FV code, second-order accurate, parallelized with MPI, in which the domain is divided in unstructured blocks of structured grids. To accomplish this, 2 cases are considered: flow between flat plates and flow over a rough surface with the presence of a model wind turbine, taking for this case the experimental data presented in [9]. In both cases the standard Smagorinsky Model (SM), the Scale Independent Dynamic Model (SIDM) and the SDDM are tested. As presented in [6][7] slight improvements are obtained with the SDDM. Nevertheless, the behavior of the scale-dependence parameter supports the generalization of the dynamic procedure proposed in the SDDM, particularly taking into account that no explicit filter is used (the implicit filter is unknown). [1] F. Porté-Agel, C. Meneveau, M.B. Parlange. "A scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric boundary layer". Journal of Fluid Mechanics, 2000, 415, 261-284. [2] E. Bou-Zeid, C. Meneveau, M. Parlante. "A scale-dependent Lagrangian dynamic model for large eddy simulation of complex turbulent flows". Physics of Fluids, 2005, 17, 025105 (18p). [3] R. Stoll, F. Porté-Agel. "Dynamic subgrid-scale models for momentum and scalar fluxes in large-eddy simulations of

Relationships between Th1 or Th2 iNKT cell activity and structures of CD1d-antigen complexes: meta-analysis of CD1d-glycolipids dynamics simulations.

Directory of Open Access Journals (Sweden)

Xavier Laurent

2014-11-01

Full Text Available A number of potentially bioactive molecules can be found in nature. In particular, marine organisms are a valuable source of bioactive compounds. The activity of an α-galactosylceramide was first discovered in 1993 via screening of a Japanese marine sponge (Agelas mauritanius. Very rapidly, a synthetic glycololipid analogue of this natural molecule was discovered, called KRN7000. Associated with the CD1d protein, this α-galactosylceramide 1 (KRN7000 interacts with the T-cell antigen receptor to form a ternary complex that yields T helper (Th 1 and Th2 responses with opposing effects. In our work, we carried out molecular dynamics simulations (11.5 µs in total involving eight different ligands (conducted in triplicate in an effort to find out correlation at the molecular level, if any, between chemical modulation of 1 and the orientation of the known biological response, Th1 or Th2. Comparative investigations of human versus mouse and Th1 versus Th2 data have been carried out. A large set of analysis tools was employed including free energy landscapes. One major result is the identification of a specific conformational state of the sugar polar head, which could be correlated, in the present study, to the biological Th2 biased response. These theoretical tools provide a structural basis for predicting the very different dynamical behaviors of α-glycosphingolipids in CD1d and might aid in the future design of new analogues of 1.

1-D hybrid code for FRM dynamics

International Nuclear Information System (INIS)

Stark, R.A.; Miley, G.H.

1985-01-01

A 1-D radial hybrid code has been written to study the start-up of the FRM via neutral-beam injection. This code, named FROST (Field Reversed One-dimensional STart-up), models the plasma as azimuthal symmetric with no axial dependence. A multi-group method in energy and canonical angular momentum describes the large-orbit ions from the beam. This method is designed to be more efficient than those employing particle tracking, since the characteristic timescale of the simulation is the ion slowing down time, rather than the much shorter cyclotron period. A time-differentiated Grad-Shafranov equation couples the ion current to massless fluid equations describing electrons and low energy ions. Flux coordinates are used in this fluid model, in preference to an Eulerian framework, so that coupling of plasma at the two different radii of a closed flux surface may be treated with ease. Since a fluid treatment for electrons is invalid near a field null, a separate model for the electron current has been included for this region, a unique feature. Results of simulation of injection into a 2XIIB-like plasma are discussed. Electron currents are found to retard, but not prevent reversal of the magnetic field at the plasma center

A model with isospin doublet U(1)D gauge symmetry

Science.gov (United States)

Nomura, Takaaki; Okada, Hiroshi

2018-05-01

We propose a model with an extra isospin doublet U(1)D gauge symmetry, in which we introduce several extra fermions with odd parity under a discrete Z2 symmetry in order to cancel the gauge anomalies out. A remarkable issue is that we impose nonzero U(1)D charge to the Standard Model Higgs, and it gives the most stringent constraint to the vacuum expectation value of a scalar field breaking the U(1)D symmetry that is severer than the LEP bound. We then explore relic density of a Majorana dark matter candidate without conflict of constraints from lepton flavor violating processes. A global analysis is carried out to search for parameters which can accommodate with the observed data.

Control oriented 1D electrochemical model of lithium ion battery

International Nuclear Information System (INIS)

Smith, Kandler A.; Rahn, Christopher D.; Wang, Chao-Yang

2007-01-01

Lithium ion (Li-ion) batteries provide high energy and power density energy storage for diverse applications ranging from cell phones to hybrid electric vehicles (HEVs). For efficient and reliable systems integration, low order dynamic battery models are needed. This paper introduces a general method to generate numerically a fully observable/controllable state variable model from electrochemical kinetic, species and charge partial differential equations that govern the discharge/charge behavior of a Li-ion battery. Validated against a 313th order nonlinear CFD model of a 6 Ah HEV cell, a 12th order state variable model predicts terminal voltage to within 1% for pulse and constant current profiles at rates up to 50 C. The state equation is constructed in modal form with constant negative real eigenvalues distributed in frequency space from 0 to 10 Hz. Open circuit potential, electrode surface concentration/reaction distribution coupling and electrolyte concentration/ionic conductivity nonlinearities are explicitly approximated in the model output equation on a local, electrode-averaged and distributed basis, respectively. The balanced realization controllability/observability gramian indicates that the fast electrode surface concentration dynamics are more observable/controllable than the electrode bulk concentration dynamics (i.e. state of charge)

3D Digital Modelling

DEFF Research Database (Denmark)

Hundebøl, Jesper

wave of new building information modelling tools demands further investigation, not least because of industry representatives' somewhat coarse parlance: Now the word is spreading -3D digital modelling is nothing less than a revolution, a shift of paradigm, a new alphabet... Research qeustions. Based...... on empirical probes (interviews, observations, written inscriptions) within the Danish construction industry this paper explores the organizational and managerial dynamics of 3D Digital Modelling. The paper intends to - Illustrate how the network of (non-)human actors engaged in the promotion (and arrest) of 3...... important to appreciate the analysis. Before turning to the presentation of preliminary findings and a discussion of 3D digital modelling, it begins, however, with an outline of industry specific ICT strategic issues. Paper type. Multi-site field study...

Sphingomyelinase D activity in model membranes: structural effects of in situ generation of ceramide-1-phosphate

DEFF Research Database (Denmark)

Stock, Roberto; Brewer, Jonathan R.; Wagner, Kerstin

2012-01-01

The toxicity of Loxosceles spider venom has been attributed to a rare enzyme, sphingomyelinase D, which transforms sphingomyelin to ceramide-1-phosphate. The bases of its inflammatory and dermonecrotic activity, however, remain unclear. In this work the effects of ceramide-1-phosphate on model...... membranes were studied both by in situ generation of this lipid using a recombinant sphingomyelinase D from the spider Loxosceles laeta and by pre-mixing it with sphingomyelin and cholesterol. The systems of choice were large unilamellar vesicles for bulk studies (enzyme kinetics, fluorescence spectroscopy...... and dynamic light scattering) and giant unilamellar vesicles for fluorescence microscopy examination using a variety of fluorescent probes. The influence of membrane lateral structure on the kinetics of enzyme activity and the consequences of enzyme activity on the structure of target membranes containing...

Using an Active-Optical Sensor to Develop an Optimal NDVI Dynamic Model for High-Yield Rice Production (Yangtze, China).

Science.gov (United States)

Liu, Xiaojun; Ferguson, Richard B; Zheng, Hengbiao; Cao, Qiang; Tian, Yongchao; Cao, Weixing; Zhu, Yan

2017-03-24

The successful development of an optimal canopy vegetation index dynamic model for obtaining higher yield can offer a technical approach for real-time and nondestructive diagnosis of rice (Oryza sativa L) growth and nitrogen (N) nutrition status. In this study, multiple rice cultivars and N treatments of experimental plots were carried out to obtain: normalized difference vegetation index (NDVI), leaf area index (LAI), above-ground dry matter (DM), and grain yield (GY) data. The quantitative relationships between NDVI and these growth indices (e.g., LAI, DM and GY) were analyzed, showing positive correlations. Using the normalized modeling method, an appropriate NDVI simulation model of rice was established based on the normalized NDVI (RNDVI) and relative accumulative growing degree days (RAGDD). The NDVI dynamic model for high-yield production in rice can be expressed by a double logistic model: RNDVI = ( 1 + e - 15.2829 × ( R A G D D i - 0.1944 ) ) - 1 - ( 1 + e - 11.6517 × ( R A G D D i - 1.0267 ) ) - 1 (R2 = 0.8577**), which can be used to accurately predict canopy NDVI dynamic changes during the entire growth period. Considering variation among rice cultivars, we constructed two relative NDVI (RNDVI) dynamic models for Japonica and Indica rice types, with R2 reaching 0.8764** and 0.8874**, respectively. Furthermore, independent experimental data were used to validate the RNDVI dynamic models. The results showed that during the entire growth period, the accuracy (k), precision (R2), and standard deviation of RNDVI dynamic models for the Japonica and Indica cultivars were 0.9991, 1.0170; 0.9084**, 0.8030**; and 0.0232, 0.0170, respectively. These results indicated that RNDVI dynamic models could accurately reflect crop growth and predict dynamic changes in high-yield crop populations, providing a rapid approach for monitoring rice growth status.

Real-Time Dynamics in the (1+1)-D Abelian Higgs Model with Fermions

NARCIS (Netherlands)

Aarts, G.; Smit, J.

2000-01-01

In approximate dynamical equations, inhomogenous classical (mean) gauge and Higgs fields are coupled to quantized fermions. The equations are solved numerically on a spacetime lattice. The fermions appear to equilibrate according to the Fermi-Dirac distribution with time-dependent temperature and

CFTSIM-ITER dynamic fuel cycle model

International Nuclear Information System (INIS)

Busigin, A.; Gierszewski, P.

1998-01-01

Dynamic system models have been developed for specific tritium systems with considerable detail and for integrated fuel cycles with lesser detail (e.g. D. Holland, B. Merrill, Analysis of tritium migration and deposition in fusion reactor systems, Proceedings of the Ninth Symposium Eng. Problems of Fusion Research (1981); M.A. Abdou, E. Vold, C. Gung, M. Youssef, K. Shin, DT fuel self-sufficiency in fusion reactors, Fusion Technol. (1986); G. Spannagel, P. Gierszewski, Dynamic tritium inventory of a NET/ITER fuel cycle with lithium salt solution blanket, Fusion Eng. Des. (1991); W. Kuan, M.A. Abdou, R.S. Willms, Dynamic simulation of a proposed ITER tritium processing system, Fusion Technol. (1995)). In order to provide a tool to understand and optimize the behavior of the ITER fuel cycle, a dynamic fuel cycle model called CFTSIM is under development. The CFTSIM code incorporates more detailed ITER models, specifically for the important isotope separation system, and also has an easier-to-use graphical interface. This paper provides an overview of CFTSIM Version 1.0. The models included are those with significant and varying tritium inventories over a test campaign: fueling, plasma and first wall, pumping, fuel cleanup, isotope separation and storage. An illustration of the results is shown. (orig.)

Na-ion dynamics in Quasi-1D compound NaV2O4

International Nuclear Information System (INIS)

Månsson, M; Umegaki, I; Nozaki, H; Higuchi, Y; Sugiyama, J; Kawasaki, I; Watanabe, I; Sakurai, H

2014-01-01

We have used the pulsed muon source at ISIS to study high-temperature Na-ion dynamics in the quasi-one-dimensional (Q1D) metallic antiferromagnet NaV 2 O 4 . By performing systematic zero-field and longitudinal-field measurements as a function of temperature we clearly distinguish that the hopping rate increases exponentially above T diff ≈ 250 K. The data is well fitted to an Arrhenius type equation typical for a diffusion process, showing that the Na-ions starts to be mobile above T diff . Such results make this compound very interesting for the tuning of Q1D magnetism using atomic-scale ion-texturing through the periodic potential from ordered Na-vacancies. Further, it also opens the door to possible use of NaV 2 O 4 and related compounds in energy related applications

Computational fluid dynamics analysis of an innovative start-up method of high temperature fuel cells using dynamic 3d model

Directory of Open Access Journals (Sweden)

Kupecki Jakub

2017-03-01

Full Text Available The article presents a numerical analysis of an innovative method for starting systems based on high temperature fuel cells. The possibility of preheating the fuel cell stacks from the cold state to the nominal working conditions encounters several limitations related to heat transfer and stability of materials. The lack of rapid and safe start-up methods limits the proliferation of MCFCs and SOFCs. For that reason, an innovative method was developed and verified using the numerical analysis presented in the paper. A dynamic 3D model was developed that enables thermo-fluidic investigations and determination of measures for shortening the preheating time of the high temperature fuel cell stacks. The model was implemented in ANSYS Fluent computational fluid dynamic (CFD software and was used for verification of the proposed start-up method. The SOFC was chosen as a reference fuel cell technology for the study. Results obtained from the study are presented and discussed.

1.5D quasilinear model and its application on beams interacting with Alfven eigenmodes in DIII-D

Energy Technology Data Exchange (ETDEWEB)

Ghantous, K.; Gorelenkov, N. N. [Princeton Plasma Physics Laboratory, PO Box 451, Princeton, New Jersey 08543-0451 (United States); Berk, H. L. [Institute for Fusion Studies, University of Texas, 2100 San Jacinto Blvd., Austin, Texas 78712-1047 (United States); Heidbrink, W. W. [Department of Physics and Astronomy, University of California Irvine, Irvine, California 92697 (United States); Van Zeeland, M. A. [General Atomics, PO Box 85608, San Diego, California 92186-560 (United States)

2012-09-15

We propose a model, denoted here by 1.5D, to study energetic particle (EP) interaction with toroidal Alfvenic eigenmodes (TAE) in the case where the local EP drive for TAE exceeds the stability limit. Based on quasilinear theory, the proposed 1.5D model assumes that the particles diffuse in phase space, flattening the pressure profile until its gradient reaches a critical value where the modes stabilize. Using local theories and NOVA-K simulations of TAE damping and growth rates, the 1.5D model calculates the critical gradient and reconstructs the relaxed EP pressure profile. Local theory is improved from previous study by including more sophisticated damping and drive mechanisms such as the numerical computation of the effect of the EP finite orbit width on the growth rate. The 1.5D model is applied on the well-diagnosed DIII-D discharges no. 142111 [M. A. Van Zeeland et al., Phys. Plasmas 18, 135001 (2011)] and no. 127112 [W. W. Heidbrink et al., Nucl. Fusion. 48, 084001 (2008)]. We achieved a very satisfactory agreement with the experimental results on the EP pressure profiles redistribution and measured losses. This agreement of the 1.5D model with experimental results allows the use of this code as a guide for ITER plasma operation where it is desired to have no more than 5% loss of fusion alpha particles as limited by the design.

SiGe HBT linear-in-dB high dynamic range RF envelope detectors and wideband high linearity amplifiers

OpenAIRE

Pan, Hsuan-yu

2010-01-01

This research work aims on exploiting SiGe HBT technologies in high dynamic range wideband RF linear-in- dB envelope detectors and linear amplifiers. First, an improved all-npn broadband highly linear SiGe HBT differential amplifier is presented based on a variation of Caprio's Quad. A broadband linear amplifier with 46dBm OIP₃ at 20MHz, 34dBm OIP₃ at 1GHz, 6dB noise figure and 10.3dBm P₁dB is demonstrated. Second, an improved exact dynamic model of a fast-settling linear-in-dB Automatic Gain...

Comparison between a 1D and a 2D numerical model of an active magnetic regenerative refrigerator

DEFF Research Database (Denmark)

Petersen, Thomas Frank; Engelbrecht, Kurt; Bahl, Christian Robert Haffenden

2008-01-01

The active magnetic regenerator (AMR) refrigeration system represents an environmentally attractive alternative to vapour-compression refrigeration. This paper compares the results of two numerical AMR models: (1) a 1D finite difference model and (2) a 2D finite element model. Both models simulate...... a reciprocating AMR and can determine the cyclical steady-state temperature profile of the system as well as performance parameters such as the refrigeration capacity, the work input and the coefficient of performance (COP). The models are used to analyse an AMR with a regenerator made of flat parallel plates...... results of overall results such as the refrigeration capacity but that a 2D model is required for a detailed analysis of the phenomena occurring inside the AMR....

Dynamic PET and Optical Imaging and Compartment Modeling using a Dual-labeled Cyclic RGD Peptide Probe

Directory of Open Access Journals (Sweden)

Lei Zhu, Ning Guo, Quanzheng Li, Ying Ma, Orit Jacboson, Seulki Lee, Hak Soo Choi, James R. Mansfield, Gang Niu, Xiaoyuan Chen

2012-01-01

Full Text Available Purpose: The aim of this study is to determine if dynamic optical imaging could provide comparable kinetic parameters to that of dynamic PET imaging by a near-infrared dye/64Cu dual-labeled cyclic RGD peptide.Methods: The integrin αvβ3 binding RGD peptide was conjugated with a macrocyclic chelator 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA for copper labeling and PET imaging and a near-infrared dye ZW-1 for optical imaging. The in vitro biological activity of RGD-C(DOTA-ZW-1 was characterized by cell staining and receptor binding assay. Sixty-min dynamic PET and optical imaging were acquired on a MDA-MB-435 tumor model. Singular value decomposition (SVD method was applied to compute the dynamic optical signal from the two-dimensional optical projection images. Compartment models were used to quantitatively analyze and compare the dynamic optical and PET data.Results: The dual-labeled probe 64Cu-RGD-C(DOTA-ZW-1 showed integrin specific binding in vitro and in vivo. The binding potential (Bp derived from dynamic optical imaging (1.762 ± 0.020 is comparable to that from dynamic PET (1.752 ± 0.026.Conclusion: The signal un-mixing process using SVD improved the accuracy of kinetic modeling of 2D dynamic optical data. Our results demonstrate that 2D dynamic optical imaging with SVD analysis could achieve comparable quantitative results as dynamic PET imaging in preclinical xenograft models.

Dynamic PET and Optical Imaging and Compartment Modeling using a Dual-labeled Cyclic RGD Peptide Probe.

Science.gov (United States)

Zhu, Lei; Guo, Ning; Li, Quanzheng; Ma, Ying; Jacboson, Orit; Lee, Seulki; Choi, Hak Soo; Mansfield, James R; Niu, Gang; Chen, Xiaoyuan

2012-01-01

The aim of this study is to determine if dynamic optical imaging could provide comparable kinetic parameters to that of dynamic PET imaging by a near-infrared dye/(64)Cu dual-labeled cyclic RGD peptide. The integrin α(v)β(3) binding RGD peptide was conjugated with a macrocyclic chelator 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) for copper labeling and PET imaging and a near-infrared dye ZW-1 for optical imaging. The in vitro biological activity of RGD-C(DOTA)-ZW-1 was characterized by cell staining and receptor binding assay. Sixty-min dynamic PET and optical imaging were acquired on a MDA-MB-435 tumor model. Singular value decomposition (SVD) method was applied to compute the dynamic optical signal from the two-dimensional optical projection images. Compartment models were used to quantitatively analyze and compare the dynamic optical and PET data. The dual-labeled probe (64)Cu-RGD-C(DOTA)-ZW-1 showed integrin specific binding in vitro and in vivo. The binding potential (Bp) derived from dynamic optical imaging (1.762 ± 0.020) is comparable to that from dynamic PET (1.752 ± 0.026). The signal un-mixing process using SVD improved the accuracy of kinetic modeling of 2D dynamic optical data. Our results demonstrate that 2D dynamic optical imaging with SVD analysis could achieve comparable quantitative results as dynamic PET imaging in preclinical xenograft models.

Homology modeling, molecular dynamics and inhibitor binding study on MurD ligase of Mycobacterium tuberculosis.

Science.gov (United States)

Arvind, Akanksha; Kumar, Vivek; Saravanan, Parameswaran; Mohan, C Gopi

2012-09-01

The cell wall of mycobacterium offers well validated targets which can be exploited for discovery of new lead compounds. MurC-MurF ligases catalyze a series of irreversible steps in the biosynthesis of peptidoglycan precursor, i.e. MurD catalyzes the ligation of D-glutamate to the nucleotide precursor UMA. The three dimensional structure of Mtb-MurD is not known and was predicted by us for the first time using comparative homology modeling technique. The accuracy and stability of the predicted Mtb-MurD structure was validated using Procheck and molecular dynamics simulation. Key interactions in Mtb-MurD were studied using docking analysis of available transition state inhibitors of E.coli-MurD. The docking analysis revealed that analogues of both L and D forms of glutamic acid have similar interaction profiles with Mtb-MurD. Further, residues His192, Arg382, Ser463, and Tyr470 are proposed to be important for inhibitor-(Mtb-MurD) interactions. We also identified few pharmacophoric features essential for Mtb-MurD ligase inhibitory activity and which can further been utilized for the discovery of putative antitubercular chemotherapy.

4D Dynamic Required Navigation Performance Final Report

Science.gov (United States)

Finkelsztein, Daniel M.; Sturdy, James L.; Alaverdi, Omeed; Hochwarth, Joachim K.

2011-01-01

New advanced four dimensional trajectory (4DT) procedures under consideration for the Next Generation Air Transportation System (NextGen) require an aircraft to precisely navigate relative to a moving reference such as another aircraft. Examples are Self-Separation for enroute operations and Interval Management for in-trail and merging operations. The current construct of Required Navigation Performance (RNP), defined for fixed-reference-frame navigation, is not sufficiently specified to be applicable to defining performance levels of such air-to-air procedures. An extension of RNP to air-to-air navigation would enable these advanced procedures to be implemented with a specified level of performance. The objective of this research effort was to propose new 4D Dynamic RNP constructs that account for the dynamic spatial and temporal nature of Interval Management and Self-Separation, develop mathematical models of the Dynamic RNP constructs, "Required Self-Separation Performance" and "Required Interval Management Performance," and to analyze the performance characteristics of these air-to-air procedures using the newly developed models. This final report summarizes the activities led by Raytheon, in collaboration with GE Aviation and SAIC, and presents the results from this research effort to expand the RNP concept to a dynamic 4D frame of reference.

Agent-Based Modelling applied to 5D model of the HIV infection

Directory of Open Access Journals (Sweden)

Toufik Laroum

2016-12-01

The simplest model was the 3D mathematical model. But the complexity of this phenomenon and the diversity of cells and actors which affect its evolution requires the use of new approaches such as multi-agents approach that we have applied in this paper. The results of our simulator on the 5D model are promising because they are consistent with biological knowledge’s. Therefore, the proposed approach is well appropriate to the study of population dynamics in general and could help to understand and predict the dynamics of HIV infection.

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

Science.gov (United States)

Way, M. J.; Aleinov, I.; Amundsen, David S.; Chandler, M. A.; Clune, T. L.; Del Genio, A.; Fujii, Y.; Kelley, M.; Kiang, N. Y.; Sohl, L.;

SAHM - Simplification of one-dimensional hydraulic networks by automated processes evaluated on 1D/2D deterministic flood models

DEFF Research Database (Denmark)

Löwe, Roland; Davidsen, Steffen; Thrysøe, Cecilie

We present an algorithm for automated simplification of 1D pipe network models. The impact of the simplifications on the flooding simulated by coupled 1D-2D models is evaluated in an Australian case study. Significant reductions of the simulation time of the coupled model are achieved by reducing...... the 1D network model. The simplifications lead to an underestimation of flooded area because interaction points between network and surface are removed and because water is transported downstream faster. These effects can be mitigated by maintaining nodes in flood-prone areas in the simplification...... and by adjusting pipe roughness to increase transport times....

Predictive assessment of models for dynamic functional connectivity

DEFF Research Database (Denmark)

Nielsen, Søren Føns Vind; Schmidt, Mikkel Nørgaard; Madsen, Kristoffer Hougaard

2018-01-01

represent functional brain networks as a meta-stable process with a discrete number of states; however, there is a lack of consensus on how to perform model selection and learn the number of states, as well as a lack of understanding of how different modeling assumptions influence the estimated state......In neuroimaging, it has become evident that models of dynamic functional connectivity (dFC), which characterize how intrinsic brain organization changes over time, can provide a more detailed representation of brain function than traditional static analyses. Many dFC models in the literature...... dynamics. To address these issues, we consider a predictive likelihood approach to model assessment, where models are evaluated based on their predictive performance on held-out test data. Examining several prominent models of dFC (in their probabilistic formulations) we demonstrate our framework...

Applications of 2D helical vortex dynamics

DEFF Research Database (Denmark)

Okulov, Valery; Sørensen, Jens Nørkær

2010-01-01

In the paper, we show how the assumption of helical symmetry in the context of 2D helical vortices can be exploited to analyse and to model various cases of rotating flows. From theory, examples of three basic applications of 2D dynamics of helical vortices embedded in flows with helical symmetry...... of the vorticity field are addressed. These included some of the problems related to vortex breakdown, instability of far wakes behind rotors and vortex theory of ideal rotors....

Modeling of light dynamic cone penetration test - Panda 3 ® in granular material by using 3D Discrete element method

Science.gov (United States)

Tran, Quoc Anh; Chevalier, Bastien; Benz, Miguel; Breul, Pierre; Gourvès, Roland

2017-06-01

The recent technological developments made on the light dynamic penetration test Panda 3 ® provide a dynamic load-penetration curve σp - sp for each impact. This curve is influenced by the mechanical and physical properties of the investigated granular media. In order to analyze and exploit the load-penetration curve, a numerical model of penetration test using 3D Discrete Element Method is proposed for reproducing tests in dynamic conditions in granular media. All parameters of impact used in this model have at first been calibrated by respecting mechanical and geometrical properties of the hammer and the rod. There is a good agreement between experimental results and the ones obtained from simulations in 2D or 3D. After creating a sample, we will simulate the Panda 3 ®. It is possible to measure directly the dynamic load-penetration curve occurring at the tip for each impact. Using the force and acceleration measured in the top part of the rod, it is possible to separate the incident and reflected waves and then calculate the tip's load-penetration curve. The load-penetration curve obtained is qualitatively similar with that obtained by experimental tests. In addition, the frequency analysis of the measured signals present also a good compliance with that measured in reality when the tip resistance is qualitatively similar.

Universality for 1d Random Band Matrices: Sigma-Model Approximation

Science.gov (United States)

Shcherbina, Mariya; Shcherbina, Tatyana

2018-02-01

The paper continues the development of the rigorous supersymmetric transfer matrix approach to the random band matrices started in (J Stat Phys 164:1233-1260, 2016; Commun Math Phys 351:1009-1044, 2017). We consider random Hermitian block band matrices consisting of W× W random Gaussian blocks (parametrized by j,k \\in Λ =[1,n]^d\\cap Z^d ) with a fixed entry's variance J_{jk}=δ _{j,k}W^{-1}+β Δ _{j,k}W^{-2} , β >0 in each block. Taking the limit W→ ∞ with fixed n and β , we derive the sigma-model approximation of the second correlation function similar to Efetov's one. Then, considering the limit β , n→ ∞, we prove that in the dimension d=1 the behaviour of the sigma-model approximation in the bulk of the spectrum, as β ≫ n , is determined by the classical Wigner-Dyson statistics.

A 1D-2D coupled SPH-SWE model applied to open channel flow simulations in complicated geometries

Science.gov (United States)

Chang, Kao-Hua; Sheu, Tony Wen-Hann; Chang, Tsang-Jung

2018-05-01

In this study, a one- and two-dimensional (1D-2D) coupled model is developed to solve the shallow water equations (SWEs). The solutions are obtained using a Lagrangian meshless method called smoothed particle hydrodynamics (SPH) to simulate shallow water flows in converging, diverging and curved channels. A buffer zone is introduced to exchange information between the 1D and 2D SPH-SWE models. Interpolated water discharge values and water surface levels at the internal boundaries are prescribed as the inflow/outflow boundary conditions in the two SPH-SWE models. In addition, instead of using the SPH summation operator, we directly solve the continuity equation by introducing a diffusive term to suppress oscillations in the predicted water depth. The performance of the two approaches in calculating the water depth is comprehensively compared through a case study of a straight channel. Additionally, three benchmark cases involving converging, diverging and curved channels are adopted to demonstrate the ability of the proposed 1D and 2D coupled SPH-SWE model through comparisons with measured data and predicted mesh-based numerical results. The proposed model provides satisfactory accuracy and guaranteed convergence.

Cell Migration in 1D and 2D Nanofiber Microenvironments.

Science.gov (United States)

Estabridis, Horacio M; Jana, Aniket; Nain, Amrinder; Odde, David J

2018-03-01

Understanding how cells migrate in fibrous environments is important in wound healing, immune function, and cancer progression. A key question is how fiber orientation and network geometry influence cell movement. Here we describe a quantitative, modeling-based approach toward identifying the mechanisms by which cells migrate in fibrous geometries having well controlled orientation. Specifically, U251 glioblastoma cells were seeded onto non-electrospinning Spinneret based tunable engineering parameters fiber substrates that consist of networks of suspended 400 nm diameter nanofibers. Cells were classified based on the local fiber geometry and cell migration dynamics observed by light microscopy. Cells were found in three distinct geometries: adhering two a single fiber, adhering to two parallel fibers, and adhering to a network of orthogonal fibers. Cells adhering to a single fiber or two parallel fibers can only move in one dimension along the fiber axis, whereas cells on a network of orthogonal fibers can move in two dimensions. We found that cells move faster and more persistently in 1D geometries than in 2D, with cell migration being faster on parallel fibers than on single fibers. To explain these behaviors mechanistically, we simulated cell migration in the three different geometries using a motor-clutch based model for cell traction forces. Using nearly identical parameter sets for each of the three cases, we found that the simulated cells naturally replicated the reduced migration in 2D relative to 1D geometries. In addition, the modestly faster 1D migration on parallel fibers relative to single fibers was captured using a correspondingly modest increase in the number of clutches to reflect increased surface area of adhesion on parallel fibers. Overall, the integrated modeling and experimental analysis shows that cell migration in response to varying fibrous geometries can be explained by a simple mechanical readout of geometry via a motor-clutch mechanism.

Hydrodynamic modelling of a tidal delta wetland using an enhanced quasi-2D model

Science.gov (United States)

Wester, Sjoerd J.; Grimson, Rafael; Minotti, Priscilla G.; Booija, Martijn J.; Brugnach, Marcela

2018-04-01

Knowledge about the hydrological regime of wetlands is key to understand their physical and biological properties. Modelling hydrological and hydrodynamic processes within a wetland is therefore becoming increasingly important. 3D models have successfully modelled wetland dynamics but depend on very detailed bathymetry and land topography. Many 1D and 2D models of river deltas highly simplify the interaction between the river and wetland area or simply neglect the wetland area. This study proposes an enhanced quasi-2D modelling strategy that captures the interaction between river discharge and moon tides and the resulting hydrodynamics, while using the scarce data available. The water flow equations are discretised with an interconnected irregular cell scheme, in which a simplification of the 1D Saint-Venant equations is used to define the water flow between cells. The spatial structure of wetlands is based on the ecogeomorphology in complex estuarine deltas. The islands within the delta are modelled with levee cells, creek cells and an interior cell representing a shallow marsh wetland. The model is calibrated for an average year and the model performance is evaluated for another average year and additionally an extreme dry three-month period and an extreme wet three-month period. The calibration and evaluation are done based on two water level measurement stations and two discharge measurement stations, all located in the main rivers. Additional calibration is carried out with field water level measurements in a wetland area. Accurate simulations are obtained for both calibration and evaluation with high correlations between observed and simulated water levels and simulated discharges in the same order of magnitude as observed discharges. Calibration against field measurements showed that the model can successfully simulate the overflow mechanism in wetland areas. A sensitivity analysis for several wetland parameters showed that these parameters are all

Investigation of reflood models by coupling REFLA-1D and multi-loop system model

International Nuclear Information System (INIS)

Sugimoto, Jun; Murao, Yoshio

1983-09-01

A system analysis code REFLA-1DS was developed by coupling reflood analysis code REFLA-1D and a multi-loop primary system model. The reflood models in the code were investigated for the development of the integral system analysis code. The REFLA-1D, which was developed with the small scale reflood experiment at JAERI, consists of one-dimensional core model and a primary system model with a constant loop resistance. The multi-loop primary system model was developed with the Cylindrical Core Test Facility of JAERI's large scale reflood tests. The components modeled in the code are the upper plenum, the steam generator, the coolant pump, the ECC injection port, the downcomer and the broken cold leg nozzle. The coupling between the two models in REFLA-1DS is accomplished by applying the equivalent flow resistance calculated with the multiloop model to the REFLA-1D. The characteristics of the code is its simplicity of the system model and the solution method which enables the fast running and the easy reflood analysis for the further model development. A fairly good agreement was obtained with the results of the Cylindrical Core Test Facility for the calculated water levels in the downcomer, the core and the upper plenum. A qualitatively good agreement was obtained concerning the parametric effects of the system pressure, the ECC flow rate and the initial clad temperature. Needs for further code improvements of the models, however, were pointed out. These include the problem concerning the generation rate of the steam and water droplets in the core in an early period, the effect of the flow oscillation on the core cooling, the heat release from the downcomer wall, and the stable system calculation. (author)

Dynamic phase transition in the kinetic spin-2 Blume-Emery-Griffiths model in an oscillating field

International Nuclear Information System (INIS)

Ertas, Mehmet; Canko, Osman; Keskin, Mustafa

2008-01-01

We extend our recent paper [M. Keskin, O. Canko, M. Ertas, J. Exp. Theor. Phys. (Sov. Phys. JETP) 105 (2007) 1190.] to present a study, within a mean-field approach, the stationary states of the kinetic spin-2 Blume-Emery-Griffiths model in the presence of a time-dependent oscillating magnetic field by using the Glauber-type of stochastic dynamics. We found 20 fundamental types of dynamic phase diagrams where exhibit more complex and richer phase diagrams than our recent paper. Especially, the obtained dynamic phase diagrams show the dynamic triple, quadruple and dynamic double critical end points besides dynamic tricritical points that depending on interaction parameters. The phase diagrams also exhibit a disordered (d) and the ferromagnetic-2 (f 2 ) phases, and the f 2 +d, f 2 +fq, fq+d, f 2 +f 1 +fq and f 2 +fq+d, where f 1 are fq the ferromagnetic-1 and ferroquadrupolar or simply quadrupolar phases respectively, coexistence phase regions that strongly depend on interaction parameters

Dynamic phase transition in the kinetic spin-2 Blume-Emery-Griffiths model in an oscillating field

Science.gov (United States)

Ertaş, Mehmet; Canko, Osman; Keskin, Mustafa

We extend our recent paper [M. Keskin, O. Canko, M. Ertaş, J. Exp. Theor. Phys. (Sov. Phys. JETP) 105 (2007) 1190.] to present a study, within a mean-field approach, the stationary states of the kinetic spin-2 Blume-Emery-Griffiths model in the presence of a time-dependent oscillating magnetic field by using the Glauber-type of stochastic dynamics. We found 20 fundamental types of dynamic phase diagrams where exhibit more complex and richer phase diagrams than our recent paper. Especially, the obtained dynamic phase diagrams show the dynamic triple, quadruple and dynamic double critical end points besides dynamic tricritical points that depending on interaction parameters. The phase diagrams also exhibit a disordered ( d) and the ferromagnetic-2 ( f2) phases, and the f2+ d, f2+ fq, fq+ d, f2+ f1+ fq and f2+ fq+ d, where f1 are fq the ferromagnetic-1 and ferroquadrupolar or simply quadrupolar phases respectively, coexistence phase regions that strongly depend on interaction parameters.

Minimal unitary representation of D(2,1;λ) and its SU(2) deformations and d=1, N=4 superconformal models

International Nuclear Information System (INIS)

Govil, Karan; Gunaydin, Murat

2013-01-01

Quantization of the geometric quasiconformal realizations of noncompact groups and supergroups leads directly to their minimal unitary representations (minreps). Using quasiconformal methods massless unitary supermultiplets of superconformal groups SU(2,2|N) and OSp(8 ⁎ |2n) in four and six dimensions were constructed as minreps and their U(1) and SU(2) deformations, respectively. In this paper we extend these results to SU(2) deformations of the minrep of N=4 superconformal algebra D(2,1;λ) in one dimension. We find that SU(2) deformations can be achieved using n pair of bosons and m pairs of fermions simultaneously. The generators of deformed minimal representations of D(2,1;λ) commute with the generators of a dual superalgebra OSp(2n ⁎ |2m) realized in terms of these bosons and fermions. We show that there exists a precise mapping between symmetry generators of N=4 superconformal models in harmonic superspace studied recently and minimal unitary supermultiplets of D(2,1;λ) deformed by a pair of bosons. This can be understood as a particular case of a general mapping between the spectra of quantum mechanical quaternionic Kähler sigma models with eight super symmetries and minreps of their isometry groups that descends from the precise mapping established between the 4d, N=2 sigma models coupled to supergravity and minreps of their isometry groups.

Modelling turbulent vertical mixing sensitivity using a 1-D version of NEMO

Science.gov (United States)

Reffray, G.; Bourdalle-Badie, R.; Calone, C.

2015-01-01

Through two numerical experiments, a 1-D vertical model called NEMO1D was used to investigate physical and numerical turbulent-mixing behaviour. The results show that all the turbulent closures tested (k+l from Blanke and Delecluse, 1993, and two equation models: generic length scale closures from Umlauf and Burchard, 2003) are able to correctly reproduce the classical test of Kato and Phillips (1969) under favourable numerical conditions while some solutions may diverge depending on the degradation of the spatial and time discretization. The performances of turbulence models were then compared with data measured over a 1-year period (mid-2010 to mid-2011) at the PAPA station, located in the North Pacific Ocean. The modelled temperature and salinity were in good agreement with the observations, with a maximum temperature error between -2 and 2 °C during the stratified period (June to October). However, the results also depend on the numerical conditions. The vertical RMSE varied, for different turbulent closures, from 0.1 to 0.3 °C during the stratified period and from 0.03 to 0.15 °C during the homogeneous period. This 1-D configuration at the PAPA station (called PAPA1D) is now available in NEMO as a reference configuration including the input files and atmospheric forcing set described in this paper. Thus, all the results described can be recovered by downloading and launching PAPA1D. The configuration is described on the NEMO site (PAPA">http://www.nemo-ocean.eu/Using-NEMO/Configurations/C1D_PAPA). This package is a good starting point for further investigation of vertical processes.

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

International Nuclear Information System (INIS)

Way, M. J.; Aleinov, I.; Amundsen, David S.; Chandler, M. A.; Genio, A. D. Del; Fujii, Y.; Kelley, M.; Kiang, N. Y.; Sohl, L.; Tsigaridis, K.; Clune, T. L.

2017-01-01

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a three-dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of solar system and exoplanetary terrestrial planets. Its parent model, known as ModelE2, is used to simulate modern Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions, including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (from slower to more rapid than modern Earth’s, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the solar system such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn’s moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents, we can then further expand its capabilities to those exoplanetary rocky worlds that have been discovered in the past, as well as those to be discovered in the future. We also discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

Energy Technology Data Exchange (ETDEWEB)

Way, M. J.; Aleinov, I.; Amundsen, David S.; Chandler, M. A.; Genio, A. D. Del; Fujii, Y.; Kelley, M.; Kiang, N. Y.; Sohl, L.; Tsigaridis, K. [NASA Goddard Institute for Space Studies, New York, NY 10025 (United States); Clune, T. L. [Global Modeling and Assimilation Office, NASA Goddard Space Flight Center (United States)

2017-07-01

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a three-dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of solar system and exoplanetary terrestrial planets. Its parent model, known as ModelE2, is used to simulate modern Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions, including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (from slower to more rapid than modern Earth’s, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the solar system such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn’s moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents, we can then further expand its capabilities to those exoplanetary rocky worlds that have been discovered in the past, as well as those to be discovered in the future. We also discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.

Coupled iterated map models of action potential dynamics in a one-dimensional cable of cardiac cells

International Nuclear Information System (INIS)

Wang Shihong; Xie Yuanfang; Qu Zhilin

2008-01-01

Low-dimensional iterated map models have been widely used to study action potential dynamics in isolated cardiac cells. Coupled iterated map models have also been widely used to investigate action potential propagation dynamics in one-dimensional (1D) coupled cardiac cells, however, these models are usually empirical and not carefully validated. In this study, we first developed two coupled iterated map models which are the standard forms of diffusively coupled maps and overcome the limitations of the previous models. We then determined the coupling strength and space constant by quantitatively comparing the 1D action potential duration profile from the coupled cardiac cell model described by differential equations with that of the coupled iterated map models. To further validate the coupled iterated map models, we compared the stability conditions of the spatially uniform state of the coupled iterated maps and those of the 1D ionic model and showed that the coupled iterated map model could well recapitulate the stability conditions, i.e. the spatially uniform state is stable unless the state is chaotic. Finally, we combined conduction into the developed coupled iterated map model to study the effects of coupling strength on wave stabilities and showed that the diffusive coupling between cardiac cells tends to suppress instabilities during reentry in a 1D ring and the onset of discordant alternans in a periodically paced 1D cable

Sphingomyelinase D activity in model membranes: structural effects of in situ generation of ceramide-1-phosphate.

Directory of Open Access Journals (Sweden)

Roberto P Stock

Full Text Available The toxicity of Loxosceles spider venom has been attributed to a rare enzyme, sphingomyelinase D, which transforms sphingomyelin to ceramide-1-phosphate. The bases of its inflammatory and dermonecrotic activity, however, remain unclear. In this work the effects of ceramide-1-phosphate on model membranes were studied both by in situ generation of this lipid using a recombinant sphingomyelinase D from the spider Loxosceles laeta and by pre-mixing it with sphingomyelin and cholesterol. The systems of choice were large unilamellar vesicles for bulk studies (enzyme kinetics, fluorescence spectroscopy and dynamic light scattering and giant unilamellar vesicles for fluorescence microscopy examination using a variety of fluorescent probes. The influence of membrane lateral structure on the kinetics of enzyme activity and the consequences of enzyme activity on the structure of target membranes containing sphingomyelin were examined. The findings indicate that: 1 ceramide-1-phosphate (particularly lauroyl ceramide-1-phosphate can be incorporated into sphingomyelin bilayers in a concentration-dependent manner and generates coexistence of liquid disordered/solid ordered domains, 2 the activity of sphingomyelinase D is clearly influenced by the supramolecular organization of its substrate in membranes and, 3 in situ ceramide-1-phosphate generation by enzymatic activity profoundly alters the lateral structure and morphology of the target membranes.

Modeling Microstructural Evolution During Dynamic Recrystallization of Alloy D9 Using Artificial Neural Network

Science.gov (United States)

Mandal, Sumantra; Sivaprasad, P. V.; Dube, R. K.

2007-12-01

An artificial neural network (ANN) model was developed to predict the microstructural evolution of a 15Cr-15Ni-2.2Mo-Ti modified austenitic stainless steel (Alloy D9) during dynamic recrystallization (DRX). The input parameters were strain, strain rate, and temperature whereas microstructural features namely, %DRX and average grain size were the output parameters. The ANN was trained with the database obtained from various industrial scale metal-forming operations like forge hammer, hydraulic press, and rolling carried out in the temperature range 1173-1473 K to various strain levels. The performance of the model was evaluated using a wide variety of statistical indices and the predictability of the model was found to be good. The combined influence of temperature and strain on microstructural features has been simulated employing the developed model. The results were found to be consistent with the relevant fundamental metallurgical phenomena.

Stochastic modeling for dynamics of HIV-1 infection using cellular automata: A review.

Science.gov (United States)

Precharattana, Monamorn

2016-02-01

Recently, the description of immune response by discrete models has emerged to play an important role to study the problems in the area of human immunodeficiency virus type 1 (HIV-1) infection, leading to AIDS. As infection of target immune cells by HIV-1 mainly takes place in the lymphoid tissue, cellular automata (CA) models thus represent a significant step in understanding when the infected population is dispersed. Motivated by these, the studies of the dynamics of HIV-1 infection using CA in memory have been presented to recognize how CA have been developed for HIV-1 dynamics, which issues have been studied already and which issues still are objectives in future studies.

Coupled 3D neutronics/thermal hydraulics modeling of the SAFARI-1 MTR

International Nuclear Information System (INIS)

Rosenkrantz, Adam; Avramova, Maria; Ivanov, Kostadin; Prinsloo, Rian; Botes, Danniëll; Elsakhawy, Khalid

2014-01-01

Highlights: • Development of 3D coupled neutronics/thermal–hydraulic model of SAFARI-1. • Verification of 3D steady-state NEM based neutronics model for SAFARI-1. • Verification of 3D COBRA-TF based thermal–hydraulic model of SAFARI-1. • Quantification of the effect of correct modeling of thermal–hydraulic feedback. - Abstract: The purpose of this study was to develop a coupled accurate multi-physics model of the SAFARI-1 Material Testing Reactor (MTR), a facility that is used for both research and the production of medical isotopes. The model was developed as part of the SAFARI-1 benchmarking project as a cooperative effort between the Pennsylvania State University (PSU) and the South African Nuclear Energy Corporation (Necsa). It was created using a multi-physics coupling of state of the art nuclear reactor simulation tools, consisting of a neutronics code and a thermal hydraulics code. The neutronics tool used was the PSU code NEM, and the results from this component were verified using the Necsa neutronics code OSCAR-4, which is utilized for SAFARI-1 core design and fuel management. On average, the multiplication factors of the neutronics models agreed to within 5 pcm and the radial assembly-averaged powers agreed to within 0.2%. The thermal hydraulics tool used was the PSU version of COBRA-TF (CTF) sub-channel code, and the results of this component were verified against another thermal hydraulics code, the RELAP5-3D system code, used at Necsa for thermal–hydraulics analysis of SAFARI-1. Although only assembly-averaged results from RELAP5-3D were available, they fell within the range of values for the corresponding assemblies in the comprehensive CTF solution. This comparison allows for the first time to perform a quantification of steady-state errors for a low-powered MTR with an advanced thermal–hydraulic code such as CTF on a per-channel basis as compared to simpler and coarser-mesh RELAP5-3D modeling. Additionally, a new cross section

Limb darkening laws for two exoplanet host stars derived from 3D stellar model atmospheres. Comparison with 1D models and HST light curve observations

Science.gov (United States)

Hayek, W.; Sing, D.; Pont, F.; Asplund, M.

2012-03-01

We compare limb darkening laws derived from 3D hydrodynamical model atmospheres and 1D hydrostatic MARCS models for the host stars of two well-studied transiting exoplanet systems, the late-type dwarfs HD 209458 and HD 189733. The surface brightness distribution of the stellar disks is calculated for a wide spectral range using 3D LTE spectrum formation and opacity sampling⋆. We test our theoretical predictions using least-squares fits of model light curves to wavelength-integrated primary eclipses that were observed with the Hubble Space Telescope (HST). The limb darkening law derived from the 3D model of HD 209458 in the spectral region between 2900 Å and 5700 Å produces significantly better fits to the HST data, removing systematic residuals that were previously observed for model light curves based on 1D limb darkening predictions. This difference arises mainly from the shallower mean temperature structure of the 3D model, which is a consequence of the explicit simulation of stellar surface granulation where 1D models need to rely on simplified recipes. In the case of HD 189733, the model atmospheres produce practically equivalent limb darkening curves between 2900 Å and 5700 Å, partly due to obstruction by spectral lines, and the data are not sufficient to distinguish between the light curves. We also analyze HST observations between 5350 Å and 10 500 Å for this star; the 3D model leads to a better fit compared to 1D limb darkening predictions. The significant improvement of fit quality for the HD 209458 system demonstrates the higher degree of realism of 3D hydrodynamical models and the importance of surface granulation for the formation of the atmospheric radiation field of late-type stars. This result agrees well with recent investigations of limb darkening in the solar continuum and other observational tests of the 3D models. The case of HD 189733 is no contradiction as the model light curves are less sensitive to the temperature stratification of

1.5D quasilinear model and its application on beams interacting with Alfvén eigenmodes in DIII-D

Science.gov (United States)

Ghantous, K.; Gorelenkov, N. N.; Berk, H. L.; Heidbrink, W. W.; Van Zeeland, M. A.

2012-09-01

We propose a model, denoted here by 1.5D, to study energetic particle (EP) interaction with toroidal Alfvenic eigenmodes (TAE) in the case where the local EP drive for TAE exceeds the stability limit. Based on quasilinear theory, the proposed 1.5D model assumes that the particles diffuse in phase space, flattening the pressure profile until its gradient reaches a critical value where the modes stabilize. Using local theories and NOVA-K simulations of TAE damping and growth rates, the 1.5D model calculates the critical gradient and reconstructs the relaxed EP pressure profile. Local theory is improved from previous study by including more sophisticated damping and drive mechanisms such as the numerical computation of the effect of the EP finite orbit width on the growth rate. The 1.5D model is applied on the well-diagnosed DIII-D discharges #142111 [M. A. Van Zeeland et al., Phys. Plasmas 18, 135001 (2011)] and #127112 [W. W. Heidbrink et al., Nucl. Fusion. 48, 084001 (2008)]. We achieved a very satisfactory agreement with the experimental results on the EP pressure profiles redistribution and measured losses. This agreement of the 1.5D model with experimental results allows the use of this code as a guide for ITER plasma operation where it is desired to have no more than 5% loss of fusion alpha particles as limited by the design.

Hypercrater Bifurcations, Attractor Coexistence, and Unfolding in a 5D Model of Economic Dynamics

Directory of Open Access Journals (Sweden)

Toichiro Asada

2011-01-01

Full Text Available Complex dynamical features are explored in a discrete interregional macrodynamic model proposed by Asada et al., using numerical methods. The model is five-dimensional with four parameters. The results demonstrate patterns of dynamical behaviour, such as bifurcation processes and coexistence of attractors, generated by high-dimensional discrete systems. In three cases of two-dimensional parameter subspaces the stability of equilibrium region is determined and its boundaries, the flip and Neimark-Hopf bifurcation curves, are identified by means of necessary coefficient criteria. In the first case closed invariant curves (CICs are found to occur through 5D-crater-type bifurcations, and for certain ranges of parameter values a stable equilibrium coexists with an unstable CIC associated with the subcritical bifurcation, as well as with an outer stable CIC. A remarkable feature of the second case is the coexistence of two attracting CICs outside the stability region. In both these cases the related hysteresis effects are illustrated by numerical simulations. In the third case a remarkable feature is the apparent unfolding of an attracting CIC before it evolves to a chaotic attractor. Examples of CICs and chaotic attractors are given in subspaces of phase space.

An Evaluative Review of Simulated Dynamic Smart 3d Objects

Science.gov (United States)

Romeijn, H.; Sheth, F.; Pettit, C. J.

2012-07-01

Three-dimensional (3D) modelling of plants can be an asset for creating agricultural based visualisation products. The continuum of 3D plants models ranges from static to dynamic objects, also known as smart 3D objects. There is an increasing requirement for smarter simulated 3D objects that are attributed mathematically and/or from biological inputs. A systematic approach to plant simulation offers significant advantages to applications in agricultural research, particularly in simulating plant behaviour and the influences of external environmental factors. This approach of 3D plant object visualisation is primarily evident from the visualisation of plants using photographed billboarded images, to more advanced procedural models that come closer to simulating realistic virtual plants. However, few programs model physical reactions of plants to external factors and even fewer are able to grow plants based on mathematical and/or biological parameters. In this paper, we undertake an evaluation of plant-based object simulation programs currently available, with a focus upon the components and techniques involved in producing these objects. Through an analytical review process we consider the strengths and weaknesses of several program packages, the features and use of these programs and the possible opportunities in deploying these for creating smart 3D plant-based objects to support agricultural research and natural resource management. In creating smart 3D objects the model needs to be informed by both plant physiology and phenology. Expert knowledge will frame the parameters and procedures that will attribute the object and allow the simulation of dynamic virtual plants. Ultimately, biologically smart 3D virtual plants that react to changes within an environment could be an effective medium to visually represent landscapes and communicate land management scenarios and practices to planners and decision-makers.

Impact of time-of-flight on indirect 3D and direct 4D parametric image reconstruction in the presence of inconsistent dynamic PET data

Science.gov (United States)

Kotasidis, F. A.; Mehranian, A.; Zaidi, H.

2016-05-01

Kinetic parameter estimation in dynamic PET suffers from reduced accuracy and precision when parametric maps are estimated using kinetic modelling following image reconstruction of the dynamic data. Direct approaches to parameter estimation attempt to directly estimate the kinetic parameters from the measured dynamic data within a unified framework. Such image reconstruction methods have been shown to generate parametric maps of improved precision and accuracy in dynamic PET. However, due to the interleaving between the tomographic and kinetic modelling steps, any tomographic or kinetic modelling errors in certain regions or frames, tend to spatially or temporally propagate. This results in biased kinetic parameters and thus limits the benefits of such direct methods. Kinetic modelling errors originate from the inability to construct a common single kinetic model for the entire field-of-view, and such errors in erroneously modelled regions could spatially propagate. Adaptive models have been used within 4D image reconstruction to mitigate the problem, though they are complex and difficult to optimize. Tomographic errors in dynamic imaging on the other hand, can originate from involuntary patient motion between dynamic frames, as well as from emission/transmission mismatch. Motion correction schemes can be used, however, if residual errors exist or motion correction is not included in the study protocol, errors in the affected dynamic frames could potentially propagate either temporally, to other frames during the kinetic modelling step or spatially, during the tomographic step. In this work, we demonstrate a new strategy to minimize such error propagation in direct 4D image reconstruction, focusing on the tomographic step rather than the kinetic modelling step, by incorporating time-of-flight (TOF) within a direct 4D reconstruction framework. Using ever improving TOF resolutions (580 ps, 440 ps, 300 ps and 160 ps), we demonstrate that direct 4D TOF image

The string difference equation of the D = 1 matrix model and W1+∞ symmetry of the KP hierarchy

International Nuclear Information System (INIS)

Awada, M.A.; Sin, S.J.

1992-01-01

In this paper, the authors give a connection between the D = 1 matrix model and the generalized KP hierarchy. First, the authors find a difference equation satisfied by F, the Legendre transformation of the free energy of the D = 1 matrix model on a circle of radius R. Then the authors show that it is a special case of the difference equation of the generalized KP hierarchy with its zero mode identified with the scaling variable of the D = 1 string theory. The authors propose that the massive D = 1 matrix model is described by the generalized KP hierarchy, which implies the manifest integrability of D = 1 string theory. The authors also show that the (generalized) KP hierarchy has an underlying W 1 + ∞ symmetry. By reduction, we prove that the generalized KdV hierarchy has a subalgebra of the above symmetry which again forms a W 1+ ∞ . The authors argue that there are no W constraints in D = 1 string theory, which is in contrast to D 1 + ∞ constraints

AxiSEM3D: broadband seismic wavefields in 3-D aspherical Earth models

Science.gov (United States)

Leng, K.; Nissen-Meyer, T.; Zad, K. H.; van Driel, M.; Al-Attar, D.

2017-12-01

Seismology is the primary tool for data-informed inference of Earth structure and dynamics. Simulating seismic wave propagation at a global scale is fundamental to seismology, but remains as one of most challenging problems in scientific computing, because of both the multiscale nature of Earth's interior and the observable frequency band of seismic data. We present a novel numerical method to simulate global seismic wave propagation in realistic 3-D Earth models. Our method, named AxiSEM3D, is a hybrid of spectral element method and pseudospectral method. It reduces the azimuthal dimension of wavefields by means of a global Fourier series parameterization, of which the number of terms can be locally adapted to the inherent azimuthal smoothness of the wavefields. AxiSEM3D allows not only for material heterogeneities, such as velocity, density, anisotropy and attenuation, but also for finite undulations on radial discontinuities, both solid-solid and solid-fluid, and thereby a variety of aspherical Earth features such as ellipticity, topography, variable crustal thickness, and core-mantle boundary topography. Such interface undulations are equivalently interpreted as material perturbations of the contiguous media, based on the "particle relabelling transformation". Efficiency comparisons show that AxiSEM3D can be 1 to 3 orders of magnitude faster than conventional 3-D methods, with the speedup increasing with simulation frequency and decreasing with model complexity, but for all realistic structures the speedup remains at least one order of magnitude. The observable frequency range of global seismic data (up to 1 Hz) has been covered for wavefield modelling upon a 3-D Earth model with reasonable computing resources. We show an application of surface wave modelling within a state-of-the-art global crustal model (Crust1.0), with the synthetics compared to real data. The high-performance C++ code is released at github.com/AxiSEM3D/AxiSEM3D.

New insights into the short pitch corrugation development enigma based on 3D-FE dynamic vehicle-track coupled modelling in frictional rolling contact

NARCIS (Netherlands)

Li, S.; Li, Z.; Nunez Vicencio, Alfredo; Dollevoet, R.P.B.J.

2017-01-01

A three-dimensional (3D) finite element (FE) dynamic frictional rolling contact model is presented for the study of short pitch corrugation that considers direct and instantaneous coupling between the contact mechanics and the structural dynamics in a vehicle-track system. In this study, we examine

Modeling and Optimal Control of a Class of Warfare Hybrid Dynamic Systems Based on Lanchester (n,1 Attrition Model

Directory of Open Access Journals (Sweden)

Xiangyong Chen

2014-01-01

hybrid dynamic systems is established based on Lanchester equation in a (n,1 battle, where a heterogeneous force of n different troop types faces a homogeneous force. This model can be characterized by the interaction of continuous-time models (governed by Lanchester equation, and discrete event systems (described by variable tactics. Furthermore, an expository discussion is presented on an optimal variable tactics control problem for warfare hybrid dynamic system. The optimal control strategies are designed based on dynamic programming and differential game theory. As an example of the consequences of this optimal control problem, we take the (2, 1 case and solve the optimal strategies in a (2, 1 case. Simulation results show the feasibility of warfare hybrid system model and the effectiveness of the optimal control strategies designed.

Integrating intracellular dynamics using CompuCell3D and Bionetsolver: applications to multiscale modelling of cancer cell growth and invasion.

Directory of Open Access Journals (Sweden)

Vivi Andasari

Full Text Available In this paper we present a multiscale, individual-based simulation environment that integrates CompuCell3D for lattice-based modelling on the cellular level and Bionetsolver for intracellular modelling. CompuCell3D or CC3D provides an implementation of the lattice-based Cellular Potts Model or CPM (also known as the Glazier-Graner-Hogeweg or GGH model and a Monte Carlo method based on the metropolis algorithm for system evolution. The integration of CC3D for cellular systems with Bionetsolver for subcellular systems enables us to develop a multiscale mathematical model and to study the evolution of cell behaviour due to the dynamics inside of the cells, capturing aspects of cell behaviour and interaction that is not possible using continuum approaches. We then apply this multiscale modelling technique to a model of cancer growth and invasion, based on a previously published model of Ramis-Conde et al. (2008 where individual cell behaviour is driven by a molecular network describing the dynamics of E-cadherin and β-catenin. In this model, which we refer to as the centre-based model, an alternative individual-based modelling technique was used, namely, a lattice-free approach. In many respects, the GGH or CPM methodology and the approach of the centre-based model have the same overall goal, that is to mimic behaviours and interactions of biological cells. Although the mathematical foundations and computational implementations of the two approaches are very different, the results of the presented simulations are compatible with each other, suggesting that by using individual-based approaches we can formulate a natural way of describing complex multi-cell, multiscale models. The ability to easily reproduce results of one modelling approach using an alternative approach is also essential from a model cross-validation standpoint and also helps to identify any modelling artefacts specific to a given computational approach.

Dynamic phase transition in the kinetic spin-2 Blume-Emery-Griffiths model in an oscillating field

Energy Technology Data Exchange (ETDEWEB)

Ertas, Mehmet [Institute of Science, Erciyes University, 38039 Kayseri (Turkey); Canko, Osman [Department of Physics, Erciyes University, 38039 Kayseri (Turkey); Keskin, Mustafa [Department of Physics, Erciyes University, 38039 Kayseri (Turkey)], E-mail: keskin@erciyes.edu.tr

2008-06-15

We extend our recent paper [M. Keskin, O. Canko, M. Ertas, J. Exp. Theor. Phys. (Sov. Phys. JETP) 105 (2007) 1190.] to present a study, within a mean-field approach, the stationary states of the kinetic spin-2 Blume-Emery-Griffiths model in the presence of a time-dependent oscillating magnetic field by using the Glauber-type of stochastic dynamics. We found 20 fundamental types of dynamic phase diagrams where exhibit more complex and richer phase diagrams than our recent paper. Especially, the obtained dynamic phase diagrams show the dynamic triple, quadruple and dynamic double critical end points besides dynamic tricritical points that depending on interaction parameters. The phase diagrams also exhibit a disordered (d) and the ferromagnetic-2 (f{sub 2}) phases, and the f{sub 2}+d, f{sub 2}+fq, fq+d, f{sub 2}+f{sub 1}+fq and f{sub 2}+fq+d, where f{sub 1} are fq the ferromagnetic-1 and ferroquadrupolar or simply quadrupolar phases respectively, coexistence phase regions that strongly depend on interaction parameters.

Dynamic phase transition properties for the mixed spin-(1/2, 1) Ising model in an oscillating magnetic field

Energy Technology Data Exchange (ETDEWEB)

Ertaş, Mehmet, E-mail: mehmetertas@erciyes.edu.tr; Keskin, Mustafa

2015-08-15

Herein we study the dynamic phase transition properties for the mixed spin-(1/2, 1) Ising model on a square lattice under a time-dependent magnetic field by means of the effective-field theory (EFT) with correlations based on Glauber dynamics. We present the dynamic phase diagrams in the reduced magnetic field amplitude and reduced temperature plane and find that the phase diagrams exhibit dynamic tricitical behavior, multicritical and zero-temperature critical points as well as reentrant behavior. We also investigate the influence of frequency (ω) and observe that for small values of ω the mixed phase disappears, but for high values it appears and the system displays reentrant behavior as well as a critical end point. - Highlights: • Dynamic behaviors of a ferrimagnetic mixed spin (1/2, 1) Ising system are studied. • We examined the effects of the Hamiltonian parameters on the dynamic behaviors. • The phase diagrams are obtained in (T-h) plane. • The dynamic phase diagrams exhibit the dynamic tricritical and reentrant behaviors.

Dynamic phase transition properties for the mixed spin-(1/2, 1) Ising model in an oscillating magnetic field

International Nuclear Information System (INIS)

Ertaş, Mehmet; Keskin, Mustafa

2015-01-01

Herein we study the dynamic phase transition properties for the mixed spin-(1/2, 1) Ising model on a square lattice under a time-dependent magnetic field by means of the effective-field theory (EFT) with correlations based on Glauber dynamics. We present the dynamic phase diagrams in the reduced magnetic field amplitude and reduced temperature plane and find that the phase diagrams exhibit dynamic tricitical behavior, multicritical and zero-temperature critical points as well as reentrant behavior. We also investigate the influence of frequency (ω) and observe that for small values of ω the mixed phase disappears, but for high values it appears and the system displays reentrant behavior as well as a critical end point. - Highlights: • Dynamic behaviors of a ferrimagnetic mixed spin (1/2, 1) Ising system are studied. • We examined the effects of the Hamiltonian parameters on the dynamic behaviors. • The phase diagrams are obtained in (T-h) plane. • The dynamic phase diagrams exhibit the dynamic tricritical and reentrant behaviors

Integrated Biogeomorphological Modeling Using Delft3D

Science.gov (United States)

Ye, Q.; Jagers, B.

2011-12-01

The skill of numerical morphological models has improved significantly from the early 2D uniform, total load sediment models (with steady state or infrequent wave updates) to recent 3D hydrodynamic models with multiple suspended and bed load sediment fractions and bed stratigraphy (online coupled with waves). Although there remain many open questions within this combined field of hydro- and morphodynamics, we observe an increasing need to include biological processes in the overall dynamics. In riverine and inter-tidal environments, there is often an important influence by riparian vegetation and macrobenthos. Over the past decade more and more researchers have started to extend the simulation environment with wrapper scripts and other quick code hacks to estimate their influence on morphological development in coastal, estuarine and riverine environments. Although one can in this way quickly analyze different approaches, these research tools have generally not been designed with reuse, performance and portability in mind. We have now implemented a reusable, flexible, and efficient two-way link between the Delft3D open source framework for hydrodynamics, waves and morphology, and the water quality and ecology modules. The same link will be used for 1D, 2D and 3D modeling on networks and both structured and unstructured grids. We will describe the concepts of the overall system, and illustrate it with some first results.

Dynamics model for real time diagnostics of Triga RC-1 system

International Nuclear Information System (INIS)

Gadomski, A.M.; Nanni, V.; Meo, G.

1988-01-01

This paper presents dynamics model of TRIGA RC-1 reactor system. The model is dedicated to the real-time early fault detection during a reactor operation in one week exploitation cycle. The algorithms are specially suited for real-time, long time and also accelerated simulations with assumed diagnostic oriented accuracy. The approximations, modular structure, numerical methods and validation are discussed. The elaborated model will be build in the TRIGA Supervisor System and TRIGA Diagnostic Simulator

Dynamical structure of water in aqueous solutions of D-glucose and D-galactose by low-frequency Raman scattering

Science.gov (United States)

Wang, Yan; Tominaga, Yasunori

1994-02-01

Low-frequency depolarized Raman spectra of aqueous solutions of D-glucose and D-galactose have been investigated in the frequency region from -250 cm-1 to 250 cm-1 at 30.0 °C as a function of concentration up to 0.04 molar ratio. The dynamical structure of water in aqueous solution is analyzed by using the reduced Raman spectrum χ`(ν¯), which corresponds to the imaginary part of the dynamical susceptibility. The reduced spectrum is fitted with the superposition of one Cole-Cole type relaxation mode and two damped harmonic oscillator modes by a nonlinear least-squares fitting. The effects of D-glucose and D-galactose on the dynamical structure of water in aqueous solution are similar. The relaxation time of hydrogen bond among water molecules becomes slower with increasing sugar concentration. The characteristic frequencies of stretching-like and bending-like vibrations among water molecules do not change in both D-glucose and D-galactose aqueous solutions.

A Tensor Statistical Model for Quantifying Dynamic Functional Connectivity.

Science.gov (United States)

Zhu, Yingying; Zhu, Xiaofeng; Kim, Minjeong; Yan, Jin; Wu, Guorong

2017-06-01

Functional connectivity (FC) has been widely investigated in many imaging-based neuroscience and clinical studies. Since functional Magnetic Resonance Image (MRI) signal is just an indirect reflection of brain activity, it is difficult to accurately quantify the FC strength only based on signal correlation. To address this limitation, we propose a learning-based tensor model to derive high sensitivity and specificity connectome biomarkers at the individual level from resting-state fMRI images. First, we propose a learning-based approach to estimate the intrinsic functional connectivity. In addition to the low level region-to-region signal correlation, latent module-to-module connection is also estimated and used to provide high level heuristics for measuring connectivity strength. Furthermore, sparsity constraint is employed to automatically remove the spurious connections, thus alleviating the issue of searching for optimal threshold. Second, we integrate our learning-based approach with the sliding-window technique to further reveal the dynamics of functional connectivity. Specifically, we stack the functional connectivity matrix within each sliding window and form a 3D tensor where the third dimension denotes for time. Then we obtain dynamic functional connectivity (dFC) for each individual subject by simultaneously estimating the within-sliding-window functional connectivity and characterizing the across-sliding-window temporal dynamics. Third, in order to enhance the robustness of the connectome patterns extracted from dFC, we extend the individual-based 3D tensors to a population-based 4D tensor (with the fourth dimension stands for the training subjects) and learn the statistics of connectome patterns via 4D tensor analysis. Since our 4D tensor model jointly (1) optimizes dFC for each training subject and (2) captures the principle connectome patterns, our statistical model gains more statistical power of representing new subject than current state

3D-QSAR (CoMFA, CoMSIA), molecular docking and molecular dynamics simulations study of 6-aryl-5-cyano-pyrimidine derivatives to explore the structure requirements of LSD1 inhibitors.

Science.gov (United States)

Ding, Lina; Wang, Zhi-Zheng; Sun, Xu-Dong; Yang, Jing; Ma, Chao-Ya; Li, Wen; Liu, Hong-Min

2017-08-01

Recently, Histone Lysine Specific Demethylase 1 (LSD1) was regarded as a promising anticancer target for the novel drug discovery. And several small molecules as LSD1 inhibitors in different structures have been reported. In this work, we carried out a molecular modeling study on the 6-aryl-5-cyano-pyrimidine fragment LSD1 inhibitors using three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q 2 =0.802, r 2 ncv =0.979, and the best CoMSIA model has q 2 =0.799, r 2 ncv =0.982. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300K. All the results can provide us more useful information for our further drug design. Copyright © 2017. Published by Elsevier Ltd.

Swimming attenuates d-galactose-induced brain aging via suppressing miR-34a-mediated autophagy impairment and abnormal mitochondrial dynamics.

Science.gov (United States)

Kou, Xianjuan; Li, Jie; Liu, Xingran; Chang, Jingru; Zhao, Qingxia; Jia, Shaohui; Fan, Jingjing; Chen, Ning

2017-06-01

microRNAs (miRNAs) have been reported to be involved in many neurodegenerative diseases. To explore the regulatory role of miR-34a in aging-related diseases such as Alzheimer's disease (AD) during exercise intervention, we constructed a rat model with d-galactose (d-gal)-induced oxidative stress and cognitive impairment coupled with dysfunctional autophagy and abnormal mitochondrial dynamics, determined the mitigation of cognitive impairment of d-gal-induced aging rats during swimming intervention, and evaluated miR-34a-mediated functional status of autophagy and abnormal mitochondrial dynamics. Meanwhile, whether the upregulation of miR-34a can lead to dysfunctional autophagy and abnormal mitochondrial dynamics was confirmed in human SH-SY5Y cells with silenced miR-34a by the transfection of a miR-34a inhibitor. Results indicated that swimming intervention could significantly attenuate cognitive impairment, prevent the upregulation of miR-34a, mitigate the dysfunctional autophagy, and inhibit the increase of dynamin-related protein 1 (DRP1) in d-gal-induced aging model rats. In contrast, the miR-34a inhibitor in cell model not only attenuated D-gal-induced the impairment of autophagy but also decreased the expression of DRP1 and mitofusin 2 (MFN2). Therefore, swimming training can delay brain aging of d-gal-induced aging rats through attenuating the impairment of miR-34a-mediated autophagy and abnormal mitochondrial dynamics, and miR-34a could be the novel therapeutic target for aging-related diseases such as AD. NEW & NOTEWORTHY In the present study, we have found that the upregulation of miR-34a is the hallmark of aging or aging-related diseases, which can result in dysfunctional autophagy and abnormal mitochondrial dynamics. In contrast, swimming intervention can delay the aging process by rescuing the impaired functional status of autophagy and abnormal mitochondrial dynamics via the suppression of miR-34a. Copyright © 2017 the American Physiological Society.

DYNSTALL: Subroutine package with a dynamic stall model

Energy Technology Data Exchange (ETDEWEB)

Bjoerck, Anders [Aeronautical Research Inst. of Sweden, Bromma (Sweden)

2001-03-01

A subroutine package, called DYNSTALL, for the calculation of 2D unsteady airfoil aerodynamics is described. The subroutines are written in FORTRAN. DYNSTALL is basically an implementation of the Beddoes-Leishman dynamic stall model. This model is a semi-empirical model for dynamic stall. It includes, however, also models for attached flow unsteady aerodynamics. It is complete in the sense that it treats attached flow as well as separated flow. Semi-empirical means that the model relies on empirically determined constants. Semi because the constants are constants in equations with some physical interpretation. It requires the input of 2D airfoil aerodynamic data via tables as function of angle of attack. The method is intended for use in an aeroelastic code with the aerodynamics solved by blade/element method. DYNSTALL was written to work for any 2D angles of attack relative to the airfoil, e.g. flow from the rear of an airfoil.

Inhibitory effects of aspirin-triggered resolvin D1 on spinal nociceptive processing in rat pain models.

Science.gov (United States)

Meesawatsom, Pongsatorn; Burston, James; Hathway, Gareth; Bennett, Andrew; Chapman, Victoria

2016-09-02

Harnessing the actions of the resolvin pathways has the potential for the treatment of a wide range of conditions associated with overt inflammatory signalling. Aspirin-triggered resolvin D1 (AT-RvD1) has robust analgesic effects in behavioural models of pain; however, the potential underlying spinal neurophysiological mechanisms contributing to these inhibitory effects in vivo are yet to be determined. This study investigated the acute effects of spinal AT-RvD1 on evoked responses of spinal neurones in vivo in a model of acute inflammatory pain and chronic osteoarthritic (OA) pain and the relevance of alterations in spinal gene expression to these neurophysiological effects. Pain behaviour was assessed in rats with established carrageenan-induced inflammatory or monosodium iodoacetate (MIA)-induced OA pain, and changes in spinal gene expression of resolvin receptors and relevant enzymatic pathways were examined. At timepoints of established pain behaviour, responses of deep dorsal horn wide dynamic range (WDR) neurones to transcutaneous electrical stimulation of the hind paw were recorded pre- and post direct spinal administration of AT-RvD1 (15 and 150 ng/50 μl). AT-RvD1 (15 ng/50 μl) significantly inhibited WDR neurone responses to electrical stimuli at C- (29 % inhibition) and Aδ-fibre (27 % inhibition) intensities. Both wind-up (53 %) and post-discharge (46 %) responses of WDR neurones in carrageenan-treated animals were significantly inhibited by AT-RvD1, compared to pre-drug response (p < 0.05). These effects were abolished by spinal pre-administration of a formyl peptide receptor 2 (FPR2/ALX) antagonist, butoxy carbonyl-Phe-Leu-Phe-Leu-Phe (BOC-2) (50 μg/50 μl). AT-RvD1 did not alter evoked WDR neurone responses in non-inflamed or MIA-treated rats. Electrophysiological effects in carrageenan-inflamed rats were accompanied by a significant increase in messenger RNA (mRNA) for chemerin (ChemR23) receptor and 5-lipoxygenase

Using Parameters of Dynamic Pulse Function for 3d Modeling in LOD3 Based on Random Textures

Science.gov (United States)

Alizadehashrafi, B.

2015-12-01

The pulse function (PF) is a technique based on procedural preprocessing system to generate a computerized virtual photo of the façade with in a fixed size square(Alizadehashrafi et al., 2009, Musliman et al., 2010). Dynamic Pulse Function (DPF) is an enhanced version of PF which can create the final photo, proportional to real geometry. This can avoid distortion while projecting the computerized photo on the generated 3D model(Alizadehashrafi and Rahman, 2013). The challenging issue that might be handled for having 3D model in LoD3 rather than LOD2, is the final aim that have been achieved in this paper. In the technique based DPF the geometries of the windows and doors are saved in an XML file schema which does not have any connections with the 3D model in LoD2 and CityGML format. In this research the parameters of Dynamic Pulse Functions are utilized via Ruby programming language in SketchUp Trimble to generate (exact position and deepness) the windows and doors automatically in LoD3 based on the same concept of DPF. The advantage of this technique is automatic generation of huge number of similar geometries e.g. windows by utilizing parameters of DPF along with defining entities and window layers. In case of converting the SKP file to CityGML via FME software or CityGML plugins the 3D model contains the semantic database about the entities and window layers which can connect the CityGML to MySQL(Alizadehashrafi and Baig, 2014). The concept behind DPF, is to use logical operations to project the texture on the background image which is dynamically proportional to real geometry. The process of projection is based on two vertical and horizontal dynamic pulses starting from upper-left corner of the background wall in down and right directions respectively based on image coordinate system. The logical one/zero on the intersections of two vertical and horizontal dynamic pulses projects/does not project the texture on the background image. It is possible to define

Monitoring conformational dynamics with solid-state R1ρ experiments

International Nuclear Information System (INIS)

Quinn, Caitlin M.; McDermott, Ann E.

2009-01-01

A new application of solid-state rotating frame (R 1ρ ) relaxation experiments to observe conformational dynamics is presented. Studies on a model compound, dimethyl sulfone (DMS), show that R 1ρ relaxation due to reorientation of a chemical shift anisotropy (CSA) tensor undergoing chemical exchange can be used to monitor slow-to-intermediate timescale conformational exchange processes. Control experiments used d 6 -DMS and alanine to confirm that the technique is monitoring reorientation of the CSA tensor rather than dipolar interactions or methyl group rotation. The application of this method to proteins could represent a new site-specific probe of conformational dynamics

Modelling atmospheric transport of α-hexachlorocyclohexane in the Northern Hemispherewith a 3-D dynamical model: DEHM-POP

Directory of Open Access Journals (Sweden)

K. M. Hansen

2004-01-01

Full Text Available The Danish Eulerian Hemispheric Model (DEHM is a 3-D dynamical atmospheric transport model originally developed to describe the atmospheric transport of sulphur into the Arctic. A new version of the model, DEHM-POP, developed to study the atmospheric transport and environmental fate of persistent organic pollutants (POPs is presented. During environmental cycling, POPs can be deposited and re-emitted several times before reaching a final destination. A description of the exchange processes between the land/ocean surfaces and the atmosphere is included in the model to account for this multi-hop transport. The α-isomer of the pesticide hexachlorocyclohexane (α-HCH is used as tracer in the model development. The structure of the model and processes included are described in detail. The results from a model simulation showing the atmospheric transport for the years 1991 to 1998 are presented and evaluated against measurements. The annual averaged atmospheric concentration of α-HCH for the 1990s is well described by the model; however, the shorter-term average concentration for most of the stations is not well captured. This indicates that the present simple surface description needs to be refined to get a better description of the air-surface exchange processes of POPs.

Modelling atmospheric transport of α-hexachlorocyclohexane in the Northern Hemispherewith a 3-D dynamical model: DEHM-POP

Science.gov (United States)

Hansen, K. M.; Christensen, J. H.; Brandt, J.; Frohn, L. M.; Geels, C.

2004-07-01

The Danish Eulerian Hemispheric Model (DEHM) is a 3-D dynamical atmospheric transport model originally developed to describe the atmospheric transport of sulphur into the Arctic. A new version of the model, DEHM-POP, developed to study the atmospheric transport and environmental fate of persistent organic pollutants (POPs) is presented. During environmental cycling, POPs can be deposited and re-emitted several times before reaching a final destination. A description of the exchange processes between the land/ocean surfaces and the atmosphere is included in the model to account for this multi-hop transport. The α-isomer of the pesticide hexachlorocyclohexane (α-HCH) is used as tracer in the model development. The structure of the model and processes included are described in detail. The results from a model simulation showing the atmospheric transport for the years 1991 to 1998 are presented and evaluated against measurements. The annual averaged atmospheric concentration of α-HCH for the 1990s is well described by the model; however, the shorter-term average concentration for most of the stations is not well captured. This indicates that the present simple surface description needs to be refined to get a better description of the air-surface exchange processes of POPs.

Quantifying Key Climate Parameter Uncertainties Using an Earth System Model with a Dynamic 3D Ocean

Science.gov (United States)

Olson, R.; Sriver, R. L.; Goes, M. P.; Urban, N.; Matthews, D.; Haran, M.; Keller, K.

2011-12-01

Climate projections hinge critically on uncertain climate model parameters such as climate sensitivity, vertical ocean diffusivity and anthropogenic sulfate aerosol forcings. Climate sensitivity is defined as the equilibrium global mean temperature response to a doubling of atmospheric CO2 concentrations. Vertical ocean diffusivity parameterizes sub-grid scale ocean vertical mixing processes. These parameters are typically estimated using Intermediate Complexity Earth System Models (EMICs) that lack a full 3D representation of the oceans, thereby neglecting the effects of mixing on ocean dynamics and meridional overturning. We improve on these studies by employing an EMIC with a dynamic 3D ocean model to estimate these parameters. We carry out historical climate simulations with the University of Victoria Earth System Climate Model (UVic ESCM) varying parameters that affect climate sensitivity, vertical ocean mixing, and effects of anthropogenic sulfate aerosols. We use a Bayesian approach whereby the likelihood of each parameter combination depends on how well the model simulates surface air temperature and upper ocean heat content. We use a Gaussian process emulator to interpolate the model output to an arbitrary parameter setting. We use Markov Chain Monte Carlo method to estimate the posterior probability distribution function (pdf) of these parameters. We explore the sensitivity of the results to prior assumptions about the parameters. In addition, we estimate the relative skill of different observations to constrain the parameters. We quantify the uncertainty in parameter estimates stemming from climate variability, model and observational errors. We explore the sensitivity of key decision-relevant climate projections to these parameters. We find that climate sensitivity and vertical ocean diffusivity estimates are consistent with previously published results. The climate sensitivity pdf is strongly affected by the prior assumptions, and by the scaling

Complex Price Dynamics in the Modified Kaldorian Model

Czech Academy of Sciences Publication Activity Database

Kodera, Jan; Van Tran, Q.; Vošvrda, Miloslav

2013-01-01

Roč. 22, č. 3 (2013), s. 358-384 ISSN 1210-0455 R&D Projects: GA ČR(CZ) GBP402/12/G097 Institutional support: RVO:67985556 Keywords : Priice dynamics, * numerical examples * two-equation model * four-equation model * nonlinear time series analysis Subject RIV: AH - Economics Impact factor: 0.208, year: 2013 http://library.utia.cas.cz/separaty/2013/E/kodera-model of price dynamics and chaos.pdf

Development of a hybrid 3-D hydrological model to simulate hillslopes and the regional unconfined aquifer system in Earth system models

Science.gov (United States)

Hazenberg, P.; Broxton, P. D.; Brunke, M.; Gochis, D.; Niu, G. Y.; Pelletier, J. D.; Troch, P. A. A.; Zeng, X.

2015-12-01

The terrestrial hydrological system, including surface and subsurface water, is an essential component of the Earth's climate system. Over the past few decades, land surface modelers have built one-dimensional (1D) models resolving the vertical flow of water through the soil column for use in Earth system models (ESMs). These models generally have a relatively coarse model grid size (~25-100 km) and only account for sub-grid lateral hydrological variations using simple parameterization schemes. At the same time, hydrologists have developed detailed high-resolution (~0.1-10 km grid size) three dimensional (3D) models and showed the importance of accounting for the vertical and lateral redistribution of surface and subsurface water on soil moisture, the surface energy balance and ecosystem dynamics on these smaller scales. However, computational constraints have limited the implementation of the high-resolution models for continental and global scale applications. The current work presents a hybrid-3D hydrological approach is presented, where the 1D vertical soil column model (available in many ESMs) is coupled with a high-resolution lateral flow model (h2D) to simulate subsurface flow and overland flow. H2D accounts for both local-scale hillslope and regional-scale unconfined aquifer responses (i.e. riparian zone and wetlands). This approach was shown to give comparable results as those obtained by an explicit 3D Richards model for the subsurface, but improves runtime efficiency considerably. The h3D approach is implemented for the Delaware river basin, where Noah-MP land surface model (LSM) is used to calculated vertical energy and water exchanges with the atmosphere using a 10km grid resolution. Noah-MP was coupled within the WRF-Hydro infrastructure with the lateral 1km grid resolution h2D model, for which the average depth-to-bedrock, hillslope width function and soil parameters were estimated from digital datasets. The ability of this h3D approach to simulate

Flood hazard assessment using 1D and 2D approaches

Science.gov (United States)

Petaccia, Gabriella; Costabile, Pierfranco; Macchione, Francesco; Natale, Luigi

2013-04-01

The EU flood risk Directive (Directive 2007/60/EC) prescribes risk assessment and mapping to develop flood risk management plans. Flood hazard mapping may be carried out with mathematical models able to determine flood-prone areas once realistic conditions (in terms of discharge or water levels) are imposed at the boundaries of the case study. The deterministic models are mainly based on shallow water equations expressed in their 1D or 2D formulation. The 1D approach is widely used, especially in technical studies, due to its relative simplicity, its computational efficiency and also because it requires topographical data not as expensive as the ones needed by 2D models. Even if in a great number of practical situations, such as modeling in-channel flows and not too wide floodplains, the 1D approach may provide results close to the prediction of a more sophisticated 2D model, it must be pointed out that the correct use of a 1D model in practical situations is more complex than it may seem. The main issues to be correctly modeled in a 1D approach are the definition of hydraulic structures such as bridges and buildings interacting with the flow and the treatment of the tributaries. Clearly all these aspects have to be taken into account also in the 2D modeling, but with fewer difficulties. The purpose of this paper is to show how the above cited issues can be described using a 1D or 2D unsteady flow modeling. In particular the Authors will show the devices that have to be implemented in 1D modeling to get reliable predictions of water levels and discharges comparable to the ones obtained using a 2D model. Attention will be focused on an actual river (Crati river) located in the South of Italy. This case study is quite complicated since it deals with the simulation of channeled flows, overbank flows, interactions with buildings, bridges and tributaries. Accurate techniques, intentionally developed by the Authors to take into account all these peculiarities in 1D and 2

Dynamic RSA: Examining parasympathetic regulatory dynamics via vector-autoregressive modeling of time-varying RSA and heart period.

Science.gov (United States)

Fisher, Aaron J; Reeves, Jonathan W; Chi, Cyrus

2016-07-01

Expanding on recently published methods, the current study presents an approach to estimating the dynamic, regulatory effect of the parasympathetic nervous system on heart period on a moment-to-moment basis. We estimated second-to-second variation in respiratory sinus arrhythmia (RSA) in order to estimate the contemporaneous and time-lagged relationships among RSA, interbeat interval (IBI), and respiration rate via vector autoregression. Moreover, we modeled these relationships at lags of 1 s to 10 s, in order to evaluate the optimal latency for estimating dynamic RSA effects. The IBI (t) on RSA (t-n) regression parameter was extracted from individual models as an operationalization of the regulatory effect of RSA on IBI-referred to as dynamic RSA (dRSA). Dynamic RSA positively correlated with standard averages of heart rate and negatively correlated with standard averages of RSA. We propose that dRSA reflects the active downregulation of heart period by the parasympathetic nervous system and thus represents a novel metric that provides incremental validity in the measurement of autonomic cardiac control-specifically, a method by which parasympathetic regulatory effects can be measured in process. © 2016 Society for Psychophysiological Research.

Dynamics model for real time diagnostics of TRIGA RC-1 system

International Nuclear Information System (INIS)

Gadomski, A.M.; Nanni, V.; Meo, G.B.

1986-01-01

This paper presents dynamics model of TRIGA RC-1 reactor system. The model is dedicated to the real-time early fault detection during a reactor operation in one week exploitation cycle. The algorithms are specially suited for real-time, long time and also accelerated simulations with assumed diagnostic oriented accuracy. The approximations, modular structure, numerical methods and validation are discussed. The elaborated model will be build in the TRIGA Supervisory System and TRIGA Diagnostic Simulator. (author)

Interacting gaps model, dynamics of order book, and stock-market fluctuations

Czech Academy of Sciences Publication Activity Database

Svorenčík, A.; Slanina, František

2007-01-01

Roč. 57, - (2007), s. 453-462 ISSN 1434-6028 R&D Projects: GA MŠk 1P04OCP10.001 Institutional research plan: CEZ:AV0Z10100520 Keywords : interacting gaps model * dynamics of order book * stock - market fluctuations Subject RIV: BE - Theoretical Physics Impact factor: 1.356, year: 2007

Dynamical Cognitive Models of Social Issues in Russia

Science.gov (United States)

Mitina, Olga; Abraham, Fred; Petrenko, Victor

We examine and model dynamics in three areas of social cognition: (1) political transformations within Russia, (2) evaluation of political trends in other countries by Russians, and (3) evaluation of Russian stereotypes concerning women. We try to represent consciousness as vectorfields and trajectories in a cognitive state space. We use psychosemantic techniques that allow definition of the state space and the systematic construction of these vectorfields and trajectories and their portrait from research data. Then we construct models to fit them, using multiple regression methods to obtain linear differential equations. These dynamical models of social cognition fit the data quite well. (1) The political transformations were modeled by a spiral repellor in a two-dimensional space of a democratic-totalitarian factor and social depression-optimism factor. (2) The evaluation of alien political trends included a flow away from a saddle toward more stable and moderate political regimes in a 2D space, of democratic-totalitarian and unstable-stable cognitive dimensions. (3) The gender study showed expectations (attractors) for more liberated, emancipated roles for women in the future.

Assessing 1D Atmospheric Solar Radiative Transfer Models: Interpretation and Handling of Unresolved Clouds.

Science.gov (United States)

Barker, H. W.; Stephens, G. L.; Partain, P. T.; Bergman, J. W.; Bonnel, B.; Campana, K.; Clothiaux, E. E.; Clough, S.; Cusack, S.; Delamere, J.; Edwards, J.; Evans, K. F.; Fouquart, Y.; Freidenreich, S.; Galin, V.; Hou, Y.; Kato, S.; Li, J.;  Mlawer, E.;  Morcrette, J.-J.;  O'Hirok, W.;  Räisänen, P.;  Ramaswamy, V.;  Ritter, B.;  Rozanov, E.;  Schlesinger, M.;  Shibata, K.;  Sporyshev, P.;  Sun, Z.;  Wendisch, M.;  Wood, N.;  Yang, F.

2003-08-01

The primary purpose of this study is to assess the performance of 1D solar radiative transfer codes that are used currently both for research and in weather and climate models. Emphasis is on interpretation and handling of unresolved clouds. Answers are sought to the following questions: (i) How well do 1D solar codes interpret and handle columns of information pertaining to partly cloudy atmospheres? (ii) Regardless of the adequacy of their assumptions about unresolved clouds, do 1D solar codes perform as intended?One clear-sky and two plane-parallel, homogeneous (PPH) overcast cloud cases serve to elucidate 1D model differences due to varying treatments of gaseous transmittances, cloud optical properties, and basic radiative transfer. The remaining four cases involve 3D distributions of cloud water and water vapor as simulated by cloud-resolving models. Results for 25 1D codes, which included two line-by-line (LBL) models (clear and overcast only) and four 3D Monte Carlo (MC) photon transport algorithms, were submitted by 22 groups. Benchmark, domain-averaged irradiance profiles were computed by the MC codes. For the clear and overcast cases, all MC estimates of top-of-atmosphere albedo, atmospheric absorptance, and surface absorptance agree with one of the LBL codes to within ±2%. Most 1D codes underestimate atmospheric absorptance by typically 15-25 W m-2 at overhead sun for the standard tropical atmosphere regardless of clouds.Depending on assumptions about unresolved clouds, the 1D codes were partitioned into four genres: (i) horizontal variability, (ii) exact overlap of PPH clouds, (iii) maximum/random overlap of PPH clouds, and (iv) random overlap of PPH clouds. A single MC code was used to establish conditional benchmarks applicable to each genre, and all MC codes were used to establish the full 3D benchmarks. There is a tendency for 1D codes to cluster near their respective conditional benchmarks, though intragenre variances typically exceed those for

Dynamic Frames Based Generation of 3D Scenes and Applications

Directory of Open Access Journals (Sweden)

Danijel Radošević

2015-05-01

Full Text Available Modern graphic/programming tools like Unity enables the possibility of creating 3D scenes as well as making 3D scene based program applications, including full physical model, motion, sounds, lightning effects etc. This paper deals with the usage of dynamic frames based generator in the automatic generation of 3D scene and related source code. The suggested model enables the possibility to specify features of the 3D scene in a form of textual specification, as well as exporting such features from a 3D tool. This approach enables higher level of code generation flexibility and the reusability of the main code and scene artifacts in a form of textual templates. An example of the generated application is presented and discussed.

Dynamic compensation temperature in the kinetic spin-1 Ising model in an oscillating external magnetic field on alternate layers of a hexagonal lattice

International Nuclear Information System (INIS)

Temizer, Umuet; Keskin, Mustafa; Canko, Osman

2009-01-01

The dynamic behavior of a two-sublattice spin-1 Ising model with a crystal-field interaction (D) in the presence of a time-varying magnetic field on a hexagonal lattice is studied by using the Glauber-type stochastic dynamics. The lattice is formed by alternate layers of spins σ=1 and S=1. For this spin arrangement, any spin at one lattice site has two nearest-neighbor spins on the same sublattice, and four on the other sublattice. The intersublattice interaction is antiferromagnetic. We employ the Glauber transition rates to construct the mean-field dynamical equations. Firstly, we study time variations of the average magnetizations in order to find the phases in the system, and the temperature dependence of the average magnetizations in a period, which is also called the dynamic magnetizations, to obtain the dynamic phase transition (DPT) points as well as to characterize the nature (continuous and discontinuous) of transitions. Then, the behavior of the total dynamic magnetization as a function of the temperature is investigated to find the types of the compensation behavior. Dynamic phase diagrams are calculated for both DPT points and dynamic compensation effect. Phase diagrams contain the paramagnetic (p) and antiferromagnetic (af) phases, the p+af and nm+p mixed phases, nm is the non-magnetic phase, and the compensation temperature or the L-type behavior that strongly depend on the interaction parameters. For D 0 >3.8275, H 0 is the magnetic field amplitude, the compensation effect does not appear in the system.

Dynamic compensation temperature in the kinetic spin-1 Ising model in an oscillating external magnetic field on alternate layers of a hexagonal lattice

Energy Technology Data Exchange (ETDEWEB)

Temizer, Umuet [Department of Physics, Bozok University, 66100 Yozgat (Turkey); Keskin, Mustafa [Department of Physics, Erciyes University, 38039 Kayseri (Turkey)], E-mail: keskin@erciyes.edu.tr; Canko, Osman [Department of Physics, Erciyes University, 38039 Kayseri (Turkey)

2009-10-15

The dynamic behavior of a two-sublattice spin-1 Ising model with a crystal-field interaction (D) in the presence of a time-varying magnetic field on a hexagonal lattice is studied by using the Glauber-type stochastic dynamics. The lattice is formed by alternate layers of spins {sigma}=1 and S=1. For this spin arrangement, any spin at one lattice site has two nearest-neighbor spins on the same sublattice, and four on the other sublattice. The intersublattice interaction is antiferromagnetic. We employ the Glauber transition rates to construct the mean-field dynamical equations. Firstly, we study time variations of the average magnetizations in order to find the phases in the system, and the temperature dependence of the average magnetizations in a period, which is also called the dynamic magnetizations, to obtain the dynamic phase transition (DPT) points as well as to characterize the nature (continuous and discontinuous) of transitions. Then, the behavior of the total dynamic magnetization as a function of the temperature is investigated to find the types of the compensation behavior. Dynamic phase diagrams are calculated for both DPT points and dynamic compensation effect. Phase diagrams contain the paramagnetic (p) and antiferromagnetic (af) phases, the p+af and nm+p mixed phases, nm is the non-magnetic phase, and the compensation temperature or the L-type behavior that strongly depend on the interaction parameters. For D<2.835 and H{sub 0}>3.8275, H{sub 0} is the magnetic field amplitude, the compensation effect does not appear in the system.

An Idealized Model of Plant and Soil Dynamics

Science.gov (United States)

Burg, David; Malkinson, Dan; Wittenberg, Lea

2014-05-01

Following wildfire events the landscape commonly becomes denuded of vegetation cover, resulting in systems prone to soil loss and degradation. In this context soil dynamics are an intricate process balanced between pedogenesis, which is a relatively slow process and erosion which depends on many inert (e.g. soil texture, slope, precipitation and wind) and biological factors such as vegetation properties, grazing intensity, and human disturbance. We develop a simple homogenous, spatially implicit, theoretical model of the global dynamics of the interactions between vegetation and soil using a system of two nonlinear differential equations describing this interdependence, assuming a double feedback between them - plants control erosion and soil availability facilitates plants growth: ( ) dV- -K-- dt = rV K - 1+ aS - V (1) dS-= σ - ɛSe-cT dt (2) where V and S represent vegetation cover and soil availability, respectively. Vegetation growth is similar to the classical logistic model with a growth rate of r(yr1), however, the "carrying capacity" (K) is dependent on soil availability (a1 is the amount of soil where V is reduced by half). Soil influxes at a constant rate σ(mm×yr1) and is eroded at a constant rategɛ (yr-1), while vegetation abates this process modeled as a decreasing exponent as the effectiveness of vegetation in reducing soil erosion (c). Parameter values were chosen from a variable range found in the literature: r=0.01 yr1, K=75%, a1=1, σ=1 mm×yr1, ɛ=0.1 yr1, c=0.08. Complex properties emerge from this model. At certain parameter values (cK≤4) the model predicts one of two steady states - full recovery of vegetation cover or a degraded barren system. However, at certain boundary conditions (cK>4 and Λ1 ≤ σ/ɛ ≤ Λ2, see Article for terms of Λ1 and Λ2) bistability may be observed. We also show that erosion seems to be the determining factor in this system, and we identify the threshold values from which beyond the systems become unstable

Modelling a 3D structure for EgDf1 from shape Echinococcus granulosus: putative epitopes, phosphorylation motifs and ligand

Science.gov (United States)

Paulino, M.; Esteves, A.; Vega, M.; Tabares, G.; Ehrlich, R.; Tapia, O.

1998-07-01

EgDf1 is a developmentally regulated protein from the parasite Echinococcus granulosus related to a family of hydrophobic ligand binding proteins. This protein could play a crucial role during the parasite life cycle development since this organism is unable to synthetize most of their own lipids de novo. Furthermore, it has been shown that two related protein from other parasitic platyhelminths (Fh15 from Fasciola hepatica and Sm14 from Schistosoma mansoni) are able to confer protective inmunity against experimental infection in animal models. A three-dimensional structure would help establishing structure/function relationships on a knowledge based manner. 3D structures for EgDf1 protein were modelled by using myelin P2 (mP2) and intestine fatty acid binding protein (I-FABP) as templates. Molecular dynamics techniques were used to validate the models. Template mP2 yielded the best 3D structure for EgDf1. Palmitic and oleic acids were docked inside EgDf1. The present theoretical results suggest definite location in the secondary structure of the epitopic regions, consensus phosphorylation motifs and oleic acid as a good ligand candidate to EgDf1. This protein might well be involved in the process of supplying hydrophobic metabolites for membrane biosynthesis and for signaling pathways.

PLUME-MoM 1.0: a new 1-D model of volcanic plumes based on the method of moments

Science.gov (United States)

de'Michieli Vitturi, M.; Neri, A.; Barsotti, S.

2015-05-01

In this paper a new mathematical model for volcanic plumes, named PlumeMoM, is presented. The model describes the steady-state 1-D dynamics of the plume in a 3-D coordinate system, accounting for continuous variability in particle distribution of the pyroclastic mixture ejected at the vent. Volcanic plumes are composed of pyroclastic particles of many different sizes ranging from a few microns up to several centimeters and more. Proper description of such a multiparticle nature is crucial when quantifying changes in grain-size distribution along the plume and, therefore, for better characterization of source conditions of ash dispersal models. The new model is based on the method of moments, which allows description of the pyroclastic mixture dynamics not only in the spatial domain but also in the space of properties of the continuous size-distribution of the particles. This is achieved by formulation of fundamental transport equations for the multiparticle mixture with respect to the different moments of the grain-size distribution. Different formulations, in terms of the distribution of the particle number, as well as of the mass distribution expressed in terms of the Krumbein log scale, are also derived. Comparison between the new moments-based formulation and the classical approach, based on the discretization of the mixture in N discrete phases, shows that the new model allows the same results to be obtained with a significantly lower computational cost (particularly when a large number of discrete phases is adopted). Application of the new model, coupled with uncertainty quantification and global sensitivity analyses, enables investigation of the response of four key output variables (mean and standard deviation (SD) of the grain-size distribution at the top of the plume, plume height and amount of mass lost by the plume during the ascent) to changes in the main input parameters (mean and SD) characterizing the pyroclastic mixture at the base of the plume

Compressive Loads on the Lumbar Spine During Lifting: 4D WATBAK versus Inverse Dynamics Calculations

Directory of Open Access Journals (Sweden)

M. H. Cole

2005-01-01

Full Text Available Numerous two- and three-dimensional biomechanical models exist for the purpose of assessing the stresses placed on the lumbar spine during the performance of a manual material handling task. More recently, researchers have utilised their knowledge to develop specific computer-based models that can be applied in an occupational setting; an example of which is 4D WATBAK. The model used by 4D WATBAK bases its predications on static calculations and it is assumed that these static loads reasonably depict the actual dynamic loads acting on the lumbar spine. Consequently, it was the purpose of this research to assess the agreement between the static predictions made by 4D WATBAK and those from a comparable dynamic model. Six individuals were asked to perform a series of five lifting tasks, which ranged from lifting 2.5 kg to 22.5 kg and were designed to replicate the lifting component of the Work Capacity Assessment Test used within Australia. A single perpendicularly placed video camera was used to film each performance in the sagittal plane. The resultant two-dimensional kinematic data were input into the 4D WATBAK software and a dynamic biomechanical model to quantify the compression forces acting at the L4/L5 intervertebral joint. Results of this study indicated that as the mass of the load increased from 2.5 kg to 22.5 kg, the static compression forces calculated by 4D WATBAK became increasingly less than those calculated using the dynamic model (mean difference ranged from 22.0% for 2.5 kg to 42.9% for 22.5 kg. This study suggested that, for research purposes, a validated three-dimensional dynamic model should be employed when a task becomes complex and when a more accurate indication of spinal compression or shear force is required. Additionally, although it is clear that 4D WATBAK is particularly suited to industrial applications, it is suggested that the limitations of such modelling tools be carefully considered when task-risk and employee

Technique development of 3D dynamic CS-EPSI for hyperpolarized 13 C pyruvate MR molecular imaging of human prostate cancer.

Science.gov (United States)

Chen, Hsin-Yu; Larson, Peder E Z; Gordon, Jeremy W; Bok, Robert A; Ferrone, Marcus; van Criekinge, Mark; Carvajal, Lucas; Cao, Peng; Pauly, John M; Kerr, Adam B; Park, Ilwoo; Slater, James B; Nelson, Sarah J; Munster, Pamela N; Aggarwal, Rahul; Kurhanewicz, John; Vigneron, Daniel B

2018-03-25

The purpose of this study was to develop a new 3D dynamic carbon-13 compressed sensing echoplanar spectroscopic imaging (EPSI) MR sequence and test it in phantoms, animal models, and then in prostate cancer patients to image the metabolic conversion of hyperpolarized [1- 13 C]pyruvate to [1- 13 C]lactate with whole gland coverage at high spatial and temporal resolution. A 3D dynamic compressed sensing (CS)-EPSI sequence with spectral-spatial excitation was designed to meet the required spatial coverage, time and spatial resolution, and RF limitations of the 3T MR scanner for its clinical translation for prostate cancer patient imaging. After phantom testing, animal studies were performed in rats and transgenic mice with prostate cancers. For patient studies, a GE SPINlab polarizer (GE Healthcare, Waukesha, WI) was used to produce hyperpolarized sterile GMP [1- 13 C]pyruvate. 3D dynamic 13 C CS-EPSI data were acquired starting 5 s after injection throughout the gland with a spatial resolution of 0.5 cm 3 , 18 time frames, 2-s temporal resolution, and 36 s total acquisition time. Through preclinical testing, the 3D CS-EPSI sequence developed in this project was shown to provide the desired spectral, temporal, and spatial 5D HP 13 C MR data. In human studies, the 3D dynamic HP CS-EPSI approach provided first-ever simultaneously volumetric and dynamic images of the LDH-catalyzed conversion of [1- 13 C]pyruvate to [1- 13 C]lactate in a biopsy-proven prostate cancer patient with full gland coverage. The results demonstrate the feasibility to characterize prostate cancer metabolism in animals, and now patients using this new 3D dynamic HP MR technique to measure k PL , the kinetic rate constant of [1- 13 C]pyruvate to [1- 13 C]lactate conversion. © 2018 International Society for Magnetic Resonance in Medicine.

Emplacement dynamics of phonolite magma into maar-diatreme structures - Correlation of field, thermal modeling and AMS analogue modeling data

Czech Academy of Sciences Publication Activity Database

Závada, Prokop; Dědeček, Petr; Mach, K.; Lexa, O.; Potužák, M.

2011-01-01

Roč. 201, č. 1-4 (2011), s. 210-226 ISSN 0377-0273 R&D Projects: GA MŠk ME09011; GA AV ČR KJB301110703 Institutional research plan: CEZ:AV0Z30120515 Keywords : phonolite * maar-diatreme * magma emplacement dynamics * thermal modeling * AMS analogue modeling * cavitation Subject RIV: DB - Geology ; Mineralogy Impact factor: 1.978, year: 2011

EFDC1D - A ONE DIMENSIONAL HYDRODYNAMIC AND SEDIMENT TRANSPORT MODEL FOR RIVER AND STREAM NETWORKS: MODEL THEORY AND USERS GUIDE

Science.gov (United States)

This technical report describes the new one-dimensional (1D) hydrodynamic and sediment transport model EFDC1D. This model that can be applied to stream networks. The model code and two sample data sets are included on the distribution CD. EFDC1D can simulate bi-directional unstea...

Fluid dynamics of moving fish in a two-dimensional multiparticle collision dynamics model in 2D

NARCIS (Netherlands)

Reid, D.A.P.; Hildenbrandt, H.; Padding, J.T.; Hemelrijk, C.K.

2012-01-01

The fluid dynamics of animal locomotion, such as that of an undulating fish, are of great interest to both biologists and engineers. However, experimentally studying these fluid dynamics is difficult and time consuming. Model studies can be of great help because of their simpler and more detailed

Modeling 3D Dynamic Rupture on Arbitrarily-Shaped faults by Boundary-Conforming Finite Difference Method

Science.gov (United States)

Zhu, D.; Zhu, H.; Luo, Y.; Chen, X.

2008-12-01

We use a new finite difference method (FDM) and the slip-weakening law to model the rupture dynamics of a non-planar fault embedded in a 3-D elastic media with free surface. The new FDM, based on boundary- conforming grid, sets up the mapping equations between the curvilinear coordinate and the Cartesian coordinate and transforms irregular physical space to regular computational space; it also employs a higher- order non-staggered DRP/opt MacCormack scheme which is of low dispersion and low dissipation so that the high accuracy and stability of our rupture modeling are guaranteed. Compared with the previous methods, not only we can compute the spontaneous rupture of an arbitrarily shaped fault, but also can model the influence of the surface topography on the rupture process of earthquake. In order to verify the feasibility of this method, we compared our results and other previous results, and found out they matched perfectly. Thanks to the boundary-conforming FDM, problems such as dynamic rupture with arbitrary dip, strike and rake over an arbitrary curved plane can be handled; and supershear or subshear rupture can be simulated with different parameters such as the initial stresses and the critical slip displacement Dc. Besides, our rupture modeling is economical to be implemented owing to its high efficiency and does not suffer from displacement leakage. With the help of inversion data of rupture by field observations, this method is convenient to model rupture processes and seismograms of natural earthquakes.

Comparison of hemodynamics of intracranial aneurysms between MR fluid dynamics using 3D cine phase-contrast MRI and MR-based computational fluid dynamics

International Nuclear Information System (INIS)

Isoda, Haruo; Sakahara, Harumi; Ohkura, Yasuhide; Kosugi, Takashi; Hirano, Masaya; Alley, Marcus T.; Bammer, Roland; Pelc, Norbert J.; Namba, Hiroki

2010-01-01

Hemodynamics is thought to play a very important role in the initiation, growth, and rupture of intracranial aneurysms. The purpose of our study was to compare hemodynamics of intracranial aneurysms of MR fluid dynamics (MRFD) using 3D cine PC MR imaging (4D-Flow) at 1.5 T and MR-based computational fluid dynamics (CFD). 4D-Flow was performed for five intracranial aneurysms by a 1.5 T MR scanner. 3D TOF MR angiography was performed for geometric information. The blood flow in the aneurysms was modeled using CFD simulation based on the finite element method. We used MR angiographic data as the vascular models and MR flow information as boundary conditions in CFD. 3D velocity vector fields, 3D streamlines, shearing velocity maps, wall shear stress (WSS) distribution maps and oscillatory shear index (OSI) distribution maps were obtained by MRFD and CFD and were compared. There was a moderate to high degree of correlation in 3D velocity vector fields and a low to moderate degree of correlation in WSS of aneurysms between MRFD and CFD using regression analysis. The patterns of 3D streamlines were similar between MRFD and CFD. The small and rotating shearing velocities and higher OSI were observed at the top of the spiral flow in the aneurysms. The pattern and location of shearing velocity in MRFD and CFD were similar. The location of high oscillatory shear index obtained by MRFD was near to that obtained by CFD. MRFD and CFD of intracranial aneurysms correlated fairly well. (orig.)

A system dynamics model for stock and flow of tritium in fusion power plant

Energy Technology Data Exchange (ETDEWEB)

Kasada, Ryuta, E-mail: r-kasada@iae.kyoto-u.ac.jp [Institute of Advanced Energy, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan); Kwon, Saerom [Institute of Advanced Energy, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan); Japan Atomic Energy Agency, Tokai-mura, Naka-gun, Ibaraki-ken 319-1195 (Japan); Konishi, Satoshi [Institute of Advanced Energy, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan); Sakamoto, Yoshiteru; Yamanishi, Toshihiko; Tobita, Kenji [Japan Atomic Energy Agency, Rokkasho-mura, Kamikita-gun, Aomori-ken 039-3212 (Japan)

2015-10-15

Highlights: • System dynamics model of tritium fuel cycle was developed for analyzing stock and flow of tritium in fusion power plants. • Sensitivity of tritium build-up to breeding ratio parameters has been assessed to two plant concepts having 3 GW and 1.5 GW fusion power. • D-D start-up absolutely without initial loading of tritium is possible for both of the 3 GW and 1.5 GW fusion power plant concepts. • Excess stock of tritium is generated by the steady state operation with the value of tritium breeding ratio over unity. - Abstract: In order to analyze self-efficiency of tritium fuel cycle (TFC) and share the systems thinking of TFC among researchers and engineers in the vast area of fusion reactor technology, we develop a system dynamics (SD) TFC model using a commercial software STELLA. The SD-TFC model is illustrated as a pipe diagram which consists of tritium stocks, such as plasma, fuel clean up, isotope separation, fueling with storage and blanket, and pipes connecting among them. By using this model, we survey a possibility of D-D start-up without initial loading of tritium on two kinds of fusion plant having different plasma parameters. The D-D start-up scenario can reduce the necessity of initial loading of tritium through the production in plasma by D-D reaction and in breeding blanket by D-D neutron. The model is also used for considering operation scenario to avoid excess stock of tritium which must be produced at tritium breeding ratio over unity.

2D and 3D Models of Convective Turbulence and Oscillations in Intermediate-Mass Main-Sequence Stars

Science.gov (United States)

Guzik, Joyce Ann; Morgan, Taylor H.; Nelson, Nicholas J.; Lovekin, Catherine; Kitiashvili, Irina N.; Mansour, Nagi N.; Kosovichev, Alexander

2015-08-01

We present multidimensional modeling of convection and oscillations in main-sequence stars somewhat more massive than the sun, using three separate approaches: 1) Applying the spherical 3D MHD ASH (Anelastic Spherical Harmonics) code to simulate the core convection and radiative zone. Our goal is to determine whether core convection can excite low-frequency gravity modes, and thereby explain the presence of low frequencies for some hybrid gamma Dor/delta Sct variables for which the envelope convection zone is too shallow for the convective blocking mechanism to drive g modes; 2) Using the 3D planar ‘StellarBox’ radiation hydrodynamics code to model the envelope convection zone and part of the radiative zone. Our goals are to examine the interaction of stellar pulsations with turbulent convection in the envelope, excitation of acoustic modes, and the role of convective overshooting; 3) Applying the ROTORC 2D stellar evolution and dynamics code to calculate evolution with a variety of initial rotation rates and extents of core convective overshooting. The nonradial adiabatic pulsation frequencies of these nonspherical models will be calculated using the 2D pulsation code NRO of Clement. We will present new insights into gamma Dor and delta Sct pulsations gained by multidimensional modeling compared to 1D model expectations.

Analytic Intermodel Consistent Modeling of Volumetric Human Lung Dynamics.

Science.gov (United States)

Ilegbusi, Olusegun; Seyfi, Behnaz; Neylon, John; Santhanam, Anand P

2015-10-01

Human lung undergoes breathing-induced deformation in the form of inhalation and exhalation. Modeling the dynamics is numerically complicated by the lack of information on lung elastic behavior and fluid-structure interactions between air and the tissue. A mathematical method is developed to integrate deformation results from a deformable image registration (DIR) and physics-based modeling approaches in order to represent consistent volumetric lung dynamics. The computational fluid dynamics (CFD) simulation assumes the lung is a poro-elastic medium with spatially distributed elastic property. Simulation is performed on a 3D lung geometry reconstructed from four-dimensional computed tomography (4DCT) dataset of a human subject. The heterogeneous Young's modulus (YM) is estimated from a linear elastic deformation model with the same lung geometry and 4D lung DIR. The deformation obtained from the CFD is then coupled with the displacement obtained from the 4D lung DIR by means of the Tikhonov regularization (TR) algorithm. The numerical results include 4DCT registration, CFD, and optimal displacement data which collectively provide consistent estimate of the volumetric lung dynamics. The fusion method is validated by comparing the optimal displacement with the results obtained from the 4DCT registration.

Dynamic modelling of n-of-1 data: powerful and flexible data analytics applied to individualised studies.

Science.gov (United States)

Vieira, Rute; McDonald, Suzanne; Araújo-Soares, Vera; Sniehotta, Falko F; Henderson, Robin

2017-09-01

N-of-1 studies are based on repeated observations within an individual or unit over time and are acknowledged as an important research method for generating scientific evidence about the health or behaviour of an individual. Statistical analyses of n-of-1 data require accurate modelling of the outcome while accounting for its distribution, time-related trend and error structures (e.g., autocorrelation) as well as reporting readily usable contextualised effect sizes for decision-making. A number of statistical approaches have been documented but no consensus exists on which method is most appropriate for which type of n-of-1 design. We discuss the statistical considerations for analysing n-of-1 studies and briefly review some currently used methodologies. We describe dynamic regression modelling as a flexible and powerful approach, adaptable to different types of outcomes and capable of dealing with the different challenges inherent to n-of-1 statistical modelling. Dynamic modelling borrows ideas from longitudinal and event history methodologies which explicitly incorporate the role of time and the influence of past on future. We also present an illustrative example of the use of dynamic regression on monitoring physical activity during the retirement transition. Dynamic modelling has the potential to expand researchers' access to robust and user-friendly statistical methods for individualised studies.

Co-targeting of Tiam1/Rac1 and Notch ameliorates chemoresistance against doxorubicin in a biomimetic 3D lymphoma model

Science.gov (United States)

Ikram, Muhammad; Lim, Yeseon; Baek, Sun-Yong; Jin, Songwan; Jeong, Young Hun; Kwak, Jong-Young; Yoon, Sik

2018-01-01

Lymphoma is a heterogeneous disease with a highly variable clinical course and prognosis. Improving the prognosis for patients with relapsed and treatment-resistant lymphoma remains challenging. Current in vitro drug testing models based on 2D cell culture lack natural tissue-like structural organization and result in disappointing clinical outcomes. The development of efficient drug testing models using 3D cell culture that more accurately reflects in vivo behaviors is vital. Our aim was to establish an in vitro 3D lymphoma model that can imitate the in vivo 3D lymphoma microenvironment. Using this model, we explored strategies to enhance chemosensitivity to doxorubicin, an important chemotherapeutic drug widely used for the treatment of hematological malignancies. Lymphoma cells grown in this model exhibited excellent biomimetic properties compared to conventional 2D culture including (1) enhanced chemotherapy resistance, (2) suppressed rate of apoptosis, (3) upregulated expression of drug resistance genes (MDR1, MRP1, BCRP and HIF-1α), (4) elevated levels of tumor aggressiveness factors including Notch (Notch-1, -2, -3, and -4) and its downstream molecules (Hes-1 and Hey-1), VEGF and MMPs (MMP-2 and MMP-9), and (5) enrichment of a lymphoma stem cell population. Tiam1, a potential biomarker of tumor progression, metastasis, and chemoresistance, was activated in our 3D lymphoma model. Remarkably, we identified two synergistic therapeutic oncotargets, Tiam1 and Notch, as a strategy to combat resistance against doxorubicin in EL4 T and A20 B lymphoma. Therefore, our data suggest that our 3D lymphoma model is a promising in vitro research platform for studying lymphoma biology and therapeutic approaches. PMID:29416753

Modeling dynamic stall on wind turbine blades under rotationally augmented flow fields

Energy Technology Data Exchange (ETDEWEB)

Guntur, S. [National Renewable Energy Laboratory (NREL), Golden, CO (United States); Schreck, S. [National Renewable Energy Laboratory (NREL), Golden, CO (United States); Sorensen, N. N. [Technical Univ. of Denmark, Lyngby (Denmark); Bergami, L. [Technical Univ. of Denmark, Lyngby (Denmark)

2015-04-22

It is well known that airfoils under unsteady flow conditions with a periodically varying angle of attack exhibit aerodynamic characteristics different from those under steady flow conditions, a phenomenon commonly known as dynamic stall. It is also well known that the steady aerodynamic characteristics of airfoils in the inboard region of a rotating blade differ from those under steady two-dimensional (2D) flow conditions, a phenomenon commonly known as rotational augmentation. This paper presents an investigation of these two phenomena together in the inboard parts of wind turbine blades. This analysis is carried out using data from three sources: (1) the National Renewable Energy Laboratory’s Unsteady Aerodynamics Experiment Phase VI experimental data, including constant as well as continuously pitching blade conditions during axial operation, (2) data from unsteady Delayed Detached Eddy Simulations (DDES) carried out using the Technical University of Denmark’s in-house flow solver Ellipsys3D, and (3) data from a simplified model based on the blade element momentum method with a dynamic stall subroutine that uses rotationally augmented steady-state polars obtained from steady Phase VI experimental sequences, instead of the traditional 2D nonrotating data. The aim of this work is twofold. First, the blade loads estimated by the DDES simulations are compared to three select cases of the N sequence experimental data, which serves as a validation of the DDES method. Results show reasonable agreement between the two data in two out of three cases studied. Second, the dynamic time series of the lift and the moment polars obtained from the experiments are compared to those from the dynamic stall subroutine that uses the rotationally augmented steady polars. This allowed the differences between the stall phenomenon on the inboard parts of harmonically pitching blades on a rotating wind turbine and the classic dynamic stall representation in 2D flow to be

Prediction of car cabin environment by means of 1D and 3D cabin model

Science.gov (United States)

Fišer, J.; Pokorný, J.; Jícha, M.

2012-04-01

Thermal comfort and also reduction of energy requirements of air-conditioning system in vehicle cabins are currently very intensively investigated and up-to-date issues. The article deals with two approaches of modelling of car cabin environment; the first model was created in simulation language Modelica (typical 1D approach without cabin geometry) and the second one was created in specialized software Theseus-FE (3D approach with cabin geometry). Performance and capabilities of this tools are demonstrated on the example of the car cabin and the results from simulations are compared with the results from the real car cabin climate chamber measurements.

On D-brane dynamics and moduli stabilization

Science.gov (United States)

Kitazawa, Noriaki

2017-09-01

We discuss the effect of the dynamics of D-branes on moduli stabilization in type IIB string theory compactifications, with reference to a concrete toy model of T6/Z 3 orientifold compactification with fractional D3-branes and anti-D3-branes at orbifold fixed points. The resulting attractive forces between anti-D3-branes and D3-branes, together with the repulsive forces between anti-D3-branes and O3-planes, can affect the stability of the compact space. There are no complex structure moduli in T6/Z 3 orientifold, which should thus capture some generic features of more general settings where all complex structure moduli are stabilized by three-form fluxes. The simultaneous presence of branes and anti-branes brings along the breaking of supersymmetry. Non-BPS combinations of this type are typical of “brane supersymmetry breaking” and are a necessary ingredient in the KKLT scenario for stabilizing the remaining Kähler moduli. The conclusion of our analysis is that, while mutual D-brane interactions sometimes help Kähler moduli stabilization, this is not always the case.

BETHSY 9.1b Test Calculation with TRACE Using 3D Vessel Component

International Nuclear Information System (INIS)

Berar, O.; Prosek, A.

2012-01-01

Recently, several advanced multidimensional computational tools for simulating reactor system behaviour during real and hypothetical transient scenarios were developed. One of such advanced, best-estimate reactor systems codes is TRAC/RELAP Advanced Computational Engine (TRACE), developed by the U.S. Nuclear Regulatory Commission. The advanced TRACE comes with a graphical user interface called SNAP (Symbolic Nuclear Analysis Package). It is intended for pre- and post-processing, running codes, RELAP5 to TRACE input deck conversion, input deck database generation etc. The TRACE code is still not fully development and it will have all the capabilities of RELAP5. The purpose of the present study was therefore to assess the 3D capability of the TRACE on BETHSY 9.1b test. The TRACE input deck was semi-converted (using SNAP and manual corrections) from the RELAP5 input deck. The 3D fluid dynamics within reactor vessel was modelled and compared to 1D fluid dynamics. The 3D calculation was compared both to TRACE 1D calculation and RELAP5 calculation. Namely, the geometry used in TRACE is basically the same, what gives very good basis for the comparison of the codes. The only exception is 3D reactor vessel model in case of TRACE 3D calculation. The TRACE V5.0 Patch 1 and RELAP5/MOD3.3 Patch 4 were used for calculations. The BETHSY 9.1b test (International Standard Problem no. 27 or ISP-27) was 5.08 cm equivalent diameter cold leg break without high pressure safety injection and with delayed ultimate procedure. BETHSY facility was a 3-loop replica of a 900 MWe FRAMATOME pressurized water reactor. For better presentation of the calculated physical phenomena and processes, an animation model using SNAP was developed. In general, the TRACE 3D code calculation is in good agreement with the BETHSY 9.1b test. The TRACE 3D calculation results are as good as or better than the RELAP5 calculated results. Also, the TRACE 3D calculation is not significantly different from TRACE 1D

Critical behavior in dome D = 1 large-N matrix models

International Nuclear Information System (INIS)

Das, S.R.; Dhar, A.; Sengupta, A.M.; Wadia, D.R.

1990-01-01

The authors study the critical behavior in D = 1 large-N matrix models. The authors also look at the subleading terms in susceptibility in order to find out the dimensions of some of the operators in the theory

Activated sludge models ASM1, ASM2, ASM2d and ASM3

DEFF Research Database (Denmark)

Henze, Mogens; Gujer, W.; Mino, T.

This book has been produced to give a total overview of the Activated Sludge Model (ASM) family at the start of 2000 and to give the reader easy access to the different models in their original versions. It thus presents ASM1, ASM2, ASM2d and ASM3 together for the first time.Modelling of activated...... sludge processes has become a common part of the design and operation of wastewater treatment plants. Today models are being used in design, control, teaching and research.ContentsASM3: Introduction, Comparison of ASM1 and ASM3, ASM3: Definition of compounds in the model, ASM3: Definition of processes...... in the Model, ASM3: Stoichiometry, ASM3: Kinetics, Limitations of ASM3, Aspects of application of ASM3, ASM3C: A Carbon based model, Conclusion ASM 2d: Introduction, Conceptual Approach, ASM 2d, Typical Wastewater Characteristics and Kinetic and Stoichiometric Constants, Limitations, Conclusion ASM 2...

Improving 1D Stellar Models with 3D Atmospheres

Science.gov (United States)

Rørsted Mosumgaard, Jakob; Silva Aguirre, Víctor; Weiss, Achim; Christensen-Dalsgaard, Jørgen; Trampedach, Regner

2017-10-01

Stellar evolution codes play a major role in present-day astrophysics, yet they share common issues. In this work we seek to remedy some of those by the use of results from realistic and highly detailed 3D hydrodynamical simulations of stellar atmospheres. We have implemented a new temperature stratification extracted directly from the 3D simulations into the Garching Stellar Evolution Code to replace the simplified atmosphere normally used. Secondly, we have implemented the use of a variable mixing-length parameter, which changes as a function of the stellar surface gravity and temperature - also derived from the 3D simulations. Furthermore, to make our models consistent, we have calculated new opacity tables to match the atmospheric simulations. Here, we present the modified code and initial results on stellar evolution using it.

Weyl and transverse diffeomorphism invariant spin-2 models in D = 2 + 1

International Nuclear Information System (INIS)

Dalmazi, Denis; Mendonca, E.L.; Santos, A.L.R. dos; Ghosh, Subir

2017-01-01

There are two covariant descriptions of massless spin-2 particles in D = 3 + 1 via a symmetric rank-2 tensor: the linearized Einstein-Hilbert (LEH) theory and the Weyl plus transverse diffeomorphism (WTDIFF) invariant model. From the LEH theory one can obtain the linearized new massive gravity (NMG) in D = 2 + 1 via Kaluza-Klein dimensional reduction followed by a dual master action. Here we show that a similar route takes us from the WTDIFF model to a linearized scalar-tensor NMG which belongs to a larger class of consistent spin-0 modifications of NMG. We also show that a traceless master action applied to a parity singlet furnishes two new spin-2 self-dual models. Moreover, we examine the singular replacement h_μ_ν → h_μ_ν - η_μ_νh/D and prove that it leads to consistent massive spin-2 models in D = 2 + 1. They include linearized versions of unimodular topologically massive gravity (TMG) and unimodular NMG. Although the free part of those unimodular theories are Weyl invariant, we do not expect any improvement in the renormalizability. Both the linearized K-term (in NMG) and the linearized gravitational Chern-Simons term (in TMG) are invariant under longitudinal reparametrizations δh_μ_ν = ∂_μ∂_νζ, which is not a symmetry of the WTDIFF Einstein-Hilbert term. Therefore, we still have one degree of freedom whose propagator behaves like 1/p"2 for large momentum. (orig.)

Weyl and transverse diffeomorphism invariant spin-2 models in D=2+1

Science.gov (United States)

Dalmazi, Denis; dos Santos, A. L. R.; Ghosh, Subir; Mendonça, E. L.

2017-09-01

There are two covariant descriptions of massless spin-2 particles in D=3+1 via a symmetric rank-2 tensor: the linearized Einstein-Hilbert (LEH) theory and the Weyl plus transverse diffeomorphism (WTDIFF) invariant model. From the LEH theory one can obtain the linearized new massive gravity (NMG) in D=2+1 via Kaluza-Klein dimensional reduction followed by a dual master action. Here we show that a similar route takes us from the WTDIFF model to a linearized scalar-tensor NMG which belongs to a larger class of consistent spin-0 modifications of NMG. We also show that a traceless master action applied to a parity singlet furnishes two new spin-2 self-dual models. Moreover, we examine the singular replacement h_{μ ν } → h_{μ ν } - η _{μ ν }h/D and prove that it leads to consistent massive spin-2 models in D=2+1. They include linearized versions of unimodular topologically massive gravity (TMG) and unimodular NMG. Although the free part of those unimodular theories are Weyl invariant, we do not expect any improvement in the renormalizability. Both the linearized K-term (in NMG) and the linearized gravitational Chern-Simons term (in TMG) are invariant under longitudinal reparametrizations δ h_{μ ν } = partial _{μ }partial _{ν }ζ , which is not a symmetry of the WTDIFF Einstein-Hilbert term. Therefore, we still have one degree of freedom whose propagator behaves like 1/p^2 for large momentum.

Dynamic gene expression analysis in a H1N1 influenza virus mouse pneumonia model.

Science.gov (United States)

Bao, Yanyan; Gao, Yingjie; Shi, Yujing; Cui, Xiaolan

2017-06-01

H1N1, a major pathogenic subtype of influenza A virus, causes a respiratory infection in humans and livestock that can range from a mild infection to more severe pneumonia associated with acute respiratory distress syndrome. Understanding the dynamic changes in the genome and the related functional changes induced by H1N1 influenza virus infection is essential to elucidating the pathogenesis of this virus and thereby determining strategies to prevent future outbreaks. In this study, we filtered the significantly expressed genes in mouse pneumonia using mRNA microarray analysis. Using STC analysis, seven significant gene clusters were revealed, and using STC-GO analysis, we explored the significant functions of these seven gene clusters. The results revealed GOs related to H1N1 virus-induced inflammatory and immune functions, including innate immune response, inflammatory response, specific immune response, and cellular response to interferon-beta. Furthermore, the dynamic regulation relationships of the key genes in mouse pneumonia were revealed by dynamic gene network analysis, and the most important genes were filtered, including Dhx58, Cxcl10, Cxcl11, Zbp1, Ifit1, Ifih1, Trim25, Mx2, Oas2, Cd274, Irgm1, and Irf7. These results suggested that during mouse pneumonia, changes in the expression of gene clusters and the complex interactions among genes lead to significant changes in function. Dynamic gene expression analysis revealed key genes that performed important functions. These results are a prelude to advancements in mouse H1N1 influenza virus infection biology, as well as the use of mice as a model organism for human H1N1 influenza virus infection studies.

Dynamic stall and 3D effects

Energy Technology Data Exchange (ETDEWEB)

Bjoerck, A.; Thor, S.E. [Aeronautical Research Inst. of Sweden, Bromma (Sweden)

1996-12-01

The JOULE II project `Dynamic stall and 3D effects` started in January 1994 and was completed in September 1995. The objective of the project has been to increase the understanding of the three-dimensional and unsteady aerodynamics of stall controlled HAWT`s. The objectives have also been to develop `engineering models` suitable for inclusion into aero-elastic codes. The project included the participation of 13 parties within Europe. This paper describes an overview of the work carried out within the project and key results. 3 refs, 4 figs

2D IR spectra of cyanide in water investigated by molecular dynamics simulations

Science.gov (United States)

Lee, Myung Won; Carr, Joshua K.; Göllner, Michael; Hamm, Peter; Meuwly, Markus

2013-01-01

Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN− solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN− molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN− and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T 1 times are sensitive to the van der Waals ranges on the CN− is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm−1 vs. 14.9 cm−1) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.

Non-cooperative Brownian donkeys: A solvable 1D model

Science.gov (United States)

Jiménez de Cisneros, B.; Reimann, P.; Parrondo, J. M. R.

2003-12-01

A paradigmatic 1D model for Brownian motion in a spatially symmetric, periodic system is tackled analytically. Upon application of an external static force F the system's response is an average current which is positive for F 0 (absolute negative mobility). Under suitable conditions, the system approaches 100% efficiency when working against the external force F.

Constitutive modeling and control of 1D smart composite structures

Science.gov (United States)

Briggs, Jonathan P.; Ostrowski, James P.; Ponte-Castaneda, Pedro

1998-07-01

Homogenization techniques for determining effective properties of composite materials may provide advantages for control of stiffness and strain in systems using hysteretic smart actuators embedded in a soft matrix. In this paper, a homogenized model of a 1D composite structure comprised of shape memory alloys and a rubber-like matrix is presented. With proportional and proportional/integral feedback, using current as the input state and global strain as an error state, implementation scenarios include the use of tractions on the boundaries and a nonlinear constitutive law for the matrix. The result is a simple model which captures the nonlinear behavior of the smart composite material system and is amenable to experiments with various control paradigms. The success of this approach in the context of the 1D model suggests that the homogenization method may prove useful in investigating control of more general smart structures. Applications of such materials could include active rehabilitation aids, e.g. wrist braces, as well as swimming/undulating robots, or adaptive molds for manufacturing processes.

Hybrid dynamics for currency modeling

OpenAIRE

Theodosopoulos, Ted; Trifunovic, Alex

2006-01-01

We present a simple hybrid dynamical model as a tool to investigate behavioral strategies based on trend following. The multiplicative symbolic dynamics are generated using a lognormal diffusion model for the at-the-money implied volatility term structure. Thus, are model exploits information from derivative markets to obtain qualititative properties of the return distribution for the underlier. We apply our model to the JPY-USD exchange rate and the corresponding 1mo., 3mo., 6mo. and 1yr. im...

Dynamic accelerator modeling

International Nuclear Information System (INIS)

Nishimura, Hiroshi.

1993-05-01

Object-Oriented Programming has been used extensively to model the LBL Advanced Light Source 1.5 GeV electron storage ring. This paper is on the present status of the class library construction with emphasis on a dynamic modeling

Prediction of car cabin environment by means of 1D and 3D cabin model

Directory of Open Access Journals (Sweden)

Jícha M.

2012-04-01

Full Text Available Thermal comfort and also reduction of energy requirements of air-conditioning system in vehicle cabins are currently very intensively investigated and up-to-date issues. The article deals with two approaches of modelling of car cabin environment; the first model was created in simulation language Modelica (typical 1D approach without cabin geometry and the second one was created in specialized software Theseus-FE (3D approach with cabin geometry. Performance and capabilities of this tools are demonstrated on the example of the car cabin and the results from simulations are compared with the results from the real car cabin climate chamber measurements.

Prediction Models for Dynamic Demand Response

Energy Technology Data Exchange (ETDEWEB)

Aman, Saima; Frincu, Marc; Chelmis, Charalampos; Noor, Muhammad; Simmhan, Yogesh; Prasanna, Viktor K.

2015-11-02

As Smart Grids move closer to dynamic curtailment programs, Demand Response (DR) events will become necessary not only on fixed time intervals and weekdays predetermined by static policies, but also during changing decision periods and weekends to react to real-time demand signals. Unique challenges arise in this context vis-a-vis demand prediction and curtailment estimation and the transformation of such tasks into an automated, efficient dynamic demand response (D²R) process. While existing work has concentrated on increasing the accuracy of prediction models for DR, there is a lack of studies for prediction models for D²R, which we address in this paper. Our first contribution is the formal definition of D²R, and the description of its challenges and requirements. Our second contribution is a feasibility analysis of very-short-term prediction of electricity consumption for D²R over a diverse, large-scale dataset that includes both small residential customers and large buildings. Our third, and major contribution is a set of insights into the predictability of electricity consumption in the context of D²R. Specifically, we focus on prediction models that can operate at a very small data granularity (here 15-min intervals), for both weekdays and weekends - all conditions that characterize scenarios for D²R. We find that short-term time series and simple averaging models used by Independent Service Operators and utilities achieve superior prediction accuracy. We also observe that workdays are more predictable than weekends and holiday. Also, smaller customers have large variation in consumption and are less predictable than larger buildings. Key implications of our findings are that better models are required for small customers and for non-workdays, both of which are critical for D²R. Also, prediction models require just few days’ worth of data indicating that small amounts of

3D Realistic Radiative Hydrodynamic Modeling of a Moderate-Mass Star: Effects of Rotation

Science.gov (United States)

Kitiashvili, Irina; Kosovichev, Alexander G.; Mansour, Nagi N.; Wray, Alan A.

2018-01-01

Recent progress in stellar observations opens new perspectives in understanding stellar evolution and structure. However, complex interactions in the turbulent radiating plasma together with effects of magnetic fields and rotation make inferences of stellar properties uncertain. The standard 1D mixing-length-based evolutionary models are not able to capture many physical processes of stellar interior dynamics, but they provide an initial approximation of the stellar structure that can be used to initialize 3D time-dependent radiative hydrodynamics simulations, based on first physical principles, that take into account the effects of turbulence, radiation, and others. In this presentation we will show simulation results from a 3D realistic modeling of an F-type main-sequence star with mass 1.47 Msun, in which the computational domain includes the upper layers of the radiation zone, the entire convection zone, and the photosphere. The simulation results provide new insight into the formation and properties of the convective overshoot region, the dynamics of the near-surface, highly turbulent layer, the structure and dynamics of granulation, and the excitation of acoustic and gravity oscillations. We will discuss the thermodynamic structure, oscillations, and effects of rotation on the dynamics of the star across these layers.

Weyl and transverse diffeomorphism invariant spin-2 models in D = 2 + 1

Energy Technology Data Exchange (ETDEWEB)

Dalmazi, Denis; Mendonca, E.L. [UNESP-Campus de Guaratingueta-DFQ, Guaratingueta, SP (Brazil); ICTP South American Institute for Fundamental Research, IFT-UNESP, Sao Paulo, SP (Brazil); Santos, A.L.R. dos [UNESP-Campus de Guaratingueta-DFQ, Guaratingueta, SP (Brazil); Ghosh, Subir [ICTP South American Institute for Fundamental Research, IFT-UNESP, Sao Paulo, SP (Brazil); Indian Statistical Institute, Physics and Applied Mathematics Unit, Kolkata (India)

2017-09-15

There are two covariant descriptions of massless spin-2 particles in D = 3 + 1 via a symmetric rank-2 tensor: the linearized Einstein-Hilbert (LEH) theory and the Weyl plus transverse diffeomorphism (WTDIFF) invariant model. From the LEH theory one can obtain the linearized new massive gravity (NMG) in D = 2 + 1 via Kaluza-Klein dimensional reduction followed by a dual master action. Here we show that a similar route takes us from the WTDIFF model to a linearized scalar-tensor NMG which belongs to a larger class of consistent spin-0 modifications of NMG. We also show that a traceless master action applied to a parity singlet furnishes two new spin-2 self-dual models. Moreover, we examine the singular replacement h{sub μν} → h{sub μν} - η{sub μν}h/D and prove that it leads to consistent massive spin-2 models in D = 2 + 1. They include linearized versions of unimodular topologically massive gravity (TMG) and unimodular NMG. Although the free part of those unimodular theories are Weyl invariant, we do not expect any improvement in the renormalizability. Both the linearized K-term (in NMG) and the linearized gravitational Chern-Simons term (in TMG) are invariant under longitudinal reparametrizations δh{sub μν} = ∂{sub μ}∂{sub ν}ζ, which is not a symmetry of the WTDIFF Einstein-Hilbert term. Therefore, we still have one degree of freedom whose propagator behaves like 1/p{sup 2} for large momentum. (orig.)

Biomechanical Analysis and Evaluation Technology Using Human Multi-Body Dynamic Model

Energy Technology Data Exchange (ETDEWEB)

Kim, Yoon Hyuk; Shin, June Ho; Khurelbaatar, Tsolmonbaatar [Kyung Hee University, Yongin (Korea, Republic of)

2011-10-15

This paper presents the biomechanical analysis and evaluation technology of musculoskeletal system by multi-body human dynamic model and 3-D motion capture data. First, medical image based geometric model and material properties of tissue were used to develop the human dynamic model and 3-D motion capture data based motion analysis techniques were develop to quantify the in-vivo joint kinematics, joint moment, joint force, and muscle force. Walking and push-up motion was investigated using the developed model. The present model and technologies would be useful to apply the biomechanical analysis and evaluation of human activities.

Understanding and Modeling Teams As Dynamical Systems

Science.gov (United States)

Gorman, Jamie C.; Dunbar, Terri A.; Grimm, David; Gipson, Christina L.

2017-01-01

By its very nature, much of teamwork is distributed across, and not stored within, interdependent people working toward a common goal. In this light, we advocate a systems perspective on teamwork that is based on general coordination principles that are not limited to cognitive, motor, and physiological levels of explanation within the individual. In this article, we present a framework for understanding and modeling teams as dynamical systems and review our empirical findings on teams as dynamical systems. We proceed by (a) considering the question of why study teams as dynamical systems, (b) considering the meaning of dynamical systems concepts (attractors; perturbation; synchronization; fractals) in the context of teams, (c) describe empirical studies of team coordination dynamics at the perceptual-motor, cognitive-behavioral, and cognitive-neurophysiological levels of analysis, and (d) consider the theoretical and practical implications of this approach, including new kinds of explanations of human performance and real-time analysis and performance modeling. Throughout our discussion of the topics we consider how to describe teamwork using equations and/or modeling techniques that describe the dynamics. Finally, we consider what dynamical equations and models do and do not tell us about human performance in teams and suggest future research directions in this area. PMID:28744231

Online dynamical downscaling of temperature and precipitation within the iLOVECLIM model (version 1.1)

Science.gov (United States)

Quiquet, Aurélien; Roche, Didier M.; Dumas, Christophe; Paillard, Didier

2018-02-01

This paper presents the inclusion of an online dynamical downscaling of temperature and precipitation within the model of intermediate complexity iLOVECLIM v1.1. We describe the following methodology to generate temperature and precipitation fields on a 40 km × 40 km Cartesian grid of the Northern Hemisphere from the T21 native atmospheric model grid. Our scheme is not grid specific and conserves energy and moisture in the same way as the original climate model. We show that we are able to generate a high-resolution field which presents a spatial variability in better agreement with the observations compared to the standard model. Although the large-scale model biases are not corrected, for selected model parameters, the downscaling can induce a better overall performance compared to the standard version on both the high-resolution grid and on the native grid. Foreseen applications of this new model feature include the improvement of ice sheet model coupling and high-resolution land surface models.

Constructing the reduced dynamical models of interannual climate variability from spatial-distributed time series

Science.gov (United States)

Mukhin, Dmitry; Gavrilov, Andrey; Loskutov, Evgeny; Feigin, Alexander

2016-04-01

We suggest a method for empirical forecast of climate dynamics basing on the reconstruction of reduced dynamical models in a form of random dynamical systems [1,2] derived from observational time series. The construction of proper embedding - the set of variables determining the phase space the model works in - is no doubt the most important step in such a modeling, but this task is non-trivial due to huge dimension of time series of typical climatic fields. Actually, an appropriate expansion of observational time series is needed yielding the number of principal components considered as phase variables, which are to be efficient for the construction of low-dimensional evolution operator. We emphasize two main features the reduced models should have for capturing the main dynamical properties of the system: (i) taking into account time-lagged teleconnections in the atmosphere-ocean system and (ii) reflecting the nonlinear nature of these teleconnections. In accordance to these principles, in this report we present the methodology which includes the combination of a new way for the construction of an embedding by the spatio-temporal data expansion and nonlinear model construction on the basis of artificial neural networks. The methodology is aplied to NCEP/NCAR reanalysis data including fields of sea level pressure, geopotential height, and wind speed, covering Northern Hemisphere. Its efficiency for the interannual forecast of various climate phenomena including ENSO, PDO, NAO and strong blocking event condition over the mid latitudes, is demonstrated. Also, we investigate the ability of the models to reproduce and predict the evolution of qualitative features of the dynamics, such as spectral peaks, critical transitions and statistics of extremes. This research was supported by the Government of the Russian Federation (Agreement No. 14.Z50.31.0033 with the Institute of Applied Physics RAS) [1] Y. I. Molkov, E. M. Loskutov, D. N. Mukhin, and A. M. Feigin, "Random

A 1D constitutive model for shape memory alloy using strain and temperature as control variables and including martensite reorientation and asymmetric behaviors

International Nuclear Information System (INIS)

Jaber, M Ben; Mehrez, S; Ghazouani, O

2014-01-01

In this paper, a new 1D constitutive model for shape memory alloy using strain and temperature as control variables is presented. The new formulation is restricted to the 1D stress case and takes into account the martensite reorientation and the asymmetry of the SMA behavior in tension and compression. Numerical implementation of the new model in a finite element code was conducted. The numerical results for superelastic behavior in tension and compression tests are presented and were compared to experimental data taken from the literature. Other numerical tests are presented, showing the model’s ability to reproduce the main aspects of SMA behavior such as the shape memory effect and the martensite reorientation under cyclic loading. Finally, to demonstrate the utility of the new constitutive model, a dynamic test of a bi-clamped SMA bending beam under forced oscillation is described. (paper)

USING AFFORDABLE DATA CAPTURING DEVICES FOR AUTOMATIC 3D CITY MODELLING

Directory of Open Access Journals (Sweden)

B. Alizadehashrafi

2017-11-01

Full Text Available In this research project, many movies from UTM Kolej 9, Skudai, Johor Bahru (See Figure 1 were taken by AR. Drone 2. Since the AR drone 2.0 has liquid lens, while flying there were significant distortions and deformations on the converted pictures of the movies. Passive remote sensing (RS applications based on image matching and Epipolar lines such as Agisoft PhotoScan have been tested to create the point clouds and mesh along with 3D models and textures. As the result was not acceptable (See Figure 2, the previous Dynamic Pulse Function based on Ruby programming language were enhanced and utilized to create the 3D models automatically in LoD3. The accuracy of the final 3D model is almost 10 to 20 cm. After rectification and parallel projection of the photos based on some tie points and targets, all the parameters were measured and utilized as an input to the system to create the 3D model automatically in LoD3 in a very high accuracy.

Using Affordable Data Capturing Devices for Automatic 3d City Modelling

Science.gov (United States)

Alizadehashrafi, B.; Abdul-Rahman, A.

2017-11-01

In this research project, many movies from UTM Kolej 9, Skudai, Johor Bahru (See Figure 1) were taken by AR. Drone 2. Since the AR drone 2.0 has liquid lens, while flying there were significant distortions and deformations on the converted pictures of the movies. Passive remote sensing (RS) applications based on image matching and Epipolar lines such as Agisoft PhotoScan have been tested to create the point clouds and mesh along with 3D models and textures. As the result was not acceptable (See Figure 2), the previous Dynamic Pulse Function based on Ruby programming language were enhanced and utilized to create the 3D models automatically in LoD3. The accuracy of the final 3D model is almost 10 to 20 cm. After rectification and parallel projection of the photos based on some tie points and targets, all the parameters were measured and utilized as an input to the system to create the 3D model automatically in LoD3 in a very high accuracy.

2D-3D Registration of CT Vertebra Volume to Fluoroscopy Projection: A Calibration Model Assessment

Directory of Open Access Journals (Sweden)

P. Bifulco

2010-01-01

Full Text Available 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.

Three-field modeling for MARS 1-D code

International Nuclear Information System (INIS)

Hwang, Moonkyu; Lim, Ho-Gon; Jeong, Jae-Jun; Chung, Bub-Dong

2006-01-01

In this study, the three-field modeling of the two-phase mixture is developed. The finite difference equations for the three-field equations thereafter are devised. The solution scheme has been implemented into the MARS 1-D code. The three-field formulations adopted are similar to those for MARS 3-D module, in a sense that the mass and momentum are treated separately for the entrained liquid and continuous liquid. As in the MARS-3D module, the entrained liquid and continuous liquid are combined into one for the energy equation, assuming thermal equilibrium between the two. All the non-linear terms are linearized to arrange the finite difference equation set into a linear matrix form with respect to the unknown arguments. The problems chosen for the assessment of the newly added entrained field consist of basic conceptual tests. Among the tests are gas-only test, liquid-only test, gas-only with supplied entrained liquid test, Edwards pipe problem, and GE level swell problem. The conceptual tests performed confirm the sound integrity of the three-field solver

Modelling reinforced concrete structures in DYNA3D

International Nuclear Information System (INIS)

Broadhouse, B.J.; Neilson, A.J.

1987-10-01

A material model for reinforced concrete has been implemented in the transient structural dynamics code DYNA3D. This paper outlines the constitutive material model, and presents comparisons of DYNA3D calculations and experiments on impulsively loaded panels, covering the full range of panel damage states from light cracking through to panel collapse or perforation. The results are presented using the post-processor code TAURUS, which has also been modified to provide mesh diagrams with superimposed crack patterns from the DYNA3D predictions. (author)

Multiobjective Traffic Signal Control Model for Intersection Based on Dynamic Turning Movements Estimation

Directory of Open Access Journals (Sweden)

Pengpeng Jiao

2014-01-01

Full Text Available The real-time traffic signal control for intersection requires dynamic turning movements as the basic input data. It is impossible to detect dynamic turning movements directly through current traffic surveillance systems, but dynamic origin-destination (O-D estimation can obtain it. However, the combined models of dynamic O-D estimation and real-time traffic signal control are rare in the literature. A framework for the multiobjective traffic signal control model for intersection based on dynamic O-D estimation (MSC-DODE is presented. A state-space model using Kalman filtering is first formulated to estimate the dynamic turning movements; then a revised sequential Kalman filtering algorithm is designed to solve the model, and the root mean square error and mean percentage error are used to evaluate the accuracy of estimated dynamic turning proportions. Furthermore, a multiobjective traffic signal control model is put forward to achieve real-time signal control parameters and evaluation indices. Finally, based on practical survey data, the evaluation indices from MSC-DODE are compared with those from Webster method. The actual and estimated turning movements are further input into MSC-DODE, respectively, and results are also compared. Case studies show that results of MSC-DODE are better than those of Webster method and are very close to unavailable actual values.

STREAM PROCESSING ALGORITHMS FOR DYNAMIC 3D SCENE ANALYSIS

Science.gov (United States)

2018-02-15

PROCESSING ALGORITHMS FOR DYNAMIC 3D SCENE ANALYSIS 5a. CONTRACT NUMBER FA8750-14-2-0072 5b. GRANT NUMBER N/A 5c. PROGRAM ELEMENT NUMBER 62788F 6...of Figures 1 The 3D processing pipeline flowchart showing key modules. . . . . . . . . . . . . . . . . 12 2 Overall view (data flow) of the proposed...pipeline flowchart showing key modules. from motion and bundle adjustment algorithm. By fusion of depth masks of the scene obtained from 3D

Dynamics modeling for parallel haptic interfaces with force sensing and control.

Science.gov (United States)

Bernstein, Nicholas; Lawrence, Dale; Pao, Lucy

2013-01-01

Closed-loop force control can be used on haptic interfaces (HIs) to mitigate the effects of mechanism dynamics. A single multidimensional force-torque sensor is often employed to measure the interaction force between the haptic device and the user's hand. The parallel haptic interface at the University of Colorado (CU) instead employs smaller 1D force sensors oriented along each of the five actuating rods to build up a 5D force vector. This paper shows that a particular manipulandum/hand partition in the system dynamics is induced by the placement and type of force sensing, and discusses the implications on force and impedance control for parallel haptic interfaces. The details of a "squaring down" process are also discussed, showing how to obtain reduced degree-of-freedom models from the general six degree-of-freedom dynamics formulation.

Dynamical compactification of D-dimensional gravity coupled to antisymmetric tensors in a 1/D expansion

International Nuclear Information System (INIS)

Foda, O.

1984-12-01

The effective potential of components of the curl of an antisymmetric tensor coupled to gravity in D dimensions is evaluated in a 1/D expansion. For large D, only highest-rank propagators contribute to leading order, while multiloop diagrams are suppressed by phase-space factors. Divergences are regulated by a cut-off LAMBDA, that we interpret as the mass-breaking scale of a larger theory that is finite. As an application we consider the bosonic sector of D=11, N=1 supergravity. If the full theory is finite, then LAMBDA is msub(SUSY): the scale below which the fermion sector decouples. For m 9 sub(SUSY)>1/akappa 2 , (kappa 2 : the D=11 Newton's coupling, a approx.= O(1)) the 11-dimensional symmetric vacuum is unstable under compactification. For m 9 sub(SUSY) 2 , it is metastable. To leading order in 1/D, all gauge dependence cancels identically, while ghosts as well as the graviton decouple. (author)

Modeling and Optimal Control of a Class of Warfare Hybrid Dynamic Systems Based on Lanchester (n,1) Attrition Model

OpenAIRE

Chen, Xiangyong; Zhang, Ancai

2014-01-01

For the particularity of warfare hybrid dynamic process, a class of warfare hybrid dynamic systems is established based on Lanchester equation in a (n,1) battle, where a heterogeneous force of n different troop types faces a homogeneous force. This model can be characterized by the interaction of continuous-time models (governed by Lanchester equation), and discrete event systems (described by variable tactics). Furthermore, an expository discussion is presented on an optimal variable tact...

ABSTRACTION OF INFORMATION FROM 2- AND 3-DIMENSIONAL PORFLOW MODELS INTO A 1-D GOLDSIM MODEL - 11404

International Nuclear Information System (INIS)

Taylor, G.; Hiergesell, R.

2010-01-01

The Savannah River National Laboratory has developed a 'hybrid' approach to Performance Assessment modeling which has been used for a number of Performance Assessments. This hybrid approach uses a multi-dimensional modeling platform (PorFlow) to develop deterministic flow fields and perform contaminant transport. The GoldSim modeling platform is used to develop the Sensitivity and Uncertainty analyses. Because these codes are performing complementary tasks, it is incumbent upon them that for the deterministic cases they produce very similar results. This paper discusses two very different waste forms, one with no engineered barriers and one with engineered barriers, each of which present different challenges to the abstraction of data. The hybrid approach to Performance Assessment modeling used at the SRNL uses a 2-D unsaturated zone (UZ) and a 3-D saturated zone (SZ) model in the PorFlow modeling platform. The UZ model consists of the waste zone and the unsaturated zoned between the waste zone and the water table. The SZ model consists of source cells beneath the waste form to the points of interest. Both models contain 'buffer' cells so that modeling domain boundaries do not adversely affect the calculation. The information pipeline between the two models is the contaminant flux. The domain contaminant flux, typically in units of moles (or Curies) per year from the UZ model is used as a boundary condition for the source cells in the SZ. The GoldSim modeling component of the hybrid approach is an integrated UZ-SZ model. The model is a 1-D representation of the SZ, typically 1-D in the UZ, but as discussed below, depending on the waste form being analyzed may contain pseudo-2-D elements. A waste form at the Savannah River Site (SRS) which has no engineered barriers is commonly referred to as a slit trench. A slit trench, as its name implies, is an unlined trench, typically 6 m deep, 6 m wide, and 200 m long. Low level waste consisting of soil, debris, rubble, wood

Generalized model development for a cryo-adsorber and 1-D results for the isobaric refueling period

Energy Technology Data Exchange (ETDEWEB)

Kumar, V. Senthil [India Science Lab, General Motors Global R and D, Creator Building, International Technology Park, Bangalore 560066 (India); Kumar, Sudarshan [Chemical Sciences and Material Systems Lab, General Motors Global R and D, Warren Technical Center Campus, 30500 Mound Road, Warren, MI 48090 (United States)

2010-04-15

We have developed 3-D model equations for a cryo-adsorption hydrogen storage tank, where the energy balance accommodates the temperature and pressure variation of all the thermodynamic properties. We then reduce the 3-D model to the 1-D isobaric system and study the isobaric refueling period, for simplified geometry and charging conditions. The hydrogen capacity evolution predicted by the 1-D axial bed model is significantly different than that predicted by the lumped-parameter model because of the presence of sharp temperature gradients during refueling. The 1-D model predicts a higher hydrogen capacity than the lumped-parameter model. This observation can be rationalized by the fact that a bed with temperature gradients on equilibration should desorb gas, whenever the adsorbed phase entropy is lower than the gas phase entropy. The 1-D analysis of the isobaric refueling period does not show any significant difference in hydrogen capacity evolution among the axial, single and multicartridge annular bed designs. Hence, a multicartridge annular design, though giving a slightly lower pressure drop, does not offer any heat and mass transfer enhancement over the single cartridge design. And, the single cartridge annular design appears to be optimal. (author)

Influence of neuromuscular noise and walking speed on fall risk and dynamic stability in a 3D dynamic walking model.

Science.gov (United States)

Roos, Paulien E; Dingwell, Jonathan B

2013-06-21

Older adults and those with increased fall risk tend to walk slower. They may do this voluntarily to reduce their fall risk. However, both slower and faster walking speeds can predict increased risk of different types of falls. The mechanisms that contribute to fall risk across speeds are not well known. Faster walking requires greater forward propulsion, generated by larger muscle forces. However, greater muscle activation induces increased signal-dependent neuromuscular noise. These speed-related increases in neuromuscular noise may contribute to the increased fall risk observed at faster walking speeds. Using a 3D dynamic walking model, we systematically varied walking speed without and with physiologically-appropriate neuromuscular noise. We quantified how actual fall risk changed with gait speed, how neuromuscular noise affected speed-related changes in fall risk, and how well orbital and local dynamic stability measures predicted changes in fall risk across speeds. When we included physiologically-appropriate noise to the 'push-off' force in our model, fall risk increased with increasing walking speed. Changes in kinematic variability, orbital, and local dynamic stability did not predict these speed-related changes in fall risk. Thus, the increased neuromuscular variability that results from increased signal-dependent noise that is necessitated by the greater muscular force requirements of faster walking may contribute to the increased fall risk observed at faster walking speeds. The lower fall risk observed at slower speeds supports experimental evidence that slowing down can be an effective strategy to reduce fall risk. This may help explain the slower walking speeds observed in older adults and others. Copyright © 2013 Elsevier Ltd. All rights reserved.

Dopamine D1 signaling organizes network dynamics underlying working memory.

Science.gov (United States)

Roffman, Joshua L; Tanner, Alexandra S; Eryilmaz, Hamdi; Rodriguez-Thompson, Anais; Silverstein, Noah J; Ho, New Fei; Nitenson, Adam Z; Chonde, Daniel B; Greve, Douglas N; Abi-Dargham, Anissa; Buckner, Randy L; Manoach, Dara S; Rosen, Bruce R; Hooker, Jacob M; Catana, Ciprian

2016-06-01

Local prefrontal dopamine signaling supports working memory by tuning pyramidal neurons to task-relevant stimuli. Enabled by simultaneous positron emission tomography-magnetic resonance imaging (PET-MRI), we determined whether neuromodulatory effects of dopamine scale to the level of cortical networks and coordinate their interplay during working memory. Among network territories, mean cortical D1 receptor densities differed substantially but were strongly interrelated, suggesting cross-network regulation. Indeed, mean cortical D1 density predicted working memory-emergent decoupling of the frontoparietal and default networks, which respectively manage task-related and internal stimuli. In contrast, striatal D1 predicted opposing effects within these two networks but no between-network effects. These findings specifically link cortical dopamine signaling to network crosstalk that redirects cognitive resources to working memory, echoing neuromodulatory effects of D1 signaling on the level of cortical microcircuits.

Registration of dynamic dopamine D2receptor images using principal component analysis

International Nuclear Information System (INIS)

Acton, P.D.; Ell, P.J.; Pilowsky, L.S.; Brammer, M.J.; Suckling, J.

1997-01-01

This paper describes a novel technique for registering a dynamic sequence of single-photon emission tomography (SPET) dopamine D 2 receptor images, using principal component analysis (PCA). Conventional methods for registering images, such as count difference and correlation coefficient algorithms, fail to take into account the dynamic nature of the data, resulting in large systematic errors when registering time-varying images. However, by using principal component analysis to extract the temporal structure of the image sequence, misregistration can be quantified by examining the distribution of eigenvalues. The registration procedures were tested using a computer-generated dynamic phantom derived from a high-resolution magnetic resonance image of a realistic brain phantom. Each method was also applied to clinical SPET images of dopamine D 2 receptors, using the ligands iodine-123 iodobenzamide and iodine-123 epidepride, to investigate the influence of misregistration on kinetic modelling parameters and the binding potential. The PCA technique gave highly significant (P 123 I-epidepride scans. The PCA method produced data of much greater quality for subsequent kinetic modelling, with an improvement of nearly 50% in the χ 2 of the fit to the compartmental model, and provided superior quality registration of particularly difficult dynamic sequences. (orig.)

Development and application of a dynamic stall model for rotating wind turbine blades

International Nuclear Information System (INIS)

Xu, B F; Yuan, Y; Wang, T G

2014-01-01

In unsteady conditions of wind turbines, both the dynamic stall phenomenon and the three-dimensional (3D) rotational effect affect the rotor aerodynamics. The dynamic stall mechanism for rotating wind turbine blades is first investigated. Through the comparison of the aerodynamic data between the rotating blade and the two-dimensional (2D) airfoil, the normal force slope in the attached flow and the separation point expression in the separated flow are modified in the Beddoes-Leishman (B-L) dynamic stall model for rotating NREL wind turbine blades. The modified model is validated by the comparison between the calculation results and the experimental results of the lift and drag coefficients at different radial positions. Both the hysteresis loop shapes and the calculation values are closer to the experiment than the 2D dynamic stall model. The present dynamic stall model is then coupled to a free vortex wake model. The coupled model is used to calculate the unsteady blade aerodynamic loads and the low speed shaft torque of the NREL wind turbine in a yawed condition. The accuracy is greatly improved by the corrections presented in the paper

Computational Fluid Dynamics Modeling of a Supersonic Nozzle and Integration into a Variable Cycle Engine Model

Science.gov (United States)

Connolly, Joseph W.; Friedlander, David; Kopasakis, George

2015-01-01

This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.

Expression of the dopaminergic D1 and D2 receptors in the anterior cingulate cortex in a model of neuropathic pain

Directory of Open Access Journals (Sweden)

Ortega-Legaspi J Manuel

2011-12-01

Full Text Available Abstract Background The anterior cingulate cortex (ACC has been related to the affective component of pain. Dopaminergic mesocortical circuits, including the ACC, are able to inhibit neuropathic nociception measured as autotomy behaviour. We determined the changes in dopamine D1 and D2 (D1R and D2R receptor expression in the ACC (cg1 and cg2 in an animal model of neuropathic pain. The neuropathic group had noxious heat applied in the right hind paw followed 30 min. later by right sciatic denervation. Autotomy score (AS was recorded for eight days and subsequently classified in low, medium and high AS groups. The control consisted of naïve animals. A semiquantitative RT-PCR procedure was done to determine mRNA levels for D1R and D2R in cg1 and cg2, and protein levels were measured by Western Blot. Results The results of D1R mRNA in cg1 showed a decrease in all groups. D2R mRNA levels in cg1 decreased in low AS and increased in medium and high AS. Regarding D1R in cg2, there was an increase in all groups. D2R expression levels in cg2 decreased in all groups. In cg1, the D2R mRNA correlated positively with autotomy behaviour. Protein levels of D2R in cg1 increased in all groups but to a higher degree in low AS. In cg2 D2R protein only decreased discretely. D1R protein was not found in either ACC region. Conclusions This is the first evidence of an increase of inhibitory dopaminergic receptor (D2R mRNA and protein in cg1 in correlation with nociceptive behaviour in a neuropathic model of pain in the rat.

Variability of dynamic source parameters inferred from kinematic models of past earthquakes

KAUST Repository

Causse, M.

2013-12-24

We analyse the scaling and distribution of average dynamic source properties (fracture energy, static, dynamic and apparent stress drops) using 31 kinematic inversion models from 21 crustal earthquakes. Shear-stress histories are computed by solving the elastodynamic equations while imposing the slip velocity of a kinematic source model as a boundary condition on the fault plane. This is achieved using a 3-D finite difference method in which the rupture kinematics are modelled with the staggered-grid-split-node fault representation method of Dalguer & Day. Dynamic parameters are then estimated from the calculated stress-slip curves and averaged over the fault plane. Our results indicate that fracture energy, static, dynamic and apparent stress drops tend to increase with magnitude. The epistemic uncertainty due to uncertainties in kinematic inversions remains small (ϕ ∼ 0.1 in log10 units), showing that kinematic source models provide robust information to analyse the distribution of average dynamic source parameters. The proposed scaling relations may be useful to constrain friction law parameters in spontaneous dynamic rupture calculations for earthquake source studies, and physics-based near-source ground-motion prediction for seismic hazard and risk mitigation.

4D rotational x-ray imaging of wrist joint dynamic motion

International Nuclear Information System (INIS)

Carelsen, Bart; Bakker, Niels H.; Strackee, Simon D.; Boon, Sjirk N.; Maas, Mario; Sabczynski, Joerg; Grimbergen, Cornelis A.; Streekstra, Geert J.

2005-01-01

Current methods for imaging joint motion are limited to either two-dimensional (2D) video fluoroscopy, or to animated motions from a series of static three-dimensional (3D) images. 3D movement patterns can be detected from biplane fluoroscopy images matched with computed tomography images. This involves several x-ray modalities and sophisticated 2D to 3D matching for the complex wrist joint. We present a method for the acquisition of dynamic 3D images of a moving joint. In our method a 3D-rotational x-ray (3D-RX) system is used to image a cyclically moving joint. The cyclic motion is synchronized to the x-ray acquisition to yield multiple sets of projection images, which are reconstructed to a series of time resolved 3D images, i.e., four-dimensional rotational x ray (4D-RX). To investigate the obtained image quality parameters the full width at half maximum (FWHM) of the point spread function (PSF) via the edge spread function and the contrast to noise ratio between air and phantom were determined on reconstructions of a bullet and rod phantom, using 4D-RX as well as stationary 3D-RX images. The CNR in volume reconstructions based on 251 projection images in the static situation and on 41 and 34 projection images of a moving phantom were 6.9, 3.0, and 2.9, respectively. The average FWHM of the PSF of these same images was, respectively, 1.1, 1.7, and 2.2 mm orthogonal to the motion and parallel to direction of motion 0.6, 0.7, and 1.0 mm. The main deterioration of 4D-RX images compared to 3D-RX images is due to the low number of projection images used and not to the motion of the object. Using 41 projection images seems the best setting for the current system. Experiments on a postmortem wrist show the feasibility of the method for imaging 3D dynamic joint motion. We expect that 4D-RX will pave the way to improved assessment of joint disorders by detection of 3D dynamic motion patterns in joints

Numerical modeling of a pitch oscillating S809 airfoil dynamic stall in 2D with application to a horizontal axis wind turbine

Energy Technology Data Exchange (ETDEWEB)

Gharali, K.; Johnson, D.A. [Waterloo Univ., ON (Canada). Dept. of Mechanical and Mechatronics Engineering, Wind Energy Group

2010-07-01

Natural wind can sometimes have a strong wind shear that causes the Dynamic Stall (DS) phenomena which may result in dynamic loads and varying lift coefficients. The DS phenomena cannot be prevented in horizontal axis wind turbines (HAWTs). Therefore, it is necessary to study the unsteady aerodynamics in order to modify common wind turbine rotor designs. This paper reported on a study that investigated the dynamic flow fields around an oscillating 2D S809 airfoil, representing the aerodynamic characteristics of HAWT airfoils for dynamic stall conditions. A computational fluid dynamic (CFD) flow solver package with Fluent was used with different turbulence models, notably the Spalart-Allmaras and Detached Eddy Simulation (DES) methods. A sliding mesh is commonly used in numerical methods for simulating an oscillating foil, but sliding meshes suffer from mesh generation complexity and increased computational time. In this study, instead of a sinusoidally pitching airfoil, the direction of the far-field flow was changed according to a user-defined function in the software to simulate a proper angle of attack for the boundary conditions in each time step. This strategy helped to decrease processing time. The simulation results were in good agreement with experimental data and the Beddoes-Leishman model results. The DES method for unsteady 2D flow was not recommended. It was concluded that the Fluent package is time efficient, reliable and economic for the wind turbine industry. 17 refs., 3 figs.

1-D hybrid code for FRM start-up

International Nuclear Information System (INIS)

Stark, R.A.; Miley, G.H.

1982-01-01

A one-D hybrid has been developed to study the start-up of the FRM via neutral-beam injection. The code uses a multi-group numerical model originally developed by J. Willenberg to describe fusion product dynamics in a solenoidal plasma. Earlier we described such a model for use in determining self-consistent ion currents and magnetic fields in FRM start-up. However, consideration of electron dynamics during start-up indicate that the electron current will oppose the injected ion current and may even foil the attempt to achieve reversal. For this reason, we have combined the multi-group ion (model) with a fluid treatment for electron dynamics to form the hybrid code FROST (Field Reversed One-dimensional STart-up). The details of this merger, along with sample results of operation of FROST, are given

A dynamical model on deposit and loan of banking: A bifurcation analysis

Science.gov (United States)

Sumarti, Novriana; Hasmi, Abrari Noor

2015-09-01

A dynamical model, which is one of sophisticated techniques using mathematical equations, can determine the observed state, for example bank profits, for all future times based on the current state. It will also show small changes in the state of the system create either small or big changes in the future depending on the model. In this research we develop a dynamical system of the form: d/D d t =f (D ,L ,rD,rL,r ), d/L d t =g (D ,L ,rD,rL,r ), Here D and rD are the volume of deposit and its rate, L and rL are the volume of loan and its rate, and r is the interbank market rate. There are parameters required in this model which give connections between two variables or between two derivative functions. In this paper we simulate the model for several parameters values. We do bifurcation analysis on the dynamics of the system in order to identify the appropriate parameters that control the stability behaviour of the system. The result shows that the system will have a limit cycle for small value of interest rate of loan, so the deposit and loan volumes are fluctuating and oscillating extremely. If the interest rate of loan is too high, the loan volume will be decreasing and vanish and the system will converge to its carrying capacity.

Modelling atmospheric transport of persistent organic pollutants in the Northern Hemisphere with a 3-D dynamical model: DEHM-POP

Science.gov (United States)

Hansen, K. M.; Christensen, J. H.; Brandt, J.; Frohn, L. M.; Geels, C.

2004-03-01

The Danish Eulerian Hemispheric Model (DEHM) is a 3-D dynamical atmospheric transport model originally developed to describe the atmospheric transport of sulphur into the Arctic. A new version of the model, DEHM-POP, developed to study the atmospheric transport and environmental fate of persistent organic pollutants (POPs) is presented. During environmental cycling, POPs can be deposited and re-emitted several times before reaching a final destination. A description of the exchange processes between the land/ocean surfaces and the atmosphere is included in the model to account for this multi-hop transport. The α-isomer of the pesticide hexachlorocyclohexane (α-HCH) is used as tracer in the model development. The structure of the model and processes included are described in detail. The results from a model simulation showing the atmospheric transport for the years 1991 to 1998 are presented and evaluated against measurements. The annual averaged atmospheric concentration of α-HCH for the 1990s is well described by the model; however, the shorter-term average concentration for most of the stations is not well captured. This indicates that the present simple surface description needs to be refined to get a better description of the air-surface exchange proceses of POPs.

Fuel temperature influence on the performance of a last generation common-rail diesel ballistic injector. Part II: 1D model development, validation and analysis

International Nuclear Information System (INIS)

Payri, R.; Salvador, F.J.; Carreres, M.; De la Morena, J.

2016-01-01

Highlights: • A 1D model of a solenoid common-rail ballistic injector is implemented in AMESim. • A detailed dimensional and a hydraulic characterization lead to a fair validation. • Fuel temperature influence on injector dynamics is assessed through 1D simulations. • Temperature impacts through changes in inlet orifice regime and viscous friction. • Cold fuel temperature leads to a slower injection opening due to high viscosity. - Abstract: A one-dimensional model of a solenoid-driven common-rail diesel injector has been developed in order to study the influence of fuel temperature on the injection process. The model has been implemented after a thorough characterization of the injector, both from the dimensional and the hydraulic point of view. In this sense, experimental tools for the determination of the geometry of the injector lines and orifices have been described in the paper, together with the hydraulic setup introduced to characterize the flow behaviour through the calibrated orifices. An extensive validation of the model has been performed by comparing the modelled mass flow rate against the experimental results introduced in the first part of the paper, which were performed for different engine-like operating conditions involving a wide range of fuel temperatures, injection pressures and energizing times. In that first part of the study, an important influence of the fuel temperature was reported, especially in terms of the dynamic behaviour of the injector, due to its ballistic nature. The results from the model have allowed to explain and further extend the findings of the experimental study by analyzing key features of the injector dynamics, such as the pressure drop established in the control volume due to the control orifices performance or the forces due to viscous friction, also assessing their influence on the needle lift laws.

Molecular Modeling Studies of 11β-Hydroxysteroid Dehydrogenase Type 1 Inhibitors through Receptor-Based 3D-QSAR and Molecular Dynamics Simulations.

Science.gov (United States)

Qian, Haiyan; Chen, Jiongjiong; Pan, Youlu; Chen, Jianzhong

2016-09-19

11β-Hydroxysteroid dehydrogenase type 1 (11β-HSD1) is a potential target for the treatment of numerous human disorders, such as diabetes, obesity, and metabolic syndrome. In this work, molecular modeling studies combining molecular docking, 3D-QSAR, MESP, MD simulations and free energy calculations were performed on pyridine amides and 1,2,4-triazolopyridines as 11β-HSD1 inhibitors to explore structure-activity relationships and structural requirement for the inhibitory activity. 3D-QSAR models, including CoMFA and CoMSIA, were developed from the conformations obtained by docking strategy. The derived pharmacophoric features were further supported by MESP and Mulliken charge analyses using density functional theory. In addition, MD simulations and free energy calculations were employed to determine the detailed binding process and to compare the binding modes of inhibitors with different bioactivities. The binding free energies calculated by MM/PBSA showed a good correlation with the experimental biological activities. Free energy analyses and per-residue energy decomposition indicated the van der Waals interaction would be the major driving force for the interactions between an inhibitor and 11β-HSD1. These unified results may provide that hydrogen bond interactions with Ser170 and Tyr183 are favorable for enhancing activity. Thr124, Ser170, Tyr177, Tyr183, Val227, and Val231 are the key amino acid residues in the binding pocket. The obtained results are expected to be valuable for the rational design of novel potent 11β-HSD1 inhibitors.

Molecular Modeling Studies of 11β-Hydroxysteroid Dehydrogenase Type 1 Inhibitors through Receptor-Based 3D-QSAR and Molecular Dynamics Simulations

Directory of Open Access Journals (Sweden)

Haiyan Qian

2016-09-01

Full Text Available 11β-Hydroxysteroid dehydrogenase type 1 (11β-HSD1 is a potential target for the treatment of numerous human disorders, such as diabetes, obesity, and metabolic syndrome. In this work, molecular modeling studies combining molecular docking, 3D-QSAR, MESP, MD simulations and free energy calculations were performed on pyridine amides and 1,2,4-triazolopyridines as 11β-HSD1 inhibitors to explore structure-activity relationships and structural requirement for the inhibitory activity. 3D-QSAR models, including CoMFA and CoMSIA, were developed from the conformations obtained by docking strategy. The derived pharmacophoric features were further supported by MESP and Mulliken charge analyses using density functional theory. In addition, MD simulations and free energy calculations were employed to determine the detailed binding process and to compare the binding modes of inhibitors with different bioactivities. The binding free energies calculated by MM/PBSA showed a good correlation with the experimental biological activities. Free energy analyses and per-residue energy decomposition indicated the van der Waals interaction would be the major driving force for the interactions between an inhibitor and 11β-HSD1. These unified results may provide that hydrogen bond interactions with Ser170 and Tyr183 are favorable for enhancing activity. Thr124, Ser170, Tyr177, Tyr183, Val227, and Val231 are the key amino acid residues in the binding pocket. The obtained results are expected to be valuable for the rational design of novel potent 11β-HSD1 inhibitors.

Online dynamical downscaling of temperature and precipitation within the iLOVECLIM model (version 1.1

Directory of Open Access Journals (Sweden)

A. Quiquet

2018-02-01

Full Text Available This paper presents the inclusion of an online dynamical downscaling of temperature and precipitation within the model of intermediate complexity iLOVECLIM v1.1. We describe the following methodology to generate temperature and precipitation fields on a 40 km  ×  40 km Cartesian grid of the Northern Hemisphere from the T21 native atmospheric model grid. Our scheme is not grid specific and conserves energy and moisture in the same way as the original climate model. We show that we are able to generate a high-resolution field which presents a spatial variability in better agreement with the observations compared to the standard model. Although the large-scale model biases are not corrected, for selected model parameters, the downscaling can induce a better overall performance compared to the standard version on both the high-resolution grid and on the native grid. Foreseen applications of this new model feature include the improvement of ice sheet model coupling and high-resolution land surface models.

Probabilistic modelling of the high-pressure arc cathode spot displacement dynamic

CERN Document Server

Coulombe, S

2003-01-01

A probabilistic modelling approach for the study of the cathode spot displacement dynamic in high-pressure arc systems is developed in an attempt to interpret the observed voltage fluctuations. The general framework of the model allows to define simple, probabilistic displacement rules, the so-called cathode spot dynamic rules, for various possible surface states (un-arced metal, arced, contaminated) and to study the resulting dynamic of the cathode spot displacements over one or several arc passages. The displacements of the type-A cathode spot (macro-spot) in a magnetically rotating arc using concentric electrodes made up of either clean or contaminated metal surfaces is considered. Experimental observations for this system revealed a 1/f sup - sup t sup i sup l sup d sup e sup 1 signature in the frequency power spectrum (FPS) of the arc voltage for anchoring arc conditions on the cathode (e.g. clean metal surface), while it shows a 'white noise' signature for conditions favouring a smooth movement (e.g. ox...

CP violation in D0- D-bar0 oscillations: general considerations and applications to the Littlest Higgs model with T-parity

International Nuclear Information System (INIS)

Bigi, Ikaros I.; Blanke, Monika; Buras, Andrzej J.; Recksiegel, Stefan

2009-01-01

The observed D 0 -D-bar 0 oscillations provide a new stage in our search for New Physics in heavy flavour dynamics. The theoretical verdict on the observed values of x D and y D remains ambiguous: while they could be totally generated by Standard Model dynamics, they could also contain a sizable or even leading contribution from New Physics. Those oscillations are likely to enhance the observability of CP violation as clear manifestations of New Physics. We present general formulae for D 0 - D-bar 0 oscillations, concentrating on the case of negligible direct CP violation. In particular we derive a general formula for the time-dependent mixing-induced CP asymmetry in decays to a CP eigenstate and its correlation with the semileptonic CP asymmetry a SL (D 0 ) in D 0 (t) → lνK. We apply our formalism to the Littlest Higgs model with T-parity, using the time-dependent CP asymmetry in D 0 → K S φ as an example. We find observable effects at a level well beyond anything possible with CKM dynamics. Comparisons with CP violation in the K and B systems offer an excellent test of this scenario and reveal the specific pattern of flavour and CP violation in the D 0 -D-bar 0 system predicted by this model. We discuss a number of charm decays that could potentially offer an insight in the dynamics of CP violation in D decays. We also apply our formalism to B s -B-bar s mixing.

Pseudo-dynamic source modelling with 1-point and 2-point statistics of earthquake source parameters

KAUST Repository

Song, S. G.

2013-12-24

Ground motion prediction is an essential element in seismic hazard and risk analysis. Empirical ground motion prediction approaches have been widely used in the community, but efficient simulation-based ground motion prediction methods are needed to complement empirical approaches, especially in the regions with limited data constraints. Recently, dynamic rupture modelling has been successfully adopted in physics-based source and ground motion modelling, but it is still computationally demanding and many input parameters are not well constrained by observational data. Pseudo-dynamic source modelling keeps the form of kinematic modelling with its computational efficiency, but also tries to emulate the physics of source process. In this paper, we develop a statistical framework that governs the finite-fault rupture process with 1-point and 2-point statistics of source parameters in order to quantify the variability of finite source models for future scenario events. We test this method by extracting 1-point and 2-point statistics from dynamically derived source models and simulating a number of rupture scenarios, given target 1-point and 2-point statistics. We propose a new rupture model generator for stochastic source modelling with the covariance matrix constructed from target 2-point statistics, that is, auto- and cross-correlations. Our sensitivity analysis of near-source ground motions to 1-point and 2-point statistics of source parameters provides insights into relations between statistical rupture properties and ground motions. We observe that larger standard deviation and stronger correlation produce stronger peak ground motions in general. The proposed new source modelling approach will contribute to understanding the effect of earthquake source on near-source ground motion characteristics in a more quantitative and systematic way.

1,25D3 enhances antitumor activity of gemcitabine and cisplatin in human bladder cancer models

Science.gov (United States)

Ma, Yingyu; Yu, Wei-Dong; Trump, Donald L.; Johnson, Candace S.

2010-01-01

Background 1,25 dihydroxyvitamin D3 (1,25D3) potentiates the cytotoxic effects of several common chemotherapeutic agents. The combination of gemcitabine and cisplatin (GC) is a current standard chemotherapy regimen for bladder cancer. We investigated whether 1,25D3 could enhance the antitumor activity of GC in bladder cancer model systems. Methods Human bladder cancer T24 and UMUC3 cells were pretreated with 1,25D3 followed by GC. Apoptosis were assessed by annexin V staining. Caspase activation was examined by immunoblot analysis and substrate-based caspase activity assay. The cytotoxic effects were examined using MTT and in vitro clonogenic assay. p73 protein levels were assessed by immunoblot analysis. Knockdown of p73 was achieved by siRNA. The in vivo antitumor activity was assessed by in vivo excision clonogenic assay and tumor regrowth delay in the T24 xenograft model. Results 1,25D3 pretreatment enhanced GC-induced apoptosis and the activities of caspases- 8, 9 and 3 in T24 and UMUC3 cells. 1,25D3 synergistically reduced GC-suppressed surviving fraction in T24 cells. 1,25D3, gemcitabine, or cisplatin induced p73 accumulation, which was enhanced by GC or 1,25D3 and GC. p73 expression was lower in human primary bladder tumor tissue compared with adjacent normal tissue. Knockdown of p73 increased clonogenic capacity of T24 cells treated with 1,25D3, GC or 1,25D3 and GC. 1,25D3 and GC combination enhanced tumor regression compared with 1,25D3 or GC alone. Conclusions 1,25D3 potentiates GC-mediated growth inhibition in human bladder cancer models in vitro and in vivo, which involves p73 induction and apoptosis. PMID:20564622

Chaos in a dynamic model of urban transportation network flow based on user equilibrium states

International Nuclear Information System (INIS)

Xu Meng; Gao Ziyou

2009-01-01

In this study, we investigate the dynamical behavior of network traffic flow. We first build a two-stage mathematical model to analyze the complex behavior of network flow, a dynamical model, which is based on the dynamical gravity model proposed by Dendrinos and Sonis [Dendrinos DS, Sonis M. Chaos and social-spatial dynamic. Berlin: Springer-Verlag; 1990] is used to estimate the number of trips. Considering the fact that the Origin-Destination (O-D) trip cost in the traffic network is hard to express as a functional form, in the second stage, the user equilibrium network assignment model was used to estimate the trip cost, which is the minimum cost of used path when user equilibrium (UE) conditions are satisfied. It is important to use UE to estimate the O-D cost, since a connection is built among link flow, path flow, and O-D flow. The dynamical model describes the variations of O-D flows over discrete time periods, such as each day and each week. It is shown that even in a system with dimensions equal to two, chaos phenomenon still exists. A 'Chaos Propagation' phenomenon is found in the given model.

3D Spatially Resolved Models of the Intracellular Dynamics of the Hepatitis C Genome Replication Cycle

KAUST Repository

Knodel, Markus

2017-10-02

Mathematical models of virus dynamics have not previously acknowledged spatial resolution at the intracellular level despite substantial arguments that favor the consideration of intracellular spatial dependence. The replication of the hepatitis C virus (HCV) viral RNA (vRNA) occurs within special replication complexes formed from membranes derived from endoplasmatic reticulum (ER). These regions, termed membranous webs, are generated primarily through specific interactions between nonstructural virus-encoded proteins (NSPs) and host cellular factors. The NSPs are responsible for the replication of the vRNA and their movement is restricted to the ER surface. Therefore, in this study we developed fully spatio-temporal resolved models of the vRNA replication cycle of HCV. Our simulations are performed upon realistic reconstructed cell structures-namely the ER surface and the membranous webs-based on data derived from immunostained cells replicating HCV vRNA. We visualized 3D simulations that reproduced dynamics resulting from interplay of the different components of our models (vRNA, NSPs, and a host factor), and we present an evaluation of the concentrations for the components within different regions of the cell. Thus far, our model is restricted to an internal portion of a hepatocyte and is qualitative more than quantitative. For a quantitative adaption to complete cells, various additional parameters will have to be determined through further in vitro cell biology experiments, which can be stimulated by the results deccribed in the present study.

A Coupled Model of the 1D River Network and 3D Estuary Based on Hydrodynamics and Suspended Sediment Simulation

Directory of Open Access Journals (Sweden)

Wei Zhang

2014-01-01

Full Text Available River networks and estuaries are very common in coastal areas. Runoff from the upper stream interacts with tidal current from open sea in these two systems, leading to a complex hydrodynamics process. Therefore, it is necessary to consider the two systems as a whole to study the flow and suspended sediment transport. Firstly, a 1D model is established in the Pearl River network and a 3D model is applied in its estuary. As sufficient mass exchanges between the river network and its estuary, a strict mathematical relationship of water level at the interfaces can be adopted to couple the 1D model with the 3D model. By doing so, the coupled model does not need to have common nested grids. The river network exchanges the suspended sediment with its estuary by adding the continuity conditions at the interfaces. The coupled model is, respectively, calibrated in the dry season and the wet season. The results demonstrate that the coupled model works excellently in simulating water level and discharge. Although there are more errors in simulating suspended sediment concentration due to some reasons, the coupled model is still good enough to evaluate the suspended sediment transport in river network and estuary systems.

Badlands: A parallel basin and landscape dynamics model

Directory of Open Access Journals (Sweden)

T. Salles

2016-01-01

Full Text Available Over more than three decades, a number of numerical landscape evolution models (LEMs have been developed to study the combined effects of climate, sea-level, tectonics and sediments on Earth surface dynamics. Most of them are written in efficient programming languages, but often cannot be used on parallel architectures. Here, I present a LEM which ports a common core of accepted physical principles governing landscape evolution into a distributed memory parallel environment. Badlands (acronym for BAsin anD LANdscape DynamicS is an open-source, flexible, TIN-based landscape evolution model, built to simulate topography development at various space and time scales.

A high sensitive 66 dB linear dynamic range receiver for 3-D laser radar

Science.gov (United States)

Ma, Rui; Zheng, Hao; Zhu, Zhangming

2017-08-01

This study presents a CMOS receiver chip realized in 0.18 μm standard CMOS technology and intended for high precision 3-D laser radar. The chip includes an adjustable gain transimpedance pre-amplifier, a post-amplifier and two timing comparators. An additional feedback is employed in the regulated cascode transimpedance amplifier to decrease the input impedance, and a variable gain transimpedance amplifier controlled by digital switches and analog multiplexer is utilized to realize four gain modes, extending the input dynamic range. The measurement shows that the highest transimpedance of the channel is 50 k {{Ω }}, the uncompensated walk error is 1.44 ns in a wide linear dynamic range of 66 dB (1:2000), and the input referred noise current is 2.3 pA/\\sqrt{{Hz}} (rms), resulting in a very low detectable input current of 1 μA with SNR = 5.

1D models for condensation induced water hammer in pipelines

International Nuclear Information System (INIS)

Bloemeling, Frank; Neuhas, Thorsten; Schaffrath, Andreas

2013-01-01

Condensation induced water hammer (CIWH) are caused by contact of steam and subcooled water. Thus, modeling the direct contact condensation is a crucial step towards the simulation of condensation induced water hammer with 1D pressure surge codes. Therefore, also the TUeV NORD SysTec GmbH and Co. KG inhouse pressure surge code DYVRO has been equipped with a new contact condensation model. The validation of DYVRO against an experiment dealing with CIWH is presented in this contribution. (orig.)

1D models for condensation induced water hammer in pipelines

Energy Technology Data Exchange (ETDEWEB)

Bloemeling, Frank; Neuhas, Thorsten; Schaffrath, Andreas [TUEV NORD SysTec GmbH und Co. KG, Hamburg (Germany)

2013-03-15

Condensation induced water hammer (CIWH) are caused by contact of steam and subcooled water. Thus, modeling the direct contact condensation is a crucial step towards the simulation of condensation induced water hammer with 1D pressure surge codes. Therefore, also the TUeV NORD SysTec GmbH and Co. KG inhouse pressure surge code DYVRO has been equipped with a new contact condensation model. The validation of DYVRO against an experiment dealing with CIWH is presented in this contribution. (orig.)

1D-Var temperature retrievals from microwave radiometer and convective scale model

Directory of Open Access Journals (Sweden)

Pauline Martinet

2015-12-01

Full Text Available This paper studies the potential of ground-based microwave radiometers (MWR for providing accurate temperature retrievals by combining convective scale numerical models and brightness temperatures (BTs. A one-dimensional variational (1D-Var retrieval technique has been tested to optimally combine MWR and 3-h forecasts from the French convective scale model AROME. A microwave profiler HATPRO (Humidity and Temperature PROfiler was operated during 6 months at the meteorological station of Bordeaux (Météo France. MWR BTs were monitored against simulations from the Atmospheric Radiative Transfer Simulator 2 radiative transfer model. An overall good agreement was found between observations and simulations for opaque V-band channels but large errors were observed for channels the most affected by liquid water and water vapour emissions (51.26 and 52.28 GHz. 1D-Var temperature retrievals are performed in clear-sky and cloudy conditions using a screening procedure based on cloud base height retrieval from ceilometer observations, infrared radiometer temperature and liquid water path derived from the MWR observations. The 1D-Var retrievals were found to improve the AROME forecasts up to 2 km with a maximum gain of approximately 50 % in root-mean-square-errors (RMSE below 500 m. They were also found to outperform neural network retrievals. A static bias correction was proposed to account for systematic instrumental errors. This correction was found to have a negligible impact on the 1D-Var retrievals. The use of low elevation angles improves the retrievals up to 12 % in RMSE in cloudy-sky in the first layers. The present implementation achieved a RMSE with respect to radiosondes within 1 K in clear-sky and 1.3 K in cloudy-sky conditions for temperature.

From Particles and Point Clouds to Voxel Models: High Resolution Modeling of Dynamic Landscapes in Open Source GIS

Science.gov (United States)

Mitasova, H.; Hardin, E. J.; Kratochvilova, A.; Landa, M.

2012-12-01

Multitemporal data acquired by modern mapping technologies provide unique insights into processes driving land surface dynamics. These high resolution data also offer an opportunity to improve the theoretical foundations and accuracy of process-based simulations of evolving landforms. We discuss development of new generation of visualization and analytics tools for GRASS GIS designed for 3D multitemporal data from repeated lidar surveys and from landscape process simulations. We focus on data and simulation methods that are based on point sampling of continuous fields and lead to representation of evolving surfaces as series of raster map layers or voxel models. For multitemporal lidar data we present workflows that combine open source point cloud processing tools with GRASS GIS and custom python scripts to model and analyze dynamics of coastal topography (Figure 1) and we outline development of coastal analysis toolbox. The simulations focus on particle sampling method for solving continuity equations and its application for geospatial modeling of landscape processes. In addition to water and sediment transport models, already implemented in GIS, the new capabilities under development combine OpenFOAM for wind shear stress simulation with a new module for aeolian sand transport and dune evolution simulations. Comparison of observed dynamics with the results of simulations is supported by a new, integrated 2D and 3D visualization interface that provides highly interactive and intuitive access to the redesigned and enhanced visualization tools. Several case studies will be used to illustrate the presented methods and tools and demonstrate the power of workflows built with FOSS and highlight their interoperability.Figure 1. Isosurfaces representing evolution of shoreline and a z=4.5m contour between the years 1997-2011at Cape Hatteras, NC extracted from a voxel model derived from series of lidar-based DEMs.

Dynamic modelling of radionuclide uptake by Fukushima coastal biota - Dynamic modelling of radionuclide uptake by marine biota: application to Fukushima assessment

Energy Technology Data Exchange (ETDEWEB)

Vives i Batlle, Jordi [Belgian Nuclear Research Centre, Boeretang 200, 2400 Mol (Belgium)

2014-07-01

Radiological assessments to non-human marine biota are usually carried out by assuming that the activity concentration in an organism is proportional to the activity concentration in an adjacent volume of water, via a concentration factor (CF). It is also assumed that radionuclides in the water are in isotopic equilibrium with the sediments via a sediment distribution coefficient (K{sub d}). These assumptions are not valid in accidental situations where the biota and the sediments react with a time delay to large variations of activity concentration in seawater. A simple dynamic model was developed to factorise the dynamics of radionuclide uptake and turnover in biota and sediments, as determined by a balance between the residence time of radionuclides in seawater/sediments and the biological half-life of elimination in the biota. The model calculates activity concentration of {sup 131}I, {sup 134}Cs, {sup 137}Cs and {sup 90}Sr in seabed sediment, fish, crustaceans, molluscs and macro-algae from surrounding activity concentrations in seawater, with which to derive internal and external dose rates. A central element of this new model is the inclusion of sediment processes in dynamic transfer modelling. The model is adapted to include depletion of radionuclides adsorbed onto suspended particulates (particle scavenging), molecular diffusion, pore water mixing and bioturbation (modelled effectively as a diffusive process) represented by a simple set of differential equations that is coupled with the biological uptake/turnover processes. In this way, the model is capable of reproducing activity concentration in sediment to give a more realistic calculation of the external dose to biota compared with the simpler approach based on CF and K{sub d} values used in previous assessments. The model is applied to the assessment of the radiological impact of the Fukushima accident on marine biota in the early phase of the accident. It is shown that previous assessment of the

Modeling an electric motor in 1-D

Science.gov (United States)

Butler, Thomas G.

1991-01-01

Quite often the dynamicist will be faced with having an electric drive motor as a link in the elastic path of a structure such that the motor's characteristics must be taken into account to properly represent the dynamics of the primary structure. He does not want to model it so accurately that he could get detailed stress and displacements in the motor proper, but just sufficiently to represent its inertia loading and elastic behavior from its mounting bolts to its drive coupling. Described here is how the rotor and stator of such a motor can be adequately modeled as a colinear pair of beams.

1.4D quasistatic profile model of transport in a field-reversed configuration (FRC)

International Nuclear Information System (INIS)

Steinhauer, L.C.

1990-01-01

Global confinement models are useful for determining how a given transport mechanism (dependent on local parameters) translates into global confinement times. Such models are also useful for inferring the overall magnitudes of transport rates, and limited information about their spatial profiles. They are especially important in a field reversed configuration (FRC) where the equilibrium and transport rates are so intimately coupled, and where the flux loss time is such an important factor. An earlier global FRC confinement model, sometimes called QUASI, was based on the assumption of a quasi-steady equilibrium. The equilibrium was assumed to have square-ends with some features of 2D equilibria: (1) equal pressure on inner and outer branches of the flux lines; and (2) average-beta relation. Models of this type have been called 1-1/4D transport models. The same general approach has now been applied in a straightforward way to an equilibrium with realistic axial structure. This might be called a 1.4D quasisteady transport model. The assumed axial structure can be that of an analytic equilibrium, or a more complicated computed equilibrium as desired. The example used here is an elongated Hill's vortex equilibrium. As will be shown later, the equilibrium is reflected by two integral quantities that appear in the quasistatic diffusion equation

Application of 6D Building Information Model (6D BIM) for Business-storage Building in Slovenia

Science.gov (United States)

Pučko, Zoran; Vincek, Dražen; Štrukelj, Andrej; Šuman, Nataša

2017-10-01

The aim of this paper is to present an application of 6D building information modelling (6D BIM) on a real business-storage building in Slovenia. First, features of building maintenance in general are described according to the current Slovenian legislation, and also a general principle of BIM is given. After that, step-by-step activities for modelling 6D BIM are exposed, namely from Element list for maintenance, determination of their lifetime and service measures, cost analysing and time analysing to 6D BIM modelling. The presented 6D BIM model is designed in a unique way in which cost analysis is performed as 5D BIM model with linked data to use BIM Construction Project Management Software (Vico Office), integrated with 3D BIM model, whereas time analysis as 4D BIM model is carried out as non-linked data with the help of Excel (without connection to 3D BIM model). The paper is intended to serve as a guide to the building owners to prepare 6D BIM and to provide an insight into the relevant dynamic information about intervals and costs for execution of maintenance works in the whole building lifecycle.

3.5D dynamic PET image reconstruction incorporating kinetics-based clusters

International Nuclear Information System (INIS)

Lu Lijun; Chen Wufan; Karakatsanis, Nicolas A; Rahmim, Arman; Tang Jing

2012-01-01

Standard 3D dynamic positron emission tomographic (PET) imaging consists of independent image reconstructions of individual frames followed by application of appropriate kinetic model to the time activity curves at the voxel or region-of-interest (ROI). The emerging field of 4D PET reconstruction, by contrast, seeks to move beyond this scheme and incorporate information from multiple frames within the image reconstruction task. Here we propose a novel reconstruction framework aiming to enhance quantitative accuracy of parametric images via introduction of priors based on voxel kinetics, as generated via clustering of preliminary reconstructed dynamic images to define clustered neighborhoods of voxels with similar kinetics. This is then followed by straightforward maximum a posteriori (MAP) 3D PET reconstruction as applied to individual frames; and as such the method is labeled ‘3.5D’ image reconstruction. The use of cluster-based priors has the advantage of further enhancing quantitative performance in dynamic PET imaging, because: (a) there are typically more voxels in clusters than in conventional local neighborhoods, and (b) neighboring voxels with distinct kinetics are less likely to be clustered together. Using realistic simulated 11 C-raclopride dynamic PET data, the quantitative performance of the proposed method was investigated. Parametric distribution-volume (DV) and DV ratio (DVR) images were estimated from dynamic image reconstructions using (a) maximum-likelihood expectation maximization (MLEM), and MAP reconstructions using (b) the quadratic prior (QP-MAP), (c) the Green prior (GP-MAP) and (d, e) two proposed cluster-based priors (CP-U-MAP and CP-W-MAP), followed by graphical modeling, and were qualitatively and quantitatively compared for 11 ROIs. Overall, the proposed dynamic PET reconstruction methodology resulted in substantial visual as well as quantitative accuracy improvements (in terms of noise versus bias performance) for parametric DV

Dynamic concision for three-dimensional reconstruction of human organ built with virtual reality modelling language (VRML)*

OpenAIRE

Yu, Zheng-yang; Zheng, Shu-sen; Chen, Lei-ting; He, Xiao-qian; Wang, Jian-jun

2005-01-01

This research studies the process of 3D reconstruction and dynamic concision based on 2D medical digital images using virtual reality modelling language (VRML) and JavaScript language, with a focus on how to realize the dynamic concision of 3D medical model with script node and sensor node in VRML. The 3D reconstruction and concision of body internal organs can be built with such high quality that they are better than those obtained from the traditional methods. With the function of dynamic c...

Dynamic Simulation of 1D Cellular Automata in the Active aTAM.

Science.gov (United States)

Jonoska, Nataša; Karpenko, Daria; Seki, Shinnosuke

2015-07-01

The Active aTAM is a tile based model for self-assembly where tiles are able to transfer signals and change identities according to the signals received. We extend Active aTAM to include deactivation signals and thereby allow detachment of tiles. We show that the model allows a dynamic simulation of cellular automata with assemblies that do not record the entire computational history but only the current updates of the states, and thus provide a way for (a) algorithmic dynamical structural changes in the assembly and (b) reusable space in self-assembly. The simulation is such that at a given location the sequence of tiles that attach and detach corresponds precisely to the sequence of states the synchronous cellular automaton generates at that location.

Proposed Robot Scheme with 5 DoF and Dynamic Modelling Using Maple Software

Directory of Open Access Journals (Sweden)

Shala Ahmet

2017-11-01

Full Text Available In this paper is represented Dynamical Modelling of robots which is commonly first important step of Modelling, Analysis and Control of robotic systems. This paper is focused on using Denavit-Hartenberg (DH convention for kinematics and Newton-Euler Formulations for dynamic modelling of 5 DoF - Degree of Freedom of 3D robot. The process of deriving of dynamical model is done using Software Maple. Derived Dynamical Model of 5 DoF robot is converted for Matlab use for future analysis, control and simulations.

Potent neutralizing anti-CD1d antibody reduces lung cytokine release in primate asthma model.

Science.gov (United States)

Nambiar, Jonathan; Clarke, Adam W; Shim, Doris; Mabon, David; Tian, Chen; Windloch, Karolina; Buhmann, Chris; Corazon, Beau; Lindgren, Matilda; Pollard, Matthew; Domagala, Teresa; Poulton, Lynn; Doyle, Anthony G

2015-01-01

CD1d is a receptor on antigen-presenting cells involved in triggering cell populations, particularly natural killer T (NKT) cells, to release high levels of cytokines. NKT cells are implicated in asthma pathology and blockade of the CD1d/NKT cell pathway may have therapeutic potential. We developed a potent anti-human CD1d antibody (NIB.2) that possesses high affinity for human and cynomolgus macaque CD1d (KD ∼100 pM) and strong neutralizing activity in human primary cell-based assays (IC50 typically <100 pM). By epitope mapping experiments, we showed that NIB.2 binds to CD1d in close proximity to the interface of CD1d and the Type 1 NKT cell receptor β-chain. Together with data showing that NIB.2 inhibited stimulation via CD1d loaded with different glycolipids, this supports a mechanism whereby NIB.2 inhibits NKT cell activation by inhibiting Type 1 NKT cell receptor β-chain interactions with CD1d, independent of the lipid antigen in the CD1d antigen-binding cleft. The strong in vitro potency of NIB.2 was reflected in vivo in an Ascaris suum cynomolgus macaque asthma model. Compared with vehicle control, NIB.2 treatment significantly reduced bronchoalveolar lavage (BAL) levels of Ascaris-induced cytokines IL-5, IL-8 and IL-1 receptor antagonist, and significantly reduced baseline levels of GM-CSF, IL-6, IL-15, IL-12/23p40, MIP-1α, MIP-1β, and VEGF. At a cellular population level NIB.2 also reduced numbers of BAL lymphocytes and macrophages, and blood eosinophils and basophils. We demonstrate that anti-CD1d antibody blockade of the CD1d/NKT pathway modulates inflammatory parameters in vivo in a primate inflammation model, with therapeutic potential for diseases where the local cytokine milieu is critical.

Method of solving conformal models in D-dimensional space 2: A family of exactly solvable models in D > 2

International Nuclear Information System (INIS)

Fradkin, E.S.; Palchik, M.Ya.

1996-02-01

We study a family of exactly solvable models of conformally-invariant quantum field theory in D-dimensional space. We demonstrate the existence of D-dimensional analogs of primary and secondary fields. Under the action of energy-momentum tensor and conserved currents, the primary fields creates an infinite set of (tensor) secondary fields of different generations. The commutators of secondary fields with zero components of current and energy-momentum tensor include anomalous operator terms. We show that the Hilbert space of conformal theory has a special sector which structure is solely defined by the Ward identities independently on the choice of dynamical model. The states of this sector are constructed from secondary fields. Definite self-consistent conditions on the states of the latter sector fix the choice of the field model uniquely. In particular, Lagrangian models do belong to this class of models. The above self-consistent conditions are formulated as follows. Special superpositions Q s , s = 1,2,... of secondary fields are constructed. Each superposition is determined by the requirement that the form of its commutators with energy-momentum tensor and current (i.e. transformation properties) should be identical to that of a primary field. Each equation Q s (x) = 0 is consistent, and defines an exactly solvable model for D ≥ 3. The structure of these models are analogous to that of well-known two dimensional conformal models. The states Q s (x) modul 0> are analogous to the null-vectors of two dimensional theory. In each of these models one can obtain a closed set of differential equations for all the higher Green functions, as well as algebraic equations relating the scale dimension of fundamental field to the D-dimensional analog of a central charge. As an example, we present a detailed discussion of a pair of exactly solvable models in even-dimensional space D ≥ 4. (author). 28 refs

Non-parametric Bayesian models of response function in dynamic image sequences

Czech Academy of Sciences Publication Activity Database

Tichý, Ondřej; Šmídl, Václav

2016-01-01

Roč. 151, č. 1 (2016), s. 90-100 ISSN 1077-3142 R&D Projects: GA ČR GA13-29225S Institutional support: RVO:67985556 Keywords : Response function * Blind source separation * Dynamic medical imaging * Probabilistic models * Bayesian methods Subject RIV: BB - Applied Statistics, Operational Research Impact factor: 2.498, year: 2016 http://library.utia.cas.cz/separaty/2016/AS/tichy-0456983.pdf

Equivalent Dynamic Models.

Science.gov (United States)

Molenaar, Peter C M

2017-01-01

Equivalences of two classes of dynamic models for weakly stationary multivariate time series are discussed: dynamic factor models and autoregressive models. It is shown that exploratory dynamic factor models can be rotated, yielding an infinite set of equivalent solutions for any observed series. It also is shown that dynamic factor models with lagged factor loadings are not equivalent to the currently popular state-space models, and that restriction of attention to the latter type of models may yield invalid results. The known equivalent vector autoregressive model types, standard and structural, are given a new interpretation in which they are conceived of as the extremes of an innovating type of hybrid vector autoregressive models. It is shown that consideration of hybrid models solves many problems, in particular with Granger causality testing.

GPU-accelerated 3-D model-based tracking

International Nuclear Information System (INIS)

Brown, J Anthony; Capson, David W

2010-01-01

Model-based approaches to tracking the pose of a 3-D object in video are effective but computationally demanding. While statistical estimation techniques, such as the particle filter, are often employed to minimize the search space, real-time performance remains unachievable on current generation CPUs. Recent advances in graphics processing units (GPUs) have brought massively parallel computational power to the desktop environment and powerful developer tools, such as NVIDIA Compute Unified Device Architecture (CUDA), have provided programmers with a mechanism to exploit it. NVIDIA GPUs' single-instruction multiple-thread (SIMT) programming model is well-suited to many computer vision tasks, particularly model-based tracking, which requires several hundred 3-D model poses to be dynamically configured, rendered, and evaluated against each frame in the video sequence. Using 6 degree-of-freedom (DOF) rigid hand tracking as an example application, this work harnesses consumer-grade GPUs to achieve real-time, 3-D model-based, markerless object tracking in monocular video.

On the Dynamics of the Fermi-Bose model

DEFF Research Database (Denmark)

Ögren, Magnus

In this talk we formulate and prove results for the exponential matrix representing the dynamics of the Fermi-Bose model in an undepleted bosonic field approximation. A recent application of this model is molecular dimmers dissociating into its atomic compounds. The problem is solved in D spatial....... In particular the results can be used for studies of threedimensional physical systems of arbitrary geometry. We illustrate the generality of our approach by giving numerical results for the dynamics of Glauber type atomic pair correlation functions for a non-isotropic three-dimensional harmonically trapped...

A Dynamic Model of Sustainment Investment

Science.gov (United States)

2015-02-01

Sustainment System Dynamics Model 11 Figure 7: Core Structure of Sustainment Work 12 Figure 8: Bandwagon Effect Loop 13 Figure 9: Limits to Growth Loop 14...Dynamics Model sustainment capacity sustainment performance gap Bandwagon Effect R1 Limits to Growth B1 S Work Smarter B3 Work Bigger B2 desired...which is of concern primarily when using the model as a vehicle for research. Figure 8 depicts a reinforcing loop called the “ Bandwagon Effect

Interim results of the sixth three-dimensional AER dynamic benchmark problem calculation. Solution of problem with DYN3D and RELAP5-3D codes

International Nuclear Information System (INIS)

Hadek, J.; Kral, P.; Macek, J.

2001-01-01

The paper gives a brief survey of the 6 th three-dimensional AER dynamic benchmark calculation results received with the codes DYN3D and RELAPS-3D at NRI Rez. This benchmark was defined at the 10 th AER Symposium. Its initiating event is a double ended break in the steam line of steam generator No. I in a WWER-440/213 plant at the end of the first fuel cycle and in hot full power conditions. Stationary and burnup calculations as well as tuning of initial state before the transient were performed with the code DYN3D. Transient calculations were made with the system code RELAPS-3D.The KASSETA library was used for the generation of reactor core neutronic parameters. The detailed six loops model of NPP Dukovany was adopted for the 6 th AER dynamic benchmark purposes. The RELAPS-3D full core neutronic model was connected with seven coolant channels thermal-hydraulic model of the core (Authors)

Crossed Beam Scattering Studies of the Dynamics of the Reaction NH.sup.+./sup.(D.sub.2./sub.; H, D) NHD.sup.+./sup., ND.sub.2./sub..sup.+./sup..

Czech Academy of Sciences Publication Activity Database

Wittemann, A.; Sadílek, Martin; Polách, Jiří; Teloy, E.; Schlier, CH.; Herman, Zdeněk

2001-01-01

Roč. 215, č. 6 (2001), s. 749-764 ISSN 0942-9352 R&D Projects: GA ČR GA203/00/0632 Institutional research plan: CEZ:AV0Z4040901 Keywords : scattering studies * Klots model * dynamics Subject RIV: CF - Physical ; Theoretical Chemistry Impact factor: 1.070, year: 2001

Event-chain algorithm for the Heisenberg model: Evidence for z≃1 dynamic scaling.

Science.gov (United States)

Nishikawa, Yoshihiko; Michel, Manon; Krauth, Werner; Hukushima, Koji

2015-12-01

We apply the event-chain Monte Carlo algorithm to the three-dimensional ferromagnetic Heisenberg model. The algorithm is rejection-free and also realizes an irreversible Markov chain that satisfies global balance. The autocorrelation functions of the magnetic susceptibility and the energy indicate a dynamical critical exponent z≈1 at the critical temperature, while that of the magnetization does not measure the performance of the algorithm. We show that the event-chain Monte Carlo algorithm substantially reduces the dynamical critical exponent from the conventional value of z≃2.

Debris Dispersion Model Using Java 3D

Science.gov (United States)

Thirumalainambi, Rajkumar; Bardina, Jorge

2004-01-01

This paper describes web based simulation of Shuttle launch operations and debris dispersion. Java 3D graphics provides geometric and visual content with suitable mathematical model and behaviors of Shuttle launch. Because the model is so heterogeneous and interrelated with various factors, 3D graphics combined with physical models provides mechanisms to understand the complexity of launch and range operations. The main focus in the modeling and simulation covers orbital dynamics and range safety. Range safety areas include destruct limit lines, telemetry and tracking and population risk near range. If there is an explosion of Shuttle during launch, debris dispersion is explained. The shuttle launch and range operations in this paper are discussed based on the operations from Kennedy Space Center, Florida, USA.

Quasi 2D hydrodynamic modelling of the flooded hinterland due to dyke breaching on the Elbe River

Directory of Open Access Journals (Sweden)

S. Huang

2007-01-01

Full Text Available In flood modeling, many 1D and 2D combination and 2D models are used to simulate diversion of water from rivers through dyke breaches into the hinterland for extreme flood events. However, these models are too demanding in data requirements and computational resources which is an important consideration when uncertainty analysis using Monte Carlo techniques is used to complement the modeling exercise. The goal of this paper is to show the development of a quasi-2D modeling approach, which still calculates the dynamic wave in 1D but the discretisation of the computational units are in 2D, allowing a better spatial representation of the flow in the hinterland due to dyke breaching without a large additional expenditure on data pre-processing and computational time. A 2D representation of the flow and velocity fields is required to model sediment and micro-pollutant transport. The model DYNHYD (1D hydrodynamics from the WASP5 modeling package was used as a basis for the simulations. The model was extended to incorporate the quasi-2D approach and a Monte-Carlo Analysis was used to conduct a flood sensitivity analysis to determine the sensitivity of parameters and boundary conditions to the resulting water flow. An extreme flood event on the Elbe River, Germany, with a possible dyke breach area was used as a test case. The results show a good similarity with those obtained from another 1D/2D modeling study.

Computational fluid dynamics modeling of mixed convection flows in buildings enclosures

Energy Technology Data Exchange (ETDEWEB)

Kayne, Alexander; Agarwal, Ramesh K. [Department of Mechanical Engineering and Materials Science, Washington University, St. Louis, MO 63130 (United States)

2013-07-01

In recent years Computational Fluid Dynamics (CFD) simulations are increasingly used to model the air circulation and temperature environment inside the rooms of residential and office buildings to gain insight into the relative energy consumptions of various HVAC systems for cooling/heating for climate control and thermal comfort. This requires accurate simulation of turbulent flow and heat transfer for various types of ventilation systems using the Reynolds-Averaged Navier-Stokes (RANS) equations of fluid dynamics. Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS) of Navier-Stokes equations is computationally intensive and expensive for simulations of this kind. As a result, vast majority of CFD simulations employ RANS equations in conjunction with a turbulence model. In order to assess the modeling requirements (mesh, numerical algorithm, turbulence model etc.) for accurate simulations, it is critical to validate the calculations against the experimental data. For this purpose, we use three well known benchmark validation cases, one for natural convection in 2D closed vertical cavity, second for forced convection in a 2D rectangular cavity and the third for mixed convection in a 2D square cavity. The simulations are performed on a number of meshes of different density using a number of turbulence models. It is found that k-epsilon two-equation turbulence model with a second-order algorithm on a reasonable mesh gives the best results. This information is then used to determine the modeling requirements (mesh, numerical algorithm, turbulence model etc.) for flows in 3D enclosures with different ventilation systems. In particular two cases are considered for which the experimental data is available. These cases are (1) air flow and heat transfer in a naturally ventilated room and (2) airflow and temperature distribution in an atrium. Good agreement with the experimental data and computations of other investigators is obtained.

Computational Modelling of Piston Ring Dynamics in 3D

Directory of Open Access Journals (Sweden)

Dlugoš Jozef

2014-12-01

Full Text Available Advanced computational models of a piston assembly based on the level of virtual prototypes require a detailed description of piston ring behaviour. Considering these requirements, the piston rings operate in regimes that cannot, in general, be simplified into an axisymmetric model. The piston and the cylinder liner do not have a perfect round shape, mainly due to machining tolerances and external thermo-mechanical loads. If the ring cannot follow the liner deformations, a local loss of contact occurs resulting in blow-by and increased consumption of lubricant oil in the engine. Current computational models are unable to implement such effects. The paper focuses on the development of a flexible 3D piston ring model based on the Timoshenko beam theory using the multibody system (MBS. The MBS model is compared to the finite element method (FEM solution.

Imaging of aortic valve dynamics in 4D OCT

Directory of Open Access Journals (Sweden)

Schnabel Christian

2015-09-01

Full Text Available The mechanical components of the heart, especially the valves and leaflets, are enormous stressed during lifetime. Therefore, those structures undergo different pathophysiological tissue transformations which affect cardiac output and in consequence living comfort of affected patients. These changes may lead to calcific aortic valve stenosis (AVS, the major heart valve disease in humans. The knowledge about changes of the dynamic behaviour during the course of this disease and the possibility of early stage diagnosis is of particular interest and could lead to the development of new treatment strategies and drug based options of prevention or therapy. 4D optical coherence tomography (OCT in combination with high-speed video microscopy were applied to characterize dynamic behaviour of the murine aortic valve and to characterize dynamic properties during artificial stimulation. We present a promising tool to investigate the aortic valve dynamics in an ex vivo disease model with a high spatial and temporal resolution using a multimodal imaging setup.

Nonlinear 2D arm dynamics in response to continuous and pulse-shaped force perturbations.

Science.gov (United States)

Happee, Riender; de Vlugt, Erwin; van Vliet, Bart

2015-01-01

Ample evidence exists regarding the nonlinearity of the neuromuscular system but linear models are widely applied to capture postural dynamics. This study quantifies the nonlinearity of human arm postural dynamics applying 2D continuous force perturbations (0.2-40 Hz) inducing three levels of hand displacement (5, 15, 45 mm RMS) followed by force-pulse perturbations inducing large hand displacements (up to 250 mm) in a position task (PT) and a relax task (RT) recording activity of eight shoulder and elbow muscles. The continuous perturbation data were used to analyze the 2D endpoint dynamics in the frequency domain and to identify reflexive and intrinsic parameters of a linear neuromuscular shoulder-elbow model. Subsequently, it was assessed to what extent the large displacements in response to force pulses could be predicted from the 'small amplitude' linear neuromuscular model. Continuous and pulse perturbation responses with varying amplitudes disclosed highly nonlinear effects. In PT, a larger continuous perturbation induced stiffening with a factor of 1.5 attributed to task adaptation evidenced by increased co-contraction and reflexive activity. This task adaptation was even more profound in the pulse responses where reflexes and displacements were strongly affected by the presence and amplitude of preceding continuous perturbations. In RT, a larger continuous perturbation resulted in yielding with a factor of 3.8 attributed to nonlinear mechanical properties as no significant reflexive activity was found. Pulse perturbations always resulted in yielding where a model fitted to the preceding 5-mm continuous perturbations predicted only 37% of the recorded peak displacements in RT and 79% in PT. This demonstrates that linear neuromuscular models, identified using continuous perturbations with small amplitudes, strongly underestimate displacements in pulse-shaped (e.g., impact) loading conditions. The data will be used to validate neuromuscular models including

Modelization and numerical simulation of atmospheric aerosols dynamics

International Nuclear Information System (INIS)

Debry, Edouard

2004-01-01

Chemical-transport models are now able to describe in a realistic way gaseous pollutants behavior in the atmosphere. Nevertheless atmospheric pollution also exists as a fine suspended particles, called aerosols which interact with gaseous phase, solar radiation, and have their own dynamic behavior. The goal of this thesis is the modelization and numerical simulation of the General Dynamic Equation of aerosols (GDE). Part I deals with some theoretical aspects of aerosol modelization. Part II is dedicated to the building of one size resolved aerosol model (SIREAM). In part III we perform the reduction of this model in order to use it in dispersion models as POLAIR3D. Several modelization issues are still opened: organic aerosol matter, externally mixed aerosols, coupling with turbulent mixing, and nano-particles. (author) [fr

CERN LEP2 constraint on 4D QED having a dynamically generated spatial dimension

International Nuclear Information System (INIS)

Cho, G.-C.; Izumi, Etsuko; Sugamoto, Akio

2002-01-01

We study 4D QED in which one spatial dimension is dynamically generated from the 3D action, following the mechanism proposed by Arkani-Hamed, Cohen, and Georgi. In this model, the generated fourth dimension is discretized by an interval parameter a. We examine the phenomenological constraint on the parameter a coming from collider experiments on the QED process e + e - →γγ. It is found that the CERN e + e - collider LEP2 experiments give the constraint of 1/a > or approx. 461 GeV. The expected bound on the same parameter a at a future e + e - linear collider is briefly discussed

Dynamics-based selective 2D 1H/1H chemical shift correlation spectroscopy under ultrafast MAS conditions

Science.gov (United States)

Zhang, Rongchun; Ramamoorthy, Ayyalusamy

2015-05-01

Dynamics plays important roles in determining the physical, chemical, and functional properties of a variety of chemical and biological materials. However, a material (such as a polymer) generally has mobile and rigid regions in order to have high strength and toughness at the same time. Therefore, it is difficult to measure the role of mobile phase without being affected by the rigid components. Herein, we propose a highly sensitive solid-state NMR approach that utilizes a dipolar-coupling based filter (composed of 12 equally spaced 90° RF pulses) to selectively measure the correlation of 1H chemical shifts from the mobile regions of a material. It is interesting to find that the rotor-synchronized dipolar filter strength decreases with increasing inter-pulse delay between the 90° pulses, whereas the dipolar filter strength increases with increasing inter-pulse delay under static conditions. In this study, we also demonstrate the unique advantages of proton-detection under ultrafast magic-angle-spinning conditions to enhance the spectral resolution and sensitivity for studies on small molecules as well as multi-phase polymers. Our results further demonstrate the use of finite-pulse radio-frequency driven recoupling pulse sequence to efficiently recouple weak proton-proton dipolar couplings in the dynamic regions of a molecule and to facilitate the fast acquisition of 1H/1H correlation spectrum compared to the traditional 2D NOESY (Nuclear Overhauser effect spectroscopy) experiment. We believe that the proposed approach is beneficial to study mobile components in multi-phase systems, such as block copolymers, polymer blends, nanocomposites, heterogeneous amyloid mixture of oligomers and fibers, and other materials.

Dynamics-based selective 2D 1H/1H chemical shift correlation spectroscopy under ultrafast MAS conditions

International Nuclear Information System (INIS)

Zhang, Rongchun; Ramamoorthy, Ayyalusamy

2015-01-01

Dynamics plays important roles in determining the physical, chemical, and functional properties of a variety of chemical and biological materials. However, a material (such as a polymer) generally has mobile and rigid regions in order to have high strength and toughness at the same time. Therefore, it is difficult to measure the role of mobile phase without being affected by the rigid components. Herein, we propose a highly sensitive solid-state NMR approach that utilizes a dipolar-coupling based filter (composed of 12 equally spaced 90° RF pulses) to selectively measure the correlation of 1 H chemical shifts from the mobile regions of a material. It is interesting to find that the rotor-synchronized dipolar filter strength decreases with increasing inter-pulse delay between the 90° pulses, whereas the dipolar filter strength increases with increasing inter-pulse delay under static conditions. In this study, we also demonstrate the unique advantages of proton-detection under ultrafast magic-angle-spinning conditions to enhance the spectral resolution and sensitivity for studies on small molecules as well as multi-phase polymers. Our results further demonstrate the use of finite-pulse radio-frequency driven recoupling pulse sequence to efficiently recouple weak proton-proton dipolar couplings in the dynamic regions of a molecule and to facilitate the fast acquisition of 1 H/ 1 H correlation spectrum compared to the traditional 2D NOESY (Nuclear Overhauser effect spectroscopy) experiment. We believe that the proposed approach is beneficial to study mobile components in multi-phase systems, such as block copolymers, polymer blends, nanocomposites, heterogeneous amyloid mixture of oligomers and fibers, and other materials

Open-end tube dynamic flow model with an oscillatory extortion

Directory of Open Access Journals (Sweden)

Tulwin Tytus

2017-01-01

Full Text Available This paper presents a derivation of dynamic 2d mathematical model for open end tube with oscillatory extortion in the region of the closed end. The aim the research is to investigate possible uses of the increased pressure in the enclosed tube chamber, especially for energy efficient lift generation. The mathematical model allows to test and predict how flow modifications impact the resultant lifting force. A derivation of the proposed mathematical model is shown. The mathematical model is then compared to the computational fluid dynamics discrete model. The results prove the accuracy of the mathematical physical model.

Dynamics of the two-spin spin-boson model with a common bath

Energy Technology Data Exchange (ETDEWEB)

Deng, Tianrui [Division of Materials Science, Nanyang Technological University, Singapore 639798 (Singapore); Centre for Optical and Electromagnetic Research, Zhejiang Provincial Key Laboratory for Sensing Technologies, Zhejiang University, Hangzhou 310058 (China); Yan, Yiying; Chen, Lipeng; Zhao, Yang, E-mail: YZhao@ntu.edu.sg [Division of Materials Science, Nanyang Technological University, Singapore 639798 (Singapore)

2016-04-14

Dynamics of the two-spin spin-boson model in the presence of Ohmic and sub-Ohmic baths is investigated by employing a multitude of the Davydov D{sub 1} trial states, also known as the multi-D{sub 1} Ansatz. Its accuracy in dynamics simulations of the two-spin SBM is improved significantly over the single D{sub 1} Ansatz, especially in the weak to moderately strong coupling regime. Validity of the multi-D{sub 1} Ansatz for various coupling strengths is also systematically examined by making use of the deviation vector which quantifies how faithfully the trial state obeys the Schrödinger equation. The time evolution of population difference and entanglement has been studied for various initial conditions and coupling strengths. Careful comparisons are carried out between our approach and three other methods, i.e., the time-dependent numerical renormalization group (TD-NRG) approach, the Bloch-Redfield theory, and a method based on a variational master equation. For strong coupling, the multi-D{sub 1} trial state yields consistent results as the TD-NRG approach in the Ohmic regime while the two disagree in the sub-Ohmic regime, where the multi-D{sub 1} trial state is shown to be more accurate. For weak coupling, the multi-D{sub 1} trial state agrees with the two master-equation methods in the presence of both Ohmic and sub-Ohmic baths, but shows considerable differences with the TD-NRG approach in the presence of a sub-Ohmic bath, calling into question the validity of the TD-NRG results at long times in the literature.

Dynamical simulation of a linear sigma model near the critical point

Energy Technology Data Exchange (ETDEWEB)

Wesp, Christian; Meistrenko, Alex; Greiner, Carsten [Institut fuer Theoretische Physik, Goethe-Universitaet Frankfurt, Max-von-Laue-Strasse 1, D-60438 Frankfurt (Germany); Hees, Hendrik van [Frankfurt Institute for Advanced Studies, Ruth-Moufang-Strasse 1, D-60438 Frankfurt (Germany)

2014-07-01

The intention of this study is the search for signatures of the chiral phase transition. To investigate the impact of fluctuations, e.g. of the baryon number, on the transition or a critical point, the linear sigma model is treated in a dynamical 3+1D numerical simulation. Chiral fields are approximated as classical fields, quarks are described by quasi particles in a Vlasov equation. Additional dynamic is implemented by quark-quark and quark-sigma-field interaction. For a consistent description of field-particle interactions, a new Monte-Carlo-Langevin-like formalism has been developed and is discussed.

Traveling waves in a continuum model of 1D schools

Science.gov (United States)

Oza, Anand; Kanso, Eva; Shelley, Michael

2017-11-01

We construct and analyze a continuum model of a 1D school of flapping swimmers. Our starting point is a delay differential equation that models the interaction between a swimmer and its upstream neighbors' wakes, which is motivated by recent experiments in the Applied Math Lab at NYU. We coarse-grain the evolution equations and derive PDEs for the swimmer density and variables describing the upstream wake. We study the equations both analytically and numerically, and find that a uniform density of swimmers destabilizes into a traveling wave. Our model makes a number of predictions about the properties of such traveling waves, and sheds light on the role of hydrodynamics in mediating the structure of swimming schools.

Coupling population dynamics with earth system models: the POPEM model.

Science.gov (United States)

Navarro, Andrés; Moreno, Raúl; Jiménez-Alcázar, Alfonso; Tapiador, Francisco J

2017-09-16

Precise modeling of CO 2 emissions is important for environmental research. This paper presents a new model of human population dynamics that can be embedded into ESMs (Earth System Models) to improve climate modeling. Through a system dynamics approach, we develop a cohort-component model that successfully simulates historical population dynamics with fine spatial resolution (about 1°×1°). The population projections are used to improve the estimates of CO 2 emissions, thus transcending the bulk approach of existing models and allowing more realistic non-linear effects to feature in the simulations. The module, dubbed POPEM (from Population Parameterization for Earth Models), is compared with current emission inventories and validated against UN aggregated data. Finally, it is shown that the module can be used to advance toward fully coupling the social and natural components of the Earth system, an emerging research path for environmental science and pollution research.

Dynamic concision for three-dimensional reconstruction of human organ built with virtual reality modelling language (VRML)*

Science.gov (United States)

Yu, Zheng-yang; Zheng, Shu-sen; Chen, Lei-ting; He, Xiao-qian; Wang, Jian-jun

2005-01-01

This research studies the process of 3D reconstruction and dynamic concision based on 2D medical digital images using virtual reality modelling language (VRML) and JavaScript language, with a focus on how to realize the dynamic concision of 3D medical model with script node and sensor node in VRML. The 3D reconstruction and concision of body internal organs can be built with such high quality that they are better than those obtained from the traditional methods. With the function of dynamic concision, the VRML browser can offer better windows for man-computer interaction in real-time environment than ever before. 3D reconstruction and dynamic concision with VRML can be used to meet the requirement for the medical observation of 3D reconstruction and have a promising prospect in the fields of medical imaging. PMID:15973760

Local and regional minimum 1D models for earthquake location and data quality assessment in complex tectonic regions: application to Switzerland

International Nuclear Information System (INIS)

Husen, S.; Clinton, J. F.; Kissling, E.

2011-01-01

One-dimensional (1D) velocity models are still widely used for computing earthquake locations at seismological centers or in regions where three-dimensional (3D) velocity models are not available due to the lack of data of sufficiently high quality. The concept of the minimum 1D model with appropriate station corrections provides a framework to compute initial hypocenter locations and seismic velocities for local earthquake tomography. Since a minimum 1D model represents a solution to the coupled hypocenter-velocity problem it also represents a suitable velocity model for earthquake location and data quality assessment, such as evaluating the consistency in assigning pre-defined weighting classes and average picking error. Nevertheless, the use of a simple 1D velocity structure in combination with station delays raises the question of how appropriate the minimum 1D model concept is when applied to complex tectonic regions with significant three-dimensional (3D) variations in seismic velocities. In this study we compute one regional minimum 1D model and three local minimum 1D models for selected subregions of the Swiss Alpine region, which exhibits a strongly varying Moho topography. We compare the regional and local minimum 1D models in terms of earthquake locations and data quality assessment to measure their performance. Our results show that the local minimum 1D models provide more realistic hypocenter locations and better data fits than a single model for the Alpine region. We attribute this to the fact that in a local minimum 1D model local and regional effects of the velocity structure can be better separated. Consequently, in tectonically complex regions, minimum 1D models should be computed in sub-regions defined by similar structure, if they are used for earthquake location and data quality assessment. (authors)

Local and regional minimum 1D models for earthquake location and data quality assessment in complex tectonic regions: application to Switzerland

Energy Technology Data Exchange (ETDEWEB)

Husen, S.; Clinton, J. F. [Swiss Seismological Service, ETH Zuerich, Zuerich (Switzerland); Kissling, E. [Institute of Geophysics, ETH Zuerich, Zuerich (Switzerland)

2011-10-15

One-dimensional (1D) velocity models are still widely used for computing earthquake locations at seismological centers or in regions where three-dimensional (3D) velocity models are not available due to the lack of data of sufficiently high quality. The concept of the minimum 1D model with appropriate station corrections provides a framework to compute initial hypocenter locations and seismic velocities for local earthquake tomography. Since a minimum 1D model represents a solution to the coupled hypocenter-velocity problem it also represents a suitable velocity model for earthquake location and data quality assessment, such as evaluating the consistency in assigning pre-defined weighting classes and average picking error. Nevertheless, the use of a simple 1D velocity structure in combination with station delays raises the question of how appropriate the minimum 1D model concept is when applied to complex tectonic regions with significant three-dimensional (3D) variations in seismic velocities. In this study we compute one regional minimum 1D model and three local minimum 1D models for selected subregions of the Swiss Alpine region, which exhibits a strongly varying Moho topography. We compare the regional and local minimum 1D models in terms of earthquake locations and data quality assessment to measure their performance. Our results show that the local minimum 1D models provide more realistic hypocenter locations and better data fits than a single model for the Alpine region. We attribute this to the fact that in a local minimum 1D model local and regional effects of the velocity structure can be better separated. Consequently, in tectonically complex regions, minimum 1D models should be computed in sub-regions defined by similar structure, if they are used for earthquake location and data quality assessment. (authors)

A 1-V 60-μW 85-dB dynamic range continuous-time third-order sigma-delta modulator

International Nuclear Information System (INIS)

Li Yuanwen; Qi Da; Dong Yifeng; Xu Jun; Ren Junyan

2009-01-01

A 1-V third order one-bit continuous-time (CT) EA modulator is presented. Designed in the SMIC mixed-signal 0.13-μm CMOS process, the modulator utilizes active RC integrators to implement the loop filter. An efficient circuit design methodology for the CT ΣΔ modulator is proposed and verified. Low power dissipation is achieved through the use of two-stage class A/AB amplifiers. The presented modulator achieves 81.4-dB SNDR and 85-dB dynamic range in a 20-kHz bandwidth with an over sampling ratio of 128. The total power consumption of the modulator is only 60 μW from a 1-V power supply and the prototype occupies an active area of 0.12 mm 2 . (semiconductor integrated circuits)

Experimental study and modelling of deuterium thermal release from Be-D co-deposited layers

Science.gov (United States)

Baldwin, M. J.; Schwarz-Selinger, T.; Doerner, R. P.

2014-07-01

A study of the thermal desorption of deuterium from 1 µm thick co-deposited Be-(0.1)D layers formed at 330 K by a magnetron sputtering technique is reported. A range of thermal desorption rates 0 ⩽ β ⩽ 1.0 K s-1 are explored with a view to studying the effectiveness of the proposed ITER wall and divertor bake procedure (β = 0 K s-1) to be carried out at 513 and 623 K. Fixed temperature bake durations up to 24 h are examined. The experimental thermal release data are used to validate a model input into the Tritium Migration and Analysis Program (TMAP-7). Good agreement with experiment is observed for a TMAP-7 model incorporating trap populations of activation energies for D release of 0.80 and 0.98 eV, and a dynamically computed surface D atomic to molecular recombination rate.

NoSQL Based 3D City Model Management System

Science.gov (United States)

Mao, B.; Harrie, L.; Cao, J.; Wu, Z.; Shen, J.

2014-04-01

To manage increasingly complicated 3D city models, a framework based on NoSQL database is proposed in this paper. The framework supports import and export of 3D city model according to international standards such as CityGML, KML/COLLADA and X3D. We also suggest and implement 3D model analysis and visualization in the framework. For city model analysis, 3D geometry data and semantic information (such as name, height, area, price and so on) are stored and processed separately. We use a Map-Reduce method to deal with the 3D geometry data since it is more complex, while the semantic analysis is mainly based on database query operation. For visualization, a multiple 3D city representation structure CityTree is implemented within the framework to support dynamic LODs based on user viewpoint. Also, the proposed framework is easily extensible and supports geoindexes to speed up the querying. Our experimental results show that the proposed 3D city management system can efficiently fulfil the analysis and visualization requirements.

Modelling the WWER-type reactor dynamics using a hybrid computer. Part 1

International Nuclear Information System (INIS)

Karpeta, C.

Results of simulation studies into reactor and steam generator dynamics of a WWER type power plant are presented. Spatial kinetics of the reactor core is described by a nodal approximation to diffusion equations, xenon poisoning equations and heat transfer equations. The simulation of the reactor model dynamics was performed on a hybrid computer. Models of both a horizontal and a vertical steam generator were developed. The dynamics was investigated over a large range of power by computing the transients on a digital computer. (author)

System dynamics and control with bond graph modeling

CERN Document Server

Kypuros, Javier

2013-01-01

Part I Dynamic System ModelingIntroduction to System DynamicsIntroductionSystem Decomposition and Model ComplexityMathematical Modeling of Dynamic SystemsAnalysis and Design of Dynamic SystemsControl of Dynamic SystemsDiagrams of Dynamic SystemsA Graph-Centered Approach to ModelingSummaryPracticeExercisesBasic Bond Graph ElementsIntroductionPower and Energy VariablesBasic 1-Port ElementsBasic 2-Ports ElementsJunction ElementsSimple Bond Graph ExamplesSummaryPracticeExercisesBond Graph Synthesis and Equation DerivationIntroductionGeneral GuidelinesMechanical TranslationMechanical RotationElectrical CircuitsHydraulic CircuitsMixed SystemsState Equation DerivationState-Space RepresentationsAlgebraic Loops and Derivative CausalitySummaryPracticeExercisesImpedance Bond GraphsIntroductionLaplace Transform of the State-Space EquationBasic 1-Port ImpedancesImpedance Bond Graph SynthesisJunctions, Transformers, and GyratorsEffort and Flow DividersSign ChangesTransfer Function DerivationAlternative Derivation of Transf...

Three-Dimensional Model Retrieval Using Dynamic Multi-Descriptor Fusion

Institute of Scientific and Technical Information of China (English)

Jau-Ling Shi; Chang-Hsing Lee; Yao-Wen Hou; Po-Ting Yeh

2017-01-01

In this paper, we propose a dynamic multi-descriptor fusion (DMDF) approach to improving the retrieval accuracy of 3-dimensional (3D) model retrieval systems. First, an independent retrieval list is generated by using each individual descriptor. Second, we propose an automatic relevant/irrelevant models selection (ARMS) approach to selecting the relevant and irrelevant 3D models automatically without any user interaction. A weighted distance, in which the weight associated with each individual descriptor is learnt by using the selected relevant and irrelevant models, is used to measure the similarity between two 3D models. Furthermore, a descriptor-dependent adaptive query point movement (AQPM) approach is employed to update every feature vector. This set of new feature vectors is used to index 3D models in the next search process. Four 3D model databases are used to compare the retrieval accuracy of our proposed DMDF approach with several descriptors as well as some well-known information fusion methods. Experimental results have shown that our proposed DMDF approach provides a promising retrieval result and always yields the best retrieval accuracy.

Modeling and measurement of the motion of the DIII-D vacuum vessel during vertical instabilities

International Nuclear Information System (INIS)

Reis, E.; Blevins, R.D.; Jensen, T.H.; Luxon, J.L.; Petersen, P.I.; Strait, E.J.

1991-11-01

The motions of the D3-D vacuum vessel during vertical instabilities of elongated plasmas have been measured and studied over the past five years. The currents flowing in the vessel wall and the plasma scrapeoff layer were also measured and correlated to a physics model. These results provide a time history load distribution on the vessel which were input to a dynamic analysis for correlation to the measured motions. The structural model of the vessel using the loads developed from the measured vessel currents showed that the calculated displacement history correlated well with the measured values. The dynamic analysis provides a good estimate of the stresses and the maximum allowable deflection of the vessel. In addition, the vessel motions produce acoustic emissions at 21 Hertz that are sufficiently loud to be felt as well as heard by the D3-D operators. Time history measurements of the sounds were correlated to the vessel displacements. An analytical model of an oscillating sphere provided a reasonable correlation to the amplitude of the measured sounds. The correlation of the theoretical and measured vessel currents, the dynamic measurements and analysis, and the acoustic measurements and analysis show that: (1) The physics model can predict vessel forces for selected values of plasma resistivity. The model also predicts poloidal and toroidal wall currents which agree with measured values; (2) The force-time history from the above model, used in conjunction with an axisymmetric structural model of the vessel, predicts vessel motions which agree well with measured values; (3) The above results, input to a simple acoustic model predicts the magnitude of sounds emitted from the vessel during disruptions which agree with acoustic measurements; (4) Correlation of measured vessel motions with structural analysis shows that a maximum vertical motion of the vessel up to 0.24 in will not overstress the vessel or its supports. 11 refs., 10 figs., 1 tab

Thermalization threshold in models of 1D fermions

Science.gov (United States)

Mukerjee, Subroto; Modak, Ranjan; Ramswamy, Sriram

2013-03-01

The question of how isolated quantum systems thermalize is an interesting and open one. In this study we equate thermalization with non-integrability to try to answer this question. In particular, we study the effect of system size on the integrability of 1D systems of interacting fermions on a lattice. We find that for a finite-sized system, a non-zero value of an integrability breaking parameter is required to make an integrable system appear non-integrable. Using exact diagonalization and diagnostics such as energy level statistics and the Drude weight, we find that the threshold value of the integrability breaking parameter scales to zero as a power law with system size. We find the exponent to be the same for different models with its value depending on the random matrix ensemble describing the non-integrable system. We also study a simple analytical model of a non-integrable system with an integrable limit to better understand how a power law emerges.

A simple 1D model with thermomechanical coupling for superelastic SMAs

International Nuclear Information System (INIS)

Zaki, W; Morin, C; Moumni, Z

2010-01-01

This paper presents an outline for a new uniaxial model for shape memory alloys that accounts for thermomechanical coupling. The coupling provides an explanation of the dependence of SMA behavior on the loading rate. 1D simulations are carried in Matlab using simple finite-difference discretization of the mechanical and thermal equations.

Amyotrophic Lateral Sclerosis Type 20 - In Silico Analysis and Molecular Dynamics Simulation of hnRNPA1.

Directory of Open Access Journals (Sweden)

Bruna Baumgarten Krebs

Full Text Available Amyotrophic Lateral Sclerosis (ALS is a fatal neurodegenerative disease that affects the upper and lower motor neurons. 5-10% of cases are genetically inherited, including ALS type 20, which is caused by mutations in the hnRNPA1 gene. The goals of this work are to analyze the effects of non-synonymous single nucleotide polymorphisms (nsSNPs on hnRNPA1 protein function, to model the complete tridimensional structure of the protein using computational methods and to assess structural and functional differences between the wild type and its variants through Molecular Dynamics simulations. nsSNP, PhD-SNP, Polyphen2, SIFT, SNAP, SNPs&GO, SNPeffect and PROVEAN were used to predict the functional effects of nsSNPs. Ab initio modeling of hnRNPA1 was made using Rosetta and refined using KoBaMIN. The structure was validated by PROCHECK, Rampage, ERRAT, Verify3D, ProSA and Qmean. TM-align was used for the structural alignment. FoldIndex, DICHOT, ELM, D2P2, Disopred and DisEMBL were used to predict disordered regions within the protein. Amino acid conservation analysis was assessed by Consurf, and the molecular dynamics simulations were performed using GROMACS. Mutations D314V and D314N were predicted to increase amyloid propensity, and predicted as deleterious by at least three algorithms, while mutation N73S was predicted as neutral by all the algorithms. D314N and D314V occur in a highly conserved amino acid. The Molecular Dynamics results indicate that all mutations increase protein stability when compared to the wild type. Mutants D314N and N319S showed higher overall dimensions and accessible surface when compared to the wild type. The flexibility level of the C-terminal residues of hnRNPA1 is affected by all mutations, which may affect protein function, especially regarding the protein ability to interact with other proteins.

GASFLOW: A Computational Fluid Dynamics Code for Gases, Aerosols, and Combustion, Volume 1: Theory and Computational Model

International Nuclear Information System (INIS)

Nichols, B.D.; Mueller, C.; Necker, G.A.; Travis, J.R.; Spore, J.W.; Lam, K.L.; Royl, P.; Redlinger, R.; Wilson, T.L.

1998-01-01

Los Alamos National Laboratory (LANL) and Forschungszentrum Karlsruhe (FzK) are developing GASFLOW, a three-dimensional (3D) fluid dynamics field code as a best-estimate tool to characterize local phenomena within a flow field. Examples of 3D phenomena include circulation patterns; flow stratification; hydrogen distribution mixing and stratification; combustion and flame propagation; effects of noncondensable gas distribution on local condensation and evaporation; and aerosol entrainment, transport, and deposition. An analysis with GASFLOW will result in a prediction of the gas composition and discrete particle distribution in space and time throughout the facility and the resulting pressure and temperature loadings on the walls and internal structures with or without combustion. A major application of GASFLOW is for predicting the transport, mixing, and combustion of hydrogen and other gases in nuclear reactor containments and other facilities. It has been applied to situations involving transporting and distributing combustible gas mixtures. It has been used to study gas dynamic behavior (1) in low-speed, buoyancy-driven flows, as well as sonic flows or diffusion dominated flows; and (2) during chemically reacting flows, including deflagrations. The effects of controlling such mixtures by safety systems can be analyzed. The code version described in this manual is designated GASFLOW 2.1, which combines previous versions of the United States Nuclear Regulatory Commission code HMS (for Hydrogen Mixing Studies) and the Department of Energy and FzK versions of GASFLOW. The code was written in standard Fortran 90. This manual comprises three volumes. Volume I describes the governing physical equations and computational model. Volume II describes how to use the code to set up a model geometry, specify gas species and material properties, define initial and boundary conditions, and specify different outputs, especially graphical displays. Sample problems are included

A mathematical model for fluid shear-sensitive 3D tissue construct development.

Science.gov (United States)

Liu, Dan; Chua, Chee-Kai; Leong, Kah-Fai

2013-01-01

This research studies dynamic culture for 3D tissue construct development with computational fluid dynamics. It proposes a mathematical model to evaluate the impact of flow rates and flow shear stress on cell growth in 3D constructs under perfusion. The modeling results show that dynamic flow, even at flow rate as low as 0.002 cm/s, can support much better mass exchange, higher cell number, and more even cell and nutrient distribution compared to static culture. Higher flow rate can further improve nutrient supply and mass exchange in the construct, promoting better nutritious environment and cell proliferation compared to lower flow rate. In addition, consideration of flow shear stress predicts much higher cell number in the construct compared to that without shear consideration. While the nutrient can dominate shear stress in influencing cell proliferation, the shear effect increases with flow rate. The proposed model helps tissue engineers better understand the cell-flow relationship at the molecular level during dynamic culture.

Testing a 1-D Analytical Salt Intrusion Model and the Predictive Equation in Malaysian Estuaries