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

NASA Astrophysics Data System (ADS)

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 troposphere mass flux show good results in particular in the Northern Hemisphere. In the Southern Hemisphere, the model tends to produce too-weak zonal-mean zonal winds and a too-narrow Hadley circulation. We discuss possible reasons for these model biases as well as planned future model improvements and applications.

Some anticipated contributions to core fluid dynamics from the GRM

NASA Technical Reports Server (NTRS)

Vanvorhies, C.

1985-01-01

It is broadly maintained that the secular variation (SV) of the large scale geomagnetic field contains information on the fluid dynamics of Earth's electrically conducting outer core. The electromagnetic theory appropriate to a simple Earth model has recently been combined with reduced geomagnetic data in order to extract some of this information and ascertain its significance. The simple Earth model consists of a rigid, electrically insulating mantle surrounding a spherical, inviscid, and perfectly conducting liquid outer core. This model was tested against seismology by using truncated spherical harmonic models of the observed geomagnetic field to locate Earth's core-mantle boundary, CMB. Further electromagnetic theory has been developed and applied to the problem of estimating the horizontal fluid motion just beneath CMB. Of particular geophysical interest are the hypotheses that these motions: (1) include appreciable surface divergence indicative of vertical motion at depth, and (2) are steady for time intervals of a decade or more. In addition to the extended testing of the basic Earth model, the proposed GRM provides a unique opportunity to test these dynamical hypotheses.

Dynamic Responses of the Earth's Outer Core to Assimilation of Observed Geomagnetic Secular Variation

NASA Technical Reports Server (NTRS)

Kuang, Weijia; Tangborn, Andrew

2014-01-01

Assimilation of surface geomagnetic observations and geodynamo models has advanced very quickly in recent years. However, compared to advanced data assimilation systems in meteorology, geomagnetic data assimilation (GDAS) is still in an early stage. Among many challenges ranging from data to models is the disparity between the short observation records and the long time scales of the core dynamics. To better utilize available observational information, we have made an effort in this study to directly assimilate the Gauss coefficients of both the core field and its secular variation (SV) obtained via global geomagnetic field modeling, aiming at understanding the dynamical responses of the core fluid to these additional observational constraints. Our studies show that the SV assimilation helps significantly to shorten the dynamo model spin-up process. The flow beneath the core-mantle boundary (CMB) responds significantly to the observed field and its SV. The strongest responses occur in the relatively small scale flow (of the degrees L is approx. 30 in spherical harmonic expansions). This part of the flow includes the axisymmetric toroidal flow (of order m = 0) and non-axisymmetric poloidal flow with m (is) greater than 5. These responses can be used to better understand the core flow and, in particular, to improve accuracies of predicting geomagnetic variability in future.

Role of core excitation in (d ,p ) transfer reactions

NASA Astrophysics Data System (ADS)

Deltuva, A.; Ross, A.; Norvaišas, E.; Nunes, F. M.

2016-10-01

Background: Recent work found that core excitations can be important in extracting structure information from (d ,p ) reactions. Purpose: Our objective is to systematically explore the role of core excitation in (d ,p ) reactions and to understand the origin of the dynamical effects. Method: Based on the particle-rotor model of n +10Be , we generate a number of models with a range of separation energies (Sn=0.1 -5.0 MeV), while maintaining a significant core excited component. We then apply the latest extension of the momentum-space-based Faddeev method, including dynamical core excitation in the reaction mechanism to all orders, to the 10Be(d ,p )11Be -like reactions, and study the excitation effects for beam energies Ed=15 -90 MeV. Results: We study the resulting angular distributions and the differences between the spectroscopic factor that would be extracted from the cross sections, when including dynamical core excitation in the reaction, and that of the original structure model. We also explore how different partial waves affect the final cross section. Conclusions: Our results show a strong beam-energy dependence of the extracted spectroscopic factors that become smaller for intermediate beam energies. This dependence increases for loosely bound systems.

A decision tree algorithm for investigation of model biases related to dynamical cores and physical parameterizations: CESM/CAM EVALUATION BY DECISION TREES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Soner Yorgun, M.; Rood, Richard B.

An object-based evaluation method using a pattern recognition algorithm (i.e., classification trees) is applied to the simulated orographic precipitation for idealized experimental setups using the National Center of Atmospheric Research (NCAR) Community Atmosphere Model (CAM) with the finite volume (FV) and the Eulerian spectral transform dynamical cores with varying resolutions. Daily simulations were analyzed and three different types of precipitation features were identified by the classification tree algorithm. The statistical characteristics of these features (i.e., maximum value, mean value, and variance) were calculated to quantify the difference between the dynamical cores and changing resolutions. Even with the simple and smoothmore » topography in the idealized setups, complexity in the precipitation fields simulated by the models develops quickly. The classification tree algorithm using objective thresholding successfully detected different types of precipitation features even as the complexity of the precipitation field increased. The results show that the complexity and the bias introduced in small-scale phenomena due to the spectral transform method of CAM Eulerian spectral dynamical core is prominent, and is an important reason for its dissimilarity from the FV dynamical core. The resolvable scales, both in horizontal and vertical dimensions, have significant effect on the simulation of precipitation. The results of this study also suggest that an efficient and informative study about the biases produced by GCMs should involve daily (or even hourly) output (rather than monthly mean) analysis over local scales.« less

A decision tree algorithm for investigation of model biases related to dynamical cores and physical parameterizations: CESM/CAM EVALUATION BY DECISION TREES

DOE PAGES

Soner Yorgun, M.; Rood, Richard B.

2016-11-11

An object-based evaluation method using a pattern recognition algorithm (i.e., classification trees) is applied to the simulated orographic precipitation for idealized experimental setups using the National Center of Atmospheric Research (NCAR) Community Atmosphere Model (CAM) with the finite volume (FV) and the Eulerian spectral transform dynamical cores with varying resolutions. Daily simulations were analyzed and three different types of precipitation features were identified by the classification tree algorithm. The statistical characteristics of these features (i.e., maximum value, mean value, and variance) were calculated to quantify the difference between the dynamical cores and changing resolutions. Even with the simple and smoothmore » topography in the idealized setups, complexity in the precipitation fields simulated by the models develops quickly. The classification tree algorithm using objective thresholding successfully detected different types of precipitation features even as the complexity of the precipitation field increased. The results show that the complexity and the bias introduced in small-scale phenomena due to the spectral transform method of CAM Eulerian spectral dynamical core is prominent, and is an important reason for its dissimilarity from the FV dynamical core. The resolvable scales, both in horizontal and vertical dimensions, have significant effect on the simulation of precipitation. The results of this study also suggest that an efficient and informative study about the biases produced by GCMs should involve daily (or even hourly) output (rather than monthly mean) analysis over local scales.« less

New quasi-geostrophic flow estimations for the Earth's core

NASA Astrophysics Data System (ADS)

Pais, M. Alexandra

2014-05-01

Quasi-geostrophic (QG) flows have been reported in numerical dynamo studies that simulate Boussinesq convection of an electrical conducting fluid inside a rapidly rotating spherical shell. In these cases, the required condition for columnar convection seems to be that inertial waves should propagate much faster in the medium than Alfvén waves. QG models are particularly appealing for studies where Earth's liquid core flows are assessed from information contained in geomagnetic data obtained at and above the Earth's surface. Here, they make the whole difference between perceiving only the core surface expression of the geodynamo or else assessing the whole interior core flow. The QG approximation has now been used in different studies to invert geomagnetic field models, providing a different kinematic interpretation of the observed geomagnetic field secular variation (SV). Under this new perspective, a large eccentric jet flowing westward under the Atlantic Hemisphere and a cyclonic column under the Pacific were pointed out as interesting features of the flow. A large eccentric jet with similar characteristics has been explained in recent numerical geodynamo simulations in terms of dynamical coupling between the solid core, the liquid core and the mantle. Nonetheless, it requires an inner core crystallization on the eastern hemisphere, contrary to what has been proposed in recent dynamical models for the inner core. Some doubts remain, as we see, concerning the dynamics that can explain the radial outward flow in the eastern core hemisphere, actually seen in inverted core flow models. This and other puzzling features justify a new assessment of core flows, taking full advantage of the recent geomagnetic field model COV-OBS and of experience, accumulated over the years, on flow inversion. Assuming the QG approximation already eliminates a large part of non-uniqueness in the inversion. Some important non-uniqueness still remains, inherent to the physical model, given our present inability to distinguish the small length scales of the internal geomagnetic field when measuring it at the Earth's surface and above. This can be dealt with in the form of a parameterization error. We recalculated flow models for the whole 1840-2010 period of COV-OBS, using the covariance matrices provided by the authors and an iterative estimation of the parameterization error. Results are compared with previous estimations. We then apply standard tools of Empirical Orthogonal Functions/ Principal Components Analysis to sort out variability modes that, hopefully, can also be identified with dynamical modes.

Role of core excitation in ( d , p ) transfer reactions

DOE PAGES

Deltuva, A.; Ross, A.; Norvaišas, E.; ...

2016-10-24

In our recent work we found that core excitations can be important in extracting structure information from (d,p) reactions. Our objective is to systematically explore the role of core excitation in (d,p) reactions and to understand the origin of the dynamical effects. Based on the particle-rotor model of n+Be 10, we generate a number of models with a range of separation energies (S n=0.1–5.0 MeV), while maintaining a significant core excited component. We then apply the latest extension of the momentum-space-based Faddeev method, including dynamical core excitation in the reaction mechanism to all orders, to the Be 10(d,p)Be 11-like reactions,more » and study the excitation effects for beam energies E d=15–90 MeV. We study the resulting angular distributions and the differences between the spectroscopic factor that would be extracted from the cross sections, when including dynamical core excitation in the reaction, and that of the original structure model. We also explore how different partial waves affect the final cross section. Our results show a strong beam-energy dependence of the extracted spectroscopic factors that become smaller for intermediate beam energies. Finally, this dependence increases for loosely bound systems.« less

Developments in the Gung Ho dynamical core

NASA Astrophysics Data System (ADS)

Melvin, Thomas

2017-04-01

Gung Ho is the new dynamical core being developed for the next generation Met Office weather and climate model, suitable for meeting the exascale challenge on emerging computer architectures. It builds upon the earlier collaborative project between the Met Office, NERC and STFC Daresbury of the same name to investigate suitable numerical methods for dynamical cores. A mixed-finite element approach is used, where different finite element spaces are used to represent various fields. This method provides a number of beneficial improvements over the current model, such a compatibility and inherent conservation on quasi-uniform unstructured meshes, whilst maintaining the accuracy and good dispersion properties of the staggered grid currently used. Furthermore, the mixed finite element approach allows a large degree of flexibility in the type of mesh, order of approximation and discretisation, providing a simple way to test alternative options to obtain the best model possible.

Dynamic analysis of gas-core reactor system

NASA Technical Reports Server (NTRS)

Turner, K. H., Jr.

1973-01-01

A heat transfer analysis was incorporated into a previously developed model CODYN to obtain a model of open-cycle gaseous core reactor dynamics which can predict the heat flux at the cavity wall. The resulting model was used to study the sensitivity of the model to the value of the reactivity coefficients and to determine the system response for twenty specified perturbations. In addition, the model was used to study the effectiveness of several control systems in controlling the reactor. It was concluded that control drums located in the moderator region capable of inserting reactivity quickly provided the best control.

Molecular Determinants and Dynamics of Hepatitis C Virus Secretion

PubMed Central

Coller, Kelly E.; Heaton, Nicholas S.; Berger, Kristi L.; Cooper, Jacob D.; Saunders, Jessica L.; Randall, Glenn

2012-01-01

The current model of hepatitis C virus (HCV) production involves the assembly of virions on or near the surface of lipid droplets, envelopment at the ER in association with components of VLDL synthesis, and egress via the secretory pathway. However, the cellular requirements for and a mechanistic understanding of HCV secretion are incomplete at best. We combined an RNA interference (RNAi) analysis of host factors for infectious HCV secretion with the development of live cell imaging of HCV core trafficking to gain a detailed understanding of HCV egress. RNAi studies identified multiple components of the secretory pathway, including ER to Golgi trafficking, lipid and protein kinases that regulate budding from the trans-Golgi network (TGN), VAMP1 vesicles and adaptor proteins, and the recycling endosome. Our results support a model wherein HCV is infectious upon envelopment at the ER and exits the cell via the secretory pathway. We next constructed infectious HCV with a tetracysteine (TC) tag insertion in core (TC-core) to monitor the dynamics of HCV core trafficking in association with its cellular cofactors. In order to isolate core protein movements associated with infectious HCV secretion, only trafficking events that required the essential HCV assembly factor NS2 were quantified. TC-core traffics to the cell periphery along microtubules and this movement can be inhibited by nocodazole. Sub-populations of TC-core localize to the Golgi and co-traffic with components of the recycling endosome. Silencing of the recycling endosome component Rab11a results in the accumulation of HCV core at the Golgi. The majority of dynamic core traffics in association with apolipoprotein E (ApoE) and VAMP1 vesicles. This study identifies many new host cofactors of HCV egress, while presenting dynamic studies of HCV core trafficking in infected cells. PMID:22241992

Acoustics and dynamics of coaxial interacting vortex rings

NASA Technical Reports Server (NTRS)

Shariff, Karim; Leonard, Anthony; Zabusky, Norman J.; Ferziger, Joel H.

1988-01-01

Using a contour dynamics method for inviscid axisymmetric flow we examine the effects of core deformation on the dynamics and acoustic signatures of coaxial interacting vortex rings. Both 'passage' and 'collision' (head-on) interactions are studied for initially identical vortices. Good correspondence with experiments is obtained. A simple model which retains only the elliptic degree of freedom in the core shape is used to explain some of the calculated features.

Density Anomalies in the Mantle and the Gravitational Core-Mantle Interaction

NASA Technical Reports Server (NTRS)

Kuang, Weijia; Liu, Lanbo

2003-01-01

Seismic studies suggest that the bulk of the mantle is heterogeneous, with density variations in depth as well as in horizontal directions (latitude and longitude). This density variation produces a three- dimensional gravity field throughout the Earth. On the other hand, the core density also varies in both time and space, due to convective core flow. Consequently, the fluid outer core and the solid mantle interact gravitationally due to the mass anomalies in both regions. This gravitational core-mantle interaction could play a significant role in exchange of angular momentum between the core and the mantle, and thus the change in Earth's rotation on time scales of decades and longer. Aiming at estimating the significance of the gravitational core-mantle interaction on Earth's rotation variation, we introduce in our MoSST core dynamics model a heterogeneous mantle, with a density distribution derived from seismic results. In this model, the core convection is driven by the buoyancy forces. And the density variation is determined dynamically with the convection. Numerical simulation is carried out with different parameter values, intending to extrapolate numerical results for geophysical implications.

Modeling analysis of pulsed magnetization process of magnetic core based on inverse Jiles-Atherton model

NASA Astrophysics Data System (ADS)

Liu, Yi; Zhang, He; Liu, Siwei; Lin, Fuchang

2018-05-01

The J-A (Jiles-Atherton) model is widely used to describe the magnetization characteristics of magnetic cores in a low-frequency alternating field. However, this model is deficient in the quantitative analysis of the eddy current loss and residual loss in a high-frequency magnetic field. Based on the decomposition of magnetization intensity, an inverse J-A model is established which uses magnetic flux density B as an input variable. Static and dynamic core losses under high frequency excitation are separated based on the inverse J-A model. Optimized parameters of the inverse J-A model are obtained based on particle swarm optimization. The platform for the pulsed magnetization characteristic test is designed and constructed. The hysteresis curves of ferrite and Fe-based nanocrystalline cores at high magnetization rates are measured. The simulated and measured hysteresis curves are presented and compared. It is found that the inverse J-A model can be used to describe the magnetization characteristics at high magnetization rates and to separate the static loss and dynamic loss accurately.

Preparing CAM-SE for Multi-Tracer Applications: CAM-SE-Cslam

NASA Astrophysics Data System (ADS)

Lauritzen, P. H.; Taylor, M.; Goldhaber, S.

2014-12-01

The NCAR-DOE spectral element (SE) dynamical core comes from the HOMME (High-Order Modeling Environment; Dennis et al., 2012) and it is available in CAM. The CAM-SE dynamical core is designed with intrinsic mimetic properties guaranteeing total energy conservation (to time-truncation errors) and mass-conservation, and has demonstrated excellent scalability on massively parallel compute platforms (Taylor, 2011). For applications involving many tracers such as chemistry and biochemistry modeling, CAM-SE has been found to be significantly more computationally costly than the current "workhorse" model CAM-FV (Finite-Volume; Lin 2004). Hence a multi-tracer efficient scheme, called the CSLAM (Conservative Semi-Lagrangian Multi-tracer; Lauritzen et al., 2011) scheme, has been implemented in the HOMME (Erath et al., 2012). The CSLAM scheme has recently been cast in flux-form in HOMME so that it can be coupled to the SE dynamical core through conventional flux-coupling methods where the SE dynamical core provides background air mass fluxes to CSLAM. Since the CSLAM scheme makes use of a finite-volume gnomonic cubed-sphere grid and hence does not operate on the SE quadrature grid, the capability of running tracer advection, the physical parameterization suite and dynamics on separate grids has been implemented in CAM-SE. The default CAM-SE-CSLAM setup is to run physics on the quasi-equal area CSLAM grid. The capability of running physics on a different grid than the SE dynamical core may provide a more consistent coupling since the physics grid option operates with quasi-equal-area cell average values rather than non-equi-distant grid-point (SE quadrature point) values. Preliminary results on the performance of CAM-SE-CSLAM will be presented.

Geodynamo Modeling of Core-Mantle Interactions

NASA Technical Reports Server (NTRS)

Kuang, Wei-Jia; Chao, Benjamin F.; Smith, David E. (Technical Monitor)

2001-01-01

Angular momentum exchange between the Earth's mantle and core influences the Earth's rotation on time scales of decades and longer, in particular in the length of day (LOD) which have been measured with progressively increasing accuracy for the last two centuries. There are four possible coupling mechanisms for transferring the axial angular momentum across the core-mantle boundary (CMB): viscous, magnetic, topography, and gravitational torques. Here we use our scalable, modularized, fully dynamic geodynamo model for the core to assess the importance of these torques. This numerical model, as an extension of the Kuang-Bloxham model that has successfully simulated the generation of the Earth's magnetic field, is used to obtain numerical results in various physical conditions in terms of specific parameterization consistent with the dynamical processes in the fluid outer core. The results show that depending on the electrical conductivity of the lower mantle and the amplitude of the boundary topography at CMB, both magnetic and topographic couplings can contribute significantly to the angular momentum exchange. This implies that the core-mantle interactions are far more complex than has been assumed and that there is unlikely a single dominant coupling mechanism for the observed decadal LOD variation.

Testing the dynamic coupling of the core-mantle and inner core boundaries

NASA Astrophysics Data System (ADS)

Driscoll, Peter E.

2015-07-01

The proposal that the seismically observed hemispherical asymmetry of Earth's inner core is controlled by the heat flux structure imposed on the outer core by the lower mantle is tested with numerical dynamo models driven by mixed thermochemical convection. We find that models driven by a single core-mantle boundary (CMB) spherical harmonic of degree and mode 2, the dominant mode in lower mantle seismic shear velocity tomography, produce a similar structure at the inner core boundary (ICB) shifted 30∘ westward. The sensitivity of the ICB to the CMB is further tested by increasing the CMB heterogeneity amplitude. In addition, two seismic tomographic models are tested: first with CMB resolution up to degree and order 4, and second with resolution up to degree and order 8. We find time-averaged ICB heat flux in these cases to be similar at large scale, with small-scale differences due to higher CMB harmonics (above degree 4). The tomographic models produce "Earth-like" magnetic fields, while similar models with twice the CMB heat flow amplitudes produce less Earth-like fields, implying that increasing CMB heterogeneity forces the model out of an Earth-like regime. The dynamic ICB heat fluxes are compared to the proposed translation mode of the inner core to test whether the CMB controls inner core growth and structure. This test indicates that, although CMB tomography is unlikely to be driving inner core translation, the ICB heat flux response is weak enough to not interfere with the most unstable translation mode, if it is occurring.

Multigrid based First-Principles Molecular Dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fattebert, Jean-Luc; Osei-Kuffuor, Daniel; Dunn, Ian

2017-06-01

MGmol ls a First-Principles Molecular Dynamics code. It relies on the Born-Oppenheimer approximation and models the electronic structure using Density Functional Theory, either LDA or PBE. Norm-conserving pseudopotentials are used to model atomic cores.

A new look at sulphur chemistry in hot cores and corinos

NASA Astrophysics Data System (ADS)

Vidal, Thomas H. G.; Wakelam, Valentine

2018-03-01

Sulphur-bearing species are often used to probe the evolution of hot cores because their abundances are particularly sensitive to physical and chemical variations. However, the chemistry of sulphur is not well understood in these regions, notably because observations of several hot cores have displayed a large variety of sulphur compositions, and because the reservoir of sulphur in dense clouds, in which hot cores form, is still poorly constrained. In order to help disentangle its complexity, we present a fresh comprehensive review of sulphur chemistry in hot cores along with a study of sulphur's sensibility to temperature and pre-collapse chemical composition. In parallel, we analyse the discrepencies that result from the use of two different types of models (static and dynamic) in order to highlight the sensitivity to the choice of model to be used in astrochemical studies. Our results show that the pre-collapse chemical composition is a critical parameter for sulphur chemistry in hot cores and that it could explain the different sulphur compositions observed. We also report that differences in abundances for a given species between the static and dynamic models can reach six orders of magnitude in the hot core, which reveals the key role of the choice of model in astrochemical studies.

DYNAMICO, an atmospheric dynamical core for high-performance climate modeling

NASA Astrophysics Data System (ADS)

Dubos, Thomas; Meurdesoif, Yann; Spiga, Aymeric; Millour, Ehouarn; Fita, Lluis; Hourdin, Frédéric; Kageyama, Masa; Traore, Abdoul-Khadre; Guerlet, Sandrine; Polcher, Jan

2017-04-01

Institut Pierre Simon Laplace has developed a very scalable atmospheric dynamical core, DYNAMICO, based on energy-conserving finite-difference/finite volume numerics on a quasi-uniform icosahedral-hexagonal mesh. Scalability is achieved by combining hybrid MPI/OpenMP parallelism to asynchronous I/O. This dynamical core has been coupled to radiative transfer physics tailored to the atmosphere of Saturn, allowing unprecedented simulations of the climate of this giant planet. For terrestrial climate studies DYNAMICO is being integrated into the IPSL Earth System Model IPSL-CM. Preliminary aquaplanet and AMIP-style simulations yield reasonable results when compared to outputs from IPSL-CM5. The observed performance suggests that an order of magnitude may be gained with respect to IPSL-CM CMIP5 simulations either on the duration of simulations or on their resolution. Longer simulations would be of interest for the study of paleoclimate, while higher resolution could improve certain aspects of the modeled climate such as extreme events, as will be explored in the HighResMIP project. Following IPSL's strategic vision of building a unified global-regional modelling system, a fully-compressible, non-hydrostatic prototype of DYNAMICO has been developed, enabling future convection-resolving simulations. Work supported by ANR project "HEAT", grant number CE23_2014_HEAT Dubos, T., Dubey, S., Tort, M., Mittal, R., Meurdesoif, Y., and Hourdin, F.: DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility, Geosci. Model Dev., 8, 3131-3150, doi:10.5194/gmd-8-3131-2015, 2015.

Prospect of Using Numerical Dynamo Model for Prediction of Geomagnetic Secular Variation

NASA Technical Reports Server (NTRS)

Kuang, Weijia; Tangborn, Andrew

2003-01-01

Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Sloan Great Wall as a complex of superclusters with collapsing cores

NASA Astrophysics Data System (ADS)

Einasto, Maret; Lietzen, Heidi; Gramann, Mirt; Tempel, Elmo; Saar, Enn; Liivamägi, Lauri Juhan; Heinämäki, Pekka; Nurmi, Pasi; Einasto, Jaan

2016-10-01

Context. The formation and evolution of the cosmic web is governed by the gravitational attraction of dark matter and antigravity of dark energy (cosmological constant). In the cosmic web, galaxy superclusters or their high-density cores are the largest objects that may collapse at present or during the future evolution. Aims: We study the dynamical state and possible future evolution of galaxy superclusters from the Sloan Great Wall (SGW), the richest galaxy system in the nearby Universe. Methods: We calculated supercluster masses using dynamical masses of galaxy groups and stellar masses of galaxies. We employed normal mixture modelling to study the structure of rich SGW superclusters and search for components (cores) in superclusters. We analysed the radial mass distribution in the high-density cores of superclusters centred approximately at rich clusters and used the spherical collapse model to study their dynamical state. Results: The lower limit of the total mass of the SGW is approximately M = 2.5 × 1016 h-1 M⊙. Different mass estimators of superclusters agree well, the main uncertainties in masses of superclusters come from missing groups and clusters. We detected three high-density cores in the richest SGW supercluster (SCl 027) and two in the second richest supercluster (SCl 019). They have masses of 1.2 - 5.9 × 1015 h-1 M⊙ and sizes of up to ≈60 h-1 Mpc. The high-density cores of superclusters are very elongated, flattened perpendicularly to the line of sight. The comparison of the radial mass distribution in the high-density cores with the predictions of spherical collapse model suggests that their central regions with radii smaller than 8 h-1 Mpc and masses of up to M = 2 × 1015 h-1 M⊙ may be collapsing. Conclusions: The rich SGW superclusters with their high-density cores represent dynamically evolving environments for studies of the properties of galaxies and galaxy systems.

Geomagnetic acceleration and rapid hydromagnetic wave dynamics in advanced numerical simulations of the geodynamo

NASA Astrophysics Data System (ADS)

Aubert, Julien

2018-04-01

Geomagnetic secular acceleration, the second temporal derivative of Earth's magnetic field, is a unique window on the dynamics taking place in Earth's core. In this study, the behaviours of the secular acceleration and underlying core dynamics are examined in new numerical simulations of the geodynamo that are dynamically closer to Earth's core conditions than earlier models. These new models reside on a theoretical path in parameter space connecting the region where most classical models are found to the natural conditions. The typical time scale for geomagnetic acceleration is found to be invariant along this path, at a value close to 10 years that matches Earth's core estimates. Despite this invariance, the spatio-temporal properties of secular acceleration show significant variability along the path, with an asymptotic regime of rapid rotation reached after 30% of this path (corresponding to a model Ekman number E = 3 - 7). In this regime, the energy of secular acceleration is entirely found at periods longer than that of planetary rotation, and the underlying flow acceleration patterns acquire a two-dimensional columnar structure representative of the rapid rotation limit. The spatial pattern of the secular acceleration at the core-mantle boundary shows significant localisation of energy within an equatorial belt. Rapid hydromagnetic wave dynamics is absent at the start of the path because of insufficient time scale separation with convective processes, weak forcing and excessive damping but can be clearly exhibited in the asymptotic regime. This study reports on ubiquitous axisymmetric geostrophic torsional waves of weak amplitude relatively to convective transport, and also stronger, laterally limited, quasi-geostrophic Alfvén waves propagating in the cylindrical radial direction from the tip of convective plumes towards the core-mantle boundary. In a system similar to Earth's core where the typical Alfvén velocity is significantly larger than the typical convective velocity, quasi-geostrophic Alfvén waves are shown to be an important carrier of flow acceleration to the core surface that links with the generation of strong, short-lived and intermittent equatorial pulses in the secular acceleration energy. The secular acceleration time scale is shown to be insensitive to magnetic signatures from torsional waves because of their weak amplitude, and from quasi-geostrophic Alfvén waves because of their intermittent character, and is therefore only indicative of convective transport phenomena that remain invariant along the parameter space path.

LFRic: Building a new Unified Model

NASA Astrophysics Data System (ADS)

Melvin, Thomas; Mullerworth, Steve; Ford, Rupert; Maynard, Chris; Hobson, Mike

2017-04-01

The LFRic project, named for Lewis Fry Richardson, aims to develop a replacement for the Met Office Unified Model in order to meet the challenges which will be presented by the next generation of exascale supercomputers. This project, a collaboration between the Met Office, STFC Daresbury and the University of Manchester, builds on the earlier GungHo project to redesign the dynamical core, in partnership with NERC. The new atmospheric model aims to retain the performance of the current ENDGame dynamical core and associated subgrid physics, while also enabling a far greater scalability and flexibility to accommodate future supercomputer architectures. Design of the model revolves around a principle of a 'separation of concerns', whereby the natural science aspects of the code can be developed without worrying about the underlying architecture, while machine dependent optimisations can be carried out at a high level. These principles are put into practice through the development of an autogenerated Parallel Systems software layer (known as the PSy layer) using a domain-specific compiler called PSyclone. The prototype model includes a re-write of the dynamical core using a mixed finite element method, in which different function spaces are used to represent the various fields. It is able to run in parallel with MPI and OpenMP and has been tested on over 200,000 cores. In this talk an overview of the both the natural science and computational science implementations of the model will be presented.

A comparison of East Asian summer monsoon simulations from CAM3.1 with three dynamic cores

NASA Astrophysics Data System (ADS)

Wei, Ting; Wang, Lanning; Dong, Wenjie; Dong, Min; Zhang, Jingyong

2011-12-01

This paper examines the sensitivity of CAM3.1 simulations of East Asian summer monsoon (EASM) to the choice of dynamic cores using three long-term simulations, one with each of the following cores: the Eulerian spectral transform method (EUL), semi-Lagrangian scheme (SLD) and finite volume approach (FV). Our results indicate that the dynamic cores significantly influence the simulated fields not only through dynamics, such as wind, but also through physical processes, such as precipitation. Generally speaking, SLD is superior to EUL and FV in simulating the climatological features of EASM and its interannual variability. The SLD version of the CAM model partially reduces its known deficiency in simulating the climatological features of East Asian summer precipitation. The strength and position of simulated western Pacific subtropical high (WPSH) and its ridge line compare more favourably with observations in SLD and FV than in EUL. They contribute to the intensification of the south-easterly along the south of WPSH and the vertical motion through the troposphere around 30° N, where the subtropical rain belt exists. Additionally, SLD simulates the scope of the westerly jet core over East Asia more realistically than the other two dynamic cores do. Considerable systematic errors of the seasonal migration of monsoon rain belt and water vapour flux exist in all of the three versions of CAM3.1 model, although it captures the broad northward shift of convection, and the simulated results share similarities. The interannual variation of EASM is found to be more accurate in SLD simulation, which reasonably reproduces the leading combined patterns of precipitation and 850-hPa winds in East Asia, as well as the 2.5- and 10-year periods of Li-Zeng EASM index. These results emphasise the importance of dynamic cores for the EASM simulation as distinct from the simulation's sensitivity to the physical parameterisations.

Models of the Earth's Core.

PubMed

Stevenson, D J

1981-11-06

Combined inferences from seismology, high-pressure experiment and theory, geomagnetism, fluid dynamics, and current views of terrestrial planetary evolution lead to models of the earth's core with the following properties. Core formation was contemporaneous with earth accretion; the core is not in chemical equilibrium with the mantle; the outer core is a fluid iron alloy containing significant quantities of lighter elements and is probably almost adiabatic and compositionally uniform; the more iron-rich inner solid core is a consequence of partial freezing of the outer core, and the energy release from this process sustains the earth's magnetic field; and the thermodynamic properties of the core are well constrained by the application of liquid-state theory to seismic and laboratory data.

Inverse stochastic-dynamic models for high-resolution Greenland ice core records

NASA Astrophysics Data System (ADS)

Boers, Niklas; Chekroun, Mickael D.; Liu, Honghu; Kondrashov, Dmitri; Rousseau, Denis-Didier; Svensson, Anders; Bigler, Matthias; Ghil, Michael

2017-12-01

Proxy records from Greenland ice cores have been studied for several decades, yet many open questions remain regarding the climate variability encoded therein. Here, we use a Bayesian framework for inferring inverse, stochastic-dynamic models from δ18O and dust records of unprecedented, subdecadal temporal resolution. The records stem from the North Greenland Ice Core Project (NGRIP), and we focus on the time interval 59-22 ka b2k. Our model reproduces the dynamical characteristics of both the δ18O and dust proxy records, including the millennial-scale Dansgaard-Oeschger variability, as well as statistical properties such as probability density functions, waiting times and power spectra, with no need for any external forcing. The crucial ingredients for capturing these properties are (i) high-resolution training data, (ii) cubic drift terms, (iii) nonlinear coupling terms between the δ18O and dust time series, and (iv) non-Markovian contributions that represent short-term memory effects.

Modelling dynamics in protein crystal structures by ensemble refinement

PubMed Central

Burnley, B Tom; Afonine, Pavel V; Adams, Paul D; Gros, Piet

2012-01-01

Single-structure models derived from X-ray data do not adequately account for the inherent, functionally important dynamics of protein molecules. We generated ensembles of structures by time-averaged refinement, where local molecular vibrations were sampled by molecular-dynamics (MD) simulation whilst global disorder was partitioned into an underlying overall translation–libration–screw (TLS) model. Modeling of 20 protein datasets at 1.1–3.1 Å resolution reduced cross-validated Rfree values by 0.3–4.9%, indicating that ensemble models fit the X-ray data better than single structures. The ensembles revealed that, while most proteins display a well-ordered core, some proteins exhibit a ‘molten core’ likely supporting functionally important dynamics in ligand binding, enzyme activity and protomer assembly. Order–disorder changes in HIV protease indicate a mechanism of entropy compensation for ordering the catalytic residues upon ligand binding by disordering specific core residues. Thus, ensemble refinement extracts dynamical details from the X-ray data that allow a more comprehensive understanding of structure–dynamics–function relationships. DOI: http://dx.doi.org/10.7554/eLife.00311.001 PMID:23251785

Theory of Auger core-valence-valence processes in simple metals. II. Dynamical and surface effects on Auger line shapes

NASA Astrophysics Data System (ADS)

Almbladh, C.-O.; Morales, A. L.

1989-02-01

Auger CVV spectra of simple metals are generally believed to be well described by one-electron-like theories in the bulk which account for matrix elements and, in some cases, also static core-hole screening effects. We present here detailed calculations on Li, Be, Na, Mg, and Al using self-consistent bulk wave functions and proper matrix elements. The resulting spectra differ markedly from experiment and peak at too low energies. To explain this discrepancy we investigate effects of the surface and dynamical effects of the sudden disappearance of the core hole in the final state. To study core-hole effects we solve Mahan-Nozières-De Dominicis (MND) model numerically over the entire band. The core-hole potential and other parameters in the MND model are determined by self-consistent calculations of the core-hole impurity. The results are compared with simpler approximations based on the final-state rule due to von Barth and Grossmann. To study surface and mean-free-path effects we perform slab calculations for Al but use a simpler infinite-barrier model in the remaining cases. The model reproduces the slab spectra for Al with very good accuracy. In all cases investigated either the effects of the surface or the effects of the core hole give important modifications and a much improved agreement with experiment.

Solvent friction effects propagate over the entire protein molecule through low-frequency collective modes.

PubMed

Moritsugu, Kei; Kidera, Akinori; Smith, Jeremy C

2014-07-24

Protein solvation dynamics has been investigated using atom-dependent Langevin friction coefficients derived directly from molecular dynamics (MD) simulations. To determine the effect of solvation on the atomic friction coefficients, solution and vacuum MD simulations were performed for lysozyme and staphylococcal nuclease and analyzed by Langevin mode analysis. The coefficients thus derived are roughly correlated with the atomic solvent-accessible surface area (ASA), as expected from the fact that friction occurs as the result of collisions with solvent molecules. However, a considerable number of atoms with higher friction coefficients are found inside the core region. Hence, the influence of solvent friction propagates into the protein core. The internal coefficients have large contributions from the low-frequency modes, yielding a simple picture of the surface-to-core long-range damping via solvation governed by collective low-frequency modes. To make use of these findings in implicit-solvent modeling, we compare the all-atom friction results with those obtained using Langevin dynamics (LD) with two empirical representations: the constant-friction and the ASA-dependent (Pastor-Karplus) friction models. The constant-friction model overestimates the core and underestimates the surface damping whereas the ASA-dependent friction model, which damps protein atoms only on the solvent-accessible surface, reproduces well the friction coefficients for both the surface and core regions observed in the explicit-solvent MD simulations. Therefore, in LD simulation, the solvent friction coefficients should be imposed only on the protein surface.

Solvent friction effects propagate over the entire protein molecule through low-frequency collective modes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Moritsugu, Kei; Kidera, Akinori; Smith, Jeremy C.

2014-06-25

Protein solvation dynamics has been investigated using atom-dependent Langevin friction coefficients derived directly from molecular dynamics (MD) simulations. To determine the effect of solvation on the atomic friction coefficients, solution and vacuum MD simulations were performed for lysozyme and staphylococcal nuclease and analyzed by Langevin mode analysis. The coefficients thus derived are roughly correlated with the atomic solvent-accessible surface area (ASA), as expected from the fact that friction occurs as the result of collisions with solvent molecules. However, a considerable number of atoms with higher friction coefficients are found inside the core region. Hence, the influence of solvent friction propagatesmore » into the protein core. The internal coefficients have large contributions from the low-frequency modes, yielding a simple picture of the surface-to-core long-range damping via solvation governed by collective low-frequency modes. To make use of these findings in implicit-solvent modeling, we compare the all-atom friction results with those obtained using Langevin dynamics (LD) with two empirical representations: the constant-friction and the ASA-dependent (Pastor Karplus) friction models. The constant-friction model overestimates the core and underestimates the surface damping whereas the ASA-dependent friction model, which damps protein atoms only on the solvent-accessible surface, reproduces well the friction coefficients for both the surface and core regions observed in the explicit-solvent MD simulations. Furthermore, in LD simulation, the solvent friction coefficients should be imposed only on the protein surface.« less

Structure, Dynamics, and Deuterium Fractionation of Massive Pre-stellar Cores

NASA Astrophysics Data System (ADS)

Goodson, Matthew D.; Kong, Shuo; Tan, Jonathan C.; Heitsch, Fabian; Caselli, Paola

2016-12-01

High levels of deuterium fraction in N2H+ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ({D}{frac}{{{N}}2{{{H}}}+}\\equiv {{{N}}}2{{{D}}}+/{{{N}}}2{{{H}}}+≳ 0.1) is longer than the free-fall time, possibly 10 times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N2H+ and N2D+. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N2D+ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H2, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios (≲ 0.01) and/or multiple free-fall times (≳ 3) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

Topology and Dynamics of the Zebrafish Segmentation Clock Core Circuit

PubMed Central

Schröter, Christian; Isakova, Alina; Hens, Korneel; Soroldoni, Daniele; Gajewski, Martin; Jülicher, Frank; Maerkl, Sebastian J.; Deplancke, Bart; Oates, Andrew C.

2012-01-01

During vertebrate embryogenesis, the rhythmic and sequential segmentation of the body axis is regulated by an oscillating genetic network termed the segmentation clock. We describe a new dynamic model for the core pace-making circuit of the zebrafish segmentation clock based on a systematic biochemical investigation of the network's topology and precise measurements of somitogenesis dynamics in novel genetic mutants. We show that the core pace-making circuit consists of two distinct negative feedback loops, one with Her1 homodimers and the other with Her7:Hes6 heterodimers, operating in parallel. To explain the observed single and double mutant phenotypes of her1, her7, and hes6 mutant embryos in our dynamic model, we postulate that the availability and effective stability of the dimers with DNA binding activity is controlled in a “dimer cloud” that contains all possible dimeric combinations between the three factors. This feature of our model predicts that Hes6 protein levels should oscillate despite constant hes6 mRNA production, which we confirm experimentally using novel Hes6 antibodies. The control of the circuit's dynamics by a population of dimers with and without DNA binding activity is a new principle for the segmentation clock and may be relevant to other biological clocks and transcriptional regulatory networks. PMID:22911291

Optimizing zonal advection of the Advanced Research WRF (ARW) dynamics for Intel MIC

NASA Astrophysics Data System (ADS)

Mielikainen, Jarno; Huang, Bormin; Huang, Allen H.

2014-10-01

The Weather Research and Forecast (WRF) model is the most widely used community weather forecast and research model in the world. There are two distinct varieties of WRF. The Advanced Research WRF (ARW) is an experimental, advanced research version featuring very high resolution. The WRF Nonhydrostatic Mesoscale Model (WRF-NMM) has been designed for forecasting operations. WRF consists of dynamics code and several physics modules. The WRF-ARW core is based on an Eulerian solver for the fully compressible nonhydrostatic equations. In the paper, we will use Intel Intel Many Integrated Core (MIC) architecture to substantially increase the performance of a zonal advection subroutine for optimization. It is of the most time consuming routines in the ARW dynamics core. Advection advances the explicit perturbation horizontal momentum equations by adding in the large-timestep tendency along with the small timestep pressure gradient tendency. We will describe the challenges we met during the development of a high-speed dynamics code subroutine for MIC architecture. Furthermore, lessons learned from the code optimization process will be discussed. The results show that the optimizations improved performance of the original code on Xeon Phi 5110P by a factor of 2.4x.

Models of the earth's core

NASA Technical Reports Server (NTRS)

Stevenson, D. J.

1981-01-01

Combined inferences from seismology, high-pressure experiment and theory, geomagnetism, fluid dynamics, and current views of terrestrial planetary evolution lead to models of the earth's core with five basic properties. These are that core formation was contemporaneous with earth accretion; the core is not in chemical equilibrium with the mantle; the outer core is a fluid iron alloy containing significant quantities of lighter elements and is probably almost adiabatic and compositionally uniform; the more iron-rich inner solid core is a consequence of partial freezing of the outer core, and the energy release from this process sustains the earth's magnetic field; and the thermodynamic properties of the core are well constrained by the application of liquid-state theory to seismic and labroatory data.

A Dynamic Model of Sustainment Investment

DTIC Science & Technology

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

Performance measurements of a pilot superconducting solenoid model core for a wind tunnel magnetic suspension and balance system

NASA Technical Reports Server (NTRS)

Goodyer, M. J.; Britcher, C. P.

1983-01-01

The results of experimental demonstrations of a superconducting solenoid model core in the Southampton University Magnetic Suspension and Balance System are detailed. Technology and techniques relevant to large-scale wind tunnel MSBSs comprise the long term goals. The magnetic moment of solenoids, difficulties peculiar to superconducting solenoid cores, lift force and pitching moment, dynamic lift calibration, and helium boil-off measurements are discussed.

Resolving Supercritical Orion Cores

NASA Astrophysics Data System (ADS)

Li, Di; Chapman, N.; Goldsmith, P.; Velusamy, T.

2009-01-01

The theoretical framework for high mass star formation (HMSF) is unclear. Observations reveal a seeming dichotomy between high- and low-mass star formation, with HMSF occurring only in Giant Molecular Clouds (GMC), mostly in clusters, and with higher star formation efficiencies than low-mass star formation. One crucial constraint to any theoretical model is the dynamical state of massive cores, in particular, whether a massive core is in supercritical collapse. Based on the mass-size relation of dust emission, we select likely unstable targets from a sample of massive cores (Li et al. 2007 ApJ 655, 351) in the nearest GMC, Orion. We have obtained N2H+ (1-0) maps using CARMA with resolution ( 2.5", 0.006 pc) significantly better than existing observations. We present observational and modeling results for ORI22. By revealing the dynamic structure down to Jeans scale, CARMA data confirms the dominance of gravity over turbulence in this cores. This work was performed by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

Outcomes of Grazing Impacts between Sub-Neptunes in Kepler  Multis

NASA Astrophysics Data System (ADS)

Hwang, Jason; Chatterjee, Sourav; Lombardi, James, Jr.; Steffen, Jason H.; Rasio, Frederic

2018-01-01

Studies of high-multiplicity, tightly packed planetary systems suggest that dynamical instabilities are common and affect both the orbits and planet structures, where the compact orbits and typically low densities make physical collisions likely outcomes. Since the structure of many of these planets is such that the mass is dominated by a rocky core, but the volume is dominated by a tenuous gas envelope, the sticky-sphere approximation, used in dynamical integrators, may be a poor model for these collisions. We perform five sets of collision calculations, including detailed hydrodynamics, sampling mass ratios, and core mass fractions typical in Kepler Multis. In our primary set of calculations, we use Kepler-36 as a nominal remnant system, as the two planets have a small dynamical separation and an extreme density ratio. We use an N-body code, Mercury 6.2, to integrate initially unstable systems and study the resultant collisions in detail. We use these collisions, focusing on grazing collisions, in combination with realistic planet models created using gas profiles from Modules for Experiments in Stellar Astrophysics and core profiles using equations of state from Seager et al. to perform hydrodynamic calculations, finding scatterings, mergers, and even a potential planet–planet binary. We dynamically integrate the remnant systems, examine the stability, and estimate the final densities, finding that the remnant densities are sensitive to the core masses, and collisions result in generally more stable systems. We provide prescriptions for predicting the outcomes and modeling the changes in mass and orbits following collisions for general use in dynamical integrators.

Geomagnetic acceleration and rapid hydromagnetic wave dynamics in advanced numerical simulations of the geodynamo

NASA Astrophysics Data System (ADS)

Aubert, Julien

2018-07-01

Geomagnetic secular acceleration, the second temporal derivative of the Earth's magnetic field, is a unique window on the dynamics taking place in the Earth's core. In this study, the behaviours of the secular acceleration and underlying core dynamics are examined in new numerical simulations of the geodynamo that are dynamically closer to the Earth's core conditions than earlier models. These new models reside on a theoretical path in parameter space connecting the region where most classical models are found to the natural conditions. The typical timescale for geomagnetic acceleration is found to be invariant along this path, at a value close to 10 yr that matches the Earth's core estimates. Despite this invariance, the spatio-temporal properties of secular acceleration show significant variability along the path, with an asymptotic regime of rapid rotation reached after 30 per cent of this path (corresponding to a model Ekman number E = 3 × 10-7). In this regime, the energy of secular acceleration is entirely found at periods longer than that of planetary rotation, and the underlying flow acceleration patterns acquire a 2-D columnar structure representative of the rapid rotation limit. The spatial pattern of the secular acceleration at the core-mantle boundary shows significant localization of energy within an equatorial belt. Rapid hydromagnetic wave dynamics is absent at the start of the path because of insufficient timescale separation with convective processes, weak forcing and excessive damping but can be clearly exhibited in the asymptotic regime. This study reports on ubiquitous axisymmetric geostrophic torsional waves of weak amplitude relatively to convective transport, and also stronger, laterally limited, quasi-geostrophic Alfvén waves propagating in the cylindrical radial direction from the tip of convective plumes towards the core-mantle boundary. In a system similar to the Earth's core where the typical Alfvén velocity is significantly larger than the typical convective velocity, quasi-geostrophic Alfvén waves are shown to be an important carrier of flow acceleration to the core surface that links with the generation of strong, short-lived and intermittent equatorial pulses in the secular acceleration energy. The secular acceleration timescale is shown to be insensitive to magnetic signatures from torsional waves because of their weak amplitude, and from quasi-geostrophic Alfvén waves because of their intermittent character, and is therefore only indicative of convective transport phenomena that remain invariant along the parameter space path.

Spontaneous polarization and dielectric relaxation dynamics of ferroelectric liquid crystals derived from 2(S)-[2(S)-ethylhexyolxy] propionic acid and its (S, R)-diastereomer

NASA Astrophysics Data System (ADS)

Huang, Lei-Ching; Fu, Chao-Ming

2015-09-01

The spontaneous polarization and molecular dynamics of four ferroelectric liquid crystals (FLCs) with two different kinds of core rings and two types of diastereomeric structures were investigated in this study. The FLCs with a biphenyl ring core structure showed higher spontaneous polarization than the FLCs with a naphthalene ring core structure. The complex dielectric spectra exhibited the Goldstone mode in the ferroelectric (SmC*) phase for all FLCs. The complex dielectric spectra of the four FLCs can be optimally fitted by the Debye model and the Cole-Cole model. Moreover, the Goldstone mode was enhanced under low DC bias fields for the FLCs with the (S, R)- diastereomeric structure, whereas the mode was suppressed for the FLCs with the (S, S)- diastereomeric structure. A microscopic molecular dynamic model is proposed to describe the underlying mechanism of the particular enhancement of the Goldstone mode. The experimental results of dielectric spectra and spontaneous polarization are explained in the discussion of the mesomorphic properties related to the FLC molecular structure.

Seismic assessment of WSDOT bridges with prestressed hollow core piles : part II.

DOT National Transportation Integrated Search

2009-12-01

This report investigates the seismic performance of a reinforced concrete : bridge with prestressed hollow core piles. Both nonlinear static and nonlinear dynamic : analyses were carried out. A three-dimensional spine model of the bridge was : ...

Methods of information geometry in computational system biology (consistency between chemical and biological evolution).

PubMed

Astakhov, Vadim

2009-01-01

Interest in simulation of large-scale metabolic networks, species development, and genesis of various diseases requires new simulation techniques to accommodate the high complexity of realistic biological networks. Information geometry and topological formalisms are proposed to analyze information processes. We analyze the complexity of large-scale biological networks as well as transition of the system functionality due to modification in the system architecture, system environment, and system components. The dynamic core model is developed. The term dynamic core is used to define a set of causally related network functions. Delocalization of dynamic core model provides a mathematical formalism to analyze migration of specific functions in biosystems which undergo structure transition induced by the environment. The term delocalization is used to describe these processes of migration. We constructed a holographic model with self-poetic dynamic cores which preserves functional properties under those transitions. Topological constraints such as Ricci flow and Pfaff dimension were found for statistical manifolds which represent biological networks. These constraints can provide insight on processes of degeneration and recovery which take place in large-scale networks. We would like to suggest that therapies which are able to effectively implement estimated constraints, will successfully adjust biological systems and recover altered functionality. Also, we mathematically formulate the hypothesis that there is a direct consistency between biological and chemical evolution. Any set of causal relations within a biological network has its dual reimplementation in the chemistry of the system environment.

Progress in the Development of a Global Quasi-3-D Multiscale Modeling Framework

NASA Astrophysics Data System (ADS)

Jung, J.; Konor, C. S.; Randall, D. A.

2017-12-01

The Quasi-3-D Multiscale Modeling Framework (Q3D MMF) is a second-generation MMF, which has following advances over the first-generation MMF: 1) The cloud-resolving models (CRMs) that replace conventional parameterizations are not confined to the large-scale dynamical-core grid cells, and are seamlessly connected to each other, 2) The CRMs sense the three-dimensional large- and cloud-scale environment, 3) Two perpendicular sets of CRM channels are used, and 4) The CRMs can resolve the steep surface topography along the channel direction. The basic design of the Q3D MMF has been developed and successfully tested in a limited-area modeling framework. Currently, global versions of the Q3D MMF are being developed for both weather and climate applications. The dynamical cores governing the large-scale circulation in the global Q3D MMF are selected from two cube-based global atmospheric models. The CRM used in the model is the 3-D nonhydrostatic anelastic Vector-Vorticity Model (VVM), which has been tested with the limited-area version for its suitability for this framework. As a first step of the development, the VVM has been reconstructed on the cubed-sphere grid so that it can be applied to global channel domains and also easily fitted to the large-scale dynamical cores. We have successfully tested the new VVM by advecting a bell-shaped passive tracer and simulating the evolutions of waves resulted from idealized barotropic and baroclinic instabilities. For improvement of the model, we also modified the tracer advection scheme to yield positive-definite results and plan to implement a new physics package that includes a double-moment microphysics and an aerosol physics. The interface for coupling the large-scale dynamical core and the VVM is under development. In this presentation, we shall describe the recent progress in the development and show some test results.

Comparison of the Tangent Linear Properties of Tracer Transport Schemes Applied to Geophysical Problems.

NASA Technical Reports Server (NTRS)

Kent, James; Holdaway, Daniel

2015-01-01

A number of geophysical applications require the use of the linearized version of the full model. One such example is in numerical weather prediction, where the tangent linear and adjoint versions of the atmospheric model are required for the 4DVAR inverse problem. The part of the model that represents the resolved scale processes of the atmosphere is known as the dynamical core. Advection, or transport, is performed by the dynamical core. It is a central process in many geophysical applications and is a process that often has a quasi-linear underlying behavior. However, over the decades since the advent of numerical modelling, significant effort has gone into developing many flavors of high-order, shape preserving, nonoscillatory, positive definite advection schemes. These schemes are excellent in terms of transporting the quantities of interest in the dynamical core, but they introduce nonlinearity through the use of nonlinear limiters. The linearity of the transport schemes used in Goddard Earth Observing System version 5 (GEOS-5), as well as a number of other schemes, is analyzed using a simple 1D setup. The linearized version of GEOS-5 is then tested using a linear third order scheme in the tangent linear version.

Early Dynamics of the Moon's Core

NASA Astrophysics Data System (ADS)

Cuk, Matija; Hamilton, Douglas; Stewart, Sarah T.

2018-04-01

The Moon has a small molten iron core (Williams et al. 2006). Remanent magnetization in lunar rocks likely derives from a past lunar dynamo (Wieczorek 2018 and references therein), which may have been powered by differential precession between the mantle and the core. The rotations of the lunar mantle and core were largely decoupled for much of lunar history, with a large mutual offset during the Cassini State Transition (Meyer and Wisdom, 2011). It is likely that the past work underestimated lunar obliquities, and therefore core offsets, during early lunar history (Cuk et al. 2016). Here we investigate the dynamics of the lunar core and mantle using a Lie-Poisson numerical integrator (Touma and Wisdom 2001) which includes interactions between triaxial core and mantle, as well as all gravitational and tidal effects included in the model of Cuk et al. (2016). Since we assume a rigid triaxial mantle, this model is applicable to the Moon only once it has acquired its current shape, which probably happened before the Moon reached 25 Earth radii. While some details of the core dynamics depend on our assumptions about the shape of the lunar core-mantle boundary, we can report some robust preliminary findings. The presence of the core does not change significantly the evolutionary scenario of Cuk et al. (2016). The core and mantle are indeed decoupled, with the core having a much smaller obliquity to the ecliptic than the mantle for almost all of the lunar history. The core was largely in an equivalent of Cassini State 2, with the vernal equinoxes (wrt the ecliptic) of the core and the mantle being anti-aligned. The core-mantle spin axis offset has been very large during the Moon's first billion years (this is true both in canonical and high-inclination tidal evolution), causing the lunar core to be sub-synchronous. If the ancient lunar magnetic dipole was rotating around the core axis that was inclined to the Moon's spin axis, then the magnetic poles would move across the lunar surface as the mantle rotates independently. This relative motion would dilute the average dipole field over much of the lunar surface, and would would restrict meaningful average fields to low lunar latitudes.

Modeled and Measured Dynamics of a Composite Beam with Periodically Varying Foam Core

NASA Technical Reports Server (NTRS)

Cabell, Randolph H.; Cano, Roberto J.; Schiller, Noah H.; Roberts Gary D.

2012-01-01

The dynamics of a sandwich beam with carbon fiber composite facesheets and foam core with periodic variations in material properties are studied. The purpose of the study is to compare finite element predictions with experimental measurements on fabricated beam specimens. For the study, three beams were fabricated: one with a compliant foam core, a second with a stiffer core, and a third with the two cores alternating down the length of the beam to create a periodic variation in properties. This periodic variation produces a bandgap in the frequency domain where vibrational energy does not readily propagate down the length of the beam. Mode shapes and natural frequencies are compared, as well as frequency responses from point force input to velocity response at the opposite end of the beam.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Helled, Ravit

We explore the change in Jupiter's normalized axial moment of inertia (NMOI) assuming that Jupiter undergoes core erosion. It is found that Jupiter's contraction combined with an erosion of 20 M{sub Circled-Plus} from a primordial core of 30 M{sub Circled-Plus} can significantly change Jupiter's NMOI over time. It is shown that Jupiter's NMOI could have changed from {approx}0.235 to {approx}0.264 throughout its evolution. We find that an NMOI value of {approx}0.235 as suggested by dynamical models could, in principle, be consistent with Jupiter's primordial internal structure. Low NMOI values, however, persist only for the first {approx}10{sup 6} years of Jupiter'smore » evolution. Re-evaluation of dynamical stability models as well as more sophisticated evolution models of Jupiter with core erosion seem to be required in order to provide more robust estimates for Jupiter's primordial NMOI.« less

Geomagnetic inverse problem and data assimilation: a progress report

NASA Astrophysics Data System (ADS)

Aubert, Julien; Fournier, Alexandre

2013-04-01

In this presentation I will present two studies recently undertaken by our group in an effort to bring the benefits of data assimilation to the study of Earth's magnetic field and the dynamics of its liquid iron core, where the geodynamo operates. In a first part I will focus on the geomagnetic inverse problem, which attempts to recover the fluid flow in the core from the temporal variation of the magnetic field (known as the secular variation). Geomagnetic data can be downward continued from the surface of the Earth down to the core-mantle boundary, but not further below, since the core is an electrical conductor. Historically, solutions to the geomagnetic inverse problem in such a sparsely observed system were thus found only for flow immediately below the core mantle boundary. We have recently shown that combining a numerical model of the geodynamo together with magnetic observations, through the use of Kalman filtering, now allows to present solutions for flow throughout the core. In a second part, I will present synthetic tests of sequential geomagnetic data assimilation aiming at evaluating the range at which the future of the geodynamo can be predicted, and our corresponding prospects to refine the current geomagnetic predictions. Fournier, Aubert, Thébault: Inference on core surface flow from observations and 3-D dynamo modelling, Geophys. J. Int. 186, 118-136, 2011, doi: 10.1111/j.1365-246X.2011.05037.x Aubert, Fournier: Inferring internal properties of Earth's core dynamics and their evolution from surface observations and a numerical geodynamo model, Nonlinear Proc. Geoph. 18, 657-674, 2011, doi:10.5194/npg-18-657-2011 Aubert: Flow throughout the Earth's core inverted from geomagnetic observations and numerical dynamo models, Geophys. J. Int., 2012, doi: 10.1093/gji/ggs051

The Tangent Linear and Adjoint of the FV3 Dynamical Core: Development and Applications

NASA Technical Reports Server (NTRS)

Holdaway, Daniel

2018-01-01

GMAO (NASA's Global Modeling and Assimilation Office) has developed a highly sophisticated adjoint modeling system based on the most recent version of the finite volume cubed sphere (FV3) dynamical core. This provides a mechanism for investigating sensitivity to initial conditions and examining observation impacts. It also allows for the computation of singular vectors and for the implementation of hybrid 4DVAR (4-Dimensional Variational Assimilation). In this work we will present the scientific assessment of the new adjoint system and show results from a number of research application of the adjoint system.

Flexible histone tails in a new mesoscopic oligonucleosome model.

PubMed

Arya, Gaurav; Zhang, Qing; Schlick, Tamar

2006-07-01

We describe a new mesoscopic model of oligonucleosomes that incorporates flexible histone tails. The nucleosome cores are modeled using the discrete surface-charge optimization model, which treats the nucleosome as an electrostatic surface represented by hundreds of point charges; the linker DNAs are treated using a discrete elastic chain model; and the histone tails are modeled using a bead/chain hydrodynamic approach as chains of connected beads where each bead represents five protein residues. Appropriate charges and force fields are assigned to each histone chain so as to reproduce the electrostatic potential, structure, and dynamics of the corresponding atomistic histone tails at different salt conditions. The dynamics of resulting oligonucleosomes at different sizes and varying salt concentrations are simulated by Brownian dynamics with complete hydrodynamic interactions. The analyses demonstrate that the new mesoscopic model reproduces experimental results better than its predecessors, which modeled histone tails as rigid entities. In particular, our model with flexible histone tails: correctly accounts for salt-dependent conformational changes in the histone tails; yields the experimentally obtained values of histone-tail mediated core/core attraction energies; and considers the partial shielding of electrostatic repulsion between DNA linkers as a result of the spatial distribution of histone tails. These effects are crucial for regulating chromatin structure but are absent or improperly treated in models with rigid histone tails. The development of this model of oligonucleosomes thus opens new avenues for studying the role of histone tails and their variants in mediating gene expression through modulation of chromatin structure.

Planetary and Primitive Object Strength Measurement and Sampling Apparatus

NASA Technical Reports Server (NTRS)

Ahrens, Thomas J.

1995-01-01

Support is requested for continuation of a program of dynamic impact (harpoon) coring of planetary, comet, or asteroid surface materials. We have previously demonstrated that good quality cores are obtainable for planetary materials with compressive strengths less than 200 MPa. Since the dynamics of penetration are observable on a Discovery class spacecraft, which images the sampling operation, these data can be used with a model developed under this project, to measure in-situ strength and frictional strength of the crust of the object. During the last year we have developed a detailed analytic model of penetrator mechanics. Progress is reported for the solid penetrators experiments, the CIT penetrator model, and the impact spall sampling apparatus.

DIFFERENT DYNAMICAL AGES FOR THE TWO YOUNG AND COEVAL LMC STAR CLUSTERS, NGC 1805 AND NGC 1818, IMPRINTED ON THEIR BINARY POPULATIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Geller, Aaron M.; Grijs, Richard de; Li, Chengyuan

2015-05-20

The two Large Magellanic Cloud star clusters, NGC 1805 and NGC 1818, are approximately the same chronological age (∼30 Myr), but show different radial trends in binary frequency. The F-type stars (1.3–2.2 M{sub ⊙}) in NGC 1818 have a binary frequency that decreases toward the core, while the binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with a peak in the core). We show here, through detailed N-body modeling, that both clusters could have formed with the same primordial binary frequency and with binary orbital elements and masses drawn from themore » same distributions (defined from observations of open clusters and the field of our Galaxy). The observed radial trends in binary frequency for both clusters are best matched with models that have initial substructure. Furthermore, both clusters may be evolving along a very similar dynamical sequence, with the key difference that NGC 1805 is dynamically older than NGC 1818. The F-type binaries in NGC 1818 still show evidence of an initial period of rapid dynamical disruptions (which occur preferentially in the core), while NGC 1805 has already begun to recover a higher core binary frequency, owing to mass segregation (which will eventually produce a distribution in binary frequency that rises only toward the core, as is observed in old Milky Way star clusters). This recovery rate increases for higher-mass binaries, and therefore even at one age in one cluster, we predict a similar dynamical sequence in the radial distribution of the binary frequency as a function of binary primary mass.« less

Gravity or turbulence? IV. Collapsing cores in out-of-virial disguise

NASA Astrophysics Data System (ADS)

Ballesteros-Paredes, Javier; Vázquez-Semadeni, Enrique; Palau, Aina; Klessen, Ralf S.

2018-06-01

We study the dynamical state of massive cores by using a simple analytical model, an observational sample, and numerical simulations of collapsing massive cores. From the analytical model, we find that cores increase their column density and velocity dispersion as they collapse, resulting in a time evolution path in the Larson velocity dispersion-size diagram from large sizes and small velocity dispersions to small sizes and large velocity dispersions, while they tend to equipartition between gravity and kinetic energy. From the observational sample, we find that: (a) cores with substantially different column densities in the sample do not follow a Larson-like linewidth-size relation. Instead, cores with higher column densities tend to be located in the upper-left corner of the Larson velocity dispersion σv, 3D-size R diagram, a result explained in the hierarchical and chaotic collapse scenario. (b) Cores appear to have overvirial values. Finally, our numerical simulations reproduce the behavior predicted by the analytical model and depicted in the observational sample: collapsing cores evolve towards larger velocity dispersions and smaller sizes as they collapse and increase their column density. More importantly, however, they exhibit overvirial states. This apparent excess is due to the assumption that the gravitational energy is given by the energy of an isolated homogeneous sphere. However, such excess disappears when the gravitational energy is correctly calculated from the actual spatial mass distribution. We conclude that the observed energy budget of cores is consistent with their non-thermal motions being driven by their self-gravity and in the process of dynamical collapse.

Numerical simulation of the geodynamo reaches Earth's core dynamical regime

NASA Astrophysics Data System (ADS)

Aubert, J.; Gastine, T.; Fournier, A.

2016-12-01

Numerical simulations of the geodynamo have been successful at reproducing a number of static (field morphology) and kinematic (secular variation patterns, core surface flows and westward drift) features of Earth's magnetic field, making them a tool of choice for the analysis and retrieval of geophysical information on Earth's core. However, classical numerical models have been run in a parameter regime far from that of the real system, prompting the question of whether we do get "the right answers for the wrong reasons", i.e. whether the agreement between models and nature simply occurs by chance and without physical relevance in the dynamics. In this presentation, we show that classical models succeed in describing the geodynamo because their large-scale spatial structure is essentially invariant as one progresses along a well-chosen path in parameter space to Earth's core conditions. This path is constrained by the need to enforce the relevant force balance (MAC or Magneto-Archimedes-Coriolis) and preserve the ratio of the convective overturn and magnetic diffusion times. Numerical simulations performed along this path are shown to be spatially invariant at scales larger than that where the magnetic energy is ohmically dissipated. This property enables the definition of large-eddy simulations that show good agreement with direct numerical simulations in the range where both are feasible, and that can be computed at unprecedented values of the control parameters, such as an Ekman number E=10-8. Combining direct and large-eddy simulations, large-scale invariance is observed over half the logarithmic distance in parameter space between classical models and Earth. The conditions reached at this mid-point of the path are furthermore shown to be representative of the rapidly-rotating, asymptotic dynamical regime in which Earth's core resides, with a MAC force balance undisturbed by viscosity or inertia, the enforcement of a Taylor state and strong-field dynamo action. We conclude that numerical modelling has advanced to a stage where it is possible to use models correctly representing the statics, kinematics and now the dynamics of the geodynamo. This opens the way to a better analysis of the geomagnetic field in the time and space domains.

A New Dynamical Core Based on the Prediction of the Curl of the Horizontal Vorticity

NASA Astrophysics Data System (ADS)

Konor, C. S.; Randall, D. A.; Heikes, R. P.

2015-12-01

The Vector-Vorticity Dynamical core (VVM) developed by Jung and Arakawa (2008) has important advantages for the use with the anelastic and unified systems of equations. The VVM predicts the horizontal vorticity vector (HVV) at each interface and the vertical vorticity at the top layer of the model. To guarantee that the three-dimensional vorticity is nondivergent, the vertical vorticity at the interior layers is diagnosed from the horizontal divergence of the HVV through a vertical integral from the top to down. To our knowledge, this is the only dynamical core that guarantees the nondivergence of the three-dimensional vorticity. The VVM uses a C-type horizontal grid, which allows a computational mode. While the computational mode does not seem to be serious in the Cartesian grid applications, it may be serious in the icosahedral grid applications because of the extra degree of freedom in such grids. Although there are special filters to minimize the effects of this computational mode, we prefer to eliminate it altogether. We have developed a new dynamical core, which uses a Z-grid to avoid the computational mode mentioned above. The dynamical core predicts the curl of the HVV and diagnoses the horizontal divergence of the HVV from the predicted vertical vorticity. The three-dimensional vorticity is guaranteed to be nondivergent as in the VVM. In this presentation, we will introduce the new dynamical core and show results obtained by using Cartesian and hexagonal grids. We will also compare the solutions to that obtained by the VVM.

Toward Connecting Core-Collapse Supernova Theory with Observations: Nucleosynthetic Yields and Distribution of Elements in a 15 M⊙ Blue Supergiant Progenitor with SN 1987A Energetics

NASA Astrophysics Data System (ADS)

Plewa, Tomasz; Handy, Timothy; Odrzywolek, Andrzej

2014-09-01

We compute and discuss the process of nucleosynthesis in a series of core-collapse explosion models of a 15 solar mass, blue supergiant progenitor. We obtain nucleosynthetic yields and study the evolution of the chemical element distribution from the moment of core bounce until young supernova remnant phase. Our models show how the process of energy deposition due to radioactive decay modifies the dynamics and the core ejecta structure on small and intermediate scales. The results are compared against observations of young supernova remnants including Cas A and the recent data obtained for SN 1987A. We compute and discuss the process of nucleosynthesis in a series of core-collapse explosion models of a 15 solar mass, blue supergiant progenitor. We obtain nucleosynthetic yields and study the evolution of the chemical element distribution from the moment of core bounce until young supernova remnant phase. Our models show how the process of energy deposition due to radioactive decay modifies the dynamics and the core ejecta structure on small and intermediate scales. The results are compared against observations of young supernova remnants including Cas A and the recent data obtained for SN 1987A. The work has been supported by the NSF grant AST-1109113 and DOE grant DE-FG52-09NA29548. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the U.S. DoE under Contract No. DE-AC02-05CH11231.

Large-scale Instability during Gravitational Collapse with Neutrino Transport and a Core-Collapse Supernova

NASA Astrophysics Data System (ADS)

Aksenov, A. G.; Chechetkin, V. M.

2018-04-01

Most of the energy released in the gravitational collapse of the cores of massive stars is carried away by neutrinos. Neutrinos play a pivotal role in explaining core-collape supernovae. Currently, mathematical models of the gravitational collapse are based on multi-dimensional gas dynamics and thermonuclear reactions, while neutrino transport is considered in a simplified way. Multidimensional gas dynamics is used with neutrino transport in the flux-limited diffusion approximation to study the role of multi-dimensional effects. The possibility of large-scale convection is discussed, which is interesting both for explaining SN II and for setting up observations to register possible high-energy (≳10MeV) neutrinos from the supernova. A new multi-dimensional, multi-temperature gas dynamics method with neutrino transport is presented.

MCore: A High-Order Finite-Volume Dynamical Core for Atmospheric General Circulation Models

NASA Astrophysics Data System (ADS)

Ullrich, P.; Jablonowski, C.

2011-12-01

The desire for increasingly accurate predictions of the atmosphere has driven numerical models to smaller and smaller resolutions, while simultaneously exponentially driving up the cost of existing numerical models. Even with the modern rapid advancement of computational performance, it is estimated that it will take more than twenty years before existing models approach the scales needed to resolve atmospheric convection. However, smarter numerical methods may allow us to glimpse the types of results we would expect from these fine-scale simulations while only requiring a fraction of the computational cost. The next generation of atmospheric models will likely need to rely on both high-order accuracy and adaptive mesh refinement in order to properly capture features of interest. We present our ongoing research on developing a set of ``smart'' numerical methods for simulating the global non-hydrostatic fluid equations which govern atmospheric motions. We have harnessed a high-order finite-volume based approach in developing an atmospheric dynamical core on the cubed-sphere. This type of method is desirable for applications involving adaptive grids, since it has been shown that spuriously reflected wave modes are intrinsically damped out under this approach. The model further makes use of an implicit-explicit Runge-Kutta-Rosenbrock (IMEX-RKR) time integrator for accurate and efficient coupling of the horizontal and vertical model components. We survey the algorithmic development of the model and present results from idealized dynamical core test cases, as well as give a glimpse at future work with our model.

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

NASA Technical Reports Server (NTRS)

Lin, Shian-Jiann; Atlas, Robert (Technical Monitor)

2002-01-01

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 km or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed-shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

NASA Technical Reports Server (NTRS)

Lin, Shian-Jiann; Atlas, Robert (Technical Monitor)

2002-01-01

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 kin or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed- shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

A moist aquaplanet variant of the Held–Suarez test for atmospheric model dynamical cores

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thatcher, Diana R.; Jablonowski, Christiane

A moist idealized test case (MITC) for atmospheric model dynamical cores is presented. The MITC is based on the Held–Suarez (HS) test that was developed for dry simulations on “a flat Earth” and replaces the full physical parameterization package with a Newtonian temperature relaxation and Rayleigh damping of the low-level winds. This new variant of the HS test includes moisture and thereby sheds light on the nonlinear dynamics–physics moisture feedbacks without the complexity of full-physics parameterization packages. In particular, it adds simplified moist processes to the HS forcing to model large-scale condensation, boundary-layer mixing, and the exchange of latent and sensible heat betweenmore » the atmospheric surface and an ocean-covered planet. Using a variety of dynamical cores of the National Center for Atmospheric Research (NCAR)'s Community Atmosphere Model (CAM), this paper demonstrates that the inclusion of the moist idealized physics package leads to climatic states that closely resemble aquaplanet simulations with complex physical parameterizations. This establishes that the MITC approach generates reasonable atmospheric circulations and can be used for a broad range of scientific investigations. This paper provides examples of two application areas. First, the test case reveals the characteristics of the physics–dynamics coupling technique and reproduces coupling issues seen in full-physics simulations. In particular, it is shown that sudden adjustments of the prognostic fields due to moist physics tendencies can trigger undesirable large-scale gravity waves, which can be remedied by a more gradual application of the physical forcing. Second, the moist idealized test case can be used to intercompare dynamical cores. These examples demonstrate the versatility of the MITC approach and suggestions are made for further application areas. Furthermore, the new moist variant of the HS test can be considered a test case of intermediate complexity.« less

A moist aquaplanet variant of the Held–Suarez test for atmospheric model dynamical cores

DOE PAGES

Thatcher, Diana R.; Jablonowski, Christiane

2016-04-04

A moist idealized test case (MITC) for atmospheric model dynamical cores is presented. The MITC is based on the Held–Suarez (HS) test that was developed for dry simulations on “a flat Earth” and replaces the full physical parameterization package with a Newtonian temperature relaxation and Rayleigh damping of the low-level winds. This new variant of the HS test includes moisture and thereby sheds light on the nonlinear dynamics–physics moisture feedbacks without the complexity of full-physics parameterization packages. In particular, it adds simplified moist processes to the HS forcing to model large-scale condensation, boundary-layer mixing, and the exchange of latent and sensible heat betweenmore » the atmospheric surface and an ocean-covered planet. Using a variety of dynamical cores of the National Center for Atmospheric Research (NCAR)'s Community Atmosphere Model (CAM), this paper demonstrates that the inclusion of the moist idealized physics package leads to climatic states that closely resemble aquaplanet simulations with complex physical parameterizations. This establishes that the MITC approach generates reasonable atmospheric circulations and can be used for a broad range of scientific investigations. This paper provides examples of two application areas. First, the test case reveals the characteristics of the physics–dynamics coupling technique and reproduces coupling issues seen in full-physics simulations. In particular, it is shown that sudden adjustments of the prognostic fields due to moist physics tendencies can trigger undesirable large-scale gravity waves, which can be remedied by a more gradual application of the physical forcing. Second, the moist idealized test case can be used to intercompare dynamical cores. These examples demonstrate the versatility of the MITC approach and suggestions are made for further application areas. Furthermore, the new moist variant of the HS test can be considered a test case of intermediate complexity.« less

A seismologically consistent compositional model of Earth's core.

PubMed

Badro, James; Côté, Alexander S; Brodholt, John P

2014-05-27

Earth's core is less dense than iron, and therefore it must contain "light elements," such as S, Si, O, or C. We use ab initio molecular dynamics to calculate the density and bulk sound velocity in liquid metal alloys at the pressure and temperature conditions of Earth's outer core. We compare the velocity and density for any composition in the (Fe-Ni, C, O, Si, S) system to radial seismological models and find a range of compositional models that fit the seismological data. We find no oxygen-free composition that fits the seismological data, and therefore our results indicate that oxygen is always required in the outer core. An oxygen-rich core is a strong indication of high-pressure and high-temperature conditions of core differentiation in a deep magma ocean with an FeO concentration (oxygen fugacity) higher than that of the present-day mantle.

A seismologically consistent compositional model of Earth’s core

PubMed Central

Badro, James; Côté, Alexander S.; Brodholt, John P.

2014-01-01

Earth’s core is less dense than iron, and therefore it must contain “light elements,” such as S, Si, O, or C. We use ab initio molecular dynamics to calculate the density and bulk sound velocity in liquid metal alloys at the pressure and temperature conditions of Earth's outer core. We compare the velocity and density for any composition in the (Fe–Ni, C, O, Si, S) system to radial seismological models and find a range of compositional models that fit the seismological data. We find no oxygen-free composition that fits the seismological data, and therefore our results indicate that oxygen is always required in the outer core. An oxygen-rich core is a strong indication of high-pressure and high-temperature conditions of core differentiation in a deep magma ocean with an FeO concentration (oxygen fugacity) higher than that of the present-day mantle. PMID:24821817

Dynamics of core accretion

DOE PAGES

Nelson, Andrew F.; Ruffert, Maximilian

2012-12-21

In this paper, we perform three-dimensional hydrodynamic simulations of gas flowing around a planetary core of mass M pl = 10M ⊕ embedded in a near Keplerian background flow, using a modified shearing box approximation. We assume an ideal gas behaviour following an equation of state with a fixed ratio of the specific heats, γ = 1.42, consistent with the conditions of a moderate-temperature background disc with solar composition. No radiative heating or cooling is included in the models. We employ a nested grid hydrodynamic code implementing the ‘Piecewise Parabolic Method’ with as many as six fixed nested grids, providingmore » spatial resolution on the finest grid comparable to the present-day diameters of Neptune and Uranus. We find that a strongly dynamically active flow develops such that no static envelope can form. The activity is not sensitive to plausible variations in the rotation curve of the underlying disc. It is sensitive to the thermodynamic treatment of the gas, as modelled by prescribed equations of state (either ‘locally isothermal’ or ‘locally isentropic’) and the temperature of the background disc material. The activity is also sensitive to the shape and depth of the core's gravitational potential, through its mass and gravitational softening coefficient. Each of these factors influences the magnitude and character of hydrodynamic feedback of the small-scale flow on the background, and we conclude that accurate modelling of such feedback is critical to a complete understanding of the core accretion process. The varying flow pattern gives rise to large, irregular eruptions of matter from the region around the core which return matter to the background flow: mass in the envelope at one time may not be found in the envelope at any later time. No net mass accretion into the envelope is observed over the course of the simulation and none is expected, due to our neglect of cooling. Except in cases of very rapid cooling however, as defined by locally isothermal or isentropic treatments, any cooling that does affect the envelope material will have limited consequences for the dynamics, since the flow quickly carries cooled material out of the core's environment entirely. The angular momentum of material in the envelope, relative to the core, varies both in magnitude and in sign on time-scales of days to months near the core and on time-scales a few years at distances comparable to the Hill radius. The dynamical activity contrasts with the largely static behaviour typically assumed within the framework of the core accretion model for Jovian planet formation. We show that material entering the dynamically active environment may suffer intense heating and cooling events the durations of which are as short as a few hours to a few days. Shorter durations are not observable in our work due to the limits of our resolution. Peak temperatures in these events range from T ~ 1000 K to as high as T ~ 3–4000 K, with densities ρ ~ 10 -9 to 10 -8 gcm -3. These time-scales, densities and temperatures span a range consistent with those required for chondrule formation in the nebular shock model. Finally, we therefore propose that dynamical activity in the Jovian planet formation environment could be responsible for the production of chondrules and other annealed silicates in the solar nebula.« less

Dynamics of core accretion

NASA Astrophysics Data System (ADS)

Nelson, Andrew F.; Ruffert, Maximilian

2013-02-01

We perform three-dimensional hydrodynamic simulations of gas flowing around a planetary core of mass Mpl = 10M⊕ embedded in a near Keplerian background flow, using a modified shearing box approximation. We assume an ideal gas behaviour following an equation of state with a fixed ratio of the specific heats, γ = 1.42, consistent with the conditions of a moderate-temperature background disc with solar composition. No radiative heating or cooling is included in the models. We employ a nested grid hydrodynamic code implementing the `Piecewise Parabolic Method' with as many as six fixed nested grids, providing spatial resolution on the finest grid comparable to the present-day diameters of Neptune and Uranus. We find that a strongly dynamically active flow develops such that no static envelope can form. The activity is not sensitive to plausible variations in the rotation curve of the underlying disc. It is sensitive to the thermodynamic treatment of the gas, as modelled by prescribed equations of state (either `locally isothermal' or `locally isentropic') and the temperature of the background disc material. The activity is also sensitive to the shape and depth of the core's gravitational potential, through its mass and gravitational softening coefficient. Each of these factors influences the magnitude and character of hydrodynamic feedback of the small-scale flow on the background, and we conclude that accurate modelling of such feedback is critical to a complete understanding of the core accretion process. The varying flow pattern gives rise to large, irregular eruptions of matter from the region around the core which return matter to the background flow: mass in the envelope at one time may not be found in the envelope at any later time. No net mass accretion into the envelope is observed over the course of the simulation and none is expected, due to our neglect of cooling. Except in cases of very rapid cooling however, as defined by locally isothermal or isentropic treatments, any cooling that does affect the envelope material will have limited consequences for the dynamics, since the flow quickly carries cooled material out of the core's environment entirely. The angular momentum of material in the envelope, relative to the core, varies both in magnitude and in sign on time-scales of days to months near the core and on time-scales a few years at distances comparable to the Hill radius. The dynamical activity contrasts with the largely static behaviour typically assumed within the framework of the core accretion model for Jovian planet formation. We show that material entering the dynamically active environment may suffer intense heating and cooling events the durations of which are as short as a few hours to a few days. Shorter durations are not observable in our work due to the limits of our resolution. Peak temperatures in these events range from T ˜ 1000 K to as high as T ˜ 3-4000 K, with densities ρ ˜ 10-9 to 10-8 g cm-3. These time-scales, densities and temperatures span a range consistent with those required for chondrule formation in the nebular shock model. We therefore propose that dynamical activity in the Jovian planet formation environment could be responsible for the production of chondrules and other annealed silicates in the solar nebula.

An accelerating high-latitude jet in Earth's core

NASA Astrophysics Data System (ADS)

Finlay, C. C.; Livermore, P. W.; Hollerbach, R.

2016-12-01

The structure of the core-generated magnetic field, and how it changes in time (its secular variation or SV), supplies an invaluable constraint on the dynamics of the outer core. At high latitude, previous studies have noted distinctive behaviour of secular change, in particular suggesting a polar vortex tied to the dynamics within the tangent cylinder region. Recent high-resolution observational models that include data from the Swarm satellites have refined the structure of observed SV, to a rapidly changing circular daisy-chain configuration centred on the north geographic pole, on or very close to the tangent cylinder itself. Motivated by theoretical considerations of the likely dynamical regime of the core, we demonstrate that this feature can be explained by a localised westwards cylindrical jet of 420 km width centred the tangent cylinder, whose amplitude appears to have increased in strength by a factor of three over the period 2000-2016 to about 40 km/yr. The current accelerating phase may be a short fragment of decadal fluctuations of the jet strength linked to both torsional wave activity and the rotation direction of the inner core.

Some conservation issues for the dynamical cores of NWP and climate models

NASA Astrophysics Data System (ADS)

Thuburn, J.

2008-03-01

The rationale for designing atmospheric numerical model dynamical cores with certain conservation properties is reviewed. The conceptual difficulties associated with the multiscale nature of realistic atmospheric flow, and its lack of time-reversibility, are highlighted. A distinction is made between robust invariants, which are conserved or nearly conserved in the adiabatic and frictionless limit, and non-robust invariants, which are not conserved in the limit even though they are conserved by exactly adiabatic frictionless flow. For non-robust invariants, a further distinction is made between processes that directly transfer some quantity from large to small scales, and processes involving a cascade through a continuous range of scales; such cascades may either be explicitly parameterized, or handled implicitly by the dynamical core numerics, accepting the implied non-conservation. An attempt is made to estimate the relative importance of different conservation laws. It is argued that satisfactory model performance requires spurious sources of a conservable quantity to be much smaller than any true physical sources; for several conservable quantities the magnitudes of the physical sources are estimated in order to provide benchmarks against which any spurious sources may be measured.

An intercomparison of GCM and RCM dynamical downscaling for characterizing the hydroclimatology of California and Nevada

NASA Astrophysics Data System (ADS)

Xu, Z.; Rhoades, A.; Johansen, H.; Ullrich, P. A.; Collins, W. D.

2017-12-01

Dynamical downscaling is widely used to properly characterize regional surface heterogeneities that shape the local hydroclimatology. However, the factors in dynamical downscaling, including the refinement of model horizontal resolution, large-scale forcing datasets and dynamical cores, have not been fully evaluated. Two cutting-edge global-to-regional downscaling methods are used to assess these, specifically the variable-resolution Community Earth System Model (VR-CESM) and the Weather Research & Forecasting (WRF) regional climate model, under different horizontal resolutions (28, 14, and 7 km). Two groups of WRF simulations are driven by either the NCEP reanalysis dataset (WRF_NCEP) or VR-CESM outputs (WRF_VRCESM) to evaluate the effects of the large-scale forcing datasets. The impacts of dynamical core are assessed by comparing the VR-CESM simulations to the coupled WRF_VRCESM simulations with the same physical parameterizations and similar grid domains. The simulated hydroclimatology (i.e., total precipitation, snow cover, snow water equivalent and surface temperature) are compared with the reference datasets. The large-scale forcing datasets are critical to the WRF simulations in more accurately simulating total precipitation, SWE and snow cover, but not surface temperature. Both the WRF and VR-CESM results highlight that no significant benefit is found in the simulated hydroclimatology by just increasing horizontal resolution refinement from 28 to 7 km. Simulated surface temperature is sensitive to the choice of dynamical core. WRF generally simulates higher temperatures than VR-CESM, alleviates the systematic cold bias of DJF temperatures over the California mountain region, but overestimates the JJA temperature in California's Central Valley.

Bioinactivation: Software for modelling dynamic microbial inactivation.

PubMed

Garre, Alberto; Fernández, Pablo S; Lindqvist, Roland; Egea, Jose A

2017-03-01

This contribution presents the bioinactivation software, which implements functions for the modelling of isothermal and non-isothermal microbial inactivation. This software offers features such as user-friendliness, modelling of dynamic conditions, possibility to choose the fitting algorithm and generation of prediction intervals. The software is offered in two different formats: Bioinactivation core and Bioinactivation SE. Bioinactivation core is a package for the R programming language, which includes features for the generation of predictions and for the fitting of models to inactivation experiments using non-linear regression or a Markov Chain Monte Carlo algorithm (MCMC). The calculations are based on inactivation models common in academia and industry (Bigelow, Peleg, Mafart and Geeraerd). Bioinactivation SE supplies a user-friendly interface to selected functions of Bioinactivation core, namely the model fitting of non-isothermal experiments and the generation of prediction intervals. The capabilities of bioinactivation are presented in this paper through a case study, modelling the non-isothermal inactivation of Bacillus sporothermodurans. This study has provided a full characterization of the response of the bacteria to dynamic temperature conditions, including confidence intervals for the model parameters and a prediction interval of the survivor curve. We conclude that the MCMC algorithm produces a better characterization of the biological uncertainty and variability than non-linear regression. The bioinactivation software can be relevant to the food and pharmaceutical industry, as well as to regulatory agencies, as part of a (quantitative) microbial risk assessment. Copyright © 2017 Elsevier Ltd. All rights reserved.

The Dynamics of Massive Starless Cores with ALMA

NASA Astrophysics Data System (ADS)

Tan, Jonathan C.; Kong, Shuo; Butler, Michael J.; Caselli, Paola; Fontani, Francesco

2013-12-01

How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The turbulent core accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their surrounding clump medium. Their internal pressure is provided by a combination of turbulence and magnetic fields. Alternatively, the competitive accretion model requires strongly sub-virial initial conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate- and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we have identified four prime examples of massive (~100 M ⊙) clumps from mid-infrared extinction mapping of infrared dark clouds. Fontani et al. found high deuteration fractions of N2H+ in these objects, which are consistent with them being starless. Here we present ALMA observations of these four clumps that probe the N2D+ (3-2) line at 2.''3 resolution. We find six N2D+ cores and determine their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the predictions of the turbulent core model of self-gravitating, magnetized (with Alfvén Mach number mA ~ 1) and virialized cores that are bounded by the high pressures of their surrounding clumps. However, in the most massive cores, with masses up to ~60 M ⊙, our results suggest that moderately enhanced magnetic fields (so that mA ~= 0.3) may be needed for the structures to be in virial and pressure equilibrium. Magnetically regulated core formation may thus be important in controlling the formation of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar IMF.

Experimental and Theoretical Investigations on Viscosity of Fe-Ni-C Liquids at High Pressures

NASA Astrophysics Data System (ADS)

Chen, B.; Lai, X.; Wang, J.; Zhu, F.; Liu, J.; Kono, Y.

2016-12-01

Understanding and modeling of Earth's core processes such as geodynamo and heat flow via convection in liquid outer cores hinges on the viscosity of candidate liquid iron alloys under core conditions. Viscosity estimates from various methods of the metallic liquid of the outer core, however, span up to 12 orders of magnitude. Due to experimental challenges, viscosity measurements of iron liquids alloyed with lighter elements are scarce and conducted at conditions far below those expected for the outer core. In this study, we adopt a synergistic approach by integrating experiments at experimentally-achievable conditions with computations up to core conditions. We performed viscosity measurements based on the modified Stokes' floating sphere viscometry method for the Fe-Ni-C liquids at high pressures in a Paris-Edinburgh press at Sector 16 of the Advanced Photon Source, Argonne National Laboratory. Our results show that the addition of 3-5 wt.% carbon to iron-nickel liquids has negligible effect on its viscosity at pressures lower than 5 GPa. The viscosity of the Fe-Ni-C liquids, however, becomes notably higher and increases by a factor of 3 at 5-8 GPa. Similarly, our first-principles molecular dynamics calculations up to Earth's core pressures show a viscosity change in Fe-Ni-C liquids at 5 GPa. The significant change in the viscosity is likely due to a liquid structural transition of the Fe-Ni-C liquids as revealed by our X-ray diffraction measurements and first-principles molecular dynamics calculations. The observed correlation between structure and physical properties of liquids permit stringent benchmark test of the computational liquid models and contribute to a more comprehensive understanding of liquid properties under high pressures. The interplay between experiments and first-principles based modeling is shown to be a practical and effective methodology for studying liquid properties under outer core conditions that are difficult to reach with the current static high-pressure capabilities. The new viscosity data from experiments and computations would provide new insights into the internal dynamics of the outer core.

Deriving dynamical models from paleoclimatic records: application to glacial millennial-scale climate variability.

PubMed

Kwasniok, Frank; Lohmann, Gerrit

2009-12-01

A method for systematically deriving simple nonlinear dynamical models from ice-core data is proposed. It offers a tool to integrate models and theories with paleoclimatic data. The method is based on the unscented Kalman filter, a nonlinear extension of the conventional Kalman filter. Here, we adopt the abstract conceptual model of stochastically driven motion in a potential that allows for two distinctly different states. The parameters of the model-the shape of the potential and the noise level-are estimated from a North Greenland ice-core record. For the glacial period from 70 to 20 ky before present, a potential is derived that is asymmetric and almost degenerate. There is a deep well corresponding to a cold stadial state and a very shallow well corresponding to a warm interstadial state.

A 3-D Finite-Volume Non-hydrostatic Icosahedral Model (NIM)

NASA Astrophysics Data System (ADS)

Lee, Jin

2014-05-01

The Nonhydrostatic Icosahedral Model (NIM) formulates the latest numerical innovation of the three-dimensional finite-volume control volume on the quasi-uniform icosahedral grid suitable for ultra-high resolution simulations. NIM's modeling goal is to improve numerical accuracy for weather and climate simulations as well as to utilize the state-of-art computing architecture such as massive parallel CPUs and GPUs to deliver routine high-resolution forecasts in timely manner. NIM dynamic corel innovations include: * A local coordinate system remapped spherical surface to plane for numerical accuracy (Lee and MacDonald, 2009), * Grid points in a table-driven horizontal loop that allow any horizontal point sequence (A.E. MacDonald, et al., 2010), * Flux-Corrected Transport formulated on finite-volume operators to maintain conservative positive definite transport (J.-L, Lee, ET. Al., 2010), *Icosahedral grid optimization (Wang and Lee, 2011), * All differentials evaluated as three-dimensional finite-volume integrals around the control volume. The three-dimensional finite-volume solver in NIM is designed to improve pressure gradient calculation and orographic precipitation over complex terrain. NIM dynamical core has been successfully verified with various non-hydrostatic benchmark test cases such as internal gravity wave, and mountain waves in Dynamical Cores Model Inter-comparisons Projects (DCMIP). Physical parameterizations suitable for NWP are incorporated into NIM dynamical core and successfully tested with multimonth aqua-planet simulations. Recently, NIM has started real data simulations using GFS initial conditions. Results from the idealized tests as well as real-data simulations will be shown in the conference.

An Intercomparison of the Dynamical Cores of Global Atmospheric Circulation Models for Mars

NASA Technical Reports Server (NTRS)

Hollingsworth, Jeffery L.; Bridger, Alison F. C.; Haberle, Robert M.

1998-01-01

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. The focus of this JRI has been to evaluate the dynamical 'cores' of two global atmospheric circulation models for Mars that are in operation at the NASA Ames Research Center. The two global circulation models in use are fundamentally different: one uses spherical harmonics in its horizontal representation of field variables; the other uses finite differences on a uniform longitude-latitude grid. Several simulations have been conducted to assess how the dynamical processors of each of these circulation models perform using identical 'simple physics' parameterizations. A variety of climate statistics (e.g., time-mean flows and eddy fields) have been compared for realistic solstitial mean basic states. Results of this research have demonstrated that the two Mars circulation models with completely different spatial representations and discretizations produce rather similar circulation statistics for first-order meteorological fields, suggestive of a tendency for convergence of numerical solutions. Second and higher-order fields can, however, vary significantly between the two models.

An Intercomparison of the Dynamical Cores of Global Atmospheric Circulation Models for Mars

NASA Technical Reports Server (NTRS)

Hollingsworth, Jeffery L.; Bridger, Alison F. C.; Haberle, Robert M.

1998-01-01

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Cen- ter and San Jose State University, Department of Meteorology. The focus of this JRI has been to evaluate the dynamical "cores" of two global atmospheric circulation models for Mars that are in operation at the NASA Ames Research Center. ne two global circulation models in use are fundamentally different: one uses spherical harmonics in its horizontal representation of field variables; the other uses finite differences on a uniform longitude-latitude grid. Several simulations have been conducted to assess how the dynamical processors of each of these circulation models perform using identical "simple physics" parameterizations. A variety of climate statistics (e.g., time-mean flows and eddy fields) have been compared for realistic solstitial mean basic states. Results of this research have demonstrated that the two Mars circulation models with completely different spatial representations and discretizations produce rather similar circulation statistics for first-order meteorological fields, suggestive of a tendency for convergence of numerical solutions. Second and higher-order fields can, however, vary significantly between the two models.

Osteocyte viability and regulation of osteoblast function in a 3D trabecular bone explant under dynamic hydrostatic pressure.

PubMed

Takai, Erica; Mauck, Robert L; Hung, Clark T; Guo, X Edward

2004-09-01

A new trabecular bone explant model was used to examine osteocyte-osteoblast interactions under DHP loading. DHP loading enhanced osteocyte viability as well as osteoblast function measured by osteoid formation. However, live osteocytes were necessary for osteoblasts to form osteoids in response to DHP, which directly show osteoblast-osteocyte interactions in this in vitro culture. A trabecular bone explant model was characterized and used to examine the effect of osteocyte and osteoblast interactions and dynamic hydrostatic pressure (DHP) loading on osteocyte viability and osteoblast function in long-term culture. Trabecular bone cores obtained from metacarpals of calves were cleaned of bone marrow and trabecular surface cells and divided into six groups, (1) live cores + dynamic hydrostatic pressure (DHP), (2) live cores + sham, (3) live cores + osteoblast + DHP, (4) live cores + osteoblast + sham, (5) devitalized cores + osteoblast + DHP, and (6) devitalized cores + osteoblast + sham, with four culture durations (2, 8, 15, and 22 days; n = 4/group). Cores from groups 3-6 were seeded with osteoblasts, and cores from groups 5 and 6 were devitalized before seeding. Groups 1, 3, and 5 were subjected to daily DHP loading. Bone histomorphometry was performed to quantify osteocyte viability based on morphology and to assess osteoblast function based on osteoid surface per bone surface (Os/Bs). TUNEL staining was performed to evaluate the mode of osteocyte death under various conditions. A portion of osteocytes remained viable for the duration of culture. DHP loading significantly enhanced osteocyte viability up to day 8, whereas the presence of seeded osteoblasts significantly decreased osteocyte viability. Cores with live osteocytes showed higher Os/Bs compared with devitalized cores, which reached significant levels over a greater range of time-points when combined with DHP loading. DHP loading did not increase Os/Bs in the absence of live osteocytes. The percentage of apoptotic cells remained the same regardless of treatment or culture duration. Enhanced osteocyte viability with DHP suggests the necessity of mechanical stimulation for osteocyte survival in vitro. Furthermore, osteocytes play a critical role in the transmission of signals from DHP loading to modulate osteoblast function. This explant culture model may be used for mechanotransduction studies in long-term cultures.

Comparing dark matter models, modified Newtonian dynamics and modified gravity in accounting for galaxy rotation curves

NASA Astrophysics Data System (ADS)

Li, Xin; Tang, Li; Lin, Hai-Nan

2017-05-01

We compare six models (including the baryonic model, two dark matter models, two modified Newtonian dynamics models and one modified gravity model) in accounting for galaxy rotation curves. For the dark matter models, we assume NFW profile and core-modified profile for the dark halo, respectively. For the modified Newtonian dynamics models, we discuss Milgrom’s MOND theory with two different interpolation functions, the standard and the simple interpolation functions. For the modified gravity, we focus on Moffat’s MSTG theory. We fit these models to the observed rotation curves of 9 high-surface brightness and 9 low-surface brightness galaxies. We apply the Bayesian Information Criterion and the Akaike Information Criterion to test the goodness-of-fit of each model. It is found that none of the six models can fit all the galaxy rotation curves well. Two galaxies can be best fitted by the baryonic model without involving nonluminous dark matter. MOND can fit the largest number of galaxies, and only one galaxy can be best fitted by the MSTG model. Core-modified model fits about half the LSB galaxies well, but no HSB galaxies, while the NFW model fits only a small fraction of HSB galaxies but no LSB galaxies. This may imply that the oversimplified NFW and core-modified profiles cannot model the postulated dark matter haloes well. Supported by Fundamental Research Funds for the Central Universities (106112016CDJCR301206), National Natural Science Fund of China (11305181, 11547305 and 11603005), and Open Project Program of State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, China (Y5KF181CJ1)

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

NASA Technical Reports Server (NTRS)

Dugan, E. T.; Kahook, S. D.; Diaz, N. J.

1996-01-01

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Statistical-dynamical modeling of the cloud-to-ground lightning activity in Portugal

NASA Astrophysics Data System (ADS)

Sousa, J. F.; Fragoso, M.; Mendes, S.; Corte-Real, J.; Santos, J. A.

2013-10-01

The present study employs a dataset of cloud-to-ground discharges over Portugal, collected by the Portuguese lightning detection network in the period of 2003-2009, to identify dynamically coherent lightning regimes in Portugal and to implement a statistical-dynamical modeling of the daily discharges over the country. For this purpose, the high-resolution MERRA reanalysis is used. Three lightning regimes are then identified for Portugal: WREG, WREM and SREG. WREG is a typical cold-core cut-off low. WREM is connected to strong frontal systems driven by remote low pressure systems at higher latitudes over the North Atlantic. SREG is a combination of an inverted trough and a mid-tropospheric cold-core nearby Portugal. The statistical-dynamical modeling is based on logistic regressions (statistical component) developed for each regime separately (dynamical component). It is shown that the strength of the lightning activity (either strong or weak) for each regime is consistently modeled by a set of suitable dynamical predictors (65-70% of efficiency). The difference of the equivalent potential temperature in the 700-500 hPa layer is the best predictor for the three regimes, while the best 4-layer lifted index is still important for all regimes, but with much weaker significance. Six other predictors are more suitable for a specific regime. For the purpose of validating the modeling approach, a regional-scale climate model simulation is carried out under a very intense lightning episode.

Adsorption of hairy particles with mobile ligands: Molecular dynamics and density functional study

NASA Astrophysics Data System (ADS)

Borówko, M.; Sokołowski, S.; Staszewski, T.; Pizio, O.

2018-01-01

We study models of hairy nanoparticles in contact with a hard wall. Each particle is built of a spherical core with a number of ligands attached to it and each ligand is composed of several spherical, tangentially jointed segments. The number of segments is the same for all ligands. Particular models differ by the numbers of ligands and of segments per ligand, but the total number of segments is constant. Moreover, our model assumes that the ligands are tethered to the core in such a manner that they can "slide" over the core surface. Using molecular dynamics simulations we investigate the differences in the structure of a system close to the wall. In order to characterize the distribution of the ligands around the core, we have calculated the end-to-end distances of the ligands and the lengths and orientation of the mass dipoles. Additionally, we also employed a density functional approach to obtain the density profiles. We have found that if the number of ligands is not too high, the proposed version of the theory is capable to predict the structure of the system with a reasonable accuracy.

Simulating the Dynamics of Earth's Core: Using NCCS Supercomputers Speeds Calculations

NASA Technical Reports Server (NTRS)

2002-01-01

If one wanted to study Earth's core directly, one would have to drill through about 1,800 miles of solid rock to reach liquid core-keeping the tunnel from collapsing under pressures that are more than 1 million atmospheres and then sink an instrument package to the bottom that could operate at 8,000 F with 10,000 tons of force crushing every square inch of its surface. Even then, several of these tunnels would probably be needed to obtain enough data. Faced with difficult or impossible tasks such as these, scientists use other available sources of information - such as seismology, mineralogy, geomagnetism, geodesy, and, above all, physical principles - to derive a model of the core and, study it by running computer simulations. One NASA researcher is doing just that on NCCS computers. Physicist and applied mathematician Weijia Kuang, of the Space Geodesy Branch, and his collaborators at Goddard have what he calls the,"second - ever" working, usable, self-consistent, fully dynamic, three-dimensional geodynamic model (see "The Geodynamic Theory"). Kuang runs his model simulations on the supercomputers at the NCCS. He and Jeremy Bloxham, of Harvard University, developed the original version, written in Fortran 77, in 1996.

Frenkel pair recombinations in UO2: Importance of explicit description of polarizability in core-shell molecular dynamics simulations

NASA Astrophysics Data System (ADS)

Devynck, Fabien; Iannuzzi, Marcella; Krack, Matthias

2012-05-01

The oxygen and uranium Frenkel pair (FP) recombination mechanisms are studied in UO2 using an empirical interatomic potential accounting for the polarizability of the ions, namely a dynamical core-shell model. The results are compared to a more conventional rigid-ion model. Both model types have been implemented into the cp2k program package and thoroughly validated. The overall picture indicates that the FP recombination mechanism is a complex process involving several phenomena. The FP recombination can happen instantaneously when the distance between the interstitial and the vacancy is small or can be thermally activated at larger separation distances. However, other criteria can prevail over the interstitial-vacancy distance. The surrounding environment of the FP defect, the mechanical stiffness of the matrix, and the orientation of the migration path are shown to be major factors acting on the FP lifetime. The core-shell and rigid-ion models provide a similar qualitative description of the FP recombination mechanism. However, the FP stabilities determined by both models significantly differ in the lower temperature range considered. Indeed, the recombination time of the oxygen and uranium FPs can be up to an order of magnitude lower in the core-shell model at T=600 K and T=1800 K, respectively. These differences highlight the importance of the explicit description of polarizability on some crucial properties such as the resistance to amorphization. This refined description of the interatomic interactions would certainly affect the description of the recrystallization process following a displacement cascade. In turn, the self-healing phase would be better accounted for in the core-shell model and the misestimate inherent to the lack of polarizability in the rigid-ion model corrected.

Ultrafast Three-Dimensional Integrated Imaging of Strain in Core/Shell Semiconductor/Metal Nanostructures

DOE PAGES

Cherukara, Mathew J.; Sasikumar, Kiran; DiChiara, Anthony; ...

2017-11-07

Visualizing the dynamical response of material heterointerfaces is increasingly important for the design of hybrid materials and structures with tailored properties for use in functional devices. In situ characterization of nanoscale heterointerfaces such as metal-semiconductor interfaces, which exhibit a complex interplay between lattice strain, electric potential, and heat transport at subnanosecond time scales, is particularly challenging. Here in this work, we use a laser pump/X-ray probe form of Bragg coherent diffraction imaging (BCDI) to visualize in three-dimension the deformation of the core of a model core/shell semiconductor-metal (ZnO/Ni) nanorod following laser heating of the shell. We observe a rich interplaymore » of radial, axial, and shear deformation modes acting at different time scales that are induced by the strain from the Ni shell. We construct experimentally informed models by directly importing the reconstructed crystal from the ultrafast experiment into a thermo-electromechanical continuum model. The model elucidates the origin of the deformation modes observed experimentally. Our integrated imaging approach represents an invaluable tool to probe strain dynamics across mixed interfaces under operando conditions.« less

Ultrafast Three-Dimensional Integrated Imaging of Strain in Core/Shell Semiconductor/Metal Nanostructures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cherukara, Mathew J.; Sasikumar, Kiran; DiChiara, Anthony

Visualizing the dynamical response of material heterointerfaces is increasingly important for the design of hybrid materials and structures with tailored properties for use in functional devices. In situ characterization of nanoscale heterointerfaces such as metal-semiconductor interfaces, which exhibit a complex interplay between lattice strain, electric potential, and heat transport at subnanosecond time scales, is particularly challenging. Here in this work, we use a laser pump/X-ray probe form of Bragg coherent diffraction imaging (BCDI) to visualize in three-dimension the deformation of the core of a model core/shell semiconductor-metal (ZnO/Ni) nanorod following laser heating of the shell. We observe a rich interplaymore » of radial, axial, and shear deformation modes acting at different time scales that are induced by the strain from the Ni shell. We construct experimentally informed models by directly importing the reconstructed crystal from the ultrafast experiment into a thermo-electromechanical continuum model. The model elucidates the origin of the deformation modes observed experimentally. Our integrated imaging approach represents an invaluable tool to probe strain dynamics across mixed interfaces under operando conditions.« less

Ultrafast Three-Dimensional Integrated Imaging of Strain in Core/Shell Semiconductor/Metal Nanostructures.

PubMed

Cherukara, Mathew J; Sasikumar, Kiran; DiChiara, Anthony; Leake, Steven J; Cha, Wonsuk; Dufresne, Eric M; Peterka, Tom; McNulty, Ian; Walko, Donald A; Wen, Haidan; Sankaranarayanan, Subramanian K R S; Harder, Ross J

2017-12-13

Visualizing the dynamical response of material heterointerfaces is increasingly important for the design of hybrid materials and structures with tailored properties for use in functional devices. In situ characterization of nanoscale heterointerfaces such as metal-semiconductor interfaces, which exhibit a complex interplay between lattice strain, electric potential, and heat transport at subnanosecond time scales, is particularly challenging. In this work, we use a laser pump/X-ray probe form of Bragg coherent diffraction imaging (BCDI) to visualize in three-dimension the deformation of the core of a model core/shell semiconductor-metal (ZnO/Ni) nanorod following laser heating of the shell. We observe a rich interplay of radial, axial, and shear deformation modes acting at different time scales that are induced by the strain from the Ni shell. We construct experimentally informed models by directly importing the reconstructed crystal from the ultrafast experiment into a thermo-electromechanical continuum model. The model elucidates the origin of the deformation modes observed experimentally. Our integrated imaging approach represents an invaluable tool to probe strain dynamics across mixed interfaces under operando conditions.

A Common Core for Active Conceptual Modeling for Learning from Surprises

NASA Astrophysics Data System (ADS)

Liddle, Stephen W.; Embley, David W.

The new field of active conceptual modeling for learning from surprises (ACM-L) may be helpful in preserving life, protecting property, and improving quality of life. The conceptual modeling community has developed sound theory and practices for conceptual modeling that, if properly applied, could help analysts model and predict more accurately. In particular, we need to associate more semantics with links, and we need fully reified high-level objects and relationships that have a clear, formal underlying semantics that follows a natural, ontological approach. We also need to capture more dynamic aspects in our conceptual models to more accurately model complex, dynamic systems. These concepts already exist, and the theory is well developed; what remains is to link them with the ideas needed to predict system evolution, thus enabling risk assessment and response planning. No single researcher or research group will be able to achieve this ambitious vision alone. As a starting point, we recommend that the nascent ACM-L community agree on a common core model that supports all aspects—static and dynamic—needed for active conceptual modeling in support of learning from surprises. A common core will more likely gain the traction needed to sustain the extended ACM-L research effort that will yield the advertised benefits of learning from surprises.

Competitive repetition suppression (CoRe) clustering: a biologically inspired learning model with application to robust clustering.

PubMed

Bacciu, Davide; Starita, Antonina

2008-11-01

Determining a compact neural coding for a set of input stimuli is an issue that encompasses several biological memory mechanisms as well as various artificial neural network models. In particular, establishing the optimal network structure is still an open problem when dealing with unsupervised learning models. In this paper, we introduce a novel learning algorithm, named competitive repetition-suppression (CoRe) learning, inspired by a cortical memory mechanism called repetition suppression (RS). We show how such a mechanism is used, at various levels of the cerebral cortex, to generate compact neural representations of the visual stimuli. From the general CoRe learning model, we derive a clustering algorithm, named CoRe clustering, that can automatically estimate the unknown cluster number from the data without using a priori information concerning the input distribution. We illustrate how CoRe clustering, besides its biological plausibility, posses strong theoretical properties in terms of robustness to noise and outliers, and we provide an error function describing CoRe learning dynamics. Such a description is used to analyze CoRe relationships with the state-of-the art clustering models and to highlight CoRe similitude with rival penalized competitive learning (RPCL), showing how CoRe extends such a model by strengthening the rival penalization estimation by means of loss functions from robust statistics.

An investigation into force-moment calibration techniques applicable to a magnetic suspension and balance system. M.S. Thesis

NASA Technical Reports Server (NTRS)

Eskins, Jonathan

1988-01-01

The problem of determining the forces and moments acting on a wind tunnel model suspended in a Magnetic Suspension and Balance System is addressed. Two calibration methods were investigated for three types of model cores, i.e., Alnico, Samarium-Cobalt, and a superconducting solenoid. Both methods involve calibrating the currents in the electromagnetic array against known forces and moments. The first is a static calibration method using calibration weights and a system of pulleys. The other method, dynamic calibration, involves oscillating the model and using its inertia to provide calibration forces and moments. Static calibration data, found to produce the most reliable results, is presented for three degrees of freedom at 0, 15, and -10 deg angle of attack. Theoretical calculations are hampered by the inability to represent iron-cored electromagnets. Dynamic calibrations, despite being quicker and easier to perform, are not as accurate as static calibrations. Data for dynamic calibrations at 0 and 15 deg is compared with the relevant static data acquired. Distortion of oscillation traces is cited as a major source of error in dynamic calibrations.

Predicting Student Academic Performance in an Engineering Dynamics Course: A Comparison of Four Types of Predictive Mathematical Models

ERIC Educational Resources Information Center

Huang, Shaobo; Fang, Ning

2013-01-01

Predicting student academic performance has long been an important research topic in many academic disciplines. The present study is the first study that develops and compares four types of mathematical models to predict student academic performance in engineering dynamics--a high-enrollment, high-impact, and core course that many engineering…

MODEL-Based Methodology for System of Systems Architecture Development with Application to the Recapitalization of the Future Towing and Salvage Platform

DTIC Science & Technology

2008-09-01

SEP) is a comprehensive , iterative and recursive problem solving process, applied sequentially top-down by integrated teams. It transforms needs...central integrated design repository. It includes a comprehensive behavior modeling notation to understand the dynamics of a design. CORE is a MBSE...37 F. DYNAMIC POSITIONING..........................................................................38 G. FIREFIGHTING

A Petascale Non-Hydrostatic Atmospheric Dynamical Core in the HOMME Framework

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tufo, Henry

The High-Order Method Modeling Environment (HOMME) is a framework for building scalable, conserva- tive atmospheric models for climate simulation and general atmospheric-modeling applications. Its spatial discretizations are based on Spectral-Element (SE) and Discontinuous Galerkin (DG) methods. These are local methods employing high-order accurate spectral basis-functions that have been shown to perform well on massively parallel supercomputers at any resolution and scale particularly well at high resolutions. HOMME provides the framework upon which the CAM-SE community atmosphere model dynamical-core is constructed. In its current incarnation, CAM-SE employs the hydrostatic primitive-equations (PE) of motion, which limits its resolution to simulations coarser thanmore » 0.1 per grid cell. The primary objective of this project is to remove this resolution limitation by providing HOMME with the capabilities needed to build nonhydrostatic models that solve the compressible Euler/Navier-Stokes equations.« less

Dynamical and Thermodynamic Elements of Modeled Climate Change at the East African Margin of Convection

NASA Astrophysics Data System (ADS)

Giannini, Alessandra; Lyon, Bradfield; Seager, Richard; Vigaud, Nicolas

2018-01-01

We propose a dynamical interpretation of model projections for an end-of-century wetting in equatorial East Africa. In the current generation of global climate models, increased atmospheric moisture content associated with warming is not the dominant process explaining the increase in rainfall, as the regional circulation is only weakly convergent even during the rainy seasons. Instead, projected wetter future conditions are generally consistent with the El Niño-like trend in tropical Pacific sea surface temperatures in climate models. In addition, a weakening in moisture convergence over the adjacent Congo Basin and Maritime Continent cores of convection results in the weakening of near-surface winds, which increases moisture advection from the Congo Basin core toward the East African margin. Overall confidence in the projections is limited by the significant biases in simulation of the regional climatology and disagreement between observed and modeled tropical Pacific sea surface temperature trends to date.

Low-Velocity Impact Response of Sandwich Beams with Functionally Graded Core

NASA Technical Reports Server (NTRS)

Apetre, N. A.; Sankar, B. V.; Ambur, D. R.

2006-01-01

The problem of low-speed impact of a one-dimensional sandwich panel by a rigid cylindrical projectile is considered. The core of the sandwich panel is functionally graded such that the density, and hence its stiffness, vary through the thickness. The problem is a combination of static contact problem and dynamic response of the sandwich panel obtained via a simple nonlinear spring-mass model (quasi-static approximation). The variation of core Young's modulus is represented by a polynomial in the thickness coordinate, but the Poisson's ratio is kept constant. The two-dimensional elasticity equations for the plane sandwich structure are solved using a combination of Fourier series and Galerkin method. The contact problem is solved using the assumed contact stress distribution method. For the impact problem we used a simple dynamic model based on quasi-static behavior of the panel - the sandwich beam was modeled as a combination of two springs, a linear spring to account for the global deflection and a nonlinear spring to represent the local indentation effects. Results indicate that the contact stiffness of thc beam with graded core Increases causing the contact stresses and other stress components in the vicinity of contact to increase. However, the values of maximum strains corresponding to the maximum impact load arc reduced considerably due to grading of thc core properties. For a better comparison, the thickness of the functionally graded cores was chosen such that the flexural stiffness was equal to that of a beam with homogeneous core. The results indicate that functionally graded cores can be used effectively to mitigate or completely prevent impact damage in sandwich composites.

Human Pyruvate Dehydrogenase Complex E2 and E3BP Core Subunits: New Models and Insights from Molecular Dynamics Simulations.

PubMed

Hezaveh, Samira; Zeng, An-Ping; Jandt, Uwe

2016-05-19

Targeted manipulation and exploitation of beneficial properties of multienzyme complexes, especially for the design of novel and efficiently structured enzymatic reaction cascades, require a solid model understanding of mechanistic principles governing the structure and functionality of the complexes. This type of system-level and quantitative knowledge has been very scarce thus far. We utilize the human pyruvate dehydrogenase complex (hPDC) as a versatile template to conduct corresponding studies. Here we present new homology models of the core subunits of the hPDC, namely E2 and E3BP, as the first time effort to elucidate the assembly of hPDC core based on molecular dynamic simulation. New models of E2 and E3BP were generated and validated at atomistic level for different properties of the proteins. The results of the wild type dimer simulations showed a strong hydrophobic interaction between the C-terminal and the hydrophobic pocket which is the main driving force in the intertrimer binding and the core self-assembly. On the contrary, the C-terminal truncated versions exhibited a drastic loss of hydrophobic interaction leading to a dimeric separation. This study represents a significant step toward a model-based understanding of structure and function of large multienzyme systems like PDC for developing highly efficient biocatalyst or bioreaction cascades.

Parameters of oscillation generation regions in open star cluster models

NASA Astrophysics Data System (ADS)

Danilov, V. M.; Putkov, S. I.

2017-07-01

We determine the masses and radii of central regions of open star cluster (OCL) models with small or zero entropy production and estimate the masses of oscillation generation regions in clustermodels based on the data of the phase-space coordinates of stars. The radii of such regions are close to the core radii of the OCL models. We develop a new method for estimating the total OCL masses based on the cluster core mass, the cluster and cluster core radii, and radial distribution of stars. This method yields estimates of dynamical masses of Pleiades, Praesepe, and M67, which agree well with the estimates of the total masses of the corresponding clusters based on proper motions and spectroscopic data for cluster stars.We construct the spectra and dispersion curves of the oscillations of the field of azimuthal velocities v φ in OCL models. Weak, low-amplitude unstable oscillations of v φ develop in cluster models near the cluster core boundary, and weak damped oscillations of v φ often develop at frequencies close to the frequencies of more powerful oscillations, which may reduce the non-stationarity degree in OCL models. We determine the number and parameters of such oscillations near the cores boundaries of cluster models. Such oscillations points to the possible role that gradient instability near the core of cluster models plays in the decrease of the mass of the oscillation generation regions and production of entropy in the cores of OCL models with massive extended cores.

Implications of Stellar Feedback for Dynamical Modeling of the Milky Way and Dwarf Galaxies

NASA Astrophysics Data System (ADS)

Wetzel, Andrew

2018-04-01

I will present recent results on dynamical modeling of stellar populations from the FIRE cosmological zoom-in baryonic simulations of Milky Way-like and dwarf galaxies. First, I will discuss the dynamical formation of the Milky Way, including the origin of thin+thick stellar disk morphology. I also will discuss the curious origin of metal-rich stars on halo-like orbits near the Sun, as recently measured by Gaia, with new insights from FIRE simulations on stellar radial migration/heating. Next, I will discuss role of stellar feedback in generating non-equilibrium fluctuations of the gravitational potential in low-mass 'dwarf' galaxies, which can explain the origin of cores in their dark-matter density profiles. In particular, we predict significant observable effects on stellar dynamics, including radial migration, size fluctuations, and population gradients, which can provide observational tests of feedback-driven core formation. Finally, this scenario can explain the formation of newly discovered 'ultra-diffuse' galaxies.

Normal Mode Derived Models of the Physical Properties of Earth's Outer Core

NASA Astrophysics Data System (ADS)

Irving, J. C. E.; Cottaar, S.; Lekic, V.; Wu, W.

2017-12-01

Earth's outer core, the largest reservoir of metal in our planet, is comprised of an iron alloy of an uncertain composition. Its dynamical behaviour is responsible for the generation of Earth's magnetic field, with convection driven both by thermal and chemical buoyancy fluxes. Existing models of the seismic velocity and density of the outer core exhibit some variation, and there are only a small number of models which aim to represent the outer core's density.It is therefore important that we develop a better understanding of the physical properties of the outer core. Though most of the outer core is likely to be well mixed, it is possible that the uppermost outer core is stably stratified: it may be enriched in light elements released during the growth of the solid, iron enriched, inner core; by elements dissolved from the mantle into the outer core; or by exsolution of compounds previously dissolved in the liquid metal which will eventually be swept into the mantle. The stratified layer may host MAC or Rossby waves and it could impede communication between the chemically differentiated mantle and outer core, including screening out some of the geodynamo's signal. We use normal mode center frequencies to estimate the physical properties of the outer core in a Bayesian framework. We estimate the mineral physical parameters needed to best produce velocity and density models of the outer core which are consistent with the normal mode observations. We require that our models satisfy realistic physical constraints. We create models of the outer core with and without a distinct uppermost layer and assess the importance of this region.Our normal mode-derived models are compared with observations of body waves which travel through the outer core. In particular, we consider SmKS waves which are especially sensitive to the uppermost outer core and are therefore an important way to understand the robustness of our models.

Next Generation Community Based Unified Global Modeling System Development and Operational Implementation Strategies at NCEP

NASA Astrophysics Data System (ADS)

Tallapragada, V.

2017-12-01

NOAA's Next Generation Global Prediction System (NGGPS) has provided the unique opportunity to develop and implement a non-hydrostatic global model based on Geophysical Fluid Dynamics Laboratory (GFDL) Finite Volume Cubed Sphere (FV3) Dynamic Core at National Centers for Environmental Prediction (NCEP), making a leap-step advancement in seamless prediction capabilities across all spatial and temporal scales. Model development efforts are centralized with unified model development in the NOAA Environmental Modeling System (NEMS) infrastructure based on Earth System Modeling Framework (ESMF). A more sophisticated coupling among various earth system components is being enabled within NEMS following National Unified Operational Prediction Capability (NUOPC) standards. The eventual goal of unifying global and regional models will enable operational global models operating at convective resolving scales. Apart from the advanced non-hydrostatic dynamic core and coupling to various earth system components, advanced physics and data assimilation techniques are essential for improved forecast skill. NGGPS is spearheading ambitious physics and data assimilation strategies, concentrating on creation of a Common Community Physics Package (CCPP) and Joint Effort for Data Assimilation Integration (JEDI). Both initiatives are expected to be community developed, with emphasis on research transitioning to operations (R2O). The unified modeling system is being built to support the needs of both operations and research. Different layers of community partners are also established with specific roles/responsibilities for researchers, core development partners, trusted super-users, and operations. Stakeholders are engaged at all stages to help drive the direction of development, resources allocations and prioritization. This talk presents the current and future plans of unified model development at NCEP for weather, sub-seasonal, and seasonal climate prediction applications with special emphasis on implementation of NCEP FV3 Global Forecast System (GFS) and Global Ensemble Forecast System (GEFS) into operations by 2019.

Composition of the low seismic velocity E' layer at the top of Earth's core

NASA Astrophysics Data System (ADS)

Badro, J.; Brodholt, J. P.

2017-12-01

Evidence for a layer (E') at the top of the outer core has been available since the '90s and while different studies suggest slightly different velocity contrasts and thicknesses, the common observation is that the layer has lower velocities than the bulk outer core (PREM). Although there are no direct measurements on the density of this layer, dynamic stability requires it to be less dense than the bulk outer core under those same pressure and temperature conditions. Using ab initio simulations on Fe-Ni-S-C-O-Si liquids we constrain the origin and composition of the low-velocity layer E' at the top of Earth's outer core. We find that increasing the concentration of any light-element always increases velocity and so a low-velocity and low-density layer (for stability) cannot be made by simply increasing light element concentration. This rules out barodiffusion or upwards sedimentation of a light phase for its origin. However, exchanging elements can—depending on the elements exchanged—produce such a layer. We evaluate three possibilities. Firstly, crystallization of a light phase from a core containing more than one light element may make such a layer, but only if the crystalizing phase is very Fe-rich, which is at odds with available phase diagrams at CMB conditions. Secondly, the E' layer may result from incomplete mixing of an early Earth core with a late impactor, depending on the light element compositions of the impactor and Earth's core, but such a primordial stratification is neither supported by dynamical models of the core nor thermodynamic models of core merger after the giant impact. The last and most plausible scenario is core-mantle chemical interaction; using thermodynamic models for metal-silicate partitioning of silicon and oxygen at CMB conditions, we show that a reaction between the core and an FeO-rich basal magma ocean can enrich the core in oxygen while depleting it in silicon, in relative amounts that produce a light and slow layer consistent with seismological observations.

Atmospheric Moisture Budget and Spatial Resolution Dependence of Precipitation Extremes in Aquaplanet Simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Qing; Leung, Lai-Yung R.; Rauscher, Sara

This study investigates the resolution dependency of precipitation extremes in an aqua-planet framework. Strong resolution dependency of precipitation extremes is seen over both tropics and extra-tropics, and the magnitude of this dependency also varies with dynamical cores. Moisture budget analyses based on aqua-planet simulations with the Community Atmosphere Model (CAM) using the Model for Prediction Across Scales (MPAS) and High Order Method Modeling Environment (HOMME) dynamical cores but the same physics parameterizations suggest that during precipitation extremes moisture supply for surface precipitation is mainly derived from advective moisture convergence. The resolution dependency of precipitation extremes mainly originates from advective moisturemore » transport in the vertical direction. At most vertical levels over the tropics and in the lower atmosphere over the subtropics, the vertical eddy transport of mean moisture field dominates the contribution to precipitation extremes and its resolution dependency. Over the subtropics, the source of moisture, its associated energy, and the resolution dependency during extremes are dominated by eddy transport of eddies moisture at the mid- and upper-troposphere. With both MPAS and HOMME dynamical cores, the resolution dependency of the vertical advective moisture convergence is mainly explained by dynamical changes (related to vertical velocity or omega), although the vertical gradients of moisture act like averaging kernels to determine the sensitivity of the overall resolution dependency to the changes in omega at different vertical levels. The natural reduction of variability with coarser resolution, represented by areal data averaging (aggregation) effect, largely explains the resolution dependency in omega. The thermodynamic changes, which likely result from non-linear feedback in response to the large dynamical changes, are small compared to the overall changes in dynamics (omega). However, after excluding the data aggregation effect in omega, thermodynamic changes become relatively significant in offsetting the effect of dynamics leading to reduce differences between the simulated and aggregated results. Compared to MPAS, the simulated stronger vertical motion with HOMME also results in larger resolution dependency. Compared to the simulation at fine resolution, the vertical motion during extremes is insufficiently resolved/parameterized at the coarser resolution even after accounting for the natural reduction in variability with coarser resolution, and this is more distinct in the simulation with HOMME. To reduce uncertainties in simulated precipitation extremes, future development in cloud parameterizations must address their sensitivity to spatial resolution as well as dynamical cores.« less

A variable resolution nonhydrostatic global atmospheric semi-implicit semi-Lagrangian model

NASA Astrophysics Data System (ADS)

Pouliot, George Antoine

2000-10-01

The objective of this project is to develop a variable-resolution finite difference adiabatic global nonhydrostatic semi-implicit semi-Lagrangian (SISL) model based on the fully compressible nonhydrostatic atmospheric equations. To achieve this goal, a three-dimensional variable resolution dynamical core was developed and tested. The main characteristics of the dynamical core can be summarized as follows: Spherical coordinates were used in a global domain. A hydrostatic/nonhydrostatic switch was incorporated into the dynamical equations to use the fully compressible atmospheric equations. A generalized horizontal variable resolution grid was developed and incorporated into the model. For a variable resolution grid, in contrast to a uniform resolution grid, the order of accuracy of finite difference approximations is formally lost but remains close to the order of accuracy associated with the uniform resolution grid provided the grid stretching is not too significant. The SISL numerical scheme was implemented for the fully compressible set of equations. In addition, the generalized minimum residual (GMRES) method with restart and preconditioner was used to solve the three-dimensional elliptic equation derived from the discretized system of equations. The three-dimensional momentum equation was integrated in vector-form to incorporate the metric terms in the calculations of the trajectories. Using global re-analysis data for a specific test case, the model was compared to similar SISL models previously developed. Reasonable agreement between the model and the other independently developed models was obtained. The Held-Suarez test for dynamical cores was used for a long integration and the model was successfully integrated for up to 1200 days. Idealized topography was used to test the variable resolution component of the model. Nonhydrostatic effects were simulated at grid spacings of 400 meters with idealized topography and uniform flow. Using a high-resolution topographic data set and the variable resolution grid, sets of experiments with increasing resolution were performed over specific regions of interest. Using realistic initial conditions derived from re-analysis fields, nonhydrostatic effects were significant for grid spacings on the order of 0.1 degrees with orographic forcing. If the model code was adapted for use in a message passing interface (MPI) on a parallel supercomputer today, it was estimated that a global grid spacing of 0.1 degrees would be achievable for a global model. In this case, nonhydrostatic effects would be significant for most areas. A variable resolution grid in a global model provides a unified and flexible approach to many climate and numerical weather prediction problems. The ability to configure the model from very fine to very coarse resolutions allows for the simulation of atmospheric phenomena at different scales using the same code. We have developed a dynamical core illustrating the feasibility of using a variable resolution in a global model.

On the coupling of fluid dynamics and electromagnetism at the top of the earth's core

NASA Technical Reports Server (NTRS)

Benton, E. R.

1985-01-01

A kinematic approach to short-term geomagnetism has recently been based upon pre-Maxwell frozen-flux electromagnetism. A complete dynamic theory requires coupling fluid dynamics to electromagnetism. A geophysically plausible simplifying assumption for the vertical vorticity balance, namely that the vertical Lorentz torque is negligible, is introduced and its consequences are developed. The simplified coupled magnetohydrodynamic system is shown to conserve a variety of magnetic and vorticity flux integrals. These provide constraints on eligible models for the geomagnetic main field, its secular variation, and the horizontal fluid motions at the top of the core, and so permit a number of tests of the underlying assumptions.

Structural response of 1/20-scale models of the Clinch River Breeder Reactor to a simulated hypothetical core-disruptive accident

DOE Office of Scientific and Technical Information (OSTI.GOV)

Romander, C M; Cagliostro, D J

Five experiments were performed to help evaluate the structural integrity of the reactor vessel and head design and to verify code predictions. In the first experiment (SM 1), a detailed model of the head was loaded statically to determine its stiffness. In the remaining four experiments (SM 2 to SM 5), models of the vessel and head were loaded dynamically under a simulated 661 MW-s hypothetical core disruptive accident (HCDA). Models SM 2 to SM 4, each of increasing complexity, systematically showed the effects of upper internals structures, a thermal liner, core support platform, and torospherical bottom on vessel response.more » Model SM 5, identical to SM 4 but more heavily instrumented, demonstrated experimental reproducibility and provided more comprehensive data. The models consisted of a Ni 200 vessel and core barrel, a head with shielding and simulated component masses, and an upper internals structure (UIS).« less

Dynamic optimization of metabolic networks coupled with gene expression.

PubMed

Waldherr, Steffen; Oyarzún, Diego A; Bockmayr, Alexander

2015-01-21

The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition. An approximation by a timescale separation leads to a coupled model of quasi-steady state constraints on the metabolic reactions, and differential equations for the substrate concentrations and biomass composition. We propose a dynamic optimization approach to determine reaction fluxes for this model, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. Discretization of the optimization problems leads to a linear program that can be efficiently solved. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle. Copyright © 2014 Elsevier Ltd. All rights reserved.

Waterlike anomalies in a two-dimensional core-softened potential

NASA Astrophysics Data System (ADS)

Bordin, José Rafael; Barbosa, Marcia C.

2018-02-01

We investigate the structural, thermodynamic, and dynamic behavior of a two-dimensional (2D) core-corona system using Langevin dynamics simulations. The particles are modeled by employing a core-softened potential which exhibits waterlike anomalies in three dimensions. In previous studies in a quasi-2D system a new region in the pressure versus temperature phase diagram of structural anomalies was observed. Here we show that for the two-dimensional case two regions in the pressure versus temperature phase diagram with structural, density, and diffusion anomalies are observed. Our findings indicate that, while the anomalous region at lower densities is due the competition between the two length scales in the potential at higher densities, the anomalous region is related to the reentrance of the melting line.

Basin-forming impacts on Mars and the coupled thermal evolution of the interior

NASA Astrophysics Data System (ADS)

Arkani-Hamed, J.; Roberts, J. H.

2015-12-01

The youngest of the Noachian giant impact basins on Mars, are either weakly magnetized or completely demagnetized, indicating that a global magnetic field was not present and that a core dynamo was not operating at the time those basins formed. Shock heating from this sequence of basin-forming impacts modified the pattern of mantle convection. The heating produced by the eight largest impacts (Acidalia, Amazonis, Ares, Chryse, Daedalia, Hellas, Scopolus, and Utopia) penetrates below the core-mantle boundary (CMB). Here, we extend previous workon coupled thermal evolution into 3D, in order to accurately model the spatial relationship between impact basins. At the time of each impact we introduce a temperature perturbation resulting from shock heating into the core and mantle. Stratification of the core occurs very quickly compared to mantle dynamics, and we horizontally average the temperature in the core.We model mantle convection using the 3D finite element code CitcomS, and the thermal evolution of the core using a 1D parameterization.Each impact alters the pattern of mantle dynamics and a significant amount of impact melt is produced in the near surface. However, only the outermost part of the core is affected; the inner core temperature is still adiabatic. Immediately following the impact, the inner core may remain convective. The top of the core will cool by conduction into the deeper core faster than across the CMB, deepening the zone of stable stratification. Further core cooling results in formation of a convecting zone at the top of the core that propagates downwards as the thermal gradient becomes adiabatic at greater depths. Our goal is to obtain a better estimate of the time scale for restoration of post-impact core dynamo activity. Because the disappearance of the magnetic field exposes the early atmosphere to solar wind activity, constraining the history of the dynamo is critical for understanding climate evolution and habitability of the surface.

Haptic Paddle Enhancements and a Formal Assessment of Student Learning in System Dynamics

ERIC Educational Resources Information Center

Gorlewicz, Jenna L.; Kratchman, Louis B.; Webster, Robert J., III

2014-01-01

The haptic paddle is a force-feedback joystick used at several universities in teaching System Dynamics, a core mechanical engineering undergraduate course where students learn to model dynamic systems in several domains. A second goal of the haptic paddle is to increase the accessibility of robotics and haptics by providing a low-cost device for…

Live-cell imaging of actin dynamics reveals mechanisms of stereocilia length regulation in the inner ear

PubMed Central

Drummond, Meghan C.; Barzik, Melanie; Bird, Jonathan E.; Zhang, Duan-Sun; Lechene, Claude P.; Corey, David P.; Cunningham, Lisa L.; Friedman, Thomas B.

2015-01-01

The maintenance of sensory hair cell stereocilia is critical for lifelong hearing; however, mechanisms of structural homeostasis remain poorly understood. Conflicting models propose that stereocilia F-actin cores are either continually renewed every 24–48 h via a treadmill or are stable, exceptionally long-lived structures. Here to distinguish between these models, we perform an unbiased survey of stereocilia actin dynamics in more than 500 utricle hair cells. Live-imaging EGFP-β-actin or dendra2-β-actin reveal stable F-actin cores with turnover and elongation restricted to stereocilia tips. Fixed-cell microscopy of wild-type and mutant β-actin demonstrates that incorporation of actin monomers into filaments is required for localization to stereocilia tips. Multi-isotope imaging mass spectrometry and live imaging of single differentiating hair cells capture stereociliogenesis and explain uniform incorporation of 15N-labelled protein and EGFP-β-actin into nascent stereocilia. Collectively, our analyses support a model in which stereocilia actin cores are stable structures that incorporate new F-actin only at the distal tips. PMID:25898120

The Blue Straggler Star Population in NGC 1261: Evidence for a Post-core-collapse Bounce State

NASA Astrophysics Data System (ADS)

Simunovic, Mirko; Puzia, Thomas H.; Sills, Alison

2014-11-01

We present a multi-passband photometric study of the Blue Straggler Star (BSS) population in the Galactic globular cluster (GC) NGC 1261, using available space- and ground-based survey data. The inner BSS population is found to have two distinct sequences in the color-magnitude diagram (CMD), similar to double BSS sequences detected in other GCs. These well defined sequences are presumably linked to single short-lived events such as core collapse, which are expected to boost the formation of BSSs. In agreement with this, we find a BSS sequence in NGC 1261 which can be well reproduced individually by a theoretical model prediction of a 2 Gyr old population of stellar collision products, which are expected to form in the denser inner regions during short-lived core contraction phases. Additionally, we report the occurrence of a group of BSSs with unusually blue colors in the CMD, which are consistent with a corresponding model of a 200 Myr old population of stellar collision products. The properties of the NGC 1261 BSS populations, including their spatial distributions, suggest an advanced dynamical evolutionary state of the cluster, but the core of this GC does not show the classical signatures of core collapse. We argue that these apparent contradictions provide evidence for a post-core-collapse bounce state seen in dynamical simulations of old GCs.

Extension of the AMBER molecular dynamics software to Intel's Many Integrated Core (MIC) architecture

NASA Astrophysics Data System (ADS)

Needham, Perri J.; Bhuiyan, Ashraf; Walker, Ross C.

2016-04-01

We present an implementation of explicit solvent particle mesh Ewald (PME) classical molecular dynamics (MD) within the PMEMD molecular dynamics engine, that forms part of the AMBER v14 MD software package, that makes use of Intel Xeon Phi coprocessors by offloading portions of the PME direct summation and neighbor list build to the coprocessor. We refer to this implementation as pmemd MIC offload and in this paper present the technical details of the algorithm, including basic models for MPI and OpenMP configuration, and analyze the resultant performance. The algorithm provides the best performance improvement for large systems (>400,000 atoms), achieving a ∼35% performance improvement for satellite tobacco mosaic virus (1,067,095 atoms) when 2 Intel E5-2697 v2 processors (2 ×12 cores, 30M cache, 2.7 GHz) are coupled to an Intel Xeon Phi coprocessor (Model 7120P-1.238/1.333 GHz, 61 cores). The implementation utilizes a two-fold decomposition strategy: spatial decomposition using an MPI library and thread-based decomposition using OpenMP. We also present compiler optimization settings that improve the performance on Intel Xeon processors, while retaining simulation accuracy.

Molecular Dynamics Simulations of Star Polymeric Molecules with Diblock Arms, a Comparative Study.

PubMed

Swope, William C; Carr, Amber C; Parker, Amanda J; Sly, Joseph; Miller, Robert D; Rice, Julia E

2012-10-09

We have performed all atom explicit solvent molecular dynamics simulations of three different star polymeric systems in water, each star molecule consisting of 16 diblock copolymer arms bound to a small adamantane core. The arms of each system consist of an inner "hydrophobic" block (either polylactide, polyvalerolactone, or polyethylene) and an outer hydrophilic block (polyethylene oxide, PEO). These models exhibit unusual structure very close to the core (clearly an artifact of our model) but which we believe becomes "normal" or bulk-like at relatively short distances from this core. We report on a number of temperature-dependent thermodynamic (structural/energetic) properties as well as kinetic properties. Our observations suggest that under physiological conditions, the hydrophobic regions of these systems may be solid and glassy, with only rare and shallow penetration by water, and that a sharp boundary exists between the hydrophobic cores and either the PEO or water. The PEO in these models is seen to be fully water-solvated at low temperatures but tends to phase separate from water as the temperature is increased, reminiscent of a lower critical solution temperature exhibited by PEO-water mixtures. Water penetration concentration and depth is composition and temperature dependent with greater water penetration for the most ester-rich star polymer.

Effect of attractive interactions on the water-like anomalies of a core-softened model potential.

PubMed

Pant, Shashank; Gera, Tarun; Choudhury, Niharendu

2013-12-28

It is now well established that water-like anomalies can be reproduced by a spherically symmetric potential with two length scales, popularly known as core-softened potential. In the present study we aim to investigate the effect of attractive interactions among the particles in a model fluid interacting with core-softened potential on the existence and location of various water-like anomalies in the temperature-pressure plane. We employ extensive molecular dynamic simulations to study anomalous nature of various order parameters and properties under isothermal compression. Order map analyses have also been done for all the potentials. We observe that all the systems with varying depth of attractive wells show structural, dynamic, and thermodynamic anomalies. As many of the previous studies involving model water and a class of core softened potentials have concluded that the structural anomaly region encloses the diffusion anomaly region, which in turn, encloses the density anomaly region, the same pattern has also been observed in the present study for the systems with less depth of attractive well. For the systems with deeper attractive well, we observe that the diffusion anomaly region shifts toward higher densities and is not always enclosed by the structural anomaly region. Also, density anomaly region is not completely enclosed by diffusion anomaly region in this case.

Toward a mineral physics reference model for the Moon's core.

PubMed

Antonangeli, Daniele; Morard, Guillaume; Schmerr, Nicholas C; Komabayashi, Tetsuya; Krisch, Michael; Fiquet, Guillaume; Fei, Yingwei

2015-03-31

The physical properties of iron (Fe) at high pressure and high temperature are crucial for understanding the chemical composition, evolution, and dynamics of planetary interiors. Indeed, the inner structures of the telluric planets all share a similar layered nature: a central metallic core composed mostly of iron, surrounded by a silicate mantle, and a thin, chemically differentiated crust. To date, most studies of iron have focused on the hexagonal closed packed (hcp, or ε) phase, as ε-Fe is likely stable across the pressure and temperature conditions of Earth's core. However, at the more moderate pressures characteristic of the cores of smaller planetary bodies, such as the Moon, Mercury, or Mars, iron takes on a face-centered cubic (fcc, or γ) structure. Here we present compressional and shear wave sound velocity and density measurements of γ-Fe at high pressures and high temperatures, which are needed to develop accurate seismic models of planetary interiors. Our results indicate that the seismic velocities proposed for the Moon's inner core by a recent reanalysis of Apollo seismic data are well below those of γ-Fe. Our dataset thus provides strong constraints to seismic models of the lunar core and cores of small telluric planets. This allows us to propose a direct compositional and velocity model for the Moon's core.

Fundamental structures of dynamic social networks.

PubMed

Sekara, Vedran; Stopczynski, Arkadiusz; Lehmann, Sune

2016-09-06

Social systems are in a constant state of flux, with dynamics spanning from minute-by-minute changes to patterns present on the timescale of years. Accurate models of social dynamics are important for understanding the spreading of influence or diseases, formation of friendships, and the productivity of teams. Although there has been much progress on understanding complex networks over the past decade, little is known about the regularities governing the microdynamics of social networks. Here, we explore the dynamic social network of a densely-connected population of ∼1,000 individuals and their interactions in the network of real-world person-to-person proximity measured via Bluetooth, as well as their telecommunication networks, online social media contacts, geolocation, and demographic data. These high-resolution data allow us to observe social groups directly, rendering community detection unnecessary. Starting from 5-min time slices, we uncover dynamic social structures expressed on multiple timescales. On the hourly timescale, we find that gatherings are fluid, with members coming and going, but organized via a stable core of individuals. Each core represents a social context. Cores exhibit a pattern of recurring meetings across weeks and months, each with varying degrees of regularity. Taken together, these findings provide a powerful simplification of the social network, where cores represent fundamental structures expressed with strong temporal and spatial regularity. Using this framework, we explore the complex interplay between social and geospatial behavior, documenting how the formation of cores is preceded by coordination behavior in the communication networks and demonstrating that social behavior can be predicted with high precision.

Fundamental structures of dynamic social networks

PubMed Central

Sekara, Vedran; Stopczynski, Arkadiusz; Lehmann, Sune

2016-01-01

Social systems are in a constant state of flux, with dynamics spanning from minute-by-minute changes to patterns present on the timescale of years. Accurate models of social dynamics are important for understanding the spreading of influence or diseases, formation of friendships, and the productivity of teams. Although there has been much progress on understanding complex networks over the past decade, little is known about the regularities governing the microdynamics of social networks. Here, we explore the dynamic social network of a densely-connected population of ∼1,000 individuals and their interactions in the network of real-world person-to-person proximity measured via Bluetooth, as well as their telecommunication networks, online social media contacts, geolocation, and demographic data. These high-resolution data allow us to observe social groups directly, rendering community detection unnecessary. Starting from 5-min time slices, we uncover dynamic social structures expressed on multiple timescales. On the hourly timescale, we find that gatherings are fluid, with members coming and going, but organized via a stable core of individuals. Each core represents a social context. Cores exhibit a pattern of recurring meetings across weeks and months, each with varying degrees of regularity. Taken together, these findings provide a powerful simplification of the social network, where cores represent fundamental structures expressed with strong temporal and spatial regularity. Using this framework, we explore the complex interplay between social and geospatial behavior, documenting how the formation of cores is preceded by coordination behavior in the communication networks and demonstrating that social behavior can be predicted with high precision. PMID:27555584

Structural response of 1/20-scale models of the Clinch River Breeder Reactor to a simulated hypothetical core disruptive accident. Technical report 4

DOE Office of Scientific and Technical Information (OSTI.GOV)

Romander, C. M.; Cagliostro, D. J.

Five experiments were performed to help evaluate the structural integrity of the reactor vessel and head design and to verify code predictions. In the first experiment (SM 1), a detailed model of the head was loaded statically to determine its stiffness. In the remaining four experiments (SM 2 to SM 5), models of the vessel and head were loaded dynamically under a simulated 661 MW-sec hypothetical core disruptive accident (HCDA). Models SM 2 to SM 4, each of increasing complexity, systematically showed the effects of upper internals structures, a thermal liner, core support platform, and torospherical bottom on vessel response.more » Model SM 5, identical to SM 4 but more heavily instrumented, demonstrated experimental reproducibility and provided more comprehensive data. The models consisted of a Ni 200 vessel and core barrel, a head with shielding and simulated component masses, an upper internals structure (UIS), and, in the more complex models SM 4 and SM 5, a Ni 200 thermal liner and core support structure. Water simulated the liquid sodium coolant and a low-density explosive simulated the HCDA loads.« less

Coupled-Circulation-Chemistry Studies with the Finite-Volume CCM: Trace Gas Transport in the Tropopause Region

NASA Technical Reports Server (NTRS)

Pawson, Steven; Lin, Shian-Jiann; Rood, Richard B.; Nebuda, Sharon; Nielsen, J. Eric; Douglass, Anne R.

2000-01-01

A joint project between the Data Assimilation Office at NASA GSFC and NCAR involves linking the physical packages from the Community Climate Model (CCM) with the flux-form semi-Lagrangian dynamical core developed by Lin and Rood in the DAO. A further development of this model includes the implementation of a chemical package developed by Douglass and colleagues in the Atmospheric Chemistry and Dynamics Branch at NASA GSFC. Results from this coupled dynamics-radiation-chemistry model will be presented, focussing on trace gas transport in the tropopause region.

Dynamical effects on the core-mantle boundary from depth-dependent thermodynamical properties of the lower mantle

NASA Technical Reports Server (NTRS)

Zhang, Shuxia; Yuen, David A.

1988-01-01

A common assumption in modeling dynamical processes in the lower mantle is that both the thermal expansivity and thermal conductivity are reasonably constant. Recent work from seismic equation of state leads to substantially higher values for the thermal conductivity and much lower thermal expansivity values in the deep mantle. The dynamical consequences of incorporating depth-dependent thermodynamic properties on the thermal-mechanical state of the lower mantle are examined with the spherical-shell mean-field equations. It is found that the thermal structure of the seismically resolved anomalous zone at the base of the mantle is strongly influenced by these variable properties and, in particular, that the convective distortion of the core-mantle boundary (CMB) is reduced with the decreasing thermal expansivity. Such a reduction of the dynamically induced topography from pure thermal convection would suggest that some other dynamical mechanism must be operating at the CMB.

Optimizing meridional advection of the Advanced Research WRF (ARW) dynamics for Intel Xeon Phi coprocessor

NASA Astrophysics Data System (ADS)

Mielikainen, Jarno; Huang, Bormin; Huang, Allen H.-L.

2015-05-01

The most widely used community weather forecast and research model in the world is the Weather Research and Forecast (WRF) model. Two distinct varieties of WRF exist. The one we are interested is the Advanced Research WRF (ARW) is an experimental, advanced research version featuring very high resolution. The WRF Nonhydrostatic Mesoscale Model (WRF-NMM) has been designed for forecasting operations. WRF consists of dynamics code and several physics modules. The WRF-ARW core is based on an Eulerian solver for the fully compressible nonhydrostatic equations. In the paper, we optimize a meridional (north-south direction) advection subroutine for Intel Xeon Phi coprocessor. Advection is of the most time consuming routines in the ARW dynamics core. It advances the explicit perturbation horizontal momentum equations by adding in the large-timestep tendency along with the small timestep pressure gradient tendency. We will describe the challenges we met during the development of a high-speed dynamics code subroutine for MIC architecture. Furthermore, lessons learned from the code optimization process will be discussed. The results show that the optimizations improved performance of the original code on Xeon Phi 7120P by a factor of 1.2x.

The Rise of China in the International Trade Network: A Community Core Detection Approach

PubMed Central

Zhu, Zhen; Cerina, Federica; Chessa, Alessandro; Caldarelli, Guido; Riccaboni, Massimo

2014-01-01

Theory of complex networks proved successful in the description of a variety of complex systems ranging from biology to computer science and to economics and finance. Here we use network models to describe the evolution of a particular economic system, namely the International Trade Network (ITN). Previous studies often assume that globalization and regionalization in international trade are contradictory to each other. We re-examine the relationship between globalization and regionalization by viewing the international trade system as an interdependent complex network. We use the modularity optimization method to detect communities and community cores in the ITN during the years 1995–2011. We find rich dynamics over time both inter- and intra-communities. In particular, the Asia-Oceania community disappeared and reemerged over time along with a switch in leadership from Japan to China. We provide a multilevel description of the evolution of the network where the global dynamics (i.e., communities disappear or reemerge) and the regional dynamics (i.e., community core changes between community members) are related. Moreover, simulation results show that the global dynamics can be generated by a simple dynamic-edge-weight mechanism. PMID:25136895

The rise of China in the International Trade Network: a community core detection approach.

PubMed

Zhu, Zhen; Cerina, Federica; Chessa, Alessandro; Caldarelli, Guido; Riccaboni, Massimo

2014-01-01

Theory of complex networks proved successful in the description of a variety of complex systems ranging from biology to computer science and to economics and finance. Here we use network models to describe the evolution of a particular economic system, namely the International Trade Network (ITN). Previous studies often assume that globalization and regionalization in international trade are contradictory to each other. We re-examine the relationship between globalization and regionalization by viewing the international trade system as an interdependent complex network. We use the modularity optimization method to detect communities and community cores in the ITN during the years 1995-2011. We find rich dynamics over time both inter- and intra-communities. In particular, the Asia-Oceania community disappeared and reemerged over time along with a switch in leadership from Japan to China. We provide a multilevel description of the evolution of the network where the global dynamics (i.e., communities disappear or reemerge) and the regional dynamics (i.e., community core changes between community members) are related. Moreover, simulation results show that the global dynamics can be generated by a simple dynamic-edge-weight mechanism.

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

NASA Technical Reports Server (NTRS)

Watson, Leela R.

2007-01-01

This report describes the work done by the Applied Meteorology Unit (AMU) in assessing the success of different model configurations in predicting warm season convection over East-Central Florida. The Weather Research and Forecasting Environmental Modeling System (WRF EMS) software allows users to choose among two dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). There are also data assimilation analysis packages available for the initialization of the WRF model - the Local Analysis and Prediction System (LAPS) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS). Besides model core and initialization options, the WRF model can be run with one- or two-way nesting. Having a series of initialization options and WRF cores, as well as many options within each core, creates challenges for local forecasters, such as determining which configuration options are best to address specific forecast concerns. This project assessed three different model intializations available to determine which configuration best predicts warm season convective initiation in East-Central Florida. The project also examined the use of one- and two-way nesting in predicting warm season convection.

Direct measurement of thermal conductivity in solid iron at planetary core conditions.

PubMed

Konôpková, Zuzana; McWilliams, R Stewart; Gómez-Pérez, Natalia; Goncharov, Alexander F

2016-06-02

The conduction of heat through minerals and melts at extreme pressures and temperatures is of central importance to the evolution and dynamics of planets. In the cooling Earth's core, the thermal conductivity of iron alloys defines the adiabatic heat flux and therefore the thermal and compositional energy available to support the production of Earth's magnetic field via dynamo action. Attempts to describe thermal transport in Earth's core have been problematic, with predictions of high thermal conductivity at odds with traditional geophysical models and direct evidence for a primordial magnetic field in the rock record. Measurements of core heat transport are needed to resolve this difference. Here we present direct measurements of the thermal conductivity of solid iron at pressure and temperature conditions relevant to the cores of Mercury-sized to Earth-sized planets, using a dynamically laser-heated diamond-anvil cell. Our measurements place the thermal conductivity of Earth's core near the low end of previous estimates, at 18-44 watts per metre per kelvin. The result is in agreement with palaeomagnetic measurements indicating that Earth's geodynamo has persisted since the beginning of Earth's history, and allows for a solid inner core as old as the dynamo.

On the Possibilities of Predicting Geomagnetic Secular Variation with Geodynamo Modeling

NASA Technical Reports Server (NTRS)

Kuang, Wei-Jia; Tangborn, Andrew; Sabaka, Terrance

2004-01-01

We use our MoSST core dynamics model and geomagnetic field at the core-mantle boundary (CMB) continued downward from surface observations to investigate possibilities of geomagnetic data assimilation, so that model results and current geomagnetic observations can be used to predict geomagnetic secular variation in future. As the first attempt, we apply data insertion technique to examine evolution of the model solution that is modified by geomagnetic input. Our study demonstrate that, with a single data insertion, large-scale poloidal magnetic field obtained from subsequent numerical simulation evolves similarly to the observed geomagnetic variation, regardless of the initial choice of the model solution (so long it is a well developed numerical solution). The model solution diverges on the time scales on the order of 60 years, similar to the time scales of the torsional oscillations in the Earth's core. Our numerical test shows that geomagnetic data assimilation is promising with our MoSST model.

Three-cluster dynamics within an ab initio framework

DOE PAGES

Quaglioni, Sofia; Romero-Redondo, Carolina; Navratil, Petr

2013-09-26

In this study, we introduce a fully antisymmetrized treatment of three-cluster dynamics within the ab initio framework of the no-core shell model/resonating-group method. Energy-independent nonlocal interactions among the three nuclear fragments are obtained from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with bound-state boundary conditions by means of the hyperspherical-harmonic method on a Lagrange mesh. We discuss the formalism in detail and give algebraic expressions for systems of two single nucleons plus a nucleus. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we apply the method to amore » 4He+n+n description of 6He and compare the results to experiment and to a six-body diagonalization of the Hamiltonian performed within the harmonic-oscillator expansions of the no-core shell model. Differences between the two calculations provide a measure of core ( 4He) polarization effects.« less

Parallel Discrete Molecular Dynamics Simulation With Speculation and In-Order Commitment*†

PubMed Central

Khan, Md. Ashfaquzzaman; Herbordt, Martin C.

2011-01-01

Discrete molecular dynamics simulation (DMD) uses simplified and discretized models enabling simulations to advance by event rather than by timestep. DMD is an instance of discrete event simulation and so is difficult to scale: even in this multi-core era, all reported DMD codes are serial. In this paper we discuss the inherent difficulties of scaling DMD and present our method of parallelizing DMD through event-based decomposition. Our method is microarchitecture inspired: speculative processing of events exposes parallelism, while in-order commitment ensures correctness. We analyze the potential of this parallelization method for shared-memory multiprocessors. Achieving scalability required extensive experimentation with scheduling and synchronization methods to mitigate serialization. The speed-up achieved for a variety of system sizes and complexities is nearly 6× on an 8-core and over 9× on a 12-core processor. We present and verify analytical models that account for the achieved performance as a function of available concurrency and architectural limitations. PMID:21822327

Parallel Discrete Molecular Dynamics Simulation With Speculation and In-Order Commitment.

PubMed

Khan, Md Ashfaquzzaman; Herbordt, Martin C

2011-07-20

Discrete molecular dynamics simulation (DMD) uses simplified and discretized models enabling simulations to advance by event rather than by timestep. DMD is an instance of discrete event simulation and so is difficult to scale: even in this multi-core era, all reported DMD codes are serial. In this paper we discuss the inherent difficulties of scaling DMD and present our method of parallelizing DMD through event-based decomposition. Our method is microarchitecture inspired: speculative processing of events exposes parallelism, while in-order commitment ensures correctness. We analyze the potential of this parallelization method for shared-memory multiprocessors. Achieving scalability required extensive experimentation with scheduling and synchronization methods to mitigate serialization. The speed-up achieved for a variety of system sizes and complexities is nearly 6× on an 8-core and over 9× on a 12-core processor. We present and verify analytical models that account for the achieved performance as a function of available concurrency and architectural limitations.

Steady induction effects in geomagnetism. Part 1A: Steady motional induction of geomagnetic chaos

NASA Technical Reports Server (NTRS)

Voorhies, Coerte V.

1992-01-01

Geomagnetic effects of magnetic induction by hypothetically steady fluid motion and steady magnetic flux diffusion near the top of Earth's core are investigated using electromagnetic theory, simple magnetic earth models, and numerical experiments with geomagnetic field models. The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation indicated by broad-scale models of the observed geomagnetic field is examined and solved. In Part 1, the steady surficial core flow estimation problem is solved in the context of the source-free mantle/frozen-flux core model. In the first paper (IA), the theory underlying such estimates is reviewed and some consequences of various kinematic and dynamic flow hypotheses are derived. For a frozen-flux core, fluid downwelling is required to change the mean square normal magnetic flux density averaged over the core-mantle boundary. For surficially geostrophic flow, downwelling implies poleward flow. The solution of the forward steady motional induction problem at the surface of a frozen-flux core is derived and found to be a fine, easily visualized example of deterministic chaos. Geomagnetic effects of statistically steady core surface flow may well dominate secular variation over several decades. Indeed, effects of persistent, if not steady, surficially geostrophic core flow are described which may help explain certain features of the present broad-scale geomagnetic field and perhaps paleomagnetic secular variation.

Analysis of the anomalous mean-field like properties of Gaussian core model in terms of entropy

NASA Astrophysics Data System (ADS)

Nandi, Manoj Kumar; Maitra Bhattacharyya, Sarika

2018-01-01

Studies of the Gaussian core model (GCM) have shown that it behaves like a mean-field model and the properties are quite different from standard glass former. In this work, we investigate the entropies, namely, the excess entropy (Sex) and the configurational entropy (Sc) and their different components to address these anomalies. Our study corroborates most of the earlier observations and also sheds new light on the high and low temperature dynamics. We find that unlike in standard glass former where high temperature dynamics is dominated by two-body correlation and low temperature by many-body correlations, in the GCM both high and low temperature dynamics are dominated by many-body correlations. We also find that the many-body entropy which is usually positive at low temperatures and is associated with activated dynamics is negative in the GCM suggesting suppression of activation. Interestingly despite the suppression of activation, the Adam-Gibbs (AG) relation that describes activated dynamics holds in the GCM, thus suggesting a non-activated contribution in AG relation. We also find an overlap between the AG relation and mode coupling power law regime leading to a power law behavior of Sc. From our analysis of this power law behavior, we predict that in the GCM the high temperature dynamics will disappear at dynamical transition temperature and below that there will be a transition to the activated regime. Our study further reveals that the activated regime in the GCM is quite narrow.

Ndarts

NASA Technical Reports Server (NTRS)

Jain, Abhinandan

2011-01-01

Ndarts software provides algorithms for computing quantities associated with the dynamics of articulated, rigid-link, multibody systems. It is designed as a general-purpose dynamics library that can be used for the modeling of robotic platforms, space vehicles, molecular dynamics, and other such applications. The architecture and algorithms in Ndarts are based on the Spatial Operator Algebra (SOA) theory for computational multibody and robot dynamics developed at JPL. It uses minimal, internal coordinate models. The algorithms are low-order, recursive scatter/ gather algorithms. In comparison with the earlier Darts++ software, this version has a more general and cleaner design needed to support a larger class of computational dynamics needs. It includes a frames infrastructure, allows algorithms to operate on subgraphs of the system, and implements lazy and deferred computation for better efficiency. Dynamics modeling modules such as Ndarts are core building blocks of control and simulation software for space, robotic, mechanism, bio-molecular, and material systems modeling.

Melting-induced stratification above the Earth's inner core due to convective translation.

PubMed

Alboussière, Thierry; Deguen, Renaud; Melzani, Mickaël

2010-08-05

In addition to its global North-South anisotropy, there are two other enigmatic seismological observations related to the Earth's inner core: asymmetry between its eastern and western hemispheres and the presence of a layer of reduced seismic velocity at the base of the outer core. This 250-km-thick layer has been interpreted as a stably stratified region of reduced composition in light elements. Here we show that this layer can be generated by simultaneous crystallization and melting at the surface of the inner core, and that a translational mode of thermal convection in the inner core can produce enough melting and crystallization on each hemisphere respectively for the dense layer to develop. The dynamical model we propose introduces a clear asymmetry between a melting and a crystallizing hemisphere which forms a basis for also explaining the East-West asymmetry. The present translation rate is found to be typically 100 million years for the inner core to be entirely renewed, which is one to two orders of magnitude faster than the growth rate of the inner core's radius. The resulting strong asymmetry of buoyancy flux caused by light elements is anticipated to have an impact on the dynamics of the outer core and on the geodynamo.

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

NASA Technical Reports Server (NTRS)

Lin, Shian-Jiann

2003-01-01

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water dynamical system. The 2D horizontal-to-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split- explicit, with large-time-step for scalar transport, and small fractional time step for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for mapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with physical parameterizations and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave growth tests and long-term integrations using the Held-Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

APPLYING THE PATUXENT LANDSCAPE UNIT MODEL TO HUMAN DOMINATED ECOSYSTEMS: THE CASE OF AGRICULTURE. (R827169)

EPA Science Inventory

Non-spatial dynamics are core to landscape simulations. Unit models simulate system interactions aggregated within one space unit of resolution used within a spatial model. For unit models to be applicable to spatial simulations they have to be formulated in a general enough w...

A Dynamic Stimulus-Driven Model of Signal Detection

ERIC Educational Resources Information Center

Turner, Brandon M.; Van Zandt, Trisha; Brown, Scott

2011-01-01

Signal detection theory forms the core of many current models of cognition, including memory, choice, and categorization. However, the classic signal detection model presumes the a priori existence of fixed stimulus representations--usually Gaussian distributions--even when the observer has no experience with the task. Furthermore, the classic…

A Nonlinear Dynamic Inversion Predictor-Based Model Reference Adaptive Controller for a Generic Transport Model

NASA Technical Reports Server (NTRS)

Campbell, Stefan F.; Kaneshige, John T.

2010-01-01

Presented here is a Predictor-Based Model Reference Adaptive Control (PMRAC) architecture for a generic transport aircraft. At its core, this architecture features a three-axis, non-linear, dynamic-inversion controller. Command inputs for this baseline controller are provided by pilot roll-rate, pitch-rate, and sideslip commands. This paper will first thoroughly present the baseline controller followed by a description of the PMRAC adaptive augmentation to this control system. Results are presented via a full-scale, nonlinear simulation of NASA s Generic Transport Model (GTM).

Toward a more efficient and scalable checkpoint/restart mechanism in the Community Atmosphere Model

NASA Astrophysics Data System (ADS)

Anantharaj, Valentine

2015-04-01

The number of cores (both CPU as well as accelerator) in large-scale systems has been increasing rapidly over the past several years. In 2008, there were only 5 systems in the Top500 list that had over 100,000 total cores (including accelerator cores) whereas the number of system with such capability has jumped to 31 in Nov 2014. This growth however has also increased the risk of hardware failure rates, necessitating the implementation of fault tolerance mechanism in applications. The checkpoint and restart (C/R) approach is commonly used to save the state of the application and restart at a later time either after failure or to continue execution of experiments. The implementation of an efficient C/R mechanism will make it more affordable to output the necessary C/R files more frequently. The availability of larger systems (more nodes, memory and cores) has also facilitated the scaling of applications. Nowadays, it is more common to conduct coupled global climate simulation experiments at 1 deg horizontal resolution (atmosphere), often requiring about 103 cores. At the same time, a few climate modeling teams that have access to a dedicated cluster and/or large scale systems are involved in modeling experiments at 0.25 deg horizontal resolution (atmosphere) and 0.1 deg resolution for the ocean. These ultrascale configurations require the order of 104 to 105 cores. It is not only necessary for the numerical algorithms to scale efficiently but the input/output (IO) mechanism must also scale accordingly. An ongoing series of ultrascale climate simulations, using the Titan supercomputer at the Oak Ridge Leadership Computing Facility (ORNL), is based on the spectral element dynamical core of the Community Atmosphere Model (CAM-SE), which is a component of the Community Earth System Model and the DOE Accelerated Climate Model for Energy (ACME). The CAM-SE dynamical core for a 0.25 deg configuration has been shown to scale efficiently across 100,000 cpu cores. At this scale, there is an increased risk that the simulation could be terminated due to hardware failures, resulting in a loss that could be as high as 105 - 106 titan core hours. Increasing the frequency of the output of C/R files could mitigate this loss but at the cost of additional C/R overhead. We are testing a more efficient C/R mechanism in CAM-SE. Our early implementation has demonstrated a nearly 3X performance improvement for a 1 deg CAM-SE (with CAM5 physics and MOZART chemistry) configuration using nearly 103 cores. We are in the process of scaling our implementation to 105 cores. This would allow us to run ultra scale simulations with more sophisticated physics and chemistry options while making better utilization of resources.

New insights into the complex regulation of the glycolytic pathway in Lactococcus lactis. I. Construction and diagnosis of a comprehensive dynamic model.

PubMed

Dolatshahi, Sepideh; Fonseca, Luis L; Voit, Eberhard O

2016-01-01

This article and the companion paper use computational systems modeling to decipher the complex coordination of regulatory signals controlling the glycolytic pathway in the dairy bacterium Lactococcus lactis. In this first article, the development of a comprehensive kinetic dynamic model is described. The model is based on in vivo NMR data that consist of concentration trends in key glycolytic metabolites and cofactors. The model structure and parameter values are identified with a customized optimization strategy that uses as its core the method of dynamic flux estimation. For the first time, a dynamic model with a single parameter set fits all available glycolytic time course data under anaerobic operation. The model captures observations that had not been addressed so far and suggests the existence of regulatory effects that had been observed in other species, but not in L. lactis. The companion paper uses this model to analyze details of the dynamic control of glycolysis under aerobic and anaerobic conditions.

Toward a mineral physics reference model for the Moon’s core

PubMed Central

Antonangeli, Daniele; Morard, Guillaume; Schmerr, Nicholas C.; Komabayashi, Tetsuya; Krisch, Michael; Fiquet, Guillaume; Fei, Yingwei

2015-01-01

The physical properties of iron (Fe) at high pressure and high temperature are crucial for understanding the chemical composition, evolution, and dynamics of planetary interiors. Indeed, the inner structures of the telluric planets all share a similar layered nature: a central metallic core composed mostly of iron, surrounded by a silicate mantle, and a thin, chemically differentiated crust. To date, most studies of iron have focused on the hexagonal closed packed (hcp, or ε) phase, as ε-Fe is likely stable across the pressure and temperature conditions of Earth’s core. However, at the more moderate pressures characteristic of the cores of smaller planetary bodies, such as the Moon, Mercury, or Mars, iron takes on a face-centered cubic (fcc, or γ) structure. Here we present compressional and shear wave sound velocity and density measurements of γ-Fe at high pressures and high temperatures, which are needed to develop accurate seismic models of planetary interiors. Our results indicate that the seismic velocities proposed for the Moon’s inner core by a recent reanalysis of Apollo seismic data are well below those of γ-Fe. Our dataset thus provides strong constraints to seismic models of the lunar core and cores of small telluric planets. This allows us to propose a direct compositional and velocity model for the Moon’s core. PMID:25775531

Development of a severe local storm prediction system: A 60-day test of a mesoscale primitive equation model

NASA Technical Reports Server (NTRS)

Paine, D. A.; Zack, J. W.; Kaplan, M. L.

1979-01-01

The progress and problems associated with the dynamical forecast system which was developed to predict severe storms are examined. The meteorological problem of severe convective storm forecasting is reviewed. The cascade hypothesis which forms the theoretical core of the nested grid dynamical numerical modelling system is described. The dynamical and numerical structure of the model used during the 1978 test period is presented and a preliminary description of a proposed multigrid system for future experiments and tests is provided. Six cases from the spring of 1978 are discussed to illustrate the model's performance and its problems. Potential solutions to the problems are examined.

Representing and reasoning about program in situation calculus

NASA Astrophysics Data System (ADS)

Yang, Bo; Zhang, Ming-yi; Wu, Mao-nian; Xie, Gang

2011-12-01

Situation calculus is an expressive tool for modeling dynamical system in artificial intelligence, changes in a dynamical world is represented naturally by the notions of action, situation and fluent in situation calculus. Program can be viewed as a discrete dynamical system, so it is possible to model program with situation calculus. To model program written in a smaller core programming language CL, notion of fluent is expanded for representing value of expression. Together with some functions returning concerned objects from expressions, a basic action theory of CL programming is constructed. Under such a theory, some properties of program, such as correctness and termination can be reasoned about.

Porosity and Permeability Evolution in Cemented Rock Cores under Reactive Flowing Conditions: Comparative Analysis between Limestone and Sandstone Host Rocks

NASA Astrophysics Data System (ADS)

Cao, P.; Karpyn, Z.; Li, L.

2013-12-01

CO2-brine has the potential to alter wellbore cement in depleted oil and gas reservoirs under geological CO2 sequestration conditions. A better understanding of CO2-brine-cement-rock interaction is needed to evaluate the seal integrity of candidate sequestration formation in the long run. This work investigates possible alteration of wellbore cement when bonded by different host formation rock upon exposure to CO2-saturated brine. Composite cement-sandstone and cement-limestone core samples were created to perform reactive coreflood experiments. After an eight-day dynamic flow-through period, both cores had a similar extent of porosity increase, while the cement-limestone core experienced a ten-fold higher increase in permeability. With the aid of X-ray Micro-CT imaging and Scanning Electron Microscopy, it is observed that cement underwent greater degradation at the cement-sandstone interface. Degradation of cement-limestone core mainly took place on the host rock matrix. Worm holes were developed and a solution channel was formed in the limestone, creating a dominant flow path that altered both flow and reaction behavior. Limestone buffered the injected acidic brine preventing further deterioration of cement near the core outlet. Changes in fluid chemistry of limestone and sandstone coreflood effluents are compared. Results from this work are aimed at assisting the development and validation of robust reactive transport models through direct measurement of cemented rock core porosity and permeability evolution as well as the effluent aqueous chemistry change. This will subsequently improve predictive capabilities of reactive transport models associated with CO2 sequestration in geologic environments. Permeability Evolution of Cement-Rock Core Sample during Dynamic Flow of CO2-Brine

Hypervelocity Impact Behaviour of CFRP-A1/HC Sandwich Panel: Finite-Element Studies

NASA Astrophysics Data System (ADS)

Phadnis, Vaibhav A.; Roy, Anish; Silberschmidt, Vadim V.

2014-06-01

The mechanical response of CFRP-Al/HC (carbon fibre- reinforced/epoxy composite face sheets with Al honeycomb core) sandwich panels to hyper-velocity impact ( 1 km/s) is studied using a finite-element model developed in ABAQUS/Explicit. The intraply damage of CFRP face sheets is analysed by the means of a user-defined material model (VUMAT) employing a combination of Hashin and Puck criteria and delamination is modelled using cohesive-zone elements. The damage of Al/HC core is assessed on the basis of a Johnson-Cook dynamic failure model while its hydrodynamic response is captured using the Mie- Gruneisen equation of state. The results obtained with the developed finite-element model showed a reasonable correlation to experimental damage patterns. The surface peeling of both face sheets was evident, with a significant delamination around the impact location accompanied by crushing of HC core.

Finite element analysis of hypervelocity impact behaviour of CFRP-Al/HC sandwich panel

NASA Astrophysics Data System (ADS)

Phadnis, Vaibhav A.; Silberschmidt, Vadim V.

2015-09-01

The mechanical response of CFRP-Al/HC (carbon fibre-reinforced/epoxy composite face sheets with Al honeycomb core) sandwich panels to hyper-velocity impact (up to 1 km/s) is studied using a finite-element model developed in ABAQUS/Explicit. The intraply damage of CFRP face sheets is analysed by mean of a user-defined material model (VUMAT) employing a combination of Hashin and Puck criteria, delamination modelled using cohesive-zone elements. The damaged Al/HC core is assessed on the basis of a Johnson Cook dynamic failure model while its hydrodynamic response is captured using the Mie-Gruneisen equation of state. The results obtained with the developed finite-element model showed a reasonable correlation to experimental damage patterns. The surface peeling of both face sheets was evident, with a significant delamination around the impact location accompanied by crushing HC core.

Computational modeling of lava domes using particle dynamics to investigate the effect of conduit flow mechanics on flow patterns

NASA Astrophysics Data System (ADS)

Husain, Taha Murtuza

Large (1--4 x 106 m3) to major (> 4 x 106 m3) dome collapses for andesitic lava domes such as Soufriere Hills Volcano, Montserrat are observed for elevated magma discharge rates (6--13 m3/s). The gas rich magma pulses lead to pressure build up in the lava dome that result in structural failure of the over steepened canyon-like walls which may lead to rockfall or pyroclastic flow. This indicates that dome collapse intimately related to magma extrusion rate. Variation in magma extrusion rate for open-system magma chambers is observed to follow alternating periods of high and low activity. Periodic behavior of magma exhibits a rich diversity in the nature of its eruptive history due to variation in magma chamber size, total crystal content, linear crystal growth rate and magma replenishment rate. Distinguished patterns of growth were observed at different magma flow rates ranging from endogenous to exogenous dome growth for magma with varying strengths. Determining the key parameters that control the transition in flow pattern of the magma during its lava dome building eruption is the main focus. This dissertation examines the mechanical effects on the morphology of the evolving lava dome on the extrusion of magma from a central vent using a 2D particle dynamics model. The particle dynamics model is coupled with a conduit flow model that incorporates the kinetics of crystallization and rheological stiffening to investigate important mechanisms during lava dome building eruptions. Chapter I of this dissertation explores lava dome growth and failure mechanics using a two-dimensional particle-dynamics model. The model follows the evolution of fractured lava, with solidification driven by degassing induced crystallization of magma. The particle-dynamics model emulates the natural development of dome growth and rearrangement of the lava dome which is difficult in mesh-based analyses due to mesh entanglement effects. The deformable talus evolves naturally as a frictional carapace that caps a ductile magma core. Extrusion rate and magma rheology together with crystallization temperature and volatile content govern the distribution of strength in the composite structure. This new model is calibrated against existing observational models of lava dome growth. Chapter II of this dissertation explores the effects of a spectrum of different rheological regimes, on eruptive style and morphologic evolution of lava domes, using a two-dimensional (2D) particle-dynamics model for a spreading viscoplastic (Bingham) fluid. We assume that the ductile magma core of a 2-D synthetic lava dome develops finite yield strength, and that deformable frictional talus evolves from a carapace that caps the magma core. Our new model is calibrated against an existing analytical model for a spreading viscoplastic lava dome and is further compared against observational data of lava dome growth. Chapter III of this dissertation explores different lava-dome styles by developing a two-dimensional particle-dynamics model. These growth patterns range from endogenous lava dome growth comprising expansion of a ductile dome core to the exogenous extrusion of a degassed lava plug resulting in generation of a lava spine. We couple conduit flow dynamics with surface growth of the evolving lava dome, fueled by an open-system magma chamber undergoing continuous replenishment. The conduit flow model accounts for the variation in rheology of ascending magma that results from degassing-induced crystallization. Chapter IV of this dissertation explores the Variation in the extruding lava flow patterns range from endogenous dome growth with a ductile core to the exogenous extrusion of a degassed lava plug that results in the generation of a spine. The variations are a manifestation of the changes in the magma rheology which is governed by magma composition and rate of decompression of the ascending magma. We simulate using a two-dimensional particle-dynamics model, the cyclic behavior of lava dome growth with endogenous growth at high discharge rates followed by exogenous extrusion of rheologically stiffened lava due to degassing induced crystallization at low discharge rates. We couple conduit flow dynamics with surface growth of the evolving lava dome which is fueled by an overpressured reservoir undergoing constant replenishment. The periodic behavior between magma chamber pressure and discharge rate is reproduced as a result of the temporal and spatial change in magma viscosity controlled by crystallization kinetics. Dimensionless numbers are used to map the flow behaviors with the changing extrusion regime. A dimensionless plot identifying the flow transition region during the growth cycle of an evolving lava dome in its lava dome eruptive period is presented. The plot provides a the threshold value of a dimensionless strength parameter (pi 2 < 3.31 x 10-4) below which the transition in flow pattern occurs from endogenously evolving lava dome with a ductile core to the development of a shear lobe for short or long lived periodic episode of the extrusion of magma. (Abstract shortened by UMI.).

Investigating the source of anomalous PKP travel-times on south-Sandwich to Alaska paths

NASA Astrophysics Data System (ADS)

Frost, D. A.; Romanowicz, B. A.

2017-12-01

Inner core anisotropy was proposed thirty years ago to explain differences in travel times of inner core phases (PKIKP) on polar and equatorial paths (Morelli et al., 1986). Over time, models of inner core anisotropy have become very complex, with evidence for depth dependence, hemispherical variations, and other localised features. Some models propose the strength of anisotropy to be in excess of 4% in the western hemisphere of the inner core. This is difficult to reconcile with predictions from mineral physics and dynamical models of inner core growth. The strong anisotropy is confined to anomalous paths between earthquakes in the south Sandwich Islands and stations in Alaska. In contrast, the strength of inner core anisotropy obtained from measurements of PKPPKP travel times on polar paths does not exceed 1-2% (Bréger et al., 2000; Frost and Romanowicz, 2017). We re-examine the trends of PKIKP travel times on polar paths, in order to reconcile the different measurements and to determine whether discrepancies can be explained by structure in the mantle, the outer core, or localized strong anisotropy in the inner core. For this, we combine existing and new measurements, taking advantage of recent deployments of broadband arrays in Alaska and Antarctica.

Glacial-interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

NASA Astrophysics Data System (ADS)

Capron, E.; Landais, A.; Buiron, D.; Cauquoin, A.; Chappellaz, J.; Debret, M.; Jouzel, J.; Leuenberger, M.; Martinerie, P.; Masson-Delmotte, V.; Mulvaney, R.; Parrenin, F.; Prié, F.

2012-12-01

Correct estimate of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice cores studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: output of a firn densification model and measurements of δ15N of N2 in air trapped in ice core. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available δ15N measurements performed from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rate and temperature conditions. While firn densification simulations are able to correctly represent most of the δ15N trends over the last deglaciation measured in the EDC, BI, TALDICE and EDML ice cores, they systematically fail to capture BI and EDML δ15N glacial levels, a mismatch previously seen for Central East Antarctic ice cores. Using empirical constraints of the EDML gas-ice depth offset during the Laschamp event (~ 41 ka), we can rule out the existence of a large convective zone as the explanation of the glacial firn model-δ15N data mismatch for this site. The good match between modelled and measured δ15N at TALDICE as well as the lack of any clear correlation between insoluble dust concentration in snow and δ15N records in the different ice cores suggest that past changes in loads of impurities are not the only main driver of glacial-interglacial changes in firn lock-in depth. We conclude that firn densification dynamics may instead be driven mostly by accumulation rate changes. The mismatch between modelled and measured δ15N may be due to inaccurate reconstruction of past accumulation rate or underestimated influence of accumulation rate in firnification models.

Dynamic characteristics of a novel damped outrigger system

NASA Astrophysics Data System (ADS)

Tan, Ping; Fang, Chuangjie; Zhou, Fulin

2014-06-01

This paper presents exact analytical solutions for a novel damped outrigger system, in which viscous dampers are vertically installed between perimeter columns and the core of a high-rise building. An improved analytical model is developed by modeling the effect of the damped outrigger as a general rotational spring acting on a Bernoulli-Euler beam. The equivalent rotational spring stiffness incorporating the combined effects of dampers and axial stiffness of perimeter columns is derived. The dynamic stiffness method (DSM) is applied to formulate the governing equation of the damped outrigger system. The accuracy and efficiency are verified in comparison with those obtained from compatibility equations and boundary equations. Parametric analysis of three non-dimensional factors is conducted to evaluate the influences of various factors, such as the stiffness ratio of the core to the beam, position of the damped outrigger, and the installed damping coefficient. Results show that the modal damping ratio is significantly influenced by the stiffness ratio of the core to the column, and is more sensitive to damping than the position of the damped outrigger. The proposed analytical model in combination with DSM can be extended to the study of structures with more outriggers.

The dynamic two-fluid model OLGA; Theory and application

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bendiksen, K.H.; Maines, D.; Moe, R.

1991-05-01

Dynamic two-fluid models have found a wide range of application in the simulation of two-phase-flow systems, particularly for the analysis of steam/water flow in the core of a nuclear reactor. Until quite recently, however, very few attempts have been made to use such models in the simulation of two-phase oil and gas flow in pipelines. This paper presents a dynamic two-fluid model, OLGA, in detail, stressing the basic equations and the two-fluid models applied. Predictions of steady-state pressure drop, liquid hold-up, and flow-regime transitions are compared with data from the SINTEF Two-Phase Flow Laboratory and from the literature. Comparisons withmore » evaluated field data are also presented.« less

Decoherence in models for hard-core bosons coupled to optical phonons

NASA Astrophysics Data System (ADS)

Dey, A.; Lone, M. Q.; Yarlagadda, S.

2015-09-01

Understanding coherent dynamics of excitons, spins, or hard-core bosons (HCBs) has tremendous scientific and technological implications for quantum computation. Here, we study decay of excited-state population and decoherence in two models for HCBs, namely, a two-site HCB model with site-dependent strong potentials and subject to non-Markovian dynamics and an infinite-range HCB model governed by Markovian dynamics. Both models are investigated in the regimes of antiadiabaticity and strong HCB-phonon coupling with each site providing a different local optical phonon environment; furthermore, the HCB systems in both models are taken to be initially uncorrelated with the environment in the polaronic frame of reference. In the case of the two-site HCB model, we show clearly that the degree of decoherence and decay of excited state are enhanced by the proximity of the site-energy difference to the eigenenergy of phonons and are most pronounced when the site-energy difference is at resonance with twice the polaronic energy; additionally, the decoherence and the decay effects are reduced when the strength of HCB-phonon coupling is increased. For the infinite-range model, when the site energies are the same, we derive an effective many-body Hamiltonian that commutes with the long-range system Hamiltonian and thus has the same set of eigenstates; consequently, a quantum-master-equation approach shows that the quantum states of the system do not decohere.

Heat Pipe Reactor Dynamic Response Tests: SAFE-100 Reactor Core Prototype

NASA Technical Reports Server (NTRS)

Bragg-Sitton, Shannon M.

2005-01-01

The SAFE-I00a test article at the NASA Marshall Space Flight Center was used to simulate a variety of potential reactor transients; the SAFEl00a is a resistively heated, stainless-steel heat-pipe (HP)-reactor core segment, coupled to a gas-flow heat exchanger (HX). For these transients the core power was controlled by a point kinetics model with reactivity feedback based on core average temperature; the neutron generation time and the temperature feedback coefficient are provided as model inputs. This type of non-nuclear test is expected to provide reasonable approximation of reactor transient behavior because reactivity feedback is very simple in a compact fast reactor (simple, negative, and relatively monotonic temperature feedback, caused mostly by thermal expansion) and calculations show there are no significant reactivity effects associated with fluid in the HP (the worth of the entire inventory of Na in the core is .

Computational Models of Relational Processes in Cognitive Development

ERIC Educational Resources Information Center

Halford, Graeme S.; Andrews, Glenda; Wilson, William H.; Phillips, Steven

2012-01-01

Acquisition of relational knowledge is a core process in cognitive development. Relational knowledge is dynamic and flexible, entails structure-consistent mappings between representations, has properties of compositionality and systematicity, and depends on binding in working memory. We review three types of computational models relevant to…

Recent changes of the Earth's core derived from satellite observations of magnetic and gravity fields.

PubMed

Mandea, Mioara; Panet, Isabelle; Lesur, Vincent; de Viron, Olivier; Diament, Michel; Le Mouël, Jean-Louis

2012-11-20

To understand the dynamics of the Earth's fluid, iron-rich outer core, only indirect observations are available. The Earth's magnetic field, originating mainly within the core, and its temporal variations can be used to infer the fluid motion at the top of the core, on a decadal and subdecadal time-scale. Gravity variations resulting from changes in the mass distribution within the Earth may also occur on the same time-scales. Such variations include the signature of the flow inside the core, though they are largely dominated by the water cycle contributions. Our study is based on 8 y of high-resolution, high-accuracy magnetic and gravity satellite data, provided by the CHAMP and GRACE missions. From the newly derived geomagnetic models we have computed the core magnetic field, its temporal variations, and the core flow evolution. From the GRACE CNES/GRGS series of time variable geoid models, we have obtained interannual gravity models by using specifically designed postprocessing techniques. A correlation analysis between the magnetic and gravity series has demonstrated that the interannual changes in the second time derivative of the core magnetic field under a region from the Atlantic to Indian Ocean coincide in phase with changes in the gravity field. The order of magnitude of these changes and proposed correlation are plausible, compatible with a core origin; however, a complete theoretical model remains to be built. Our new results and their broad geophysical significance could be considered when planning new Earth observation space missions and devising more sophisticated Earth's interior models.

Development of the Joint NASA/NCAR General Circulation Model

NASA Technical Reports Server (NTRS)

Lin, S.-J.; Rood, R. B.

1999-01-01

The Data Assimilation Office at NASA/Goddard Space Flight Center is collaborating with NCAR/CGD in an ambitious proposal for the development of a unified climate, numerical weather prediction, and chemistry transport model which is suitable for global data assimilation of the physical and chemical state of the Earth's atmosphere. A prototype model based on the NCAR CCM3 physics and the NASA finite-volume dynamical core has been built. A unique feature of the NASA finite-volume dynamical core is its advanced tracer transport algorithm on the floating Lagrangian control-volume coordinate. The model currently has a highly idealized ozone production/loss chemistry derived from the observed 2D (latitude-height) climatology of the recent decades. Nevertheless, the simulated horizontal wave structure of the total ozone is in good qualitative agreement with the observed (TOMS). Long term climate simulations and NWP experiments have been carried out. Current up to date status and futur! e plan will be discussed in the conference.

Implementing Molecular Dynamics for Hybrid High Performance Computers - 1. Short Range Forces

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brown, W Michael; Wang, Peng; Plimpton, Steven J

The use of accelerators such as general-purpose graphics processing units (GPGPUs) have become popular in scientific computing applications due to their low cost, impressive floating-point capabilities, high memory bandwidth, and low electrical power requirements. Hybrid high performance computers, machines with more than one type of floating-point processor, are now becoming more prevalent due to these advantages. In this work, we discuss several important issues in porting a large molecular dynamics code for use on parallel hybrid machines - 1) choosing a hybrid parallel decomposition that works on central processing units (CPUs) with distributed memory and accelerator cores with shared memory,more » 2) minimizing the amount of code that must be ported for efficient acceleration, 3) utilizing the available processing power from both many-core CPUs and accelerators, and 4) choosing a programming model for acceleration. We present our solution to each of these issues for short-range force calculation in the molecular dynamics package LAMMPS. We describe algorithms for efficient short range force calculation on hybrid high performance machines. We describe a new approach for dynamic load balancing of work between CPU and accelerator cores. We describe the Geryon library that allows a single code to compile with both CUDA and OpenCL for use on a variety of accelerators. Finally, we present results on a parallel test cluster containing 32 Fermi GPGPUs and 180 CPU cores.« less

Density-Dependent Regulation of Brook Trout Population Dynamics along a Core-Periphery Distribution Gradient in a Central Appalachian Watershed

PubMed Central

Huntsman, Brock M.; Petty, J. Todd

2014-01-01

Spatial population models predict strong density-dependence and relatively stable population dynamics near the core of a species' distribution with increasing variance and importance of density-independent processes operating towards the population periphery. Using a 10-year data set and an information-theoretic approach, we tested a series of candidate models considering density-dependent and density-independent controls on brook trout population dynamics across a core-periphery distribution gradient within a central Appalachian watershed. We sampled seven sub-populations with study sites ranging in drainage area from 1.3–60 km2 and long-term average densities ranging from 0.335–0.006 trout/m. Modeled response variables included per capita population growth rate of young-of-the-year, adult, and total brook trout. We also quantified a stock-recruitment relationship for the headwater population and coefficients of variability in mean trout density for all sub-populations over time. Density-dependent regulation was prevalent throughout the study area regardless of stream size. However, density-independent temperature models carried substantial weight and likely reflect the effect of year-to-year variability in water temperature on trout dispersal between cold tributaries and warm main stems. Estimated adult carrying capacities decreased exponentially with increasing stream size from 0.24 trout/m in headwaters to 0.005 trout/m in the main stem. Finally, temporal variance in brook trout population size was lowest in the high-density headwater population, tended to peak in mid-sized streams and declined slightly in the largest streams with the lowest densities. Our results provide support for the hypothesis that local density-dependent processes have a strong control on brook trout dynamics across the entire distribution gradient. However, the mechanisms of regulation likely shift from competition for limited food and space in headwater streams to competition for thermal refugia in larger main stems. It also is likely that source-sink dynamics and dispersal from small headwater habitats may partially influence brook trout population dynamics in the main stem. PMID:24618602

Design of Broadband High Dynamic-Range Fiber Optic Links

NASA Astrophysics Data System (ADS)

Monsurrò, P.; Tommasino, P.; Trifiletti, A.; Vannucci, A.

2018-04-01

An analytic design-oriented model of microwave optical links has been developed. The core of the model is the non-linear and noise model of a Mach-Zehnder LiNbO3 interferometer. Both a 100 MHz-20 GHz link and a linearized microwave link, comprising an auxiliary modulator, have been designed and prototyped by using the model.

Earth system modelling on system-level heterogeneous architectures: EMAC (version 2.42) on the Dynamical Exascale Entry Platform (DEEP)

NASA Astrophysics Data System (ADS)

Christou, Michalis; Christoudias, Theodoros; Morillo, Julián; Alvarez, Damian; Merx, Hendrik

2016-09-01

We examine an alternative approach to heterogeneous cluster-computing in the many-core era for Earth system models, using the European Centre for Medium-Range Weather Forecasts Hamburg (ECHAM)/Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model as a pilot application on the Dynamical Exascale Entry Platform (DEEP). A set of autonomous coprocessors interconnected together, called Booster, complements a conventional HPC Cluster and increases its computing performance, offering extra flexibility to expose multiple levels of parallelism and achieve better scalability. The EMAC model atmospheric chemistry code (Module Efficiently Calculating the Chemistry of the Atmosphere (MECCA)) was taskified with an offload mechanism implemented using OmpSs directives. The model was ported to the MareNostrum 3 supercomputer to allow testing with Intel Xeon Phi accelerators on a production-size machine. The changes proposed in this paper are expected to contribute to the eventual adoption of Cluster-Booster division and Many Integrated Core (MIC) accelerated architectures in presently available implementations of Earth system models, towards exploiting the potential of a fully Exascale-capable platform.

Cooperative effects in spherical spasers: Ab initio analytical model

NASA Astrophysics Data System (ADS)

Bordo, V. G.

2017-06-01

A fully analytical semiclassical theory of cooperative optical processes which occur in an ensemble of molecules embedded in a spherical core-shell nanoparticle is developed from first principles. Both the plasmonic Dicke effect and spaser generation are investigated for the designs in which a shell/core contains an arbitrarily large number of active molecules in the vicinity of a metallic core/shell. An essential aspect of the theory is an ab initio account of the feedback from the core/shell boundaries which significantly modifies the molecular dynamics. The theory provides rigorous, albeit simple and physically transparent, criteria for both plasmonic superradiance and surface plasmon generation.

Modifications to WRFs dynamical core to improve the treatment of moisture for large-eddy simulations

DOE PAGES

Xiao, Heng; Endo, Satoshi; Wong, May; ...

2015-10-29

Yamaguchi and Feingold (2012) note that the cloud fields in their large-eddy simulations (LESs) of marine stratocumulus using the Weather Research and Forecasting (WRF) model exhibit a strong sensitivity to time stepping choices. In this study, we reproduce and analyze this sensitivity issue using two stratocumulus cases, one marine and one continental. Results show that (1) the sensitivity is associated with spurious motions near the moisture jump between the boundary layer and the free atmosphere, and (2) these spurious motions appear to arise from neglecting small variations in water vapor mixing ratio (qv) in the pressure gradient calculation in themore » acoustic sub-stepping portion of the integration procedure. We show that this issue is remedied in the WRF dynamical core by replacing the prognostic equation for the potential temperature θ with one for the moist potential temperature θm=θ(1+1.61qv), which allows consistent treatment of moisture in the calculation of pressure during the acoustic sub-steps. With this modification, the spurious motions and the sensitivity to the time stepping settings (i.e., the dynamic time step length and number of acoustic sub-steps) are eliminated in both of the example stratocumulus cases. In conclusion, this modification improves the applicability of WRF for LES applications, and possibly other models using similar dynamical core formulations, and also permits the use of longer time steps than in the original code.« less

The Effect of Surface Topography on the Nonlinear Dynamics of Rossby Waves

NASA Technical Reports Server (NTRS)

Abarzhi, S. I.; Desjardins, O.; Pitsch, H.

2003-01-01

Boussinesq convection in rotating systems attracts a sustained attention of the fluid dynamics community, because it has intricate non-linear dynamics (Cross & Hohenberg 1993) and plays an important role in geophysical and astrophysical applications, such as the motion of the liquid outer core of Earth, the Red Spot in Jupiter, the giant cells in the Sun etc. (Alridge et al. 1990). A fundamental distinction between the real geo- and astrophysical problems and the idealized laboratory studies is that natural systems are inhomogeneous (Alridge et al. 1990). Heterogeneities modulate the flow and influence significantly the dynamics of convective patterns (Alridge et al. 1990; Hide 1971). The effect of modulations on pattern formation and transition to turbulence in Boussinesq convection is far from being completely understood (Cross & Hohenberg 1993; Aranson & Kramer 2002). It is generally accepted that in the liquid outer core of the Earth the transport of the angular momentum and internal heat occurs via thermal Rossby waves (Zhang et al. 2001; Kuang & Bloxham 1999). These waves been visualized in laboratory experiments in rotating liquid-filled spheres and concentric spherical shells (Zhang et al. 2001; Kuang & Bloxham 1999). The basic dynamical features of Rossby waves have been reproduced in a cylindrical annulus, a system much simpler than the spherical ones (Busse & Or 1986; Or & Busse 1987). For convection in a cylindrical annulus, the fluid motion is two-dimensional, and gravity is replaced by a centrifugal force, (Busse & Or 1986; Or & Busse 1987). Hide (1971) has suggested that the momentum and heat transport in the core might be influenced significantly by so-called bumps, which are heterogeneities on the mantle-core boundary. To model the effect of surface topography on the transport of momentum and energy in the liquid outer core of the Earth, Bell & Soward (1996), Herrmann & Busse (1998) and Westerburg & Busse (2001) have studied the nonlinear dynamics of thermal Rossby waves in a cylindrical annulus with azimuthally modulated height.

Analysis of the uncertainties in the physical calculations of water-moderated power reactors of the VVER type by the parameters of models of preparing few-group constants

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bryukhin, V. V., E-mail: bryuhin@yandex.ru; Kurakin, K. Yu.; Uvakin, M. A.

The article covers the uncertainty analysis of the physical calculations of the VVER reactor core for different meshes of the reference values of the feedback parameters (FBP). Various numbers of nodes of the parametric axes of FBPs and different ranges between them are investigated. The uncertainties of the dynamic calculations are analyzed using RTS RCCA ejection as an example within the framework of the model with the boundary conditions at the core inlet and outlet.

Rock Physical Interpretation of the Relationship between Dynamic and Static Young's Moduli of Sedimentary Rocks

NASA Astrophysics Data System (ADS)

Takahashi, T.

2017-12-01

The static Young's modulus (deformability) of a rock is indispensable for designing and constructing tunnels, dams and underground caverns in civil engineering. Static Young's modulus which is an elastic modulus at large strain level is usually obtained with the laboratory tests of rock cores sampled in boreholes drilled in a rock mass. A deformability model of the entire rock mass is then built by extrapolating the measurements based on a rock mass classification obtained in geological site characterization. However, model-building using data obtained from a limited number of boreholes in the rock mass, especially a complex rock mass, may cause problems in the accuracy and reliability of the model. On the other hand, dynamic Young's modulus which is the modulus at small strain level can be obtained from seismic velocity. If dynamic Young's modulus can be rationally converted to static one, a seismic velocity model by the seismic method can be effectively used to build a deformability model of the rock mass. In this study, we have, therefore, developed a rock physics model (Mavko et al., 2009) to estimate static Young's modulus from dynamic one for sedimentary rocks. The rock physics model has been generally applied to seismic properties at small strain level. In the proposed model, however, the sandy shale model, one of rock physics models, is extended for modeling the static Young's modulus at large strain level by incorporating the mixture of frictional and frictionless grain contacts into the Hertz-Mindlin model. The proposed model is verified through its application to the dynamic Young's moduli derived from well log velocities and static Young's moduli measured in the tri-axial compression tests of rock cores sampled in the same borehole as the logs were acquired. This application proves that the proposed rock physics model can be possibly used to estimate static Young's modulus (deformability) which is required in many types of civil engineering applications from seismically derived dynamic Young's modulus. References:Mavko, G., Mukerji, T. and Dvorkin, J., 2009, The Rock Physics Handbook, 2nd Edition, Cambridge University Press, Cambridge.

Convective radiation fluid-dynamics: formation and early evolution of ultra low-mass objects

NASA Astrophysics Data System (ADS)

Wuchterl, G.

2005-12-01

The formation process of ultra low-mass objects is some kind of extension of the star formation process. The physical changes towards lower mass are discussed by investigating the collapse of cloud cores that are modelled as Bonnor-Ebert spheres. Their collapse is followed by solving the equations of fluid dynamics with radiation and a model of time-dependent convection that has been calibrated to the Sun. For a sequence of cloud-cores with 1 to 0.01 solar masses, evolutionary tracks and isochrones are shown in the mass-radius diagram, the Hertzsprung-Russel diagram and the effective temperature-surface gravity or Kiel diagram. The collapse and the early hydrostatic evolution to ages of few Ma are briefly discussed and compared to observations of objects in Upper Scorpius and the low-mass components of GG Tau.

The Density Jump at the Inner Core Boundary in the Eastern and Western Hemispheres

NASA Astrophysics Data System (ADS)

Krasnoshchekov, D. N.; Ovtchinnikov, V. M.

2018-02-01

The results of analysis of more than 1300 new PKiKP/PcP amplitude ratios measured in Southeast Asia and South America at the epicentral distances of 3.2°-35.2° are presented. The density jump in the Eastern Hemisphere of the Earth's inner core (IC) is 0.3 g/cm3, and it is 0.9 g/cm3 in the Western one. Taking the large discrepancy in the obtained estimates into consideration, maintenance of such large lateral variations in the mosaic properties of the IC reflecting surface requires considerable variations in the thermodynamic parameters (mostly temperature) of the inner-outer core transition. However, if the observed asymmetry in the density jump distribution is of a global character, the data presented support the translation model of the IC dynamics. This model implies IC crystallization in the Western Hemisphere and melting in the Eastern one, not vice versa, as suggested by another geodynamic model based on thermochemical convection in the outer core and the thermal balance of the core-mantle system.

Intracellular hepatitis C modeling predicts infection dynamics and viral protein mechanisms

DOE PAGES

Aunins, Thomas R.; Marsh, Katherine M.; Subramanya, Gitanjali; ...

2018-03-21

Hepatitis C virus infection is a global health problem, with nearly 2 million new infections occurring every year and up to 85% of these becoming chronic infections that pose serious long-term health risks. To effectively reduce the prevalence of HCV infection and associated diseases, it is important to understand the intracellular dynamics of the viral lifecycle. Here, we present a detailed mathematical model that represents the full hepatitis C lifecycle. It is the first full HCV model to be fit to acute intracellular infection data and the first to explore the functions of distinct viral proteins, probing multiple hypotheses ofmore » cis- and trans-acting mechanisms to provide insights for drug targeting. Model parameters were derived from the literature, experiments, and fitting to experimental intracellular viral RNA, extracellular viral titer, and HCV core and NS3 protein kinetic data from viral inoculation to steady-state. Our model predicts faster rates for protein translation and polyprotein cleavage than previous replicon models and demonstrates that the processes of translation and synthesis of viral RNA have the most influence on the levels of the species we tracked in experiments. Overall, our experimental data and the resulting mathematical infection model reveal information about the regulation of core protein during infection, produce specific insights into the roles of the viral core, NS5A, and NS5B proteins, and demonstrate the sensitivities of viral proteins and RNA to distinct reactions within the lifecycle.« less

Intracellular hepatitis C modeling predicts infection dynamics and viral protein mechanisms

DOE Office of Scientific and Technical Information (OSTI.GOV)

Aunins, Thomas R.; Marsh, Katherine M.; Subramanya, Gitanjali

Hepatitis C virus infection is a global health problem, with nearly 2 million new infections occurring every year and up to 85% of these becoming chronic infections that pose serious long-term health risks. To effectively reduce the prevalence of HCV infection and associated diseases, it is important to understand the intracellular dynamics of the viral lifecycle. Here, we present a detailed mathematical model that represents the full hepatitis C lifecycle. It is the first full HCV model to be fit to acute intracellular infection data and the first to explore the functions of distinct viral proteins, probing multiple hypotheses ofmore » cis- and trans-acting mechanisms to provide insights for drug targeting. Model parameters were derived from the literature, experiments, and fitting to experimental intracellular viral RNA, extracellular viral titer, and HCV core and NS3 protein kinetic data from viral inoculation to steady-state. Our model predicts faster rates for protein translation and polyprotein cleavage than previous replicon models and demonstrates that the processes of translation and synthesis of viral RNA have the most influence on the levels of the species we tracked in experiments. Overall, our experimental data and the resulting mathematical infection model reveal information about the regulation of core protein during infection, produce specific insights into the roles of the viral core, NS5A, and NS5B proteins, and demonstrate the sensitivities of viral proteins and RNA to distinct reactions within the lifecycle.« less

Symbolic Model of Perception in Dynamic 3D Environments

DTIC Science & Technology

2006-11-01

can retrieve memories , work on goals, recognize visual or aural percepts, and perform actions. ACT-R has been selected for the current...types of memory . Procedural memory is the store of condition- action productions that are selected and executed by the core production system...a declarative memory chunk that is made available to the core production system through the vision module . 4 The vision module has been

Research on Shock Responses of Three Types of Honeycomb Cores

NASA Astrophysics Data System (ADS)

Peng, Fei; Yang, Zhiguang; Jiang, Liangliang; Ren, Yanting

2018-03-01

The shock responses of three kinds of honeycomb cores have been investigated and analyzed based on explicit dynamics analysis. According to the real geometric configuration and the current main manufacturing methods of aluminum alloy honeycomb cores, the finite element models of honeycomb cores with three different cellular configurations (conventional hexagon honeycomb core, rectangle honeycomb core and auxetic honeycomb core with negative Poisson’s ratio) have been established through FEM parametric modeling method based on Python and Abaqus. In order to highlight the impact response characteristics of the above three honeycomb cores, a 5 mm thick panel with the same mass and material was taken as contrast. The analysis results showed that the peak values of longitudinal acceleration history curves of the three honeycomb cores were lower than those of the aluminum alloy panel in all three reference points under the loading of a longitudinal pulse pressure load with the peak value of 1 MPa and the pulse width of 1 μs. It could be concluded that due to the complex reflection and diffraction of stress wave induced by shock in honeycomb structures, the impact energy was redistributed which led to a decrease in the peak values of the longitudinal acceleration at the measuring points of honeycomb cores relative to the panel.

Modelling the optical properties of aerosols in a chemical transport model

NASA Astrophysics Data System (ADS)

Andersson, E.; Kahnert, M.

2015-12-01

According to the IPCC fifth assessment report (2013), clouds and aerosols still contribute to the largest uncertainty when estimating and interpreting changes to the Earth's energy budget. Therefore, understanding the interaction between radiation and aerosols is both crucial for remote sensing observations and modelling the climate forcing arising from aerosols. Carbon particles are the largest contributor to the aerosol absorption of solar radiation, thereby enhancing the warming of the planet. Modelling the radiative properties of carbon particles is a hard task and involves many uncertainties arising from the difficulties of accounting for the morphologies and heterogeneous chemical composition of the particles. This study aims to compare two ways of modelling the optical properties of aerosols simulated by a chemical transport model. The first method models particle optical properties as homogeneous spheres and are externally mixed. This is a simple model that is particularly easy to use in data assimilation methods, since the optics model is linear. The second method involves a core-shell internal mixture of soot, where sulphate, nitrate, ammonia, organic carbon, sea salt, and water are contained in the shell. However, by contrast to previously used core-shell models, only part of the carbon is concentrated in the core, while the remaining part is homogeneously mixed with the shell. The chemical transport model (CTM) simulations are done regionally over Europe with the Multiple-scale Atmospheric Transport and CHemistry (MATCH) model, developed by the Swedish Meteorological and Hydrological Institute (SMHI). The MATCH model was run with both an aerosol dynamics module, called SALSA, and with a regular "bulk" approach, i.e., a mass transport model without aerosol dynamics. Two events from 2007 are used in the analysis, one with high (22/12-2007) and one with low (22/6-2007) levels of elemental carbon (EC) over Europe. The results of the study help to assess the significance of aerosol morphology for modelling radiative forcing and aerosol optical properties relevant to interpreting remote sensing observations. The uncertainties introduced by the optics model are gauged by comparing them to model uncertainties related to the inclusion or omission of aerosol dynamic processes.

Core surface magnetic field evolution 2000-2010

NASA Astrophysics Data System (ADS)

Finlay, C. C.; Jackson, A.; Gillet, N.; Olsen, N.

2012-05-01

We present new dedicated core surface field models spanning the decade from 2000.0 to 2010.0. These models, called gufm-sat, are based on CHAMP, Ørsted and SAC-C satellite observations along with annual differences of processed observatory monthly means. A spatial parametrization of spherical harmonics up to degree and order 24 and a temporal parametrization of sixth-order B-splines with 0.25 yr knot spacing is employed. Models were constructed by minimizing an absolute deviation measure of misfit along with measures of spatial and temporal complexity at the core surface. We investigate traditional quadratic or maximum entropy regularization in space, and second or third time derivative regularization in time. Entropy regularization allows the construction of models with approximately constant spectral slope at the core surface, avoiding both the divergence characteristic of the crustal field and the unrealistic rapid decay typical of quadratic regularization at degrees above 12. We describe in detail aspects of the models that are relevant to core dynamics. Secular variation and secular acceleration are found to be of lower amplitude under the Pacific hemisphere where the core field is weaker. Rapid field evolution is observed under the eastern Indian Ocean associated with the growth and drift of an intense low latitude flux patch. We also find that the present axial dipole decay arises from a combination of subtle changes in the southern hemisphere field morphology.

Axisymmetric simulations of magnetorotational core collapse: approximate inclusion of general relativistic effects

NASA Astrophysics Data System (ADS)

Obergaulinger, M.; Aloy, M. A.; Dimmelmeier, H.; Müller, E.

2006-10-01

We continue our investigations of the magnetorotational collapse of stellar cores by discussing simulations performed with a modified Newtonian gravitational potential that mimics general relativistic effects. The approximate TOV gravitational potential used in our simulations captures several basic features of fully relativistic simulations quite well. In particular, it is able to correctly reproduce the behavior of models that show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. For models where the dynamics and gravitational wave signals are already captured qualitatively correctly by a Newtonian potential, the results of the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential, allowing for a more efficient amplification of the magnetic field by differential rotation. The strength of the saturation fields (˜ 1015 ~ G at the surface of the inner core) is a factor of two to three higher than in Newtonian gravity. Due to the more efficient field amplification, the influence of magnetic fields is considerably more pronounced than in the Newtonian case for some of the models. As in the Newtonian case, sufficiently strong magnetic fields slow down the core's rotation and trigger a secular contraction phase to higher densities. More clearly than in Newtonian models, the collapsed cores of these models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce at subnuclear densities does not suffice for ejecting any mass, and the temporarily stabilized core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Therefore, in the first paper of this series we introduced a new type of gravitational wave signal, which we call type IV or “magnetic type”. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Most of our weak-field models are marginally detectable with the current LIGO interferometer for a source located at a distance of 10 kpc. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and ⪉ 100 kHz denotes the generation of a jet-like hydromagnetic outflow.

A coarse grained molecular dynamics simulation study on the structural properties of carbon nanotube-dendrimer composites.

PubMed

Kavyani, Sajjad; Dadvar, Mitra; Modarress, Hamid; Amjad-Iranagh, Sepideh

2018-04-25

By employing coarse grained (CG) molecular dynamics (MD) simulation, the effect of the size and hydrophilic/hydrophobic properties of the interior/exterior structures of the dendrimers in carbon nanotube (CNT)-dendrimer composites has been studied, to find a stable composite with high solubility in water and the capability to be used in drug delivery applications. For this purpose, composites consisting of core-shell dendrimer complexes including: [PPI{core}-PAMAM{shell}], [PAMAM{core}-polyethyleneglycol (PEG){shell}] and [PAMAM{core}-fattyacid (FTA){shell}] were constructed. A new CG model for the fatty acid (FTA) molecules as functionalized to the dendrimer was developed, which, unlike the previous models, could generate the structural conformations of the FTA properly. The obtained results indicated that the dendrimer complexes with short FTA chains can form stable composites with the CNT. Also, it was found that the pristine PAMAM and PPI-PAMAM with small PPI, and PAMAM-PEG dendrimers with short PEG chains, can distribute their chains into the water medium and interact with the CNT efficiently, to form a stable water-soluble CNT-dendrimer composite. The results demonstrated that the structural difference between the interior and exterior of a core-shell dendrimer complex can prevent the core and the interior layers of the dendrimer complex from interacting with the CNT. An overall analysis of the results manifested that the CNT-PAMAM:4-PEG:4 is the most stable composite, due to strong binding of the dendrimer with the CNT while also having high solubility in water, and its core retains its structure properly and unchanged, suitable for encapsulating drugs in the targeted delivery applications.

Recent Developments in Smart Adaptive Structures for Solar Sailcraft

NASA Technical Reports Server (NTRS)

Whorton, M. S.; Kim, Y. K.; Oakley, J.; Adetona, O.; Keel, L. H.

2007-01-01

The "Smart Adaptive Structures for Solar Sailcraft" development activity at MSFC has investigated issues associated with understanding how to model and scale the subsystem and multi-body system dynamics of a gossamer solar sailcraft with the objective of designing sailcraft attitude control systems. This research and development activity addressed three key tasks that leveraged existing facilities and core competencies of MSFC to investigate dynamics and control issues of solar sails. Key aspects of this effort included modeling and testing of a 30 m deployable boom; modeling of the multi-body system dynamics of a gossamer sailcraft; investigation of control-structures interaction for gossamer sailcraft; and development and experimental demonstration of adaptive control technologies to mitigate control-structures interaction.

ReaDDy - A Software for Particle-Based Reaction-Diffusion Dynamics in Crowded Cellular Environments

PubMed Central

Schöneberg, Johannes; Noé, Frank

2013-01-01

We introduce the software package ReaDDy for simulation of detailed spatiotemporal mechanisms of dynamical processes in the cell, based on reaction-diffusion dynamics with particle resolution. In contrast to other particle-based reaction kinetics programs, ReaDDy supports particle interaction potentials. This permits effects such as space exclusion, molecular crowding and aggregation to be modeled. The biomolecules simulated can be represented as a sphere, or as a more complex geometry such as a domain structure or polymer chain. ReaDDy bridges the gap between small-scale but highly detailed molecular dynamics or Brownian dynamics simulations and large-scale but little-detailed reaction kinetics simulations. ReaDDy has a modular design that enables the exchange of the computing core by efficient platform-specific implementations or dynamical models that are different from Brownian dynamics. PMID:24040218

Normal-mode and free-Air gravity constraints on lateral variations in velocity and density of Earth's mantle

PubMed

Ishii; Tromp

1999-08-20

With the use of a large collection of free-oscillation data and additional constraints imposed by the free-air gravity anomaly, lateral variations in shear velocity, compressional velocity, and density within the mantle; dynamic topography on the free surface; and topography on the 660-km discontinuity and the core-mantle boundary were determined. The velocity models are consistent with existing models based on travel-time and waveform inversions. In the lowermost mantle, near the core-mantle boundary, denser than average material is found beneath regions of upwellings centered on the Pacific Ocean and Africa that are characterized by slow shear velocities. These anomalies suggest the existence of compositional heterogeneity near the core-mantle boundary.

Theoretical Models of Protostellar Binary and Multiple Systems with AMR Simulations

NASA Astrophysics Data System (ADS)

Matsumoto, Tomoaki; Tokuda, Kazuki; Onishi, Toshikazu; Inutsuka, Shu-ichiro; Saigo, Kazuya; Takakuwa, Shigehisa

2017-05-01

We present theoretical models for protostellar binary and multiple systems based on the high-resolution numerical simulation with an adaptive mesh refinement (AMR) code, SFUMATO. The recent ALMA observations have revealed early phases of the binary and multiple star formation with high spatial resolutions. These observations should be compared with theoretical models with high spatial resolutions. We present two theoretical models for (1) a high density molecular cloud core, MC27/L1521F, and (2) a protobinary system, L1551 NE. For the model for MC27, we performed numerical simulations for gravitational collapse of a turbulent cloud core. The cloud core exhibits fragmentation during the collapse, and dynamical interaction between the fragments produces an arc-like structure, which is one of the prominent structures observed by ALMA. For the model for L1551 NE, we performed numerical simulations of gas accretion onto protobinary. The simulations exhibit asymmetry of a circumbinary disk. Such asymmetry has been also observed by ALMA in the circumbinary disk of L1551 NE.

Molecular Simulation Studies of Covalently and Ionically Grafted Nanoparticles

NASA Astrophysics Data System (ADS)

Hong, Bingbing

Solvent-free covalently- or ionically-grafted nanoparticles (CGNs and IGNs) are a new class of organic-inorganic hybrid composite materials exhibiting fluid-like behaviors around room temperature. With similar structures to prior systems, e.g. nanocomposites, neutral or charged colloids, ionic liquids, etc, CGNs and IGNs inherit the functionality of inorganic nanopariticles, the facile processibility of polymers, as well as conductivity and nonvolatility from their constituent materials. In spite of the extensive prior experimental research having covered synthesis and measurements of thermal and dynamic properties, little progress in understanding of these new materials at the molecular level has been achieved, because of the lack of simulation work in this new area. Atomistic and coarse-grained molecular dynamics simulations have been performed in this thesis to investigate the thermodynamics, structure, and dynamics of these systems and to seek predictive methods predictable for their properties. Starting from poly(ethylene oxide) oligomers (PEO) melts, we established atomistic models based on united-atom representations of methylene. The Green-Kubo and Einstein-Helfand formulas were used to calculate the transport properties. The simulations generate densities, viscosities, diffusivities, in good agreement with experimental data. The chain-length dependence of the transport properties suggests that neither Rouse nor reptation models are applicable in the short-chain regime investigated. Coupled with thermodynamic integration methods, the models give good predictions of pressure-composition-density relations for CO 2 + PEO oligomers. Water effects on the Henry's constant of CO 2 in PEO have also been investigated. The dependence of the calculated Henry's constants on the weight percentage of water falls on a temperature-dependent master curve, irrespective of PEO chain length. CGNs are modeled by the inclusion of solid-sphere nanoparticles into the atomistic oligomers. The calculated viscosities from Green-Kubo relationships and temperature extrapolation are of the same order of magnitude as experimental values, but show a smaller activation energy relative to real CGNs systems. Grafted systems have higher viscosities, smaller diffusion coefficients, and slower chain dynamics than the ungrafted counterparts - nanocomposites - at high temperatures. At lower temperatures, grafted systems exhibit faster dynamics for both nanoparticles and chains relative to ungrafted systems, because of lower aggregation of nanoparticles and enhanced correlations between nanoparticles and chains. This agrees with the experimental observation that the new materials have liquid-like behavior in the absence of a solvent. To lower the simulated temperatures into the experimental range, we established a coarse-grained CGNs model by matching structural distribution functions to atomistic simulation data. In contrast with linear polymer systems, for which coarse-graining always accelerate dynamics, coarse-graining of grafted nanoparticles can either accelerate or slowdown the core motions, depending on the length of the grafted chains. This can be qualitatively predicted by a simple transition-state theory. Similar atomistic models to CGNs were developed for IGNs, with ammonium counterions described by an explicit-hydrogen way; these were in turn compared with "generic" coarse-grained IGNs. The elimination of chemical details in the coarse-grained models does not bring in qualitative changes to the radial distribution functions and diffusion of atomistic IGNs, but saves considerable simulation resources and make simulations near room temperatures affordable. The chain counterions in both atomistic and coarse-grained models are mobile, moving from site to site and from nanoparticle to nanoparticle. At the same temperature and the same core volume fractions, the nanoparticle diffusivities in coarse-grained IGNs are slower by a factor ten than the cores of CGNs. The coarse-grained IGNs models are later used to investigate the system dynamics through analysis of the dependence on temperature and structural parameters of the transport properties (self-diffusion coefficients, viscosities and conductivities). Further, migration kinetics of oligomeric counterions is analyzed in a manner analogous to unimer exchange between micellar aggregates. The counterion migrations follow the "double-core" mechanism and are kinetically controlled by neighboring-core collisions. (Abstract shortened by UMI.)

Dynamic Response of Functionally Graded Carbon Nanotube Reinforced Sandwich Plate

NASA Astrophysics Data System (ADS)

Mehar, Kulmani; Panda, Subrata Kumar

2018-03-01

In this article, the dynamic response of the carbon nanotube-reinforced functionally graded sandwich composite plate has been studied numerically with the help of finite element method. The face sheets of the sandwich composite plate are made of carbon nanotube- reinforced composite for two different grading patterns whereas the core phase is taken as isotropic material. The final properties of the structure are calculated using the rule of mixture. The geometrical model of the sandwich plate is developed and discretized suitably with the help of available shell element in ANSYS library. Subsequently, the corresponding numerical dynamic responses computed via batch input technique (parametric design language code in ANSYS) of ANSYS including Newmark’s integration scheme. The stability of the sandwich structural numerical model is established through the proper convergence study. Further, the reliability of the sandwich model is checked by comparison study between present and available results from references. As a final point, some numerical problems have been solved to examine the effect of different design constraints (carbon nanotube distribution pattern, core to face thickness ratio, volume fractions of the nanotube, length to thickness ratio, aspect ratio and constraints at edges) on the time-responses of sandwich plate.

A finite-volume module for all-scale Earth-system modelling at ECMWF

NASA Astrophysics Data System (ADS)

Kühnlein, Christian; Malardel, Sylvie; Smolarkiewicz, Piotr

2017-04-01

We highlight recent advancements in the development of the finite-volume module (FVM) (Smolarkiewicz et al., 2016) for the IFS at ECMWF. FVM represents an alternative dynamical core that complements the operational spectral dynamical core of the IFS with new capabilities. Most notably, these include a compact-stencil finite-volume discretisation, flexible meshes, conservative non-oscillatory transport and all-scale governing equations. As a default, FVM solves the compressible Euler equations in a geospherical framework (Szmelter and Smolarkiewicz, 2010). The formulation incorporates a generalised terrain-following vertical coordinate. A hybrid computational mesh, fully unstructured in the horizontal and structured in the vertical, enables efficient global atmospheric modelling. Moreover, a centred two-time-level semi-implicit integration scheme is employed with 3D implicit treatment of acoustic, buoyant, and rotational modes. The associated 3D elliptic Helmholtz problem is solved using a preconditioned Generalised Conjugate Residual approach. The solution procedure employs the non-oscillatory finite-volume MPDATA advection scheme that is bespoke for the compressible dynamics on the hybrid mesh (Kühnlein and Smolarkiewicz, 2017). The recent progress of FVM is illustrated with results of benchmark simulations of intermediate complexity, and comparison to the operational spectral dynamical core of the IFS. C. Kühnlein, P.K. Smolarkiewicz: An unstructured-mesh finite-volume MPDATA for compressible atmospheric dynamics, J. Comput. Phys. (2017), in press. P.K. Smolarkiewicz, W. Deconinck, M. Hamrud, C. Kühnlein, G. Mozdzynski, J. Szmelter, N.P. Wedi: A finite-volume module for simulating global all-scale atmospheric flows, J. Comput. Phys. 314 (2016) 287-304. J. Szmelter, P.K. Smolarkiewicz: An edge-based unstructured mesh discretisation in geospherical framework, J. Comput. Phys. 229 (2010) 4980-4995.

Constraints on Inner Core Anisotropy Using Array Observations of P'P'

NASA Astrophysics Data System (ADS)

Frost, Daniel A.; Romanowicz, Barbara

2017-11-01

Recent studies of PKPdf travel times suggest strong anisotropy (4% or more) in the quasi-western inner core hemisphere. However, the availability of paths sampling at low angles to the Earth's rotation axis (the fast axis) is limited. To augment this sampling, we collected a travel time data set for the phase P'P'df (PKPPKPdf), for which at least one inner core leg is quasi-polar, at two high latitude seismic arrays. We find that the inferred anisotropy is weak (on the order of 0.5 to 1.5%), confirming previous results based on a much smaller P'P' data set. While previous models of inner core anisotropy required very strong alignment of anisotropic iron grains, our results are more easily explained by current dynamic models of inner core growth. We observe large travel time anomalies when one leg of P'P'df is along the South Sandwich to Alaska path, consistent with PKPdf observations, and warranting further investigation.

Magnetic properties and core electron binding energies of liquid water

NASA Astrophysics Data System (ADS)

Galamba, N.; Cabral, Benedito J. C.

2018-01-01

The magnetic properties and the core and inner valence electron binding energies of liquid water are investigated. The adopted methodology relies on the combination of molecular dynamics and electronic structure calculations. Born-Oppenheimer molecular dynamics with the Becke and Lee-Yang-Parr functionals for exchange and correlation, respectively, and includes an empirical correction (BLYP-D3) functional and classical molecular dynamics with the TIP4P/2005-F model were carried out. The Keal-Tozer functional was applied for predicting magnetic shielding and spin-spin coupling constants. Core and inner valence electron binding energies in liquid water were calculated with symmetry adapted cluster-configuration interaction. The relationship between the magnetic shielding constant σ(17O), the role played by the oxygen atom as a proton acceptor and donor, and the tetrahedral organisation of liquid water are investigated. The results indicate that the deshielding of the oxygen atom in water is very dependent on the order parameter (q) describing the tetrahedral organisation of the hydrogen bond network. The strong sensitivity of magnetic properties on changes of the electronic density in the nuclei environment is illustrated by a correlation between σ(17O) and the energy gap between the 1a1[O1s] (core) and the 2a1 (inner valence) orbitals of water. Although several studies discussed the eventual connection between magnetic properties and core electron binding energies, such a correlation could not be clearly established. Here, we demonstrate that for liquid water this correlation exists although involving the gap between electron binding energies of core and inner valence orbitals.

The effect of ilmenite viscosity on the dynamics and evolution of an overturned lunar cumulate mantle

NASA Astrophysics Data System (ADS)

Zhang, Nan; Dygert, Nick; Liang, Yan; Parmentier, E. M.

2017-07-01

Lunar cumulate mantle overturn and the subsequent upwelling of overturned mantle cumulates provide a potential framework for understanding the first-order thermochemical evolution of the Moon. Upwelling of ilmenite-bearing cumulates (IBCs) after the overturn has a dominant influence on the dynamics and long-term thermal evolution of the lunar mantle. An important parameter determining the stability and convective behavior of the IBC is its viscosity, which was recently constrained through rock deformation experiments. To examine the effect of IBC viscosity on the upwelling of overturned lunar cumulate mantle, here we conduct three-dimensional mantle convection models with an evolving core superposed by an IBC-rich layer, which resulted from mantle overturn after magma ocean solidification. Our modeling shows that a reduction of mantle viscosity by 1 order of magnitude, due to the presence of ilmenite, can dramatically change convective planform and long-term lunar mantle evolution. Our model results suggest a relatively stable partially molten IBC layer that has surrounded the lunar core to the present day.