Identification of damping in a bridge using a moving instrumented vehicle

NASA Astrophysics Data System (ADS)

González, A.; OBrien, E. J.; McGetrick, P. J.

2012-08-01

In recent years, there has been a significant increase in the number of bridges which are being instrumented and monitored on an ongoing basis. This is in part due to the introduction of bridge management systems designed to provide a high level of protection to the public and early warning if the bridge becomes unsafe. This paper investigates a novel alternative; a low-cost method consisting of the use of a vehicle fitted with accelerometers on its axles to monitor the dynamic behaviour of bridges. A simplified half-car vehicle-bridge interaction model is used in theoretical simulations to test the effectiveness of the approach in identifying the damping ratio of the bridge. The method is tested for a range of bridge spans and vehicle velocities using theoretical simulations and the influences of road roughness, initial vibratory condition of the vehicle, signal noise, modelling errors and frequency matching on the accuracy of the results are investigated.

Sensitivity of equatorial Pacific and Indian Ocean watermasses to the position of the Indonesian Throughflow

NASA Astrophysics Data System (ADS)

Rodgers, Keith B.; Latif, Mojib; Legutke, Stephanie

2000-09-01

The sensitivity of the thermal structure of the equatorial Pacific and Indian Ocean pycnoclines to a model's representation of the Indonesian Straits connecting the two basins is investigated. Two integrations are performed using the global HOPE ocean model. The initial conditions and surface forcing for both cases are identical; the only difference between the runs is that one has an opening for the Indonesian Straits which spans the equator on the Pacific side, and the other has an opening which lies fully north of the equator. The resulting sensitivity throughout much of the upper ocean is greater than 0.5°C for both the equatorial Indian and Pacific. A realistic simulation of net Indonesian Throughflow (ITF) transport (measured in Sverdrups) is not sufficient for an adequate simulation of equatorial watermasses. The ITF must also contain a realistic admixture of northern and southern Pacific source water.

Simulating and Forecasting Flooding Events in the City of Jeddah, Saudi Arabia

NASA Astrophysics Data System (ADS)

Ghostine, Rabih; Viswanadhapalli, Yesubabu; Hoteit, Ibrahim

2014-05-01

Metropolitan cities in the Kingdom of Saudi Arabia, as Jeddah and Riyadh, are more frequently experiencing flooding events caused by strong convective storms that produce intense precipitation over a short span of time. The flooding in the city of Jeddah in November 2009 was described by civil defense officials as the worst in 27 years. As of January 2010, 150 people were reported killed and more than 350 were missing. Another flooding event, less damaging but comparably spectacular, occurred one year later (Jan 2011) in Jeddah. Anticipating floods before they occur could minimize human and economic losses through the implementation of appropriate protection, provision and rescue plans. We have developed a coupled hydro-meteorological model for simulating and predicting flooding events in the city of Jeddah. We use the Weather Research Forecasting (WRF) model assimilating all available data in the Jeddah region for simulating the storm events in Jeddah. The resulting rain is then used on 10 minutes intervals to feed up an advanced numerical shallow water model that has been discretized on an unstructured grid using different numerical schemes based on the finite elements or finite volume techniques. The model was integrated on a high-resolution grid size varying between 0.5m within the streets of Jeddah and 500m outside the city. This contribution will present the flooding simulation system and the simulation results, focusing on the comparison of the different numerical schemes on the system performances in terms of accuracy and computational efficiency.

Workplace Simulation: An Integrated Approach to Training University Students in Professional Communication

ERIC Educational Resources Information Center

Ismail, Norhayati; Sabapathy, Chitra

2016-01-01

In the redesign of a professional communication course for real estate students, a workplace simulation was implemented, spanning the entire 12-week duration of the course. The simulation was achieved through the creation of an online company presence, the infusion of communication typically encountered in the workplace, and an intensive and…

Span efficiency of wings with leading edge protuberances

NASA Astrophysics Data System (ADS)

Custodio, Derrick; Henoch, Charles; Johari, Hamid

2013-11-01

Past work has shown that sinusoidal leading edge protuberances resembling those found on humpback whale flippers alter the lift and drag coefficients of full- and finite-span foils and wings depending on the angle of attack and leading edge geometry. Although the load characteristics of protuberance modified finite-span wings have been reported for flipper-like geometries at higher Reynolds numbers and for rectangular planforms at lower Reynolds numbers, the effects of leading edge geometry on the span efficiency, which is indicative of the deviation of the spanwise lift distribution from elliptical and the viscous effects, for a range of planforms and Reynolds numbers have not been addressed. The lift and drag coefficients of 7 rectangular, 2 swept, and 2 flipper-like planform models with aspect ratios of 4.3, 4.0, and 8.86, respectively, were used to compute the span efficiency at Reynolds numbers ranging from 0.9 to 4.5 × 105. The span efficiency, based on the data at lower angles of attack, of modified wings was compared with the unmodified models. For the cases considered, the span efficiencies of the leading edge modified models were less than those of the equivalent unmodified models. The dependence of span efficiency on the leading edge geometry, planform, and Reynolds number will be presented. Supported by the ONR-ULI program.

Optimized image acquisition for breast tomosynthesis in projection and reconstruction space.

PubMed

Chawla, Amarpreet S; Lo, Joseph Y; Baker, Jay A; Samei, Ehsan

2009-11-01

Breast tomosynthesis has been an exciting new development in the field of breast imaging. While the diagnostic improvement via tomosynthesis is notable, the full potential of tomosynthesis has not yet been realized. This may be attributed to the dependency of the diagnostic quality of tomosynthesis on multiple variables, each of which needs to be optimized. Those include dose, number of angular projections, and the total angular span of those projections. In this study, the authors investigated the effects of these acquisition parameters on the overall diagnostic image quality of breast tomosynthesis in both the projection and reconstruction space. Five mastectomy specimens were imaged using a prototype tomosynthesis system. 25 angular projections of each specimen were acquired at 6.2 times typical single-view clinical dose level. Images at lower dose levels were then simulated using a noise modification routine. Each projection image was supplemented with 84 simulated 3 mm 3D lesions embedded at the center of 84 nonoverlapping ROIs. The projection images were then reconstructed using a filtered backprojection algorithm at different combinations of acquisition parameters to investigate which of the many possible combinations maximizes the performance. Performance was evaluated in terms of a Laguerre-Gauss channelized Hotelling observer model-based measure of lesion detectability. The analysis was also performed without reconstruction by combining the model results from projection images using Bayesian decision fusion algorithm. The effect of acquisition parameters on projection images and reconstructed slices were then compared to derive an optimization rule for tomosynthesis. The results indicated that projection images yield comparable but higher performance than reconstructed images. Both modes, however, offered similar trends: Performance improved with an increase in the total acquisition dose level and the angular span. Using a constant dose level and angular span, the performance rolled off beyond a certain number of projections, indicating that simply increasing the number of projections in tomosynthesis may not necessarily improve its performance. The best performance for both projection images and tomosynthesis slices was obtained for 15-17 projections spanning an angular are of approximately 45 degrees--the maximum tested in our study, and for an acquisition dose equal to single-view mammography. The optimization framework developed in this framework is applicable to other reconstruction techniques and other multiprojection systems.

The Pilot Phase of the Global Soil Wetness Project Phase 3

NASA Astrophysics Data System (ADS)

Kim, H.; Oki, T.

2015-12-01

After the second phase of the Global Soil Wetness Project (GSWP2) as an early global continuous gridded multi-model analysis, a comprehensive set of land surface fluxes and state variables became available. It has been broadly utilized in the hydrology community, and its success has evolved to take advantages of recent scientific progress and to extend the relatively short time span (1986-1995) of the previous project. In the third phase proposed here (GSWP3), an extensive set of quantities for hydro-energy-eco systems will be produced to investigate their long-term (1901-2010) changes. The energy-water-carbon cycles and their interactions are also examined subcomponent-wise with appropriate model verifications in ensemble land simulations. In this study, the preliminary results and problems found from the first round analysis of the GSWP3 pilot study are shown. Also, it is discussed how the global offline simulation activity contributes to wider communities and a bigger scope such as Climate Model Intercomparison Project Phase 6 (CMIP6).

Continuum kinetic modeling of the tokamak plasma edge

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

Dorf, M. A.; Dorr, M. R.; Hittinger, J. A.

2016-05-15

The first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalous radial transport.

Experimental prediction of tube support interaction characteristics in steam generators: Volume 2, Westinghouse Model 51 flow entrance region: Topical report

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

Haslinger, K.H.

Tube-to-tube support interaction characterisitics were determined experimentally on a single tube, multi-span geometry, representative of the Westinghouse Model 51 steam generator economizer design. Results, in part, became input for an autoclave type wear test program on steam generator tubes, performed by Kraftwerk Union (KWU). More importantly, the test data reported here have been used to validate two analytical wear prediction codes; the WECAN code, which was developed by Westinghouse, and the ABAQUS code which has been enhanced for EPRI by Foster Wheeler to enable simulation of gap conditions (including fluid film effects) for various support geometries.

Taking the measure of a landscape: Comparing a simulated and natural landscape in the Virginia Coastal Plain

NASA Astrophysics Data System (ADS)

Howard, Alan D.; Tierney, Heather E.

2012-01-01

A landform evolution model is used to investigate the historical evolution of a fluvial landscape along the Potomac River in Virginia, USA. The landscape has developed on three terraces whose ages span 3.5 Ma. The simulation model specifies the temporal evolution of base level control by the river as having a high-frequency component of the response of the Potomac River to sea level fluctuations superimposed on a long-term epeirogenic uplift. The wave-cut benches are assumed to form instantaneously during sea level highstands. The region is underlain by relatively soft coastal plain sediments with high intrinsic erodibility. The survival of portions of these terrace surfaces, up to 3.5 Ma, is attributable to a protective cover of vegetation. The vegetation influence is parameterized as a critical shear stress to fluvial erosion whose magnitude decreases with increasing contributing area. The simulation model replicates the general pattern of dissection of the natural landscape, with decreasing degrees of dissection of the younger terrace surfaces. Channel incision and relief increase in headwater areas are most pronounced during the relatively brief periods of river lowstands. Imposition of the wave-cut terraces onto the simulated landscape triggers a strong incisional response. By qualitative and quantitative measures the model replicates, in a general way, the landform evolution and present morphology of the target region.

An elemental mercury diffusion coefficient for natural waters determined by molecular dynamics simulation.

PubMed

Kuss, Joachim; Holzmann, Jörg; Ludwig, Ralf

2009-05-01

Mercury is a priority pollutant as its mobility between the hydrosphere and the atmosphere threatens the biosphere globally. The air-water gas transfer of elemental mercury (Hg0) is controlled by its diffusion through the water-side boundary layer and thus by its diffusion coefficient, D(Hg), the value of which, however, has not been established. Here, the diffusion of Hg0 in water was modeled by molecular dynamics (MD) simulation and the diffusion coefficient subsequently determined. Therefore the movement of either Hg(0) or xenon and 1000 model water molecules (TIP4P-Ew) were traced for time spans of 50 ns. The modeled D(Xe) of the monatomic noble gas agreed well with measured data; thus, MD simulation was assumed to be a reliable approach to determine D(Hg) for monatomic Hg(0) as well. Accordingly, Hg(0) diffusion was then simulated for freshwater and seawater, and the data were well-described by the equation of Eyring. The activation energies for the diffusion of Hg0 in freshwater was 17.0 kJ mol(-1) and in seawater 17.8 kJ mol(-1). The newly determined D(Hg) is clearly lower than the one previously used for an oceanic mercury budget. Thus, its incorporation into the model should lead to lower estimates of global ocean mercury emissions.

A new vertical grid nesting capability in the Weather Research and Forecasting (WRF) Model

DOE PAGES

Daniels, Megan H.; Lundquist, Katherine A.; Mirocha, Jeffrey D.; ...

2016-09-16

Mesoscale atmospheric models are increasingly used for high-resolution (<3 km) simulations to better resolve smaller-scale flow details. Increased resolution is achieved using mesh refinement via grid nesting, a procedure where multiple computational domains are integrated either concurrently or in series. A constraint in the concurrent nesting framework offered by the Weather Research and Forecasting (WRF) Model is that mesh refinement is restricted to the horizontal dimensions. This limitation prevents control of the grid aspect ratio, leading to numerical errors due to poor grid quality and preventing grid optimization. Here, a procedure permitting vertical nesting for one-way concurrent simulation is developedmore » and validated through idealized cases. The benefits of vertical nesting are demonstrated using both mesoscale and large-eddy simulations (LES). Mesoscale simulations of the Terrain-Induced Rotor Experiment (T-REX) show that vertical grid nesting can alleviate numerical errors due to large aspect ratios on coarse grids, while allowing for higher vertical resolution on fine grids. Furthermore, the coarsening of the parent domain does not result in a significant loss of accuracy on the nested domain. LES of neutral boundary layer flow shows that, by permitting optimal grid aspect ratios on both parent and nested domains, use of vertical nesting yields improved agreement with the theoretical logarithmic velocity profile on both domains. Lastly, vertical grid nesting in WRF opens the path forward for multiscale simulations, allowing more accurate simulations spanning a wider range of scales than previously possible.« less

A new vertical grid nesting capability in the Weather Research and Forecasting (WRF) Model

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

Daniels, Megan H.; Lundquist, Katherine A.; Mirocha, Jeffrey D.

Mesoscale atmospheric models are increasingly used for high-resolution (<3 km) simulations to better resolve smaller-scale flow details. Increased resolution is achieved using mesh refinement via grid nesting, a procedure where multiple computational domains are integrated either concurrently or in series. A constraint in the concurrent nesting framework offered by the Weather Research and Forecasting (WRF) Model is that mesh refinement is restricted to the horizontal dimensions. This limitation prevents control of the grid aspect ratio, leading to numerical errors due to poor grid quality and preventing grid optimization. Here, a procedure permitting vertical nesting for one-way concurrent simulation is developedmore » and validated through idealized cases. The benefits of vertical nesting are demonstrated using both mesoscale and large-eddy simulations (LES). Mesoscale simulations of the Terrain-Induced Rotor Experiment (T-REX) show that vertical grid nesting can alleviate numerical errors due to large aspect ratios on coarse grids, while allowing for higher vertical resolution on fine grids. Furthermore, the coarsening of the parent domain does not result in a significant loss of accuracy on the nested domain. LES of neutral boundary layer flow shows that, by permitting optimal grid aspect ratios on both parent and nested domains, use of vertical nesting yields improved agreement with the theoretical logarithmic velocity profile on both domains. Lastly, vertical grid nesting in WRF opens the path forward for multiscale simulations, allowing more accurate simulations spanning a wider range of scales than previously possible.« less

A Global Lake Ecological Observatory Network (GLEON) for synthesising high-frequency sensor data for validation of deterministic ecological models

USGS Publications Warehouse

David, Hamilton P; Carey, Cayelan C.; Arvola, Lauri; Arzberger, Peter; Brewer, Carol A.; Cole, Jon J; Gaiser, Evelyn; Hanson, Paul C.; Ibelings, Bas W; Jennings, Eleanor; Kratz, Tim K; Lin, Fang-Pang; McBride, Christopher G.; de Motta Marques, David; Muraoka, Kohji; Nishri, Ami; Qin, Boqiang; Read, Jordan S.; Rose, Kevin C.; Ryder, Elizabeth; Weathers, Kathleen C.; Zhu, Guangwei; Trolle, Dennis; Brookes, Justin D

2014-01-01

A Global Lake Ecological Observatory Network (GLEON; www.gleon.org) has formed to provide a coordinated response to the need for scientific understanding of lake processes, utilising technological advances available from autonomous sensors. The organisation embraces a grassroots approach to engage researchers from varying disciplines, sites spanning geographic and ecological gradients, and novel sensor and cyberinfrastructure to synthesise high-frequency lake data at scales ranging from local to global. The high-frequency data provide a platform to rigorously validate process- based ecological models because model simulation time steps are better aligned with sensor measurements than with lower-frequency, manual samples. Two case studies from Trout Bog, Wisconsin, USA, and Lake Rotoehu, North Island, New Zealand, are presented to demonstrate that in the past, ecological model outputs (e.g., temperature, chlorophyll) have been relatively poorly validated based on a limited number of directly comparable measurements, both in time and space. The case studies demonstrate some of the difficulties of mapping sensor measurements directly to model state variable outputs as well as the opportunities to use deviations between sensor measurements and model simulations to better inform process understanding. Well-validated ecological models provide a mechanism to extrapolate high-frequency sensor data in space and time, thereby potentially creating a fully 3-dimensional simulation of key variables of interest.

Studying Turbulence Using Numerical Simulation Databases - IX: Proceedings of the 2002 Summer Program

NASA Technical Reports Server (NTRS)

Bradshaw, Peter (Editor); Rogers, Michael M. (Technical Monitor)

2002-01-01

The ninth Summer Program of the Center for Turbulence Research was held during the period July 29th - August 23rd, 2002. The increase in number of participants, noted in the Preface to the Proceedings of the 2000 Program, continues: this year there were 50 participants from ten countries, and 30 hosts from Stanford and NASA-Ames. This Proceedings volume contains 32 papers that span a wide range of topics and an enormous range of physical scales. The papers have been divided into seven groups: Acoustics, RANS modeling, Combustion, Large-eddy simulation (LES), LES Numerics, Stratified Flows, and Fundamentals, In several cases, a paper could have fitted in more than one group so the classification is somewhat arbitrary.

Assessment of dynamic properties and stiffness of composite bridges with pavement defects

NASA Astrophysics Data System (ADS)

Kartopol'tsev, Vladimir; Kartopol'tsev, Andrei; Kolmakov, Boris

2017-01-01

This paper is aimed at assessing the dynamic properties and stiffness of the reinforced concrete roadway slab under live loads that impact composite bridge girders considering pavement defects. A special attention is paid to the reinforced concrete roadway slab as a transfer member of forced oscillations. The test results obtained for bridges with different spans ranging from 24 to 110 m are presented to assess the behavior of the reinforced concrete roadway slab and the dynamic stiffness of bridge span allowed for the pavement defects. Dynamic tests are carried out under controlled and random loads that simulate live load interaction with the span and the pavement with defects. The differential equations are presented for vertical oscillations of spans, pavement defect parameter, Eigen frequency and others. As a result of the experimental research the equation is derived to ascertain the dynamic stiffness of the vehicle-span system.

An Assessment of Multimodel Simulations for the Variability of Western North Pacific Tropical Cyclones and Its Association with ENSO

NASA Technical Reports Server (NTRS)

Han, Rongqing; Wang, Hui; Hu, Zeng-Zhen; Kumar, Arun; Li, Weijing; Long, Lindsey N.; Schemm, Jae-Kyung E.; Peng, Peitao; Wang, Wanqiu; Si, Dong;

Exploratory and problem-solving consumer behavior across the life span.

PubMed

Lesser, J A; Kunkel, S R

1991-09-01

Different cognitive functioning, social, and personality changes appear to occur systematically during the adult life span. This article synthesizes research on life span changes in order to develop age-specific models of shopping behavior. The models are tested within a naturalistic field study of shoppers.

Conditional Probabilities of Large Earthquake Sequences in California from the Physics-based Rupture Simulator RSQSim

NASA Astrophysics Data System (ADS)

Gilchrist, J. J.; Jordan, T. H.; Shaw, B. E.; Milner, K. R.; Richards-Dinger, K. B.; Dieterich, J. H.

2017-12-01

Within the SCEC Collaboratory for Interseismic Simulation and Modeling (CISM), we are developing physics-based forecasting models for earthquake ruptures in California. We employ the 3D boundary element code RSQSim (Rate-State Earthquake Simulator of Dieterich & Richards-Dinger, 2010) to generate synthetic catalogs with tens of millions of events that span up to a million years each. This code models rupture nucleation by rate- and state-dependent friction and Coulomb stress transfer in complex, fully interacting fault systems. The Uniform California Earthquake Rupture Forecast Version 3 (UCERF3) fault and deformation models are used to specify the fault geometry and long-term slip rates. We have employed the Blue Waters supercomputer to generate long catalogs of simulated California seismicity from which we calculate the forecasting statistics for large events. We have performed probabilistic seismic hazard analysis with RSQSim catalogs that were calibrated with system-wide parameters and found a remarkably good agreement with UCERF3 (Milner et al., this meeting). We build on this analysis, comparing the conditional probabilities of sequences of large events from RSQSim and UCERF3. In making these comparisons, we consider the epistemic uncertainties associated with the RSQSim parameters (e.g., rate- and state-frictional parameters), as well as the effects of model-tuning (e.g., adjusting the RSQSim parameters to match UCERF3 recurrence rates). The comparisons illustrate how physics-based rupture simulators might assist forecasters in understanding the short-term hazards of large aftershocks and multi-event sequences associated with complex, multi-fault ruptures.

Connections between residence time distributions and watershed characteristics across the continental US

NASA Astrophysics Data System (ADS)

Condon, L. E.; Maxwell, R. M.; Kollet, S. J.; Maher, K.; Haggerty, R.; Forrester, M. M.

2016-12-01

Although previous studies have demonstrated fractal residence time distributions in small watersheds, analyzing residence time scaling over large spatial areas is difficult with existing observational methods. For this study we use a fully integrated groundwater surface water simulation combined with Lagrangian particle tracking to evaluate connections between residence time distributions and watershed characteristics such as geology, topography and climate. Our simulation spans more than six million square kilometers of the continental US, encompassing a broad range of watershed sizes and physiographic settings. Simulated results demonstrate power law residence time distributions with peak age rages from 1.5 to 10.5 years. These ranges agree well with previous observational work and demonstrate the feasibility of using integrated models to simulate residence times. Comparing behavior between eight major watersheds, we show spatial variability in both the peak and the variance of the residence time distributions that can be related to model inputs. Peak age is well correlated with basin averaged hydraulic conductivity and the semi-variance corresponds to aridity. While power law age distributions have previously been attributed to fractal topography, these results illustrate the importance of subsurface characteristics and macro climate as additional controls on groundwater configuration and residence times.

Aeroservoelastic Testing of Free Flying Wind Tunnel Models Part 1: A Sidewall Supported Semispan Model Tested for Gust Load Alleviation and Flutter Suppression

NASA Technical Reports Server (NTRS)

Scott, Robert C.; Vetter, Travis K.; Penning, Kevin B.; Coulson, David A.; Heeg, Jennifer.

2013-01-01

of a two part document. Part 2 is titled: "Aeroservoelastic Testing of Free Flying Wind Tunnel Models, Part 2: A Centerline Supported Fullspan Model Tested for Gust Load Alleviation." A team comprised of the Air Force Research Laboratory (AFRL), Northrop Grumman, Lockheed Martin, and the NASA Langley Research Center conducted three aeroservoelastic wind tunnel tests in the Transonic Dynamics Tunnel to demonstrate active control technologies relevant to large, flexible vehicles. In the first of these three tests, a semispan, aeroelastically scaled, wind tunnel model of a flying wing SensorCraft vehicle was mounted to a force balance to demonstrate gust load alleviation. In the second and third tests, the same wing was mated to a new, multi-degree of freedom, sidewall mount. This mount allowed the half-span model to translate vertically and pitch at the wing root, allowing better simulation of the full span vehicle's rigid body modes. Gust load alleviation (GLA) and Body freedom flutter (BFF) suppression were successfully demonstrated. The rigid body degrees-of-freedom required that the model be flown in the wind tunnel using an active control system. This risky mode of testing necessitated that a model arrestment system be integrated into the new mount. The safe and successful completion of these free flying tests required the development and integration of custom hardware and software. This paper describes the many systems, software, and procedures that were developed as part of this effort.

Simulated Wake Characteristics Data for Closely Spaced Parallel Runway Operations Analysis

NASA Technical Reports Server (NTRS)

Guerreiro, Nelson M.; Neitzke, Kurt W.

2012-01-01

A simulation experiment was performed to generate and compile wake characteristics data relevant to the evaluation and feasibility analysis of closely spaced parallel runway (CSPR) operational concepts. While the experiment in this work is not tailored to any particular operational concept, the generated data applies to the broader class of CSPR concepts, where a trailing aircraft on a CSPR approach is required to stay ahead of the wake vortices generated by a lead aircraft on an adjacent CSPR. Data for wake age, circulation strength, and wake altitude change, at various lateral offset distances from the wake-generating lead aircraft approach path were compiled for a set of nine aircraft spanning the full range of FAA and ICAO wake classifications. A total of 54 scenarios were simulated to generate data related to key parameters that determine wake behavior. Of particular interest are wake age characteristics that can be used to evaluate both time- and distance- based in-trail separation concepts for all aircraft wake-class combinations. A simple first-order difference model was developed to enable the computation of wake parameter estimates for aircraft models having weight, wingspan and speed characteristics similar to those of the nine aircraft modeled in this work.

An economical model for simulating droplet spectrum evolution in turbulent cloud chambers and wind tunnels

NASA Astrophysics Data System (ADS)

Krueger, Steven; Cantrell, W.; Niedermeier, D.; Shaw, R.; Stratmann, F.

2017-11-01

Although airborne instruments provide detailed information about the microphysical structure of clouds, the measurements provide only a few snapshots of each cloud. Deducing the droplet spectrum evolution from such measurements is next to impossible. We are using two alternative approaches: laboratory studies and numerical simulations. The former relies on a new turbulent cloud chamber (the Pi Chamber) at Michigan Technical University, as well as the first humid turbulent wind tunnel (LACIS-T) at the Leibniz Institute for Tropospheric Research. Both produce conditions for droplet growth (i.e., supersaturation) by mixing saturated vapor at different temperatures. The Pi Chamber produces turbulence by inducing Rayleigh-Bénard convection, while the wind tunnel generates turbulence with a grid. We are using the Explicit Mixing Parcel Model (EMPM) to numerically simulate droplet spectrum evolution in these flows. The EMPM explicitly links turbulent mixing and droplet spectrum evolution by representing a turbulent flow in a 1D domain with the linear eddy model. The EMPM can economically span scales from those of the smallest turbulent eddies to those of the largest. The EMPM grows or evaporates thousands of individual cloud droplets according to their local environments.

Developmental Regulation across the Life Span: Toward a New Synthesis

ERIC Educational Resources Information Center

Haase, Claudia M.; Heckhausen, Jutta; Wrosch, Carsten

2013-01-01

How can individuals regulate their own development to live happy, healthy, and productive lives? Major theories of developmental regulation across the life span have been proposed (e.g., dual-process model of assimilation and accommodation; motivational theory of life-span development; model of selection, optimization, and compensation), but they…

Streamflow response to future land-cover change at a headwaters catchment spanning the alpine-subalpine transition on the Colorado Front Range

NASA Astrophysics Data System (ADS)

Barnhart, T. B.; Vukomanovic, J.; Bourgeron, P.; Molotch, N. P.

2017-12-01

Land-cover change at the alpine-subalpine interface has the potential to change the water balance of mountainous, snow-dominated catchments due to the influence of vegetation on blowing snow, effective precipitation, evapotranspiration, and other processes. Understanding how land-cover change will impact water resources in snow-dominated regions is of critical importance as these locations produce a disproportionate amount of runoff relative to their land area. We coupled the LANdscape DIsturbance and Succession (LANDIS-II) model with a spatially explicit, physics-based, watershed process model, the Regional Hydro-Ecologic Simulation System (RHESSys), to simulate land-cover change and its impact on the water balance in a 6.6 km­2 headwater catchment that spans the alpine-subalpine transition on the Colorado Front Range. We simulated two potential futures of air temperature warming (+4 °C/century) to 2100: a) increased precipitation (+15%, MP) and b) decreased precipitation (-15%, LP). As the LANDIS-II model simulates forest succession in a stochastic manner, we use three LANDIS-II model runs each for the MP and LP future forcing conditions. For both MP and LP, the RHESSys forcing data set was updated to reflect the changes in precipitation and temperature used to generate the land-cover futures. Forest cover in the catchment increased from 72% in 2000 to 84% and 83% in 2050 and to 95% and 92% in 2100 for MP and LP, respectively. Somewhat surprisingly, this increase in forest cover led to mean increases in streamflow production of 9% for MP and 3% for LP in 2050. In 2100, mean streamflow production increased by 15% and 6% for the MP and LP scenarios, respectively. This is likely due to increases in effective precipitation as the catchment forested and blowing snow decreased. Indeed, catchment effective precipitation increased from 94% in 2000 to 97% and 99% in 2050 and 2100, respectively, for both MP and LP conditions. This result counters previous work as runoff production increased with forested area, highlighting the need to better understand the impacts of forest expansion on blowing snow and effective precipitation. Identifying the hydrologic response of mountainous areas to climate warming induced land-cover change is of critical importance due to the potential water resources impacts in downstream regions.

Modelling carbon responses of tundra ecosystems to historical and projected climate: A comparison of a plot- and a global-scale ecosystem model to identify process-based uncertainties

USGS Publications Warehouse

Clein, Joy S.; Kwiatkowski, B.L.; McGuire, A.D.; Hobbie, J.E.; Rastetter, E.B.; Melillo, J.M.; Kicklighter, D.W.

2000-01-01

We are developing a process-based modelling approach to investigate how carbon (C) storage of tundra across the entire Arctic will respond to projected climate change. To implement the approach, the processes that are least understood, and thus have the most uncertainty, need to be identified and studied. In this paper, we identified a key uncertainty by comparing the responses of C storage in tussock tundra at one site between the simulations of two models - one a global-scale ecosystem model (Terrestrial Ecosystem Model, TEM) and one a plot-scale ecosystem model (General Ecosystem Model, GEM). The simulations spanned the historical period (1921-94) and the projected period (1995-2100). In the historical period, the model simulations of net primary production (NPP) differed in their sensitivity to variability in climate. However, the long-term changes in C storage were similar in both simulations, because the dynamics of heterotrophic respiration (RH) were similar in both models. In contrast, the responses of C storage in the two model simulations diverged during the projected period. In the GEM simulation for this period, increases in RH tracked increases in NPP, whereas in the TEM simulation increases in RH lagged increases in NPP. We were able to make the long-term C dynamics of the two simulations agree by parameterizing TEM to the fast soil C pools of GEM. We concluded that the differences between the long-term C dynamics of the two simulations lay in modelling the role of the recalcitrant soil C. These differences, which reflect an incomplete understanding of soil processes, lead to quite different projections of the response of pan-Arctic C storage to global change. For example, the reference parameterization of TEM resulted in an estimate of cumulative C storage of 2032 g C m-2 for moist tundra north of 50??N, which was substantially higher than the 463 g C m-2 estimated for a parameterization of fast soil C dynamics. This uncertainty in the depiction of the role of recalcitrant soil C in long-term ecosystem C dynamics resulted from our incomplete understanding of controls over C and N transformations in Arctic soils. Mechanistic studies of these issues are needed to improve our ability to model the response of Arctic ecosystems to global change.

A multicriteria decision analysis model and risk assessment framework for carbon capture and storage.

PubMed

Humphries Choptiany, John Michael; Pelot, Ronald

2014-09-01

Multicriteria decision analysis (MCDA) has been applied to various energy problems to incorporate a variety of qualitative and quantitative criteria, usually spanning environmental, social, engineering, and economic fields. MCDA and associated methods such as life-cycle assessments and cost-benefit analysis can also include risk analysis to address uncertainties in criteria estimates. One technology now being assessed to help mitigate climate change is carbon capture and storage (CCS). CCS is a new process that captures CO2 emissions from fossil-fueled power plants and injects them into geological reservoirs for storage. It presents a unique challenge to decisionmakers (DMs) due to its technical complexity, range of environmental, social, and economic impacts, variety of stakeholders, and long time spans. The authors have developed a risk assessment model using a MCDA approach for CCS decisions such as selecting between CO2 storage locations and choosing among different mitigation actions for reducing risks. The model includes uncertainty measures for several factors, utility curve representations of all variables, Monte Carlo simulation, and sensitivity analysis. This article uses a CCS scenario example to demonstrate the development and application of the model based on data derived from published articles and publicly available sources. The model allows high-level DMs to better understand project risks and the tradeoffs inherent in modern, complex energy decisions. © 2014 Society for Risk Analysis.

Simulating Air-Entrapment in Low Permeability Mudrocks using a Macroscopic Invasion Percolation Model

NASA Astrophysics Data System (ADS)

Singh, A.; Holt, R. M.; Ramarao, B.; Clemo, T.

2011-12-01

Three radioactive waste disposal landfills at the Waste Control Specialists (WCS) facility in Andrews County, Texas are constructed below grade, within the low-permeability Dockum Group mudrocks (Cooper Canyon Formation) of Triassic age. Recent site investigations at the WCS disposal facilities indicate the presence of a trapped and compressed gas phase in the mudrocks. The Dockum is a low-permeability medium with vertical and horizontal effective hydraulic conductivities of 1.2E-9 cm/s and 2.9E-7 cm/s. The upper 300+ feet of the Dockum is in the unsaturated zone, with an average saturation of 0.87 and average capillary pressure of 2.8 MPa determined from core samples. Air entry pressures on core samples range from from 0.016 to 9.8 MPa, with a mean of 1.0 MPa. Heat dissipation sensors, thermocouple psychrometers, and advanced tensiometers installed in Dockum borehole arrays generally show capillary pressures one order of magnitude less than those measured on core samples. These differences with core data are attributed to the presence of a trapped and compressed gas phase within Dockum materials. In the vicinity of an instrumented borehole, the gas phase pressure equilibrates with atmospheric pressure, lowering the capillary pressure. We have developed a new macroscopic invasion percolation (MIP) model to illustrate the origin of the trapped gas phase in the Dockum rocks. An MIP model differs from invasion percolation (IP) through the definition of macro-scale capillarity. Individual pore throats and necks are not considered. Instead, a near pore-scale block is defined and characterized by a local threshold spanning pressure (a local block-scale breakthrough pressure) that represents the behavior of the subscale network. The model domain is discretized into an array of grid blocks with assigned spanning pressures. An invasion pressure for each block is then determined by the sum of spanning pressure, buoyance forces, and viscous forces. An IP algorithm sorts the invadable blocks, selects the block connected to the growing cluster with the lowest invasion pressure, and invades it. Our new MIP model incorporates several new features, including an efficient three-dimensional clustering algorithm; simultaneous invasion/reinvasion of water and air phases; hysteresis in water and air drainage curves; capability for distributed porosities and drainage parameters; and gas-phase compression and trapping. We apply this model in simulations representing the WCS site and illustrate the origin of the trapped and compressed gas phase in Dockum mudrocks.

Future intensification of hydro-meteorological extremes: downscaling using the weather research and forecasting model

NASA Astrophysics Data System (ADS)

El-Samra, R.; Bou-Zeid, E.; Bangalath, H. K.; Stenchikov, G.; El-Fadel, M.

2017-12-01

A set of ten downscaling simulations at high spatial resolution (3 km horizontally) were performed using the Weather Research and Forecasting (WRF) model to generate future climate projections of annual and seasonal temperature and precipitation changes over the Eastern Mediterranean (with a focus on Lebanon). The model was driven with the High Resolution Atmospheric Model (HiRAM), running over the whole globe at a resolution of 25 km, under the conditions of two Representative Concentration Pathways (RCP) (4.5 and 8.5). Each downscaling simulation spanned one year. Two past years (2003 and 2008), also forced by HiRAM without data assimilation, were simulated to evaluate the model's ability to capture the cold and wet (2003) and hot and dry (2008) extremes. The downscaled data were in the range of recent observed climatic variability, and therefore corrected for the cold bias of HiRAM. Eight future years were then selected based on an anomaly score that relies on the mean annual temperature and accumulated precipitation to identify the worst year per decade from a water resources perspective. One hot and dry year per decade, from 2011 to 2050, and per scenario was simulated and compared to the historic 2008 reference. The results indicate that hot and dry future extreme years will be exacerbated and the study area might be exposed to a significant decrease in annual precipitation (rain and snow), reaching up to 30% relative to the current extreme conditions.

Near-realtime simulations of biolelectric activity in small mammalian hearts using graphical processing units

PubMed Central

Vigmond, Edward J.; Boyle, Patrick M.; Leon, L. Joshua; Plank, Gernot

2014-01-01

Simulations of cardiac bioelectric phenomena remain a significant challenge despite continual advancements in computational machinery. Spanning large temporal and spatial ranges demands millions of nodes to accurately depict geometry, and a comparable number of timesteps to capture dynamics. This study explores a new hardware computing paradigm, the graphics processing unit (GPU), to accelerate cardiac models, and analyzes results in the context of simulating a small mammalian heart in real time. The ODEs associated with membrane ionic flow were computed on traditional CPU and compared to GPU performance, for one to four parallel processing units. The scalability of solving the PDE responsible for tissue coupling was examined on a cluster using up to 128 cores. Results indicate that the GPU implementation was between 9 and 17 times faster than the CPU implementation and scaled similarly. Solving the PDE was still 160 times slower than real time. PMID:19964295

A comparison of cylindrical and row trenched cooling holes with alignment angle of 0 degree near the combustor endwall

NASA Astrophysics Data System (ADS)

Kianpour, E.; Nor Azwadi, C. S.; Golshokouh, I.

2013-12-01

We studied the effects of cylindrical and row trenched cooling holes with alignment angle of 0° at BR=3.18 on the film cooling performance near the endwall surface of a combustor simulator. In this research, a three-dimensional presentation of gas turbine engine was simulated and analyzed with a commercial finite volume package FLUENT 6.2.26 to gain fundamental data. The current study has been performed with Reynolds-averaged Navier-Stokes turbulence model (RANS) on internal cooling passages. This combustor simulator combined the interaction of two rows of dilution jets, which were staggered in the stream wise direction and aligned in the span wise direction. The entire findings of the study declared that with using the row trenched holes near the enwall surface; film cooling effectiveness is doubled compared to the cooling performance of baseline case.

Whole Atmosphere Modeling and Data Analysis: Success Stories, Challenges and Perspectives

NASA Astrophysics Data System (ADS)

Yudin, V. A.; Akmaev, R. A.; Goncharenko, L. P.; Fuller-Rowell, T. J.; Matsuo, T.; Ortland, D. A.; Maute, A. I.; Solomon, S. C.; Smith, A. K.; Liu, H.; Wu, Q.

2015-12-01

At the end of the 20-th century Raymond Roble suggested an ambitious target of developing an atmospheric general circulation model (GCM) that spans from the surface to the thermosphere for modeling the coupled atmosphere-ionosphere with drivers from terrestrial meteorology and solar-geomagnetic inputs. He pointed out several areas of research and applications that would benefit highly from the development and improvement of whole atmosphere modeling. At present several research groups using middle and whole atmosphere models have attempted to perform coupled ionosphere-thermosphere predictions to interpret the "unexpected" anomalies in the electron content, ions and plasma drifts observed during recent stratospheric warming events. The recent whole atmosphere inter-comparison case studies also displayed striking differences in simulations of prevailing flows, planetary waves and dominant tidal modes even when the lower atmosphere domain of those models were constrained by similar meteorological analyses. We will present the possible reasons of such differences between data-constrained whole atmosphere simulations when analyses with 6-hour time resolution are used and discuss the potential model-data and model-model differences above the stratopause. The possible shortcomings of the whole atmosphere simulations associated with model physics, dynamical cores and resolutions will be discussed. With the increased confidence in the space-borne temperature, winds and ozone observations and extensive collections of ground-based upper atmosphere observational facilities, the whole atmosphere modelers will be able to quantify annual and year-to-variability of the zonal mean flows, planetary wave and tides. We will demonstrate the value of tidal and planetary wave variability deduced from the space-borne data and ground-based systems for evaluation and tune-up of whole atmosphere simulations including corrections of systematic model errors. Several success stories on the middle and whole atmosphere simulations coupled with the ionosphere models will be highlighted, and future perspectives for links of the space and terrestrial weather predictions constrained by current and scheduled ionosphere-thermosphere-mesosphere satellite missions will be presented

Revival of the fittest: exploding core-collapse supernovae from 12 to 25 M⊙

NASA Astrophysics Data System (ADS)

Vartanyan, David; Burrows, Adam; Radice, David; Skinner, M. Aaron; Dolence, Joshua

2018-07-01

We present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M⊙. Four of the models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si-O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, highlighting the Si-O interface as one key to explosion. Furthermore, we show that all of the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including reasonable changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of a few × 1050 erg at the end of our simulation, and are rising, emphasizing the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in nature. Morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. We do not find evidence for the Lepton-number Emission Self-sustained Asymmetry. Finally, we look at proto-neutron star (PNS) properties and explore the role of dimension in our simulations. We find that convection in the PNS produces larger PNS radii as well as greater `νμ' luminosities in 2D compared to 1D.

Revival of the Fittest: Exploding Core-Collapse Supernovae from 12 to 25 M⊙

NASA Astrophysics Data System (ADS)

Vartanyan, David; Burrows, Adam; Radice, David; Skinner, M. Aaron; Dolence, Joshua

2018-03-01

We present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M⊙. Four of the models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si-O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, highlighting the Si-O interface as one key to explosion. Furthermore, we show that all of the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including reasonable changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of few × 1050 ergs at the end of our simulation, and are rising, emphasizing the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in Nature. Morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. We do not find evidence for the Lepton-number Emission Self-Sustained Asymmetry. Finally, we look at PNS properties and explore the role of dimension in our simulations. We find that convection in the proto-neutron star (PNS) produces larger PNS radii as well as greater "νμ" luminosities in 2D compared to 1D.

Direct calculation of liquid-vapor phase equilibria from transition matrix Monte Carlo simulation

NASA Astrophysics Data System (ADS)

Errington, Jeffrey R.

2003-06-01

An approach for directly determining the liquid-vapor phase equilibrium of a model system at any temperature along the coexistence line is described. The method relies on transition matrix Monte Carlo ideas developed by Fitzgerald, Picard, and Silver [Europhys. Lett. 46, 282 (1999)]. During a Monte Carlo simulation attempted transitions between states along the Markov chain are monitored as opposed to tracking the number of times the chain visits a given state as is done in conventional simulations. Data collection is highly efficient and very precise results are obtained. The method is implemented in both the grand canonical and isothermal-isobaric ensemble. The main result from a simulation conducted at a given temperature is a density probability distribution for a range of densities that includes both liquid and vapor states. Vapor pressures and coexisting densities are calculated in a straightforward manner from the probability distribution. The approach is demonstrated with the Lennard-Jones fluid. Coexistence properties are directly calculated at temperatures spanning from the triple point to the critical point.

Development of metamodels for predicting aerosol dispersion in ventilated spaces

NASA Astrophysics Data System (ADS)

Hoque, Shamia; Farouk, Bakhtier; Haas, Charles N.

2011-04-01

Artificial neural network (ANN) based metamodels were developed to describe the relationship between the design variables and their effects on the dispersion of aerosols in a ventilated space. A Hammersley sequence sampling (HSS) technique was employed to efficiently explore the multi-parameter design space and to build numerical simulation scenarios. A detailed computational fluid dynamics (CFD) model was applied to simulate these scenarios. The results derived from the CFD simulations were used to train and test the metamodels. Feed forward ANN's were developed to map the relationship between the inputs and the outputs. The predictive ability of the neural network based metamodels was compared to linear and quadratic metamodels also derived from the same CFD simulation results. The ANN based metamodel performed well in predicting the independent data sets including data generated at the boundaries. Sensitivity analysis showed that particle tracking time to residence time and the location of input and output with relation to the height of the room had more impact than the other dimensionless groups on particle behavior.

Operations analysis (study 2.1): Program manual and users guide for the LOVES computer code

NASA Technical Reports Server (NTRS)

Wray, S. T., Jr.

1975-01-01

Information is provided necessary to use the LOVES Computer Program in its existing state, or to modify the program to include studies not properly handled by the basic model. The Users Guide defines the basic elements assembled together to form the model for servicing satellites in orbit. As the program is a simulation, the method of attack is to disassemble the problem into a sequence of events, each occurring instantaneously and each creating one or more other events in the future. The main driving force of the simulation is the deterministic launch schedule of satellites and the subsequent failure of the various modules which make up the satellites. The LOVES Computer Program uses a random number generator to simulate the failure of module elements and therefore operates over a long span of time typically 10 to 15 years. The sequence of events is varied by making several runs in succession with different random numbers resulting in a Monte Carlo technique to determine statistical parameters of minimum value, average value, and maximum value.

Experimental and Numerical Studies on the Formability of Materials in Hot Stamping and Cold Die Quenching Processes

NASA Astrophysics Data System (ADS)

Li, N.; Mohamed, M. S.; Cai, J.; Lin, J.; Balint, D.; Dean, T. A.

2011-05-01

Formability of steel and aluminium alloys in hot stamping and cold die quenching processes is studied in this research. Viscoplastic-damage constitutive equations are developed and determined from experimental data for the prediction of viscoplastic flow and ductility of the materials. The determined unified constitutive equations are then implemented into the commercial Finite Element code Abaqus/Explicit via a user defined subroutine, VUMAT. An FE process simulation model and numerical procedures are established for the modeling of hot stamping processes for a spherical part with a central hole. Different failure modes (failure takes place either near the central hole or in the mid span of the part) are obtained. To validate the simulation results, a test programme is developed, a test die set has been designed and manufactured, and tests have been carried out for the materials with different forming rates. It has been found that very close agreements between experimental and numerical process simulation results are obtained for the ranges of temperatures and forming rates carried out.

Experimental study of the effect on span loading on aircraft wakes

NASA Technical Reports Server (NTRS)

Corsiglia, V. R.; Rossow, V. J.; Ciffone, D. L.

1975-01-01

Measurements were made in the NASA-Ames 40- by 80-foot wind tunnel of the rolling moment induced on a following model in the wake 13.6 spans behind a subsonic transport model for a variety of trailing edge flap settings of the generator. It was found that the rolling moment on the following model was reduced substantially, compared to the conventional landing configuration, by reshaping the span loading on the generating model to approximate a span loading, found in earlier studies, which resulted in reduced wake velocities. This was accomplished by retracting the outboard trailing edge flaps. It was concluded, based on flow visualization conducted in the wind tunnel as well as in a water tow facility, that this flap arrangement redistributes the vorticity shed by the wing along the span to form three vortex pairs that interact to disperse the wake.

Evaluation of tropical channel refinement using MPAS-A aquaplanet simulations

DOE PAGES

Martini, Matus N.; Gustafson, Jr., William I.; O'Brien, Travis A.; ...

2015-09-13

Climate models with variable-resolution grids offer a computationally less expensive way to provide more detailed information at regional scales and increased accuracy for processes that cannot be resolved by a coarser grid. This study uses the Model for Prediction Across Scales–Atmosphere (MPAS22A), consisting of a nonhydrostatic dynamical core and a subset of Advanced Research Weather Research and Forecasting (ARW-WRF) model atmospheric physics that have been modified to include the Community Atmosphere Model version 5 (CAM5) cloud fraction parameterization, to investigate the potential benefits of using increased resolution in an tropical channel. The simulations are performed with an idealized aquaplanet configurationmore » using two quasi-uniform grids, with 30 km and 240 km grid spacing, and two variable-resolution grids spanning the same grid spacing range; one with a narrow (20°S–20°N) and one with a wide (30°S–30°N) tropical channel refinement. Results show that increasing resolution in the tropics impacts both the tropical and extratropical circulation. Compared to the quasi-uniform coarse grid, the narrow-channel simulation exhibits stronger updrafts in the Ferrel cell as well as in the middle of the upward branch of the Hadley cell. The wider tropical channel has a closer correspondence to the 30 km quasi-uniform simulation. However, the total atmospheric poleward energy transports are similar in all simulations. The largest differences are in the low-level cloudiness. The refined channel simulations show improved tropical and extratropical precipitation relative to the global 240 km simulation when compared to the global 30 km simulation. All simulations have a single ITCZ. Furthermore, the relatively small differences in mean global and tropical precipitation rates among the simulations are a promising result, and the evidence points to the tropical channel being an effective method for avoiding the extraneous numerical artifacts seen in earlier studies that only refined portion of the tropics.« less

Design and fabrication of brazed Rene 41 honeycomb sandwich structural panels for advanced space transportation systems

NASA Technical Reports Server (NTRS)

Hepler, A. K.; Swegle, A. R.

1981-01-01

The design and fabrication of two large brazed Rene 41 honeycomb panels, the establishment of a test plan, the design and fabrication of a test fixture to subject the panels to cyclic thermal gradients and mechanical loads equivalent to those imposed on an advanced space transportation vehicle during its boost and entry trajectories are discussed. The panels will be supported at four points, creating three spans. The outer spans are 45.7 cm (18 in.) and the center span 76.2 cm (30 in). Specimen width is 30.5 cm (12 in.). The panels were primarily designed by boost conditions simulated by subjecting the panels to liquid nitrogen, 77K (-320 F) on one side and 455K (360 F) on the other side and by mechanically imposing loads representing vehicle fuel pressure loads. Entry conditions were simulated by radiant heating to 1034K (1400 F). The test program subjected the panels to 500 boost thermal conditions. Results are presented.

Thermal stratification hinders gyrotactic micro-organism rising in free-surface turbulence

NASA Astrophysics Data System (ADS)

Lovecchio, Salvatore; Zonta, Francesco; Marchioli, Cristian; Soldati, Alfredo

2017-05-01

Thermal stratification in water bodies influences the exchange of heat, momentum, and chemical species across the air-water interface by modifying the sub-surface turbulence characteristics. Turbulence modifications may in turn prevent small motile algae (phytoplankton, in particular) from reaching the heated surface. We examine how different regimes of stable thermal stratification affect the motion of these microscopic organisms (modelled as gyrotactic self-propelling cells) in a free-surface turbulent channel flow. This archetypal setup mimics an environmentally plausible situation that can be found in lakes and oceans. Results from direct numerical simulations of turbulence coupled with Lagrangian tracking reveal that rising of bottom-heavy self-propelling cells depends strongly on the strength of stratification, especially near the thermocline where high temperature and velocity gradients occur: Here hydrodynamic shear may disrupt directional cell motility and hamper near-surface accumulation. For all gyrotactic re-orientation times considered in this study (spanning two orders of magnitude), we observe a reduction of the cell rising speed and temporary confinement under the thermocline: If re-orientation is fast, cells eventually trespass the thermocline within the simulated time span; if re-orientation is slow, confinement lasts much longer because cells align in the streamwise direction and their vertical swimming is practically annihilated.

Analysis of resonance mechanism and conditions of train bridge system

NASA Astrophysics Data System (ADS)

Xia, H.; Zhang, N.; Guo, W. W.

2006-11-01

In this paper, the resonance mechanism and conditions of train-bridge system are investigated through theoretical derivations, numerical simulations and experimental data analyses. The resonant responses of the bridge induced by moving trains are classified into three types according to different resonance mechanisms: the first is related to the periodical actions of moving load series of the vertical weights, lateral centrifugal and wind forces of vehicles; the second is induced by the loading rate of moving load series of vehicles; the third is owing to the periodically loading of the swing forces of the train vehicles excited by track irregularities and wheel hunting movements. The vehicle resonance is induced by the periodical action of regular arrangement of bridge spans and their deflections. The resonant conditions are proposed and the corresponding resonant train speeds are determined. The application scopes of resonance conditions are discussed. The resonance of the train-bridge system is affected by the span, total length, lateral and vertical stiffness of the bridge, the compositions of the train, and the axle arrangements and natural frequencies of the vehicles. The resonant train speeds for some bridges are estimated and are compared with the critical train speeds obtained from the dynamic simulation of train-bridge interaction model or from the field measurements.

Effects of Combined Stellar Feedback on Star Formation in Stellar Clusters

NASA Astrophysics Data System (ADS)

Wall, Joshua Edward; McMillan, Stephen; Pellegrino, Andrew; Mac Low, Mordecai; Klessen, Ralf; Portegies Zwart, Simon

2018-01-01

We present results of hybrid MHD+N-body simulations of star cluster formation and evolution including self consistent feedback from the stars in the form of radiation, winds, and supernovae from all stars more massive than 7 solar masses. The MHD is modeled with the adaptive mesh refinement code FLASH, while the N-body computations are done with a direct algorithm. Radiation is modeled using ray tracing along long characteristics in directions distributed using the HEALPIX algorithm, and causes ionization and momentum deposition, while winds and supernova conserve momentum and energy during injection. Stellar evolution is followed using power-law fits to evolution models in SeBa. We use a gravity bridge within the AMUSE framework to couple the N-body dynamics of the stars to the gas dynamics in FLASH. Feedback from the massive stars alters the structure of young clusters as gas ejection occurs. We diagnose this behavior by distinguishing between fractal distribution and central clustering using a Q parameter computed from the minimum spanning tree of each model cluster. Global effects of feedback in our simulations will also be discussed.

Statistics of Dark Matter Halos from Gravitational Lensing.

PubMed

Jain; Van Waerbeke L

2000-02-10

We present a new approach to measure the mass function of dark matter halos and to discriminate models with differing values of Omega through weak gravitational lensing. We measure the distribution of peaks from simulated lensing surveys and show that the lensing signal due to dark matter halos can be detected for a wide range of peak heights. Even when the signal-to-noise ratio is well below the limit for detection of individual halos, projected halo statistics can be constrained for halo masses spanning galactic to cluster halos. The use of peak statistics relies on an analytical model of the noise due to the intrinsic ellipticities of source galaxies. The noise model has been shown to accurately describe simulated data for a variety of input ellipticity distributions. We show that the measured peak distribution has distinct signatures of gravitational lensing, and its non-Gaussian shape can be used to distinguish models with different values of Omega. The use of peak statistics is complementary to the measurement of field statistics, such as the ellipticity correlation function, and is possibly not susceptible to the same systematic errors.

Critical end point in the presence of a chiral chemical potential

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

Cui, Z. -F.; Cloët, I. C.; Lu, Y.

A class of Polyakov-loop-modified Nambu-Jona-Lasinio models has been used to support a conjecture that numerical simulations of lattice-regularized QCD defined with a chiral chemical potential can provide information about the existence and location of a critical end point in the QCD phase diagram drawn in the plane spanned by baryon chemical potential and temperature. That conjecture is challenged by conflicts between the model results and analyses of the same problem using simulations of lattice-regularized QCD (lQCD) and well-constrained Dyson-Schwinger equation (DSE) studies. We find the conflict is resolved in favor of the lQCD and DSE predictions when both a physicallymore » motivated regularization is employed to suppress the contribution of high-momentum quark modes in the definition of the effective potential connected with the Polyakov-loop-modified Nambu-Jona-Lasinio models and the four-fermion coupling in those models does not react strongly to changes in the mean field that is assumed to mock-up Polyakov-loop dynamics. With the lQCD and DSE predictions thus confirmed, it seems unlikely that simulations of lQCD with mu(5) > 0 can shed any light on a critical end point in the regular QCD phase diagram.« less

QUIESCENT PROMINENCES IN THE ERA OF ALMA: SIMULATED OBSERVATIONS USING THE 3D WHOLE-PROMINENCE FINE STRUCTURE MODEL

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

Gunár, Stanislav; Heinzel, Petr; Mackay, Duncan H.

2016-12-20

We use the detailed 3D whole-prominence fine structure model to produce the first simulated high-resolution ALMA observations of a modeled quiescent solar prominence. The maps of synthetic brightness temperature and optical thickness shown in the present paper are produced using a visualization method for synthesis of the submillimeter/millimeter radio continua. We have obtained the simulated observations of both the prominence at the limb and the filament on the disk at wavelengths covering a broad range that encompasses the full potential of ALMA. We demonstrate here extent to which the small-scale and large-scale prominence and filament structures will be visible inmore » the ALMA observations spanning both the optically thin and thick regimes. We analyze the relationship between the brightness and kinetic temperature of the prominence plasma. We also illustrate the opportunities ALMA will provide for studying the thermal structure of the prominence plasma from the cores of the cool prominence fine structure to the prominence–corona transition region. In addition, we show that detailed 3D modeling of entire prominences with their numerous fine structures will be important for the correct interpretation of future ALMA observations of prominences.« less

Estimating daily time series of streamflow using hydrological model calibrated based on satellite observations of river water surface width: Toward real world applications.

PubMed

Sun, Wenchao; Ishidaira, Hiroshi; Bastola, Satish; Yu, Jingshan

2015-05-01

Lacking observation data for calibration constrains applications of hydrological models to estimate daily time series of streamflow. Recent improvements in remote sensing enable detection of river water-surface width from satellite observations, making possible the tracking of streamflow from space. In this study, a method calibrating hydrological models using river width derived from remote sensing is demonstrated through application to the ungauged Irrawaddy Basin in Myanmar. Generalized likelihood uncertainty estimation (GLUE) is selected as a tool for automatic calibration and uncertainty analysis. Of 50,000 randomly generated parameter sets, 997 are identified as behavioral, based on comparing model simulation with satellite observations. The uncertainty band of streamflow simulation can span most of 10-year average monthly observed streamflow for moderate and high flow conditions. Nash-Sutcliffe efficiency is 95.7% for the simulated streamflow at the 50% quantile. These results indicate that application to the target basin is generally successful. Beyond evaluating the method in a basin lacking streamflow data, difficulties and possible solutions for applications in the real world are addressed to promote future use of the proposed method in more ungauged basins. Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.

A numerical forecast model for road meteorology

NASA Astrophysics Data System (ADS)

Meng, Chunlei

2017-05-01

A fine-scale numerical model for road surface parameters prediction (BJ-ROME) is developed based on the Common Land Model. The model is validated using in situ observation data measured by the ROSA road weather stations of Vaisala Company, Finland. BJ-ROME not only takes into account road surface factors, such as imperviousness, relatively low albedo, high heat capacity, and high heat conductivity, but also considers the influence of urban anthropogenic heat, impervious surface evaporation, and urban land-use/land-cover changes. The forecast time span and the update interval of BJ-ROME in vocational operation are 24 and 3 h, respectively. The validation results indicate that BJ-ROME can successfully simulate the diurnal variation of road surface temperature both under clear-sky and rainfall conditions. BJ-ROME can simulate road water and snow depth well if the artificial removing was considered. Road surface energy balance in rainy days is quite different from that in clear-sky conditions. Road evaporation could not be neglected in road surface water cycle research. The results of sensitivity analysis show solar radiation correction coefficient, asphalt depth, and asphalt heat conductivity are important parameters in road interface temperatures simulation. The prediction results could be used as a reference of maintenance decision support system to mitigate the traffic jam and urban water logging especially in large cities.

Response of the Antarctic ice sheet to ocean forcing using the POPSICLES coupled ice sheet-ocean model

NASA Astrophysics Data System (ADS)

Martin, D. F.; Asay-Davis, X.; Price, S. F.; Cornford, S. L.; Maltrud, M. E.; Ng, E. G.; Collins, W.

2014-12-01

We present the response of the continental Antarctic ice sheet to sub-shelf-melt forcing derived from POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period 1990 to 2010. Simulations are performed at 0.1 degree (~5 km) ocean resolution and ice sheet resolution as fine as 500 m using adaptive mesh refinement. A comparison of fully-coupled and comparable standalone ice-sheet model results demonstrates the importance of two-way coupling between the ice sheet and the ocean. The POPSICLES model couples the POP2x ocean model, a modified version of the Parallel Ocean Program (Smith and Gent, 2002), and the BISICLES ice-sheet model (Cornford et al., 2012). BISICLES makes use of adaptive mesh refinement to fully resolve dynamically-important regions like grounding lines and employs a momentum balance similar to the vertically-integrated formulation of Schoof and Hindmarsh (2009). Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests like MISMIP3D (Pattyn et al., 2013) and realistic configurations (Favier et al. 2014). POP2x includes sub-ice-shelf circulation using partial top cells (Losch, 2008) and boundary layer physics following Holland and Jenkins (1999), Jenkins (2001), and Jenkins et al. (2010). Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP; Losch, 2008) and other continental-scale simulations and melt-rate observations (Kimura et al., 2013; Rignot et al., 2013). A companion presentation, "Present-day circum-Antarctic simulations using the POPSICLES coupled land ice-ocean model" in session C027 describes the ocean-model perspective of this work, while we focus on the response of the ice sheet and on details of the model. The figure shows the BISICLES-computed vertically-integrated ice velocity field about 1 month into a 20-year coupled Antarctic run. Groundling lines are shown in green.

Improved representation of plant functional types and physiology in the Joint UK Land Environment Simulator (JULES v4.2) using plant trait information

NASA Astrophysics Data System (ADS)

Harper, Anna B.; Cox, Peter M.; Friedlingstein, Pierre; Wiltshire, Andy J.; Jones, Chris D.; Sitch, Stephen; Mercado, Lina M.; Groenendijk, Margriet; Robertson, Eddy; Kattge, Jens; Bönisch, Gerhard; Atkin, Owen K.; Bahn, Michael; Cornelissen, Johannes; Niinemets, Ülo; Onipchenko, Vladimir; Peñuelas, Josep; Poorter, Lourens; Reich, Peter B.; Soudzilovskaia, Nadjeda A.; van Bodegom, Peter

2016-07-01

Dynamic global vegetation models are used to predict the response of vegetation to climate change. They are essential for planning ecosystem management, understanding carbon cycle-climate feedbacks, and evaluating the potential impacts of climate change on global ecosystems. JULES (the Joint UK Land Environment Simulator) represents terrestrial processes in the UK Hadley Centre family of models and in the first generation UK Earth System Model. Previously, JULES represented five plant functional types (PFTs): broadleaf trees, needle-leaf trees, C3 and C4 grasses, and shrubs. This study addresses three developments in JULES. First, trees and shrubs were split into deciduous and evergreen PFTs to better represent the range of leaf life spans and metabolic capacities that exists in nature. Second, we distinguished between temperate and tropical broadleaf evergreen trees. These first two changes result in a new set of nine PFTs: tropical and temperate broadleaf evergreen trees, broadleaf deciduous trees, needle-leaf evergreen and deciduous trees, C3 and C4 grasses, and evergreen and deciduous shrubs. Third, using data from the TRY database, we updated the relationship between leaf nitrogen and the maximum rate of carboxylation of Rubisco (Vcmax), and updated the leaf turnover and growth rates to include a trade-off between leaf life span and leaf mass per unit area.Overall, the simulation of gross and net primary productivity (GPP and NPP, respectively) is improved with the nine PFTs when compared to FLUXNET sites, a global GPP data set based on FLUXNET, and MODIS NPP. Compared to the standard five PFTs, the new nine PFTs simulate a higher GPP and NPP, with the exception of C3 grasses in cold environments and C4 grasses that were previously over-productive. On a biome scale, GPP is improved for all eight biomes evaluated and NPP is improved for most biomes - the exceptions being the tropical forests, savannahs, and extratropical mixed forests where simulated NPP is too high. With the new PFTs, the global present-day GPP and NPP are 128 and 62 Pg C year-1, respectively. We conclude that the inclusion of trait-based data and the evergreen/deciduous distinction has substantially improved productivity fluxes in JULES, in particular the representation of GPP. These developments increase the realism of JULES, enabling higher confidence in simulations of vegetation dynamics and carbon storage.

Development and Utilization of Regional Oceanic Modeling System (ROMS). Delicacy, Imprecision, and Uncertainty of Oceanic Simulations: An Investigation with the Regional Oceanic Modeling System (ROMS). Submesoscale Flows and Mixing in the Ocean Surface Layer Using the Regional Oceanic Modeling System (ROMS). Eddy Effects in General Circulation, Spanning Mean Currents, Mesoscale Eddies, and Topographic Generation, including Submesoscale Nests

DTIC Science & Technology

2012-09-30

unbalanced motions is likely to occur. Due to an rapidly expanding set of investigation on oceanic flows at submesoscales, it is increasingly clear...Uchiyama, E. M. Lane, J. M. Restrepo, & J. C. McWilliams, 2011: A vortex force analysis of the interaction of rip currents and gravity waves. J. Geophys...particular topographic features, the torque is pervasively positive (cyclonic) along the Stream, in opposition to the anticyclonic wind curl in the

A parallel algorithm for the initial screening of space debris collisions prediction using the SGP4/SDP4 models and GPU acceleration

NASA Astrophysics Data System (ADS)

Lin, Mingpei; Xu, Ming; Fu, Xiaoyu

2017-05-01

Currently, a tremendous amount of space debris in Earth's orbit imperils operational spacecraft. It is essential to undertake risk assessments of collisions and predict dangerous encounters in space. However, collision predictions for an enormous amount of space debris give rise to large-scale computations. In this paper, a parallel algorithm is established on the Compute Unified Device Architecture (CUDA) platform of NVIDIA Corporation for collision prediction. According to the parallel structure of NVIDIA graphics processors, a block decomposition strategy is adopted in the algorithm. Space debris is divided into batches, and the computation and data transfer operations of adjacent batches overlap. As a consequence, the latency to access shared memory during the entire computing process is significantly reduced, and a higher computing speed is reached. Theoretically, a simulation of collision prediction for space debris of any amount and for any time span can be executed. To verify this algorithm, a simulation example including 1382 pieces of debris, whose operational time scales vary from 1 min to 3 days, is conducted on Tesla C2075 of NVIDIA. The simulation results demonstrate that with the same computational accuracy as that of a CPU, the computing speed of the parallel algorithm on a GPU is 30 times that on a CPU. Based on this algorithm, collision prediction of over 150 Chinese spacecraft for a time span of 3 days can be completed in less than 3 h on a single computer, which meets the timeliness requirement of the initial screening task. Furthermore, the algorithm can be adapted for multiple tasks, including particle filtration, constellation design, and Monte-Carlo simulation of an orbital computation.

NMR spin-rotation relaxation and diffusion of methane

NASA Astrophysics Data System (ADS)

Singer, P. M.; Asthagiri, D.; Chapman, W. G.; Hirasaki, G. J.

2018-05-01

The translational diffusion-coefficient and the spin-rotation contribution to the 1H NMR relaxation rate for methane (CH4) are investigated using MD (molecular dynamics) simulations, over a wide range of densities and temperatures, spanning the liquid, supercritical, and gas phases. The simulated diffusion-coefficients agree well with measurements, without any adjustable parameters in the interpretation of the simulations. A minimization technique is developed to compute the angular velocity for non-rigid spherical molecules, which is used to simulate the autocorrelation function for spin-rotation interactions. With increasing diffusivity, the autocorrelation function shows increasing deviations from the single-exponential decay predicted by the Langevin theory for rigid spheres, and the deviations are quantified using inverse Laplace transforms. The 1H spin-rotation relaxation rate derived from the autocorrelation function using the "kinetic model" agrees well with measurements in the supercritical/gas phase, while the relaxation rate derived using the "diffusion model" agrees well with measurements in the liquid phase. 1H spin-rotation relaxation is shown to dominate over the MD-simulated 1H-1H dipole-dipole relaxation at high diffusivity, while the opposite is found at low diffusivity. At high diffusivity, the simulated spin-rotation correlation time agrees with the kinetic collision time for gases, which is used to derive a new expression for 1H spin-rotation relaxation, without any adjustable parameters.

The zoom lens of attention: Simulating shuffled versus normal text reading using the SWIFT model

PubMed Central

Schad, Daniel J.; Engbert, Ralf

2012-01-01

Assumptions on the allocation of attention during reading are crucial for theoretical models of eye guidance. The zoom lens model of attention postulates that attentional deployment can vary from a sharp focus to a broad window. The model is closely related to the foveal load hypothesis, i.e., the assumption that the perceptual span is modulated by the difficulty of the fixated word. However, these important theoretical concepts for cognitive research have not been tested quantitatively in eye movement models. Here we show that the zoom lens model, implemented in the SWIFT model of saccade generation, captures many important patterns of eye movements. We compared the model's performance to experimental data from normal and shuffled text reading. Our results demonstrate that the zoom lens of attention might be an important concept for eye movement control in reading. PMID:22754295

Electrophysiological models of neural processing.

PubMed

Nelson, Mark E

2011-01-01

The brain is an amazing information processing system that allows organisms to adaptively monitor and control complex dynamic interactions with their environment across multiple spatial and temporal scales. Mathematical modeling and computer simulation techniques have become essential tools in understanding diverse aspects of neural processing ranging from sub-millisecond temporal coding in the sound localization circuity of barn owls to long-term memory storage and retrieval in humans that can span decades. The processing capabilities of individual neurons lie at the core of these models, with the emphasis shifting upward and downward across different levels of biological organization depending on the nature of the questions being addressed. This review provides an introduction to the techniques for constructing biophysically based models of individual neurons and local networks. Topics include Hodgkin-Huxley-type models of macroscopic membrane currents, Markov models of individual ion-channel currents, compartmental models of neuronal morphology, and network models involving synaptic interactions among multiple neurons.

Remapping dark matter halo catalogues between cosmological simulations

NASA Astrophysics Data System (ADS)

Mead, A. J.; Peacock, J. A.

2014-05-01

We present and test a method for modifying the catalogue of dark matter haloes produced from a given cosmological simulation, so that it resembles the result of a simulation with an entirely different set of parameters. This extends the method of Angulo & White, which rescales the full particle distribution from a simulation. Working directly with the halo catalogue offers an advantage in speed, and also allows modifications of the internal structure of the haloes to account for non-linear differences between cosmologies. Our method can be used directly on a halo catalogue in a self-contained manner without any additional information about the overall density field; although the large-scale displacement field is required by the method, this can be inferred from the halo catalogue alone. We show proof of concept of our method by rescaling a matter-only simulation with no baryon acoustic oscillation (BAO) features to a more standard Λ cold dark matter model containing a cosmological constant and a BAO signal. In conjunction with the halo occupation approach, this method provides a basis for the rapid generation of mock galaxy samples spanning a wide range of cosmological parameters.

Three dimensional dynamics of a flexible Motorised Momentum Exchange Tether

NASA Astrophysics Data System (ADS)

Ismail, N. A.; Cartmell, M. P.

2016-03-01

This paper presents a new flexural model for the three dimensional dynamics of the Motorised Momentum Exchange Tether (MMET) concept. This study has uncovered the relationships between planar and nonplanar motions, and the effect of the coupling between these two parameters on pragmatic circular and elliptical orbits. The tether sub-spans are modelled as stiffened strings governed by partial differential equations of motion, with specific boundary conditions. The tether sub-spans are flexible and elastic, thereby allowing three dimensional displacements. The boundary conditions lead to a specific frequency equation and the eigenvalues from this provide the natural frequencies of the orbiting flexible motorised tether when static, accelerating in monotonic spin, and at terminal angular velocity. A rotation transformation matrix has been utilised to get the position vectors of the system's components in an assumed inertial frame. Spatio-temporal coordinates are transformed to modal coordinates before applying Lagrange's equations, and pre-selected linear modes are included to generate the equations of motion. The equations of motion contain inertial nonlinearities which are essentially of cubic order, and these show the potential for intricate intermodal coupling effects. A simulation of planar and non-planar motions has been undertaken and the differences in the modal responses, for both motions, and between the rigid body and flexible models are highlighted and discussed.

How long will my mouse live? Machine learning approaches for prediction of mouse life span.

PubMed

Swindell, William R; Harper, James M; Miller, Richard A

2008-09-01

Prediction of individual life span based on characteristics evaluated at middle-age represents a challenging objective for aging research. In this study, we used machine learning algorithms to construct models that predict life span in a stock of genetically heterogeneous mice. Life-span prediction accuracy of 22 algorithms was evaluated using a cross-validation approach, in which models were trained and tested with distinct subsets of data. Using a combination of body weight and T-cell subset measures evaluated before 2 years of age, we show that the life-span quartile to which an individual mouse belongs can be predicted with an accuracy of 35.3% (+/-0.10%). This result provides a new benchmark for the development of life-span-predictive models, but improvement can be expected through identification of new predictor variables and development of computational approaches. Future work in this direction can provide tools for aging research and will shed light on associations between phenotypic traits and longevity.

A wind turbine hybrid simulation framework considering aeroelastic effects

NASA Astrophysics Data System (ADS)

Song, Wei; Su, Weihua

2015-04-01

In performing an effective structural analysis for wind turbine, the simulation of turbine aerodynamic loads is of great importance. The interaction between the wake flow and the blades may impact turbine blades loading condition, energy yield and operational behavior. Direct experimental measurement of wind flow field and wind profiles around wind turbines is very helpful to support the wind turbine design. However, with the growth of the size of wind turbines for higher energy output, it is not convenient to obtain all the desired data in wind-tunnel and field tests. In this paper, firstly the modeling of dynamic responses of large-span wind turbine blades will consider nonlinear aeroelastic effects. A strain-based geometrically nonlinear beam formulation will be used for the basic structural dynamic modeling, which will be coupled with unsteady aerodynamic equations and rigid-body rotations of the rotor. Full wind turbines can be modeled by using the multi-connected beams. Then, a hybrid simulation experimental framework is proposed to potentially address this issue. The aerodynamic-dominant components, such as the turbine blades and rotor, are simulated as numerical components using the nonlinear aeroelastic model; while the turbine tower, where the collapse of failure may occur under high level of wind load, is simulated separately as the physical component. With the proposed framework, dynamic behavior of NREL's 5MW wind turbine blades will be studied and correlated with available numerical data. The current work will be the basis of the authors' further studies on flow control and hazard mitigation on wind turbine blades and towers.

Implementation of a Marauding Insect Module (MIM, version 1.0) in the Integrated BIosphere Simulator (IBIS, version 2.6b4) Dynamic Vegetation-Land Surface Model

NASA Astrophysics Data System (ADS)

Landry, J.-S.; Price, D. T.; Ramankutty, N.; Parrott, L.; Matthews, H. D.

2015-12-01

Insects defoliate and kill plants in many ecosystems worldwide. The consequences of these natural processes on terrestrial ecology and nutrient cycling are well established, and their potential climatic effects resulting from modified land-atmosphere exchanges of carbon, energy, and water are increasingly being recognized. We developed a Marauding Insect Module (MIM) to quantify, in the Integrated BIosphere Simulator (IBIS), the consequences of insect activity on biogeochemical and biogeophysical fluxes, also accounting for the effects of altered vegetation dynamics. MIM can simulate damage from broadleaf defoliators, needleleaf defoliators, and bark beetles, with the resulting impacts being estimated by IBIS based on the new, insect-modified state of the vegetation. MIM further accounts for the physical presence and gradual fall of insect-killed dead standing trees. The design of MIM should facilitate the addition of other insect types besides the ones already included and could guide the development of similar modules for other process-based vegetation models. After describing IBIS-MIM, we illustrate the usefulness of the model by presenting results spanning daily to centennial timescales for vegetation dynamics and cycling of carbon, energy, and water following a simulated outbreak of the mountain pine beetle. We then show that these simulated impacts agree with many previous studies based on field measurements, satellite data, or modelling. MIM and similar tools should therefore be of great value in assessing the wide array of impacts resulting from insect-induced plant damage in the Earth system.

Comparison of Cirrus Cloud Models: A Project of the GEWEX Cloud System Study (GCSS) Working Group on Cirrus Cloud Systems

NASA Technical Reports Server (NTRS)

Starr, David O'C.; Benedetti, Angela; Boehm, Matt; Brown, Philip R. A.; Gierens, Klaus M.; Girard, Eric; Giraud, Vincent; Jakob, Christian; Jensen, Eric

2000-01-01

The GEWEX Cloud System Study (GCSS, GEWEX is the Global Energy and Water Cycle Experiment) is a community activity aiming to promote development of improved cloud parameterizations for application in the large-scale general circulation models (GCMs) used for climate research and for numerical weather prediction. The GCSS strategy is founded upon the use of cloud-system models (CSMs). These are "process" models with sufficient spatial and temporal resolution to represent individual cloud elements, but spanning a wide range of space and time scales to enable statistical analysis of simulated cloud systems. GCSS also employs single-column versions of the parametric cloud models (SCMs) used in GCMs. GCSS has working groups on boundary-layer clouds, cirrus clouds, extratropical layer cloud systems, precipitating deep convective cloud systems, and polar clouds.

Continuum kinetic modeling of the tokamak plasma edge

DOE PAGES

Dorf, M. A.; Dorr, M.; Rognlien, T.; ...

2016-03-10

In this study, the first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasmatransport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalousmore » radial transport.« less

VLF Radio Wave Propagation Across the Day/Night Terminator

NASA Astrophysics Data System (ADS)

Burch, H.; Moore, R. C.

2016-12-01

In May 2016, a new array of VLF radio receivers was deployed spanning the East Coast of the United States. We present preliminary observations from the array, which was designed in part to track the propagation of the narrowband VLF transmitter signal, NAA (24.0 kHz), down the coast from Cutler, Maine. Amplitude, phase, and polarization observations are compared over multiple days and at different times of year to investigate the dependence of VLF propagation characteristics on solar zenith angle. Measurements are compared to simulations using the Long Wave Propagation Capability code (LWPC) in order to evaluate the accuracy of LWPC's built-in ionosphere model. Efforts to improve the ionosphere model based on observations are discussed.

Modeling the Influence of Hemispheric Transport on Trends in ...

EPA Pesticide Factsheets

We describe the development and application of the hemispheric version of the CMAQ to examine the influence of long-range pollutant transport on trends in surface level O3 distributions. The WRF-CMAQ model is expanded to hemispheric scales and multi-decadal model simulations were recently performed for the period spanning 1990-2010 to examine changes in hemispheric air pollution resulting from changes in emissions over this period. Simulated trends in ozone and precursor species concentrations across the U.S. and the northern hemisphere over the past two decades are compared with those inferred from available measurements during this period. Additionally, the decoupled direct method (DDM) in CMAQ is used to estimate the sensitivity of O3 to emissions from different source regions across the northern hemisphere. The seasonal variations in source region contributions to background O3 is then estimated from these sensitivity calculations and will be discussed. A reduced form model combining these source region sensitivities estimated from DDM with the multi-decadal simulations of O3 distributions and emissions trends, is then developed to characterize the changing contributions of different source regions to background O3 levels across North America. The National Exposure Research Laboratory (NERL) Computational Exposure Division (CED) develops and evaluates data, decision-support tools, and models to be applied to media-specific or receptor-specific problem areas

Influence of disturbance on temperate forest productivity

USGS Publications Warehouse

Peters, Emily B.; Wythers, Kirk R.; Bradford, John B.; Reich, Peter B.

2013-01-01

Climate, tree species traits, and soil fertility are key controls on forest productivity. However, in most forest ecosystems, natural and human disturbances, such as wind throw, fire, and harvest, can also exert important and lasting direct and indirect influence over productivity. We used an ecosystem model, PnET-CN, to examine how disturbance type, intensity, and frequency influence net primary production (NPP) across a range of forest types from Minnesota and Wisconsin, USA. We assessed the importance of past disturbances on NPP, net N mineralization, foliar N, and leaf area index at 107 forest stands of differing types (aspen, jack pine, northern hardwood, black spruce) and disturbance history (fire, harvest) by comparing model simulations with observations. The model reasonably predicted differences among forest types in productivity, foliar N, leaf area index, and net N mineralization. Model simulations that included past disturbances minimally improved predictions compared to simulations without disturbance, suggesting the legacy of past disturbances played a minor role in influencing current forest productivity rates. Modeled NPP was more sensitive to the intensity of soil removal during a disturbance than the fraction of stand mortality or wood removal. Increasing crown fire frequency resulted in lower NPP, particularly for conifer forest types with longer leaf life spans and longer recovery times. These findings suggest that, over long time periods, moderate frequency disturbances are a relatively less important control on productivity than climate, soil, and species traits.

From single fiber to macro-level mechanics: A structural finite-element model for elastomeric fibrous biomaterials.

PubMed

D'Amore, Antonio; Amoroso, Nicholas; Gottardi, Riccardo; Hobson, Christopher; Carruthers, Christopher; Watkins, Simon; Wagner, William R; Sacks, Michael S

2014-11-01

In the present work, we demonstrate that the mesoscopic in-plane mechanical behavior of membrane elastomeric scaffolds can be simulated by replication of actual quantified fibrous geometries. Elastomeric electrospun polyurethane (ES-PEUU) scaffolds, with and without particulate inclusions, were utilized. Simulations were developed from experimentally-derived fiber network geometries, based on a range of scaffold isotropic and anisotropic behaviors. These were chosen to evaluate the effects on macro-mechanics based on measurable geometric parameters such as fiber intersections, connectivity, orientation, and diameter. Simulations were conducted with only the fiber material model parameters adjusted to match the macro-level mechanical test data. Fiber model validation was performed at the microscopic level by individual fiber mechanical tests using AFM. Results demonstrated very good agreement to the experimental data, and revealed the formation of extended preferential fiber orientations spanning the entire model space. We speculate that these emergent structures may be responsible for the tissue-like macroscale behaviors observed in electrospun scaffolds. To conclude, the modeling approach has implications for (1) gaining insight on the intricate relationship between fabrication variables, structure, and mechanics to manufacture more functional devices/materials, (2) elucidating the effects of cell or particulate inclusions on global construct mechanics, and (3) fabricating better performing tissue surrogates that could recapitulate native tissue mechanics. Copyright © 2014 Elsevier Ltd. All rights reserved.

The Abacus Cosmos: A Suite of Cosmological N-body Simulations

NASA Astrophysics Data System (ADS)

Garrison, Lehman H.; Eisenstein, Daniel J.; Ferrer, Douglas; Tinker, Jeremy L.; Pinto, Philip A.; Weinberg, David H.

2018-06-01

We present a public data release of halo catalogs from a suite of 125 cosmological N-body simulations from the ABACUS project. The simulations span 40 wCDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, 4 × 1010 h ‑1 M ⊙ and 1 × 1010 h ‑1 M ⊙, in 1.1 h ‑1 Gpc and 720 h ‑1 Mpc boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning z = 1.5 to 0.1 are available for friends-of-friends and ROCKSTAR halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos.

Stability and dynamics of membrane-spanning DNA nanopores

NASA Astrophysics Data System (ADS)

Maingi, Vishal; Burns, Jonathan R.; Uusitalo, Jaakko J.; Howorka, Stefan; Marrink, Siewert J.; Sansom, Mark S. P.

2017-03-01

Recently developed DNA-based analogues of membrane proteins have advanced synthetic biology. A fundamental question is how hydrophilic nanostructures reside in the hydrophobic environment of the membrane. Here, we use multiscale molecular dynamics (MD) simulations to explore the structure, stability and dynamics of an archetypical DNA nanotube inserted via a ring of membrane anchors into a phospholipid bilayer. Coarse-grained MD reveals that the lipids reorganize locally to interact closely with the membrane-spanning section of the DNA tube. Steered simulations along the bilayer normal establish the metastable nature of the inserted pore, yielding a force profile with barriers for membrane exit due to the membrane anchors. Atomistic, equilibrium simulations at two salt concentrations confirm the close packing of lipid around of the stably inserted DNA pore and its cation selectivity, while revealing localized structural fluctuations. The wide-ranging and detailed insight informs the design of next-generation DNA pores for synthetic biology or biomedicine.

Exploring the effect of East Antarctic ice mass loss on GIA-induced horizontal bedrock motions

NASA Astrophysics Data System (ADS)

Konfal, S. A.; Whitehouse, P. L.; Hermans, T.; van der Wal, W.; Wilson, T. J.; Bevis, M. G.; Kendrick, E. C.; Dalziel, I.; Smalley, R., Jr.

2017-12-01

Ice history inputs used in Antarctic models of GIA include major centers of ice mass loss in West Antarctica. In the Transantarctic Mountains (TAM) region spanning the boundary between East and West Antarctica, horizontal crustal motions derived from GPS observations from the Antarctic Network (ANET) component of the Polar Earth Observing Network (POLENET) are towards these West Antarctic ice mass centers, opposite to the pattern of radial crustal motion expected in an unloading scenario. We investigate alternative ice history and earth structure inputs to GIA models in an attempt to reproduce observed crustal motions in the region. The W12 ice history model is altered to create scenarios including ice unloading in the Wilkes Subglacial Basin based on available glaciological records. These altered ice history models, along with the unmodified W12 ice history model, are coupled with 60 radially varying (1D) earth model combinations, including approximations of optimal earth profiles identified in published GIA models. The resulting model-predicted motions utilizing both the modified and unmodified ice history models fit ANET GPS-derived crustal motions in the northern TAM region for a suite of earth model combinations. Further south, where the influence of simulated Wilkes unloading is weakest and West Antarctic unloading is strongest, observed and predicted motions do not agree. The influence of simulated Wilkes ice unloading coupled with laterally heterogeneous earth models is also investigated. The resulting model-predicted motions do not differ significantly between the original W12 and W12 with simulated Wilkes unloading ice histories.

Ion Transport through Membrane-Spanning Nanopores Studied by Molecular Dynamics Simulations and Continuum Electrostatics Calculations

PubMed Central

Peter, Christine; Hummer, Gerhard

2005-01-01

Narrow hydrophobic regions are a common feature of biological channels, with possible roles in ion-channel gating. We study the principles that govern ion transport through narrow hydrophobic membrane pores by molecular dynamics simulation of model membranes formed of hexagonally packed carbon nanotubes. We focus on the factors that determine the energetics of ion translocation through such nonpolar nanopores and compare the resulting free-energy barriers for pores with different diameters corresponding to the gating regions in closed and open forms of potassium channels. Our model system also allows us to compare the results from molecular dynamics simulations directly to continuum electrostatics calculations. Both simulations and continuum calculations show that subnanometer wide pores pose a huge free-energy barrier for ions, but a small increase in the pore diameter to ∼1 nm nearly eliminates that barrier. We also find that in those wider channels the ion mobility is comparable to that in the bulk phase. By calculating local electrostatic potentials, we show that the long range Coulomb interactions of ions are strongly screened in the wide water-filled channels. Whereas continuum calculations capture the overall energetics reasonably well, the local water structure, which is not accounted for in this model, leads to interesting effects such as the preference of hydrated ions to move along the pore wall rather than through the center of the pore. PMID:16006629

Spatial and temporal agreement in climate model simulations of the Interdecadal Pacific Oscillation

DOE PAGES

Henley, Benjamin J.; Meehl, Gerald; Power, Scott B.; ...

2017-01-31

Accelerated warming and hiatus periods in the long-term rise of Global Mean Surface Temperature (GMST) have, in recent decades, been associated with the Interdecadal Pacific Oscillation (IPO). Critically, decadal climate prediction relies on the skill of state-of-the-art climate models to reliably represent these low-frequency climate variations. We undertake a systematic evaluation of the simulation of the IPO in the suite of Coupled Model Intercomparison Project 5 (CMIP5) models. We track the IPO in pre-industrial (control) and all-forcings (historical) experiments using the IPO tripole index (TPI). The TPI is explicitly aligned with the observed spatial pattern of the IPO, and circumventsmore » assumptions about the nature of global warming. We find that many models underestimate the ratio of decadal-to-total variance in sea surface temperatures (SSTs). However, the basin-wide spatial pattern of positive and negative phases of the IPO are simulated reasonably well, with spatial pattern correlation coefficients between observations and models spanning the range 0.4–0.8. Deficiencies are mainly in the extratropical Pacific. Models that better capture the spatial pattern of the IPO also tend to more realistically simulate the ratio of decadal to total variance. Of the 13% of model centuries that have a fractional bias in the decadal-to-total TPI variance of 0.2 or less, 84% also have a spatial pattern correlation coefficient with the observed pattern exceeding 0.5. This result is highly consistent across both IPO positive and negative phases. This is evidence that the IPO is related to one or more inherent dynamical mechanisms of the climate system.« less

Spatial and temporal agreement in climate model simulations of the Interdecadal Pacific Oscillation

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

Henley, Benjamin J.; Meehl, Gerald; Power, Scott B.

Accelerated warming and hiatus periods in the long-term rise of Global Mean Surface Temperature (GMST) have, in recent decades, been associated with the Interdecadal Pacific Oscillation (IPO). Critically, decadal climate prediction relies on the skill of state-of-the-art climate models to reliably represent these low-frequency climate variations. We undertake a systematic evaluation of the simulation of the IPO in the suite of Coupled Model Intercomparison Project 5 (CMIP5) models. We track the IPO in pre-industrial (control) and all-forcings (historical) experiments using the IPO tripole index (TPI). The TPI is explicitly aligned with the observed spatial pattern of the IPO, and circumventsmore » assumptions about the nature of global warming. We find that many models underestimate the ratio of decadal-to-total variance in sea surface temperatures (SSTs). However, the basin-wide spatial pattern of positive and negative phases of the IPO are simulated reasonably well, with spatial pattern correlation coefficients between observations and models spanning the range 0.4–0.8. Deficiencies are mainly in the extratropical Pacific. Models that better capture the spatial pattern of the IPO also tend to more realistically simulate the ratio of decadal to total variance. Of the 13% of model centuries that have a fractional bias in the decadal-to-total TPI variance of 0.2 or less, 84% also have a spatial pattern correlation coefficient with the observed pattern exceeding 0.5. This result is highly consistent across both IPO positive and negative phases. This is evidence that the IPO is related to one or more inherent dynamical mechanisms of the climate system.« less

Seamless atmospheric modeling across the hydrostatic-nonhydrostatic scales - preliminary results using an unstructured-Voronoi mesh for weather prediction.

NASA Astrophysics Data System (ADS)

Skamarock, W. C.

2015-12-01

One of the major problems in atmospheric model applications is the representation of deep convection within the models; explicit simulation of deep convection on fine meshes performs much better than sub-grid parameterized deep convection on coarse meshes. Unfortunately, the high cost of explicit convective simulation has meant it has only been used to down-scale global simulations in weather prediction and regional climate applications, typically using traditional one-way interactive nesting technology. We have been performing real-time weather forecast tests using a global non-hydrostatic atmospheric model (the Model for Prediction Across Scales, MPAS) that employs a variable-resolution unstructured Voronoi horizontal mesh (nominally hexagons) to span hydrostatic to nonhydrostatic scales. The smoothly varying Voronoi mesh eliminates many downscaling problems encountered using traditional one- or two-way grid nesting. Our test weather forecasts cover two periods - the 2015 Spring Forecast Experiment conducted at the NOAA Storm Prediction Center during the month of May in which we used a 50-3 km mesh, and the PECAN field program examining nocturnal convection over the US during the months of June and July in which we used a 15-3 km mesh. An important aspect of this modeling system is that the model physics be scale-aware, particularly the deep convection parameterization. These MPAS simulations employ the Grell-Freitas scale-aware convection scheme. Our test forecasts show that the scheme produces a gradual transition in the deep convection, from the deep unstable convection being handled entirely by the convection scheme on the coarse mesh regions (dx > 15 km), to the deep convection being almost entirely explicit on the 3 km NA region of the meshes. We will present results illustrating the performance of critical aspects of the MPAS model in these tests.

Ground-water hydrology, historical water use, and simulated ground-water flow in Cretaceous-age Coastal Plain aquifers near Charleston and Florence, South Carolina

USGS Publications Warehouse

Campbell, B.G.; van Heeswijk, Marijke

1996-01-01

A quasi-three-dimensional, transient, digital, ground-water flow model representing the Coastal Plain aquifers of South Carolina, has been constructed to assist in defining the ground- water-flow system of Cretaceous aquifers near Charleston and Florence, S.C. Both cities are near the centers of large (greater than 150 feet) potentiometric declines in the Middendorf aquifer. In 1989, the diameter of the depressions was approximately 40 miles at Charleston and 15 miles at Florence. The potentiometric decline occurred between predevelopment (1926) and 1982 near Florence, and between predevelopment (1879) and 1989 near Charleston. The city of Charleston does not withdraw water from these aquifers; however, some of the small communities in the area use these aquifers for a potable water supply. The model simulates flow in and between four aquifer systems. The model has a variable-cell-size grid, and spans the Coastal Plain from the Savannah River in the southwest to the Cape Fear Arch in the northeast, and from the Fall Line in the northwest to approximately 30 miles offshore to the southeast. Model-grid cell size is 1 by 1 mile in a 48 by 48 mile area centered in Charleston, and in a 36 by 48 mile area centered in Florence. The model cell size gradually increases to a maximum of 4 by 4 miles outside the two study areas. The entire grid consists of 115 by 127 cells and covers an area of 39,936 square miles. The model was calibrated to historical water-level data. The calibration relied on three techniques: (1) matching simulated and observed potentiometric map surfaces, (2) statistical comparison of observed and simulated heads, and (3) comparison of observed and simulated well hydrographs. Systematic changes in model parameters showed that simulated heads are most sensitive to changes in aquifer transmissivity. Eight predictive ground-water-use scenarios were simulated for the Mount Pleasant area, which presently (1993) uses the Middendorf aquifer as a sole-source of potable water. These simulations use various combinations of spatial distribution, and injection of treated wastewater effluent for existing and future Middendorf aquifer wells.

Present-day Circum-Antarctic Simulations using the POPSICLES Coupled Ice Sheet-Ocean Model

NASA Astrophysics Data System (ADS)

Asay-Davis, X.; Martin, D. F.; Price, S. F.; Maltrud, M. E.; Collins, W.

2014-12-01

We present POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period 1990 to 2010. Simulations are performed at 0.1o (~5 km) ocean resolution and with adaptive ice-sheet model resolution as fine as 500 m. We compare time-averaged melt rates below a number of major ice shelves with those reported by Rignot et al. (2013) as well as other recent studies. We also present seasonal variability and decadal trends in submarine melting from several Antarctic regions. Finally, we explore the influence on basal melting and system dynamics resulting from two different choices of climate forcing: a "normal-year" climatology and the CORE v. 2 forcing data (Large and Yeager 2008).POPSICLES couples the POP2x ocean model, a modified version of the Parallel Ocean Program (Smith and Gent, 2002), and the BISICLES ice-sheet model (Cornford et al., 2012). POP2x includes sub-ice-shelf circulation using partial top cells (Losch, 2008) and boundary layer physics following Holland and Jenkins (1999), Jenkins (2001), and Jenkins et al. (2010). Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP; Losch, 2008) and other continental-scale simulations and melt-rate observations (Kimura et al., 2013; Rignot et al., 2013). BISICLES makes use of adaptive mesh refinement and a 1st-order accurate momentum balance similar to the L1L2 model of Schoof and Hindmarsh (2009) to accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests (MISMIP-3D; Pattyn et al., 2013) and realistic configurations (Favier et al. 2014).A companion presentation, "Response of the Antarctic Ice Sheet to ocean forcing using the POPSICLES coupled ice sheet-ocean model" in session C024 covers the ice-sheet response to these melt rates in the coupled simulation. The figure shows eddy activity in the vertically integrated (barotropic) velocity nearly six years into a POPSICLES simulation of the Antarctic region.

Low-speed impact phenomena and orbital resonances in the moon- and planet-building process

NASA Technical Reports Server (NTRS)

Chapman, C. R.

1977-01-01

A simulation of collisional and gravitational interaction in the early solar system generates planets approximately 1000 km in diameter from an initial swarm of kilometer sized planetesimals. The model treats collisions according to experimental and theoretical impact results (such as rebound, cratering, and catastrophic fragmentation) for a variety of materials whose parameters span plausible values for early solid objects. The small planets form in approximately 1000 yr, during which time most of the mass of the system continues to reside in particles near the original size. The simulation is terminated when the largest objects' random motion is of smaller dimension than their collision cross-sections. The few 1000 km planets may act as seeds for the subsequent, gradual, accretional growth into full-sized planets.

Facilitated sequence counting and assembly by template mutagenesis

PubMed Central

Levy, Dan; Wigler, Michael

2014-01-01

Presently, inferring the long-range structure of the DNA templates is limited by short read lengths. Accurate template counts suffer from distortions occurring during PCR amplification. We explore the utility of introducing random mutations in identical or nearly identical templates to create distinguishable patterns that are inherited during subsequent copying. We simulate the applications of this process under assumptions of error-free sequencing and perfect mapping, using cytosine deamination as a model for mutation. The simulations demonstrate that within readily achievable conditions of nucleotide conversion and sequence coverage, we can accurately count the number of otherwise identical molecules as well as connect variants separated by long spans of identical sequence. We discuss many potential applications, such as transcript profiling, isoform assembly, haplotype phasing, and de novo genome assembly. PMID:25313059

Quartz c-axis orientation patterns in fracture cement as a measure of fracture opening rate and a validation tool for fracture pattern models

DOE PAGES

Ukar, Estibalitz; Laubach, Stephen E.; Marrett, Randall

2016-03-09

Here, we evaluate a published model for crystal growth patterns in quartz cement in sandstone fractures by comparing crystal fracture-spanning predictions to quartz c-axis orientation distributions measured by electron backscatter diffraction (EBSD) of spanning quartz deposits. Samples from eight subvertical opening-mode fractures in four sandstone formations, the Jurassic– Cretaceous Nikanassin Formation, northwestern Alberta Foothills (Canada), Cretaceous Mesaverde Group (USA; Cozzette Sandstone Member of the Iles Formation), Piceance Basin, Colorado (USA), and upper Jurassic–lower Cretaceous Cotton Valley Group (Taylor sandstone) and overlying Travis Peak Formation, east Texas, have similar quartzose composition and grain size but contain fractures with different temperature historiesmore » and opening rates based on fluid inclusion assemblages and burial history. Spherical statistical analysis shows that, in agreement with model predictions, bridging crystals have a preferred orientation with c-axis orientations at a high angle to fracture walls. The second form of validation is for spanning potential that depends on the size of cut substrate grains. Using measured cut substrate grain sizes and c-axis orientations of spanning bridges, we calculated the required orientation for the smallest cut grain to span the maximum gap size and the required orientation of the crystal with the least spanning potential to form overgrowths that span across maximum measured gap sizes. We find that within a 10° error all spanning crystals conform to model predictions. Using crystals with the lowest spanning potential based on crystallographic orientation (c-axis parallel to fracture wall) and a temperature range for fracture opening measured from fluid inclusion assemblages, we calculate maximum fracture opening rates that allow crystals to span. These rates are comparable to those derived independently from fracture temperature histories based on burial history and multiple sequential fluid inclusion assemblages. Results support the R. Lander and S. Laubach model, which predicts that for quartz deposited synchronously with fracture opening, spanning potential, or likelihood of quartz deposits that are thick enough to span between fracture walls, depends on temperature history, fracture opening rate, size of opening increments, and size, mineralogy, and crystallographic orientation of substrates in the fracture wall (transected grains). Results suggest that EBSD maps, which can be more rapidly acquired than measurement of tens to hundreds of fluid inclusion assemblages, can provide a useful measure of relative opening rates within populations of quartz-filled fractures formed under sedimentary basin conditions. Such data are useful for evaluating fracture pattern development models.« less

Quartz c-axis orientation patterns in fracture cement as a measure of fracture opening rate and a validation tool for fracture pattern models

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

Ukar, Estibalitz; Laubach, Stephen E.; Marrett, Randall

Here, we evaluate a published model for crystal growth patterns in quartz cement in sandstone fractures by comparing crystal fracture-spanning predictions to quartz c-axis orientation distributions measured by electron backscatter diffraction (EBSD) of spanning quartz deposits. Samples from eight subvertical opening-mode fractures in four sandstone formations, the Jurassic– Cretaceous Nikanassin Formation, northwestern Alberta Foothills (Canada), Cretaceous Mesaverde Group (USA; Cozzette Sandstone Member of the Iles Formation), Piceance Basin, Colorado (USA), and upper Jurassic–lower Cretaceous Cotton Valley Group (Taylor sandstone) and overlying Travis Peak Formation, east Texas, have similar quartzose composition and grain size but contain fractures with different temperature historiesmore » and opening rates based on fluid inclusion assemblages and burial history. Spherical statistical analysis shows that, in agreement with model predictions, bridging crystals have a preferred orientation with c-axis orientations at a high angle to fracture walls. The second form of validation is for spanning potential that depends on the size of cut substrate grains. Using measured cut substrate grain sizes and c-axis orientations of spanning bridges, we calculated the required orientation for the smallest cut grain to span the maximum gap size and the required orientation of the crystal with the least spanning potential to form overgrowths that span across maximum measured gap sizes. We find that within a 10° error all spanning crystals conform to model predictions. Using crystals with the lowest spanning potential based on crystallographic orientation (c-axis parallel to fracture wall) and a temperature range for fracture opening measured from fluid inclusion assemblages, we calculate maximum fracture opening rates that allow crystals to span. These rates are comparable to those derived independently from fracture temperature histories based on burial history and multiple sequential fluid inclusion assemblages. Results support the R. Lander and S. Laubach model, which predicts that for quartz deposited synchronously with fracture opening, spanning potential, or likelihood of quartz deposits that are thick enough to span between fracture walls, depends on temperature history, fracture opening rate, size of opening increments, and size, mineralogy, and crystallographic orientation of substrates in the fracture wall (transected grains). Results suggest that EBSD maps, which can be more rapidly acquired than measurement of tens to hundreds of fluid inclusion assemblages, can provide a useful measure of relative opening rates within populations of quartz-filled fractures formed under sedimentary basin conditions. Such data are useful for evaluating fracture pattern development models.« less

Active vibration suppression of self-excited structures using an adaptive LMS algorithm

NASA Astrophysics Data System (ADS)

Danda Roy, Indranil

The purpose of this investigation is to study the feasibility of an adaptive feedforward controller for active flutter suppression in representative linear wing models. The ability of the controller to suppress limit-cycle oscillations in wing models having root springs with freeplay nonlinearities has also been studied. For the purposes of numerical simulation, mathematical models of a rigid and a flexible wing structure have been developed. The rigid wing model is represented by a simple three-degree-of-freedom airfoil while the flexible wing is modelled by a multi-degree-of-freedom finite element representation with beam elements for bending and rod elements for torsion. Control action is provided by one or more flaps attached to the trailing edge and extending along the entire wing span for the rigid model and a fraction of the wing span for the flexible model. Both two-dimensional quasi-steady aerodynamics and time-domain unsteady aerodynamics have been used to generate the airforces in the wing models. An adaptive feedforward controller has been designed based on the filtered-X Least Mean Squares (LMS) algorithm. The control configuration for the rigid wing model is single-input single-output (SISO) while both SISO and multi-input multi-output (MIMO) configurations have been applied on the flexible wing model. The controller includes an on-line adaptive system identification scheme which provides the LMS controller with a reasonably accurate model of the plant. This enables the adaptive controller to track time-varying parameters in the plant and provide effective control. The wing models in closed-loop exhibit highly damped responses at airspeeds where the open-loop responses are destructive. Simulations with the rigid and the flexible wing models in a time-varying airstream show a 63% and 53% increase, respectively, over their corresponding open-loop flutter airspeeds. The ability of the LMS controller to suppress wing store flutter in the two models has also been investigated. With 10% measurement noise introduced in the flexible wing model, the controller demonstrated good robustness to the extraneous disturbances. In the examples studied it is found that adaptation is rapid enough to successfully control flutter at accelerations in the airstream of up to 15 ft/sec2 for the rigid wing model and 9 ft/sec2 for the flexible wing model.

Accessing the capability of TRMM 3B42 V7 to simulate streamflow during extreme rain events: Case study for a Himalayan River Basin

NASA Astrophysics Data System (ADS)

Kumar, Brijesh; Lakshmi, Venkat

2018-03-01

The paper examines the quality of Tropical Rainfall Monitoring Mission (TRMM) 3B42 V7 precipitation product to simulate the streamflow using Soil Water Assessment Tool (SWAT) model for various rainfall intensities over the Himalayan region. The SWAT model has been set up for Gandak River Basin with 41 sub-basins and 420 HRUs. Five stream gauge locations are used to simulate the streamflow for a time span of 10 years (2000-2010). Daily streamflow for the simulation period is collected from Central Water Commission (CWC), India and Department of Hydrology and Meteorology (DHM), Nepal. The simulation results are found good in terms of Nash-Sutcliffe efficiency (NSE) {>}0.65, coefficient of determination (R2) {>}0.67 and Percentage Bias (PBIAS){<}15%, at each stream gauge sites. Thereafter, we have calculated the PBIAS and RMSE-observations standard deviation ratio (RSR) statistics between TRMM simulated and observed streamflow for various rainfall intensity classes, viz., light ({<}7.5 mm/d), moderate (7.5 to 35.4 mm/d), heavy (35.5 to 124.4 mm/d) and extremely heavy ({>}124.4 mm/d). The PBIAS and RSR show that TRMM simulated streamflow is suitable for moderate to heavy rainfall intensities. However, it does not perform well for light- and extremely-heavy rainfall intensities. The finding of the present work is useful for the problems related to water resources management, irrigation planning and hazard analysis over the Himalayan regions.

“Modeling Trends in Air Pollutant Concentrations over the ...

EPA Pesticide Factsheets

Regional model calculations over annual cycles have pointed to the need for accurately representing impacts of long-range transport. Linking regional and global scale models have met with mixed success as biases in the global model can propagate and influence regional calculations and often confound interpretation of model results. Since transport is efficient in the free-troposphere and since simulations over Continental scales and annual cycles provide sufficient opportunity for “atmospheric turn-over”, i.e., exchange between the free-troposphere and the boundary-layer, a conceptual framework is needed wherein interactions between processes occurring at various spatial and temporal scales can be consistently examined. The coupled WRF-CMAQ model is expanded to hemispheric scales and model simulations over period spanning 1990-current are analyzed to examine changes in hemispheric air pollution resulting from changes in emissions over this period. The National Exposure Research Laboratory (NERL) Atmospheric Modeling and Analysis Division (AMAD) conducts research in support of EPA mission to protect human health and the environment. AMAD research program is engaged in developing and evaluating predictive atmospheric models on all spatial and temporal scales for forecasting the air quality and for assessing changes in air quality and air pollutant exposures, as affected by changes in ecosystem management and regulatory decisions. AMAD is responsible for pr

Quantum Mechanics/Molecular Mechanics Method Combined with Hybrid All-Atom and Coarse-Grained Model: Theory and Application on Redox Potential Calculations.

PubMed

Shen, Lin; Yang, Weitao

2016-04-12

We developed a new multiresolution method that spans three levels of resolution with quantum mechanical, atomistic molecular mechanical, and coarse-grained models. The resolution-adapted all-atom and coarse-grained water model, in which an all-atom structural description of the entire system is maintained during the simulations, is combined with the ab initio quantum mechanics and molecular mechanics method. We apply this model to calculate the redox potentials of the aqueous ruthenium and iron complexes by using the fractional number of electrons approach and thermodynamic integration simulations. The redox potentials are recovered in excellent accordance with the experimental data. The speed-up of the hybrid all-atom and coarse-grained water model renders it computationally more attractive. The accuracy depends on the hybrid all-atom and coarse-grained water model used in the combined quantum mechanical and molecular mechanical method. We have used another multiresolution model, in which an atomic-level layer of water molecules around redox center is solvated in supramolecular coarse-grained waters for the redox potential calculations. Compared with the experimental data, this alternative multilayer model leads to less accurate results when used with the coarse-grained polarizable MARTINI water or big multipole water model for the coarse-grained layer.

GWVis: A Tool for Comparative Ground-Water Data Visualization

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

Best, Daniel M.; Lewis, Robert R.

2010-11-01

The Ground-Water Visualization application (GWVis) presents ground-water data visually in order to educate the public on ground-water issues. It is also intended for presentations to government and other funding agencies. Current three dimensional models of ground-water are overly complex, while the two dimensional representations (i.e., on paper) are neither comprehensive, nor engaging. At present, GWVis operates on water head elevation data over a given time span, together with a matching (fixed) underlying geography. Two elevation scenarios are compared with each other, typically a control data set (actual field data) and a simulation. Scenario comparison can be animated for the timemore » span provided. We developed GWVis using the Python programming language, associated libraries, and pyOpenGL extension packages to improve performance and control of attributes of the mode (such as color, positioning, scale, and interpolation). GWVis bridges the gap between two dimensional and dynamic three dimensional research visualizations by providing an intuitive, interactive design that allows participants to view the model from different perspectives and to infer information about scenarios. By incorporating scientific data in an environment that can be easily understood, GWVis allows the information to be presented to a large audience base.« less

Historical (1850-2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols:methodology and application

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

Lamarque, J. F.; Bond, Tami C.; Eyring, Veronika

2010-08-11

We present and discuss a new dataset of gridded emissions covering the historical period (1850-2000) in decadal increments at a horizontal resolution of 0.5° in latitude and longitude. The primary purpose of this inventory is to provide consistent gridded emissions of reactive gases and aerosols for use in chemistry model simulations needed by climate models for the Climate Model Intercomparison Program #5 (CMIP5) in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment report. Our best estimate for the year 2000 inventory represents a combination of existing regional and global inventories to capture the best information available atmore » this point; 40 regions and 12 sectors were used to combine the various sources. The historical reconstruction of each emitted compound, for each region and sector, was then forced to agree with our 2000 estimate, ensuring continuity between past and 2000 emissions. Application of these emissions into two chemistry-climate models is used to test their ability to capture long-term changes in atmospheric ozone, carbon monoxide and aerosols distributions. The simulated long-term change in the Northern mid-latitudes surface and mid-troposphere ozone is not quite as rapid as observed. However, stations outside this latitude band show much better agreement in both present-day and long-term trend. The model simulations consistently underestimate the carbon monoxide trend, while capturing the long-term trend at the Mace Head station. The simulated sulfate and black carbon deposition over Greenland is in very good agreement with the ice-core observations spanning the simulation period. Finally, aerosol optical depth and additional aerosol diagnostics are shown to be in good agreement with previously published estimates.« less

Simplex-method based transmission performance optimization for 100G PDM-QPSK systems with non-identical spans

NASA Astrophysics Data System (ADS)

Li, Yuanyuan; Gao, Guanjun; Zhang, Jie; Zhang, Kai; Chen, Sai; Yu, Xiaosong; Gu, Wanyi

2015-06-01

A simplex-method based optimizing (SMO) strategy is proposed to improve the transmission performance for dispersion uncompensated (DU) coherent optical systems with non-identical spans. Through analytical expression of quality of transmission (QoT), this strategy improves the Q factors effectively, while minimizing the number of erbium-doped optical fiber amplifier (EDFA) that needs to be optimized. Numerical simulations are performed for 100 Gb/s polarization-division multiplexed quadrature phase shift keying (PDM-QPSK) channels over 10-span standard single mode fiber (SSMF) with randomly distributed span-lengths. Compared to the EDFA configurations with complete span loss compensation, the Q factor of the SMO strategy is improved by approximately 1 dB at the optimal transmitter launch power. Moreover, instead of adjusting the gains of all the EDFAs to their optimal value, the number of EDFA that needs to be adjusted for SMO is reduced from 8 to 2, showing much less tuning costs and almost negligible performance degradation.

Finite element modelling of FRC beams containing PVA and Basalt fibres: A comparative study

NASA Astrophysics Data System (ADS)

Ayub, Tehmina; Khan, Sadaqat Ullah

2017-09-01

The endeavour of current study is to compare the flexural behaviour and three dimensional (3D) finite element analysis (FEA) and the results of FEM are compared with the experimental results of 07 HPFRC beams. Out of seven (07), 01 beam of plain concrete without fibres was cast as a control beam. Three (03) beams containing 1, 2 and 3% volume of PVA fibres were prepared by using HPFRC mixes while, the remaining other three (03) beams were prepared using HPFRC mixes containing 1, 2 and 3% volume of Basalt fibres. In order to ensure flexural failure, three-point bending load was applied at the mid span of all beams. The maximum flexural load and corresponding deflection and strains at the mid span attained prior to the failure were obtained as flexural test results. The FEM results were obtained by simulating all beams in ATENA 3D program and verified through flexural test results. Both of the results of FEM and Experiment showed good agreement with each other.

A Comparison of Computation Span and Reading Span Working Memory Measures' Relations With Problem-Solving Criteria.

PubMed

Perlow, Richard; Jattuso, Mia

2018-06-01

Researchers have operationalized working memory in different ways and although working memory-performance relationships are well documented, there has been relatively less attention devoted to determining whether seemingly similar measures yield comparable relations with performance outcomes. Our objective is to assess whether two working memory measures deploying the same processes but different item content yield different relations with two problem-solving criteria. Participants completed a computation-based working memory measure and a reading-based measure prior to performing a computerized simulation. Results reveal differential relations with one of the two criteria and support the notion that the two working memory measures tap working memory capacity and other cognitive abilities. One implication for theory development is that researchers should consider incorporating other cognitive abilities in their working memory models and that the selection of those abilities should correspond to the criterion of interest. One practical implication is that researchers and practitioners shouldn't automatically assume that different phonological loop-based working memory scales are interchangeable.

Modeling the evolution and distribution of the frequency's second derivative and the braking index of pulsar spin

NASA Astrophysics Data System (ADS)

Xie, Yi; Zhang, Shuang-Nan; Liao, Jin-Yuan

2015-07-01

We model the evolution of the spin frequency's second derivative v̈ and the braking index n of radio pulsars with simulations within the phenomenological model of their surface magnetic field evolution, which contains a long-term power-law decay modulated by short-term oscillations. For the pulsar PSR B0329+54, a model with three oscillation components can reproduce its v̈ variation. We show that the “averaged” n is different from the instantaneous n, and its oscillation magnitude decreases abruptly as the time span increases, due to the “averaging” effect. The simulated timing residuals agree with the main features of the reported data. Our model predicts that the averaged v̈ of PSR B0329+54 will start to decrease rapidly with newer data beyond those used in Hobbs et al. We further perform Monte Carlo simulations for the distribution of the reported data in |v̈| and |n| versus characteristic age τC diagrams. It is found that the magnetic field oscillation model with decay index α = 0 can reproduce the distributions quite well. Compared with magnetic field decay due to the ambipolar diffusion (α = 0.5) and the Hall cascade (α = 1.0), the model with no long term decay (α = 0) is clearly preferred for old pulsars by the p-values of the two-dimensional Kolmogorov-Smirnov test. Supported by the National Natural Science Foundation of China.

Calibration of an Unsteady Groundwater Flow Model for a Complex, Strongly Heterogeneous Aquifer

NASA Astrophysics Data System (ADS)

Curtis, Z. K.; Liao, H.; Li, S. G.; Phanikumar, M. S.; Lusch, D.

2016-12-01

Modeling of groundwater systems characterized by complex three-dimensional structure and heterogeneity remains a significant challenge. Most of today's groundwater models are developed based on relatively simple conceptual representations in favor of model calibratibility. As more complexities are modeled, e.g., by adding more layers and/or zones, or introducing transient processes, more parameters have to be estimated and issues related to ill-posed groundwater problems and non-unique calibration arise. Here, we explore the use of an alternative conceptual representation for groundwater modeling that is fully three-dimensional and can capture complex 3D heterogeneity (both systematic and "random") without over-parameterizing the aquifer system. In particular, we apply Transition Probability (TP) geostatistics on high resolution borehole data from a water well database to characterize the complex 3D geology. Different aquifer material classes, e.g., `AQ' (aquifer material), `MAQ' (marginal aquifer material'), `PCM' (partially confining material), and `CM' (confining material), are simulated, with the hydraulic properties of each material type as tuning parameters during calibration. The TP-based approach is applied to simulate unsteady groundwater flow in a large, complex, and strongly heterogeneous glacial aquifer system in Michigan across multiple spatial and temporal scales. The resulting model is calibrated to observed static water level data over a time span of 50 years. The results show that the TP-based conceptualization enables much more accurate and robust calibration/simulation than that based on conventional deterministic layer/zone based conceptual representations.

A Motivational Theory of Life-Span Development

PubMed Central

Heckhausen, Jutta; Wrosch, Carsten; Schulz, Richard

2010-01-01

This article had four goals. First, the authors identified a set of general challenges and questions that a life-span theory of development should address. Second, they presented a comprehensive account of their Motivational Theory of Life-Span Development. They integrated the model of optimization in primary and secondary control and the action-phase model of developmental regulation with their original life-span theory of control to present a comprehensive theory of development. Third, they reviewed the relevant empirical literature testing key propositions of the Motivational Theory of Life-Span Development. Finally, because the conceptual reach of their theory goes far beyond the current empirical base, they pointed out areas that deserve further and more focused empirical inquiry. PMID:20063963

A mixing timescale model for TPDF simulations of turbulent premixed flames

DOE PAGES

Kuron, Michael; Ren, Zhuyin; Hawkes, Evatt R.; ...

2017-02-06

Transported probability density function (TPDF) methods are an attractive modeling approach for turbulent flames as chemical reactions appear in closed form. However, molecular micro-mixing needs to be modeled and this modeling is considered a primary challenge for TPDF methods. In the present study, a new algebraic mixing rate model for TPDF simulations of turbulent premixed flames is proposed, which is a key ingredient in commonly used molecular mixing models. The new model aims to properly account for the transition in reactive scalar mixing rate behavior from the limit of turbulence-dominated mixing to molecular mixing behavior in flamelets. An a priorimore » assessment of the new model is performed using direct numerical simulation (DNS) data of a lean premixed hydrogen–air jet flame. The new model accurately captures the mixing timescale behavior in the DNS and is found to be a significant improvement over the commonly used constant mechanical-to-scalar mixing timescale ratio model. An a posteriori TPDF study is then performed using the same DNS data as a numerical test bed. The DNS provides the initial conditions and time-varying input quantities, including the mean velocity, turbulent diffusion coefficient, and modeled scalar mixing rate for the TPDF simulations, thus allowing an exclusive focus on the mixing model. Here, the new mixing timescale model is compared with the constant mechanical-to-scalar mixing timescale ratio coupled with the Euclidean Minimum Spanning Tree (EMST) mixing model, as well as a laminar flamelet closure. It is found that the laminar flamelet closure is unable to properly capture the mixing behavior in the thin reaction zones regime while the constant mechanical-to-scalar mixing timescale model under-predicts the flame speed. Furthermore, the EMST model coupled with the new mixing timescale model provides the best prediction of the flame structure and flame propagation among the models tested, as the dynamics of reactive scalar mixing across different flame regimes are appropriately accounted for.« less

A mixing timescale model for TPDF simulations of turbulent premixed flames

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

Kuron, Michael; Ren, Zhuyin; Hawkes, Evatt R.

Transported probability density function (TPDF) methods are an attractive modeling approach for turbulent flames as chemical reactions appear in closed form. However, molecular micro-mixing needs to be modeled and this modeling is considered a primary challenge for TPDF methods. In the present study, a new algebraic mixing rate model for TPDF simulations of turbulent premixed flames is proposed, which is a key ingredient in commonly used molecular mixing models. The new model aims to properly account for the transition in reactive scalar mixing rate behavior from the limit of turbulence-dominated mixing to molecular mixing behavior in flamelets. An a priorimore » assessment of the new model is performed using direct numerical simulation (DNS) data of a lean premixed hydrogen–air jet flame. The new model accurately captures the mixing timescale behavior in the DNS and is found to be a significant improvement over the commonly used constant mechanical-to-scalar mixing timescale ratio model. An a posteriori TPDF study is then performed using the same DNS data as a numerical test bed. The DNS provides the initial conditions and time-varying input quantities, including the mean velocity, turbulent diffusion coefficient, and modeled scalar mixing rate for the TPDF simulations, thus allowing an exclusive focus on the mixing model. Here, the new mixing timescale model is compared with the constant mechanical-to-scalar mixing timescale ratio coupled with the Euclidean Minimum Spanning Tree (EMST) mixing model, as well as a laminar flamelet closure. It is found that the laminar flamelet closure is unable to properly capture the mixing behavior in the thin reaction zones regime while the constant mechanical-to-scalar mixing timescale model under-predicts the flame speed. Furthermore, the EMST model coupled with the new mixing timescale model provides the best prediction of the flame structure and flame propagation among the models tested, as the dynamics of reactive scalar mixing across different flame regimes are appropriately accounted for.« less

SIDM on FIRE: hydrodynamical self-interacting dark matter simulations of low-mass dwarf galaxies

NASA Astrophysics Data System (ADS)

Robles, Victor H.; Bullock, James S.; Elbert, Oliver D.; Fitts, Alex; González-Samaniego, Alejandro; Boylan-Kolchin, Michael; Hopkins, Philip F.; Faucher-Giguère, Claude-André; Kereš, Dušan; Hayward, Christopher C.

2017-12-01

We compare a suite of four simulated dwarf galaxies formed in 1010 M⊙ haloes of collisionless cold dark matter (CDM) with galaxies simulated in the same haloes with an identical galaxy formation model but a non-zero cross-section for DM self-interactions. These cosmological zoom-in simulations are part of the Feedback In Realistic Environments (FIRE) project and utilize the FIRE-2 model for hydrodynamics and galaxy formation physics. We find the stellar masses of the galaxies formed in self-interacting dark matter (SIDM) with σ/m = 1 cm2 g-1 are very similar to those in CDM (spanning M⋆ ≈ 105.7-7.0M⊙) and all runs lie on a similar stellar mass-size relation. The logarithmic DM density slope (α = d log ρ/d log r) in the central 250-500 pc remains steeper than α = -0.8 for the CDM-Hydro simulations with stellar mass M⋆ ∼ 106.6 M⊙ and core-like in the most massive galaxy. In contrast, every SIDM hydrodynamic simulation yields a flatter profile, with α > -0.4. Moreover, the central density profiles predicted in SIDM runs without baryons are similar to the SIDM runs that include FIRE-2 baryonic physics. Thus, SIDM appears to be much more robust to the inclusion of (potentially uncertain) baryonic physics than CDM on this mass scale, suggesting that SIDM will be easier to falsify than CDM using low-mass galaxies. Our FIRE simulations predict that galaxies less massive than M⋆ ≲ 3 × 106 M⊙ provide potentially ideal targets for discriminating models, with SIDM producing substantial cores in such tiny galaxies and CDM producing cusps.

Validation of a mixture-averaged thermal diffusion model for premixed lean hydrogen flames

NASA Astrophysics Data System (ADS)

Schlup, Jason; Blanquart, Guillaume

2018-03-01

The mixture-averaged thermal diffusion model originally proposed by Chapman and Cowling is validated using multiple flame configurations. Simulations using detailed hydrogen chemistry are done on one-, two-, and three-dimensional flames. The analysis spans flat and stretched, steady and unsteady, and laminar and turbulent flames. Quantitative and qualitative results using the thermal diffusion model compare very well with the more complex multicomponent diffusion model. Comparisons are made using flame speeds, surface areas, species profiles, and chemical source terms. Once validated, this model is applied to three-dimensional laminar and turbulent flames. For these cases, thermal diffusion causes an increase in the propagation speed of the flames as well as increased product chemical source terms in regions of high positive curvature. The results illustrate the necessity for including thermal diffusion, and the accuracy and computational efficiency of the mixture-averaged thermal diffusion model.

Decadal-Scale Response of the Antarctic Ice sheet to a Warming Ocean using the POPSICLES Coupled Ice Sheet-Ocean model

NASA Astrophysics Data System (ADS)

Martin, D. F.; Asay-Davis, X.; Cornford, S. L.; Price, S. F.; Ng, E. G.; Collins, W.

2015-12-01

We present POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period from 1990 to 2010. We use the CORE v. 2 interannual forcing data to force the ocean model. Simulations are performed at 0.1o(~5 km) ocean resolution with adaptive ice sheet resolution as fine as 500 m to adequately resolve the grounding line dynamics. We discuss the effect of improved ocean mixing and subshelf bathymetry (vs. the standard Bedmap2 bathymetry) on the behavior of the coupled system, comparing time-averaged melt rates below a number of major ice shelves with those reported in the literature. We also present seasonal variability and decadal melting trends from several Antarctic regions, along with the response of the ice shelves and the consequent dynamic response of the grounded ice sheet.POPSICLES couples the POP2x ocean model, a modified version of the Parallel Ocean Program, and the BISICLES ice-sheet model. POP2x includes sub-ice-shelf circulation using partial top cells and the commonly used three-equation boundary layer physics. Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP) and other continental-scale simulations and melt-rate observations. BISICLES makes use of adaptive mesh refinement and a 1st-order accurate momentum balance similar to the L1L2 model of Schoof and Hindmarsh to accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests (MISMIP-3d) and realistic configurations.The figure shows the BISICLES-computed vertically-integrated grounded ice velocity field 5 years into a 20-year coupled full-continent Antarctic-Southern-Ocean simulation. Submarine melt rates are painted onto the surface of the floating ice shelves. Grounding lines are shown in green.

Revival of the fittest: exploding core-collapse supernovae from 12 to 25 M⊙

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

Vartanyan, David; Burrows, Adam; Radice, David

Here, we present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M⊙. Four of the models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si–O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, highlighting the Si–O interface as one key to explosion. Furthermore, we show that all ofmore » the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including reasonable changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of a few × 10 50 erg at the end of our simulation, and are rising, emphasizing the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in nature. Morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. We do not find evidence for the Lepton-number Emission Self-sustained Asymmetry. Finally, we look at proto-neutron star (PNS) properties and explore the role of dimension in our simulations. We find that convection in the PNS produces larger PNS radii as well as greater ‘ν μ’ luminosities in 2D compared to 1D.« less

Revival of the fittest: exploding core-collapse supernovae from 12 to 25 M⊙

DOE PAGES

Vartanyan, David; Burrows, Adam; Radice, David; ...

2018-03-28

Here, we present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M⊙. Four of the models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si–O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, highlighting the Si–O interface as one key to explosion. Furthermore, we show that all ofmore » the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including reasonable changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of a few × 10 50 erg at the end of our simulation, and are rising, emphasizing the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in nature. Morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. We do not find evidence for the Lepton-number Emission Self-sustained Asymmetry. Finally, we look at proto-neutron star (PNS) properties and explore the role of dimension in our simulations. We find that convection in the PNS produces larger PNS radii as well as greater ‘ν μ’ luminosities in 2D compared to 1D.« less

Uncertainty Analysis of Runoff Simulations and Parameter Identifiability in the Community Land Model – Evidence from MOPEX Basins

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

Huang, Maoyi; Hou, Zhangshuan; Leung, Lai-Yung R.

2013-12-01

With the emergence of earth system models as important tools for understanding and predicting climate change and implications to mitigation and adaptation, it has become increasingly important to assess the fidelity of the land component within earth system models to capture realistic hydrological processes and their response to the changing climate and quantify the associated uncertainties. This study investigates the sensitivity of runoff simulations to major hydrologic parameters in version 4 of the Community Land Model (CLM4) by integrating CLM4 with a stochastic exploratory sensitivity analysis framework at 20 selected watersheds from the Model Parameter Estimation Experiment (MOPEX) spanning amore » wide range of climate and site conditions. We found that for runoff simulations, the most significant parameters are those related to the subsurface runoff parameterizations. Soil texture related parameters and surface runoff parameters are of secondary significance. Moreover, climate and soil conditions play important roles in the parameter sensitivity. In general, site conditions within water-limited hydrologic regimes and with finer soil texture result in stronger sensitivity of output variables, such as runoff and its surface and subsurface components, to the input parameters in CLM4. This study demonstrated the feasibility of parameter inversion for CLM4 using streamflow observations to improve runoff simulations. By ranking the significance of the input parameters, we showed that the parameter set dimensionality could be reduced for CLM4 parameter calibration under different hydrologic and climatic regimes so that the inverse problem is less ill posed.« less

Estimation of river and stream temperature trends under haphazard sampling

USGS Publications Warehouse

Gray, Brian R.; Lyubchich, Vyacheslav; Gel, Yulia R.; Rogala, James T.; Robertson, Dale M.; Wei, Xiaoqiao

2015-01-01

Long-term temporal trends in water temperature in rivers and streams are typically estimated under the assumption of evenly-spaced space-time measurements. However, sampling times and dates associated with historical water temperature datasets and some sampling designs may be haphazard. As a result, trends in temperature may be confounded with trends in time or space of sampling which, in turn, may yield biased trend estimators and thus unreliable conclusions. We address this concern using multilevel (hierarchical) linear models, where time effects are allowed to vary randomly by day and date effects by year. We evaluate the proposed approach by Monte Carlo simulations with imbalance, sparse data and confounding by trend in time and date of sampling. Simulation results indicate unbiased trend estimators while results from a case study of temperature data from the Illinois River, USA conform to river thermal assumptions. We also propose a new nonparametric bootstrap inference on multilevel models that allows for a relatively flexible and distribution-free quantification of uncertainties. The proposed multilevel modeling approach may be elaborated to accommodate nonlinearities within days and years when sampling times or dates typically span temperature extremes.

MATLAB/Simulink Framework for Modeling Complex Coolant Flow Configurations of Advanced Automotive Thermal Management Systems: Preprint

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

Titov, Eugene; Lustbader, Jason; Leighton, Daniel

The National Renewable Energy Laboratory's (NREL's) CoolSim MATLAB/Simulink modeling framework was extended by including a newly developed coolant loop solution method aimed at reducing the simulation effort for arbitrarily complex thermal management systems. The new approach does not require the user to identify specific coolant loops and their flow. The user only needs to connect the fluid network elements in a manner consistent with the desired schematic. Using the new solution method, a model of NREL's advanced combined coolant loop system for electric vehicles was created that reflected the test system architecture. This system was built using components provided bymore » the MAHLE Group and included both air conditioning and heat pump modes. Validation with test bench data and verification with the previous solution method were performed for 10 operating points spanning a range of ambient temperatures between -2 degrees C and 43 degrees C. The largest root mean square difference between pressure, temperature, energy and mass flow rate data and simulation results was less than 7%.« less

MATLAB/Simulink Framework for Modeling Complex Coolant Flow Configurations of Advanced Automotive Thermal Management Systems

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

Titov, Gene; Lustbader, Jason; Leighton, Daniel

The National Renewable Energy Laboratory's (NREL's) CoolSim MATLAB/Simulink modeling framework was extended by including a newly developed coolant loop solution method aimed at reducing the simulation effort for arbitrarily complex thermal management systems. The new approach does not require the user to identify specific coolant loops and their flow. The user only needs to connect the fluid network elements in a manner consistent with the desired schematic. Using the new solution method, a model of NREL's advanced combined coolant loop system for electric vehicles was created that reflected the test system architecture. This system was built using components provided bymore » the MAHLE Group and included both air conditioning and heat pump modes. Validation with test bench data and verification with the previous solution method were performed for 10 operating points spanning a range of ambient temperatures between -2 degrees C and 43 degrees C. The largest root mean square difference between pressure, temperature, energy and mass flow rate data and simulation results was less than 7%.« less

Actin cable distribution and dynamics arising from cross-linking, motor pulling, and filament turnover

PubMed Central

Tang, Haosu; Laporte, Damien; Vavylonis, Dimitrios

2014-01-01

The growth of fission yeast relies on the polymerization of actin filaments nucleated by formin For3p, which localizes at tip cortical sites. These actin filaments bundle to form actin cables that span the cell and guide the movement of vesicles toward the cell tips. A big challenge is to develop a quantitative understanding of these cellular actin structures. We used computer simulations to study the spatial and dynamical properties of actin cables. We simulated individual actin filaments as semiflexible polymers in three dimensions composed of beads connected with springs. Polymerization out of For3p cortical sites, bundling by cross-linkers, pulling by type V myosin, and severing by cofilin are simulated as growth, cross-linking, pulling, and turnover of the semiflexible polymers. With the foregoing mechanisms, the model generates actin cable structures and dynamics similar to those observed in live-cell experiments. Our simulations reproduce the particular actin cable structures in myoVΔ cells and predict the effect of increased myosin V pulling. Increasing cross-linking parameters generates thicker actin cables. It also leads to antiparallel and parallel phases with straight or curved cables, consistent with observations of cells overexpressing α-actinin. Finally, the model predicts that clustering of formins at cell tips promotes actin cable formation. PMID:25103242

Prolongation of SMAP to Spatiotemporally Seamless Coverage of Continental U.S. Using a Deep Learning Neural Network

NASA Astrophysics Data System (ADS)

Fang, Kuai; Shen, Chaopeng; Kifer, Daniel; Yang, Xiao

2017-11-01

The Soil Moisture Active Passive (SMAP) mission has delivered valuable sensing of surface soil moisture since 2015. However, it has a short time span and irregular revisit schedules. Utilizing a state-of-the-art time series deep learning neural network, Long Short-Term Memory (LSTM), we created a system that predicts SMAP level-3 moisture product with atmospheric forcings, model-simulated moisture, and static physiographic attributes as inputs. The system removes most of the bias with model simulations and improves predicted moisture climatology, achieving small test root-mean-square errors (<0.035) and high-correlation coefficients >0.87 for over 75% of Continental United States, including the forested southeast. As the first application of LSTM in hydrology, we show the proposed network avoids overfitting and is robust for both temporal and spatial extrapolation tests. LSTM generalizes well across regions with distinct climates and environmental settings. With high fidelity to SMAP, LSTM shows great potential for hindcasting, data assimilation, and weather forecasting.

3D ion flow measurements and simulations near a boundary at oblique incidence to a magnetic field

NASA Astrophysics Data System (ADS)

Thompson, Derek S.; Keniley, Shane; Khaziev, Rinat; Curreli, Davide; Good, Timothy N.; Henriquez, Miguel; McIlvain, Julianne; Siddiqui, M. Umair; Scime, Earl E.

2016-10-01

Boundaries at oblique incidence to magnetic fields are abundant in magnetic confinement plasmas. The ion dynamics near these boundaries has implications for applications such as tokamak divertor wall loading and Hall thruster channel erosion. We present 3D, non-perturbative measurements of ion velocity distribution functions (IVDFs), providing ion temperatures and flows upstream of a grounded stainless steel limiter plate immersed in an argon plasma, oriented obliquely to the background axial magnetic field (ψ = 74°). The spatial resolution of the measurements is sufficient to probe the kinetic details of magnetic presheath structures, which span several ion Larmor radii ( 1 cm). Furthermore, we report probe measurements of electron density and temperature, and of local electric potential. To complement these measurements, results from particle-in-cell and Boltzmann models of the same region are presented. These models allow for point-to-point comparison of simulated and measured electrostatic structures and IVDFs at high spatial resolution. NSF Award PHYS-1360278.

CFD Simulation of a Wing-In-Ground-Effect UAV

NASA Astrophysics Data System (ADS)

Lao, C. T.; Wong, E. T. T.

2018-05-01

This paper reports a numerical analysis on a wing section used for a Wing-In-Ground-Effect (WIG) unmanned aerial vehicle (UAV). The wing geometry was created by SolidWorks and the incompressible Reynolds-averaged Navier-Stokes (RANS) equations were solved with the Spalart–Allmaras turbulence model using CFD software ANSYS FLUENT. In FLUENT, the Spalart-Allmaras model has been implemented to use wall functions when the mesh resolution is not sufficiently fine. This might make it the best choice for relatively crude simulations on coarse meshes where accurate turbulent flow computations are not critical. The results show that the lift coefficient and lift-drag ratio derived excellent performance enhancement by ground effect. However, the moment coefficient shows inconsistency when the wing is operating in very low altitude - this is owing to the difficulty on the stability control of WIG vehicle. A drag polar estimation based on the analysis also indicated that the Oswald (or span) efficiency of the wing was improved by ground effect.

Technical note: Simultaneous fully dynamic characterization of multiple input–output relationships in climate models

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

Kravitz, Ben; MacMartin, Douglas G.; Rasch, Philip J.

We introduce system identification techniques to climate science wherein multiple dynamic input–output relationships can be simultaneously characterized in a single simulation. This method, involving multiple small perturbations (in space and time) of an input field while monitoring output fields to quantify responses, allows for identification of different timescales of climate response to forcing without substantially pushing the climate far away from a steady state. We use this technique to determine the steady-state responses of low cloud fraction and latent heat flux to heating perturbations over 22 regions spanning Earth's oceans. We show that the response characteristics are similar to thosemore » of step-change simulations, but in this new method the responses for 22 regions can be characterized simultaneously. Moreover, we can estimate the timescale over which the steady-state response emerges. The proposed methodology could be useful for a wide variety of purposes in climate science, including characterization of teleconnections and uncertainty quantification to identify the effects of climate model tuning parameters.« less

Technical note: Simultaneous fully dynamic characterization of multiple input–output relationships in climate models

DOE PAGES

Kravitz, Ben; MacMartin, Douglas G.; Rasch, Philip J.; ...

2017-02-17

We introduce system identification techniques to climate science wherein multiple dynamic input–output relationships can be simultaneously characterized in a single simulation. This method, involving multiple small perturbations (in space and time) of an input field while monitoring output fields to quantify responses, allows for identification of different timescales of climate response to forcing without substantially pushing the climate far away from a steady state. We use this technique to determine the steady-state responses of low cloud fraction and latent heat flux to heating perturbations over 22 regions spanning Earth's oceans. We show that the response characteristics are similar to thosemore » of step-change simulations, but in this new method the responses for 22 regions can be characterized simultaneously. Moreover, we can estimate the timescale over which the steady-state response emerges. The proposed methodology could be useful for a wide variety of purposes in climate science, including characterization of teleconnections and uncertainty quantification to identify the effects of climate model tuning parameters.« less

TTLEM: Open access tool for building numerically accurate landscape evolution models in MATLAB

NASA Astrophysics Data System (ADS)

Campforts, Benjamin; Schwanghart, Wolfgang; Govers, Gerard

2017-04-01

Despite a growing interest in LEMs, accuracy assessment of the numerical methods they are based on has received little attention. Here, we present TTLEM which is an open access landscape evolution package designed to develop and test your own scenarios and hypothesises. TTLEM uses a higher order flux-limiting finite-volume method to simulate river incision and tectonic displacement. We show that this scheme significantly influences the evolution of simulated landscapes and the spatial and temporal variability of erosion rates. Moreover, it allows the simulation of lateral tectonic displacement on a fixed grid. Through the use of a simple GUI the software produces visible output of evolving landscapes through model run time. In this contribution, we illustrate numerical landscape evolution through a set of movies spanning different spatial and temporal scales. We focus on the erosional domain and use both spatially constant and variable input values for uplift, lateral tectonic shortening, erodibility and precipitation. Moreover, we illustrate the relevance of a stochastic approach for realistic hillslope response modelling. TTLEM is a fully open source software package, written in MATLAB and based on the TopoToolbox platform (topotoolbox.wordpress.com). Installation instructions can be found on this website and the therefore designed GitHub repository.

Explicit Convection over the Western Pacific Warm Pool in the Community Atmospheric Model.

NASA Astrophysics Data System (ADS)

Ziemiaski, Micha Z.; Grabowski, Wojciech W.; Moncrieff, Mitchell W.

2005-05-01

This paper reports on the application of the cloud-resolving convection parameterization (CRCP) to the Community Atmospheric Model (CAM), the atmospheric component of the Community Climate System Model (CCSM). The cornerstone of CRCP is the use of a two-dimensional zonally oriented cloud-system-resolving model to represent processes on mesoscales at the subgrid scale of a climate model. Herein, CRCP is applied at each climate model column over the tropical western Pacific warm pool, in a domain spanning 10°S-10°N, 150°-170°E. Results from the CRCP simulation are compared with CAM in its standard configuration.The CRCP simulation shows significant improvements of the warm pool climate. The cloud condensate distribution is much improved as well as the bias of the tropopause height. More realistic structure of the intertropical convergence zone (ITCZ) during the boreal winter and better representation of the variability of convection are evident. In particular, the diurnal cycle of precipitation has phase and amplitude in good agreement with observations. Also improved is the large-scale organization of the tropical convection, especially superclusters associated with Madden-Julian oscillation (MJO)-like systems. Location and propagation characteristics, as well as lower-tropospheric cyclonic and upper-tropospheric anticyclonic gyres, are more realistic than in the standard CAM. Finally, the simulations support an analytic theory of dynamical coupling between organized convection and equatorial beta-plane vorticity dynamics associated with MJO-like systems.

Assessing uncertainties in crop and pasture ensemble model simulations of productivity and N2 O emissions.

PubMed

Ehrhardt, Fiona; Soussana, Jean-François; Bellocchi, Gianni; Grace, Peter; McAuliffe, Russel; Recous, Sylvie; Sándor, Renáta; Smith, Pete; Snow, Val; de Antoni Migliorati, Massimiliano; Basso, Bruno; Bhatia, Arti; Brilli, Lorenzo; Doltra, Jordi; Dorich, Christopher D; Doro, Luca; Fitton, Nuala; Giacomini, Sandro J; Grant, Brian; Harrison, Matthew T; Jones, Stephanie K; Kirschbaum, Miko U F; Klumpp, Katja; Laville, Patricia; Léonard, Joël; Liebig, Mark; Lieffering, Mark; Martin, Raphaël; Massad, Raia S; Meier, Elizabeth; Merbold, Lutz; Moore, Andrew D; Myrgiotis, Vasileios; Newton, Paul; Pattey, Elizabeth; Rolinski, Susanne; Sharp, Joanna; Smith, Ward N; Wu, Lianhai; Zhang, Qing

2018-02-01

Simulation models are extensively used to predict agricultural productivity and greenhouse gas emissions. However, the uncertainties of (reduced) model ensemble simulations have not been assessed systematically for variables affecting food security and climate change mitigation, within multi-species agricultural contexts. We report an international model comparison and benchmarking exercise, showing the potential of multi-model ensembles to predict productivity and nitrous oxide (N 2 O) emissions for wheat, maize, rice and temperate grasslands. Using a multi-stage modelling protocol, from blind simulations (stage 1) to partial (stages 2-4) and full calibration (stage 5), 24 process-based biogeochemical models were assessed individually or as an ensemble against long-term experimental data from four temperate grassland and five arable crop rotation sites spanning four continents. Comparisons were performed by reference to the experimental uncertainties of observed yields and N 2 O emissions. Results showed that across sites and crop/grassland types, 23%-40% of the uncalibrated individual models were within two standard deviations (SD) of observed yields, while 42 (rice) to 96% (grasslands) of the models were within 1 SD of observed N 2 O emissions. At stage 1, ensembles formed by the three lowest prediction model errors predicted both yields and N 2 O emissions within experimental uncertainties for 44% and 33% of the crop and grassland growth cycles, respectively. Partial model calibration (stages 2-4) markedly reduced prediction errors of the full model ensemble E-median for crop grain yields (from 36% at stage 1 down to 4% on average) and grassland productivity (from 44% to 27%) and to a lesser and more variable extent for N 2 O emissions. Yield-scaled N 2 O emissions (N 2 O emissions divided by crop yields) were ranked accurately by three-model ensembles across crop species and field sites. The potential of using process-based model ensembles to predict jointly productivity and N 2 O emissions at field scale is discussed. © 2017 John Wiley & Sons Ltd.

Improved representations of coupled soil–canopy processes in the CABLE land surface model (Subversion revision 3432)

DOE PAGES

Haverd, Vanessa; Cuntz, Matthias; Nieradzik, Lars P.; ...

2016-09-07

CABLE is a global land surface model, which has been used extensively in offline and coupled simulations. While CABLE performs well in comparison with other land surface models, results are impacted by decoupling of transpiration and photosynthesis fluxes under drying soil conditions, often leading to implausibly high water use efficiencies. Here, we present a solution to this problem, ensuring that modelled transpiration is always consistent with modelled photosynthesis, while introducing a parsimonious single-parameter drought response function which is coupled to root water uptake. We further improve CABLE's simulation of coupled soil–canopy processes by introducing an alternative hydrology model with amore » physically accurate representation of coupled energy and water fluxes at the soil–air interface, including a more realistic formulation of transfer under atmospherically stable conditions within the canopy and in the presence of leaf litter. The effects of these model developments are assessed using data from 18 stations from the global eddy covariance FLUXNET database, selected to span a large climatic range. Here, marked improvements are demonstrated, with root mean squared errors for monthly latent heat fluxes and water use efficiencies being reduced by 40 %. Results highlight the important roles of deep soil moisture in mediating drought response and litter in dampening soil evaporation.« less

A Free Wake Numerical Simulation for Darrieus Vertical Axis Wind Turbine Performance Prediction

NASA Astrophysics Data System (ADS)

Belu, Radian

2010-11-01

In the last four decades, several aerodynamic prediction models have been formulated for the Darrieus wind turbine performances and characteristics. We can identified two families: stream-tube and vortex. The paper presents a simplified numerical techniques for simulating vertical axis wind turbine flow, based on the lifting line theory and a free vortex wake model, including dynamic stall effects for predicting the performances of a 3-D vertical axis wind turbine. A vortex model is used in which the wake is composed of trailing stream-wise and shedding span-wise vortices, whose strengths are equal to the change in the bound vortex strength as required by the Helmholz and Kelvin theorems. Performance parameters are computed by application of the Biot-Savart law along with the Kutta-Jukowski theorem and a semi-empirical stall model. We tested the developed model with an adaptation of the earlier multiple stream-tube performance prediction model for the Darrieus turbines. Predictions by using our method are shown to compare favorably with existing experimental data and the outputs of other numerical models. The method can predict accurately the local and global performances of a vertical axis wind turbine, and can be used in the design and optimization of wind turbines for built environment applications.

Improved representations of coupled soil–canopy processes in the CABLE land surface model (Subversion revision 3432)

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

Haverd, Vanessa; Cuntz, Matthias; Nieradzik, Lars P.

CABLE is a global land surface model, which has been used extensively in offline and coupled simulations. While CABLE performs well in comparison with other land surface models, results are impacted by decoupling of transpiration and photosynthesis fluxes under drying soil conditions, often leading to implausibly high water use efficiencies. Here, we present a solution to this problem, ensuring that modelled transpiration is always consistent with modelled photosynthesis, while introducing a parsimonious single-parameter drought response function which is coupled to root water uptake. We further improve CABLE's simulation of coupled soil–canopy processes by introducing an alternative hydrology model with amore » physically accurate representation of coupled energy and water fluxes at the soil–air interface, including a more realistic formulation of transfer under atmospherically stable conditions within the canopy and in the presence of leaf litter. The effects of these model developments are assessed using data from 18 stations from the global eddy covariance FLUXNET database, selected to span a large climatic range. Here, marked improvements are demonstrated, with root mean squared errors for monthly latent heat fluxes and water use efficiencies being reduced by 40 %. Results highlight the important roles of deep soil moisture in mediating drought response and litter in dampening soil evaporation.« less

Improved representations of coupled soil-canopy processes in the CABLE land surface model (Subversion revision 3432)

NASA Astrophysics Data System (ADS)

Haverd, Vanessa; Cuntz, Matthias; Nieradzik, Lars P.; Harman, Ian N.

2016-09-01

CABLE is a global land surface model, which has been used extensively in offline and coupled simulations. While CABLE performs well in comparison with other land surface models, results are impacted by decoupling of transpiration and photosynthesis fluxes under drying soil conditions, often leading to implausibly high water use efficiencies. Here, we present a solution to this problem, ensuring that modelled transpiration is always consistent with modelled photosynthesis, while introducing a parsimonious single-parameter drought response function which is coupled to root water uptake. We further improve CABLE's simulation of coupled soil-canopy processes by introducing an alternative hydrology model with a physically accurate representation of coupled energy and water fluxes at the soil-air interface, including a more realistic formulation of transfer under atmospherically stable conditions within the canopy and in the presence of leaf litter. The effects of these model developments are assessed using data from 18 stations from the global eddy covariance FLUXNET database, selected to span a large climatic range. Marked improvements are demonstrated, with root mean squared errors for monthly latent heat fluxes and water use efficiencies being reduced by 40 %. Results highlight the important roles of deep soil moisture in mediating drought response and litter in dampening soil evaporation.

Predicted thermal response of a cryogenic fuel tank exposed to simulated aerodynamic heating profiles with different cryogens and fill levels

NASA Technical Reports Server (NTRS)

Hanna, Gregory J.; Stephens, Craig A.

1991-01-01

A two dimensional finite difference thermal model was developed to predict the effects of heating profile, fill level, and cryogen type prior to experimental testing the Generic Research Cryogenic Tank (GRCT). These numerical predictions will assist in defining test scenarios, sensor locations, and venting requirements for the GRCT experimental tests. Boiloff rates, tank-wall and fluid temperatures, and wall heat fluxes were determined for 20 computational test cases. The test cases spanned three discrete fill levels and three heating profiles for hydrogen and nitrogen.

The U.S. Geological Survey Monthly Water Balance Model Futures Portal

USGS Publications Warehouse

Bock, Andrew R.; Hay, Lauren E.; Markstrom, Steven L.; Emmerich, Christopher; Talbert, Marian

2017-05-03

The U.S. Geological Survey Monthly Water Balance Model Futures Portal (https://my.usgs.gov/mows/) is a user-friendly interface that summarizes monthly historical and simulated future conditions for seven hydrologic and meteorological variables (actual evapotranspiration, potential evapotranspiration, precipitation, runoff, snow water equivalent, atmospheric temperature, and streamflow) at locations across the conterminous United States (CONUS).The estimates of these hydrologic and meteorological variables were derived using a Monthly Water Balance Model (MWBM), a modular system that simulates monthly estimates of components of the hydrologic cycle using monthly precipitation and atmospheric temperature inputs. Precipitation and atmospheric temperature from 222 climate datasets spanning historical conditions (1952 through 2005) and simulated future conditions (2020 through 2099) were summarized for hydrographic features and used to drive the MWBM for the CONUS. The MWBM input and output variables were organized into an open-access database. An Open Geospatial Consortium, Inc., Web Feature Service allows the querying and identification of hydrographic features across the CONUS. To connect the Web Feature Service to the open-access database, a user interface—the Monthly Water Balance Model Futures Portal—was developed to allow the dynamic generation of summary files and plots  based on plot type, geographic location, specific climate datasets, period of record, MWBM variable, and other options. Both the plots and the data files are made available to the user for download 

Spatial Cognition and Map Interpretation

DTIC Science & Technology

1987-09-01

Terrain association Spatial cognition Map reading Videogames aa mldm II naeaaaaiy and Hontlty by block numbor) Spatial memory span Orientation...ability. Finally, field and classroom performance was compared to wayfinding in a simulated ( videogame ) environment in which position coordinates were...a simulated ( videogame ) environment. Findings: MITAC instruction significantly improved the experimental group’s ability to perform terrain

raaSAFT: A framework enabling coarse-grained molecular dynamics simulations based on the SAFT- γ Mie force field

NASA Astrophysics Data System (ADS)

Ervik, Åsmund; Serratos, Guadalupe Jiménez; Müller, Erich A.

2017-03-01

We describe here raaSAFT, a Python code that enables the setup and running of coarse-grained molecular dynamics simulations in a systematic and efficient manner. The code is built on top of the popular HOOMD-blue code, and as such harnesses the computational power of GPUs. The methodology makes use of the SAFT- γ Mie force field, so the resulting coarse grained pair potentials are both closely linked to and consistent with the macroscopic thermodynamic properties of the simulated fluid. In raaSAFT both homonuclear and heteronuclear models are implemented for a wide range of compounds spanning from linear alkanes, to more complicated fluids such as water and alcohols, all the way up to nonionic surfactants and models of asphaltenes and resins. Adding new compounds as well as new features is made straightforward by the modularity of the code. To demonstrate the ease-of-use of raaSAFT, we give a detailed walkthrough of how to simulate liquid-liquid equilibrium of a hydrocarbon with water. We describe in detail how both homonuclear and heteronuclear compounds are implemented. To demonstrate the performance and versatility of raaSAFT, we simulate a large polymer-solvent mixture with 300 polystyrene molecules dissolved in 42 700 molecules of heptane, reproducing the experimentally observed temperature-dependent solubility of polystyrene. For this case we obtain a speedup of more than three orders of magnitude as compared to atomistically-detailed simulations.

Controls on surface soil drying rates observed by SMAP and simulated by the Noah land surface model

NASA Astrophysics Data System (ADS)

Shellito, Peter J.; Small, Eric E.; Livneh, Ben

2018-03-01

Drydown periods that follow precipitation events provide an opportunity to assess controls on soil evaporation on a continental scale. We use SMAP (Soil Moisture Active Passive) observations and Noah simulations from drydown periods to quantify the role of soil moisture, potential evaporation, vegetation cover, and soil texture on soil drying rates. Rates are determined using finite differences over intervals of 1 to 3 days. In the Noah model, the drying rates are a good approximation of direct soil evaporation rates, and our work suggests that SMAP-observed drying is also predominantly affected by direct soil evaporation. Data cover the domain of the North American Land Data Assimilation System Phase 2 and span the first 1.8 years of SMAP's operation. Drying of surface soil moisture observed by SMAP is faster than that simulated by Noah. SMAP drying is fastest when surface soil moisture levels are high, potential evaporation is high, and when vegetation cover is low. Soil texture plays a minor role in SMAP drying rates. Noah simulations show similar responses to soil moisture and potential evaporation, but vegetation has a minimal effect and soil texture has a much larger effect compared to SMAP. When drying rates are normalized by potential evaporation, SMAP observations and Noah simulations both show that increases in vegetation cover lead to decreases in evaporative efficiency from the surface soil. However, the magnitude of this effect simulated by Noah is much weaker than that determined from SMAP observations.

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

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

Martínez-Sykora, Juan; Pontieu, Bart De; Hansteen, Viggo H.

2017-09-20

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolarmore » diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.« less

The Sherwood simulation suite: overview and data comparisons with the Lyman α forest at redshifts 2 ≤ z ≤ 5

NASA Astrophysics Data System (ADS)

Bolton, James S.; Puchwein, Ewald; Sijacki, Debora; Haehnelt, Martin G.; Kim, Tae-Sun; Meiksin, Avery; Regan, John A.; Viel, Matteo

2017-01-01

We introduce a new set of large-scale, high-resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes of 103-1603h-3 comoving Mpc3, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high-resolution, high signal-to-noise observations of the Ly α forest over the redshift range 2 ≤ z ≤ 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Ly α forest is a natural consequence of the web-like distribution of matter arising in Λcold dark matter cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Ly α absorption lines with column densities N_{H I}>10^{14.5} cm^{-2} at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z = 2.7 and a flatter temperature-density relation is required at z = 2.4, consistent with the expected effects of non-equilibrium ionization during He II reionization.

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

NASA Astrophysics Data System (ADS)

Martínez-Sykora, Juan; De Pontieu, Bart; Carlsson, Mats; Hansteen, Viggo H.; Nóbrega-Siverio, Daniel; Gudiksen, Boris V.

2017-09-01

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion-neutral interaction effects using the generalized Ohm’s law, I.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolar diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

A strategy for analysis of (molecular) equilibrium simulations: Configuration space density estimation, clustering, and visualization

NASA Astrophysics Data System (ADS)

Hamprecht, Fred A.; Peter, Christine; Daura, Xavier; Thiel, Walter; van Gunsteren, Wilfred F.

2001-02-01

We propose an approach for summarizing the output of long simulations of complex systems, affording a rapid overview and interpretation. First, multidimensional scaling techniques are used in conjunction with dimension reduction methods to obtain a low-dimensional representation of the configuration space explored by the system. A nonparametric estimate of the density of states in this subspace is then obtained using kernel methods. The free energy surface is calculated from that density, and the configurations produced in the simulation are then clustered according to the topography of that surface, such that all configurations belonging to one local free energy minimum form one class. This topographical cluster analysis is performed using basin spanning trees which we introduce as subgraphs of Delaunay triangulations. Free energy surfaces obtained in dimensions lower than four can be visualized directly using iso-contours and -surfaces. Basin spanning trees also afford a glimpse of higher-dimensional topographies. The procedure is illustrated using molecular dynamics simulations on the reversible folding of peptide analoga. Finally, we emphasize the intimate relation of density estimation techniques to modern enhanced sampling algorithms.

Comparison of Cirrus Cloud Models: A Project of the GEWEX Cloud System Study (GCSS) Working Group on Cirrus Cloud Systems

NASA Technical Reports Server (NTRS)

Starr, David OC.; Benedetti, Angela; Boehm, Matt; Brown, Philip R. A.; Gierens, Klaus M.; Girard, Eric; Giraud, Vincent; Jakob, Christian; Jensen, Eric; Khvorostyanov, Vitaly;

Modeling of earthquake ground motion in the frequency domain

NASA Astrophysics Data System (ADS)

Thrainsson, Hjortur

In recent years, the utilization of time histories of earthquake ground motion has grown considerably in the design and analysis of civil structures. It is very unlikely, however, that recordings of earthquake ground motion will be available for all sites and conditions of interest. Hence, there is a need for efficient methods for the simulation and spatial interpolation of earthquake ground motion. In addition to providing estimates of the ground motion at a site using data from adjacent recording stations, spatially interpolated ground motions can also be used in design and analysis of long-span structures, such as bridges and pipelines, where differential movement is important. The objective of this research is to develop a methodology for rapid generation of horizontal earthquake ground motion at any site for a given region, based on readily available source, path and site characteristics, or (sparse) recordings. The research includes two main topics: (i) the simulation of earthquake ground motion at a given site, and (ii) the spatial interpolation of earthquake ground motion. In topic (i), models are developed to simulate acceleration time histories using the inverse discrete Fourier transform. The Fourier phase differences, defined as the difference in phase angle between adjacent frequency components, are simulated conditional on the Fourier amplitude. Uniformly processed recordings from recent California earthquakes are used to validate the simulation models, as well as to develop prediction formulas for the model parameters. The models developed in this research provide rapid simulation of earthquake ground motion over a wide range of magnitudes and distances, but they are not intended to replace more robust geophysical models. In topic (ii), a model is developed in which Fourier amplitudes and Fourier phase angles are interpolated separately. A simple dispersion relationship is included in the phase angle interpolation. The accuracy of the interpolation model is assessed using data from the SMART-1 array in Taiwan. The interpolation model provides an effective method to estimate ground motion at a site using recordings from stations located up to several kilometers away. Reliable estimates of differential ground motion are restricted to relatively limited ranges of frequencies and inter-station spacings.

Modeling Vegetation Growth Impact on Groundwater Recharge

NASA Astrophysics Data System (ADS)

Anurag, H.; Ng, G. H. C.; Tipping, R.

2017-12-01

Vegetation growth is affected by variability in climate and land-cover / land-use over a range of temporal and spatial scales. Vegetation also modifies water budget through interception and evapotranspiration and thus has a significant impact on groundwater recharge. Most groundwater recharge assessments represent vegetation using specified, static parameter, such as for leaf-area-index, but this neglects the effect of vegetation dynamics on recharge estimates. Our study addresses this gap by including vegetation growth in model simulations of recharge. We use NCAR's Community Land Model v4.5 with its BGC module (BGC is the new CLM4.5 biogeochemistry). It integrates prognostic vegetation growth with land-surface and subsurface hydrological processes and can thus capture the effect of vegetation on groundwater. A challenge, however, is the need to resolve uncertainties in model inputs ranging from vegetation growth parameters all the way down to the water table. We have compiled diverse data spanning meteorological inputs to subsurface geology and use these to implement ensemble model simulations to evaluate the possible effects of dynamic vegetation growth (versus specified, static vegetation parameterizations) on estimating groundwater recharge. We present preliminary results for select data-intensive test locations throughout the state of Minnesota (USA), which has a sharp east-west precipitation gradient that makes it an apt testbed for examining ecohydrologic relationships across different temperate climatic settings and ecosystems. Using the ensemble simulations, we examine the effect of seasonal to interannual variability of vegetation growth on recharge and water table depths, which has implications for predicting the combined impact of climate, vegetation, and geology on groundwater resources. Future work will include distributed model simulations over the entire state, as well as conditioning uncertain vegetation and subsurface parameters on remote sensing data and statewide water table records using data assimilation.

Evaluation of a Cloud Resolving Model Using TRMM Observations for Multiscale Modeling Applications

NASA Technical Reports Server (NTRS)

Posselt, Derek J.; L'Ecuyer, Tristan; Tao, Wei-Kuo; Hou, Arthur Y.; Stephens, Graeme L.

2007-01-01

The climate change simulation community is moving toward use of global cloud resolving models (CRMs), however, current computational resources are not sufficient to run global CRMs over the hundreds of years necessary to produce climate change estimates. As an intermediate step between conventional general circulation models (GCMs) and global CRMs, many climate analysis centers are embedding a CRM in each grid cell of a conventional GCM. These Multiscale Modeling Frameworks (MMFs) represent a theoretical advance over the use of conventional GCM cloud and convection parameterizations, but have been shown to exhibit an overproduction of precipitation in the tropics during the northern hemisphere summer. In this study, simulations of clouds, precipitation, and radiation over the South China Sea using the CRM component of the NASA Goddard MMF are evaluated using retrievals derived from the instruments aboard the Tropical Rainfall Measuring Mission (TRMM) satellite platform for a 46-day time period that spans 5 May - 20 June 1998. The NASA Goddard Cumulus Ensemble (GCE) model is forced with observed largescale forcing derived from soundings taken during the intensive observing period of the South China Sea Monsoon Experiment. It is found that the GCE configuration used in the NASA Goddard MMF responds too vigorously to the imposed large-scale forcing, accumulating too much moisture and producing too much cloud cover during convective phases, and overdrying the atmosphere and suppressing clouds during monsoon break periods. Sensitivity experiments reveal that changes to ice cloud microphysical parameters have a relatively large effect on simulated clouds, precipitation, and radiation, while changes to grid spacing and domain length have little effect on simulation results. The results motivate a more detailed and quantitative exploration of the sources and magnitude of the uncertainty associated with specified cloud microphysical parameters in the CRM components of MMFs.

Risk assessment of watershed erosion at Naesung Stream, South Korea.

PubMed

Ji, Un; Velleux, Mark; Julien, Pierre Y; Hwang, Manha

2014-04-01

A three-tiered approach was used to assess erosion risks within the Nakdong River Basin in South Korea and included: (1) a screening based on topography and land use; (2) a lumped parameter analysis using RUSLE; and (3) a detailed analysis using TREX, a fully distributed watershed model. These tiers span a range of spatial and temporal scales, with each tier providing increasing detail and resolution. The first two tiers were applied to the entire Nakdong River Basin and the Naesung Stream watershed was identified as having the highest soil erosion risk and potential for sedimentation problems. For the third tier, the TREX watershed model simulated runoff, channel flow, soil erosion, and stream sediment transport in the Naesung Stream watershed at very high resolution. TREX was calibrated for surface flows and sediment transport, and was used to simulate conditions for a large design storm. Highly erosive areas were identified along ridgelines in several headwater areas, with the northeast area of Songriwon having a particularly high erosion potential. Design storm simulations also indicated that sediment deposition of up to 55 cm could occur. Copyright © 2014 Elsevier Ltd. All rights reserved.

Comparison of 2D simulations of detached divertor plasmas with divertor Thomson measurements in the DIII-D tokamak

DOE PAGES

Rognlien, Thomas D.; McLean, Adam G.; Fenstermacher, Max E.; ...

2017-01-27

A modeling study is reported using new 2D data from DIII-D tokamak divertor plasmas and improved 2D transport model that includes large cross-field drifts for the numerically difficult H-mode regime. The data set, which spans a range of plasmas densities for both forward and reverse toroidal magnetic field (B t) over a range of plasma densities, is provided by divertor Thomson scattering (DTS). Measurements utilizing X-point sweeping give corresponding 2D profiles of electron temperature (T e) and density (n e) across both divertor legs for individual discharges. The calculations show the same features of in/out plasma asymmetries as measured inmore » the experiment, with the normal B t direction (ion ∇B drift toward the X-point) having higher n e and lower T e in the inner divertor leg than outer. Corresponding emission data for total radiated power shows a strong inner-divertor/outer-divertor asymmetry that is reproduced by the simulations. Furthermore, these 2D UEDGE transport simulations are enabled for steep-gradient H-mode conditions by newly implemented algorithms to control isolated grid-scale irregularities.« less

Meso- to micro-scale coupled simulations of flow over complex terrain at the Perdigao site

NASA Astrophysics Data System (ADS)

Neher, J.; van Veen, L.; Chow, F. K.; Mirocha, J. D.; Lundquist, J. K.

2017-12-01

In this work, the site of the 2017 Perdigao field campaign is analyzed with high resolution large-eddy simulations generated using the Weather Research and Forecasting (WRF) model as a coupled mesoscale to microscale model. The fine topographic features of the site, with its ridgelines a mere 1.2 km apart, the occurrence of intermittent turbulence at night, and the presence of a wind turbine on one of the ridgelines pose a challenge for many current numerical models. Key test cases in the observational data that demonstrate these modelling difficulties are identified, and advanced modeling techniques for overcoming these issues in the WRF model are presented. These techniques include vertical grid nesting for control of the grid aspect ratio, the cell perturbation method for accelerating the generation of turbulence at the boundary, the dynamic reconstruction model as a closure model that allows for backscatter of turbulence, and the actuator disk model for representing the turbine wake. Multiple nesting configurations are considered, with special consideration given to spanning the `grey zone' where neither PBL nor LES closures are effective. Comparisons between model results and measured sounding, meteorological tower, and Lidar data are used to evaluate the effectiveness of these techniques, and the model results are evaluated to provide a broader view of the flow field and the turbine wake interactions at the site.

GCM simulations of volcanic aerosol forcing. I - Climate changes induced by steady-state perturbations

NASA Technical Reports Server (NTRS)

Pollack, James B.; Rind, David; Lacis, Andrew; Hansen, James E.; Sato, Makiko; Ruedy, Reto

1993-01-01

The response of the climate system to a temporally and spatially constant amount of volcanic particles is simulated using a general circulation model (GCM). The optical depth of the aerosols is chosen so as to produce approximately the same amount of forcing as results from doubling the present CO2 content of the atmosphere and from the boundary conditions associated with the peak of the last ice age. The climate changes produced by long-term volcanic aerosol forcing are obtained by differencing this simulation and one made for the present climate with no volcanic aerosol forcing. The simulations indicate that a significant cooling of the troposphere and surface can occur at times of closely spaced multiple sulfur-rich volcanic explosions that span time scales of decades to centuries. The steady-state climate response to volcanic forcing includes a large expansion of sea ice, especially in the Southern Hemisphere; a resultant large increase in surface and planetary albedo at high latitudes; and sizable changes in the annually and zonally averaged air temperature.

Modeling and simulating the neuromuscular mechanisms regulating ankle and knee joint stiffness during human locomotion.

PubMed

Sartori, Massimo; Maculan, Marco; Pizzolato, Claudio; Reggiani, Monica; Farina, Dario

2015-10-01

This work presents an electrophysiologically and dynamically consistent musculoskeletal model to predict stiffness in the human ankle and knee joints as derived from the joints constituent biological tissues (i.e., the spanning musculotendon units). The modeling method we propose uses electromyography (EMG) recordings from 13 muscle groups to drive forward dynamic simulations of the human leg in five healthy subjects during overground walking and running. The EMG-driven musculoskeletal model estimates musculotendon and resulting joint stiffness that is consistent with experimental EMG data as well as with the experimental joint moments. This provides a framework that allows for the first time observing 1) the elastic interplay between the knee and ankle joints, 2) the individual muscle contribution to joint stiffness, and 3) the underlying co-contraction strategies. It provides a theoretical description of how stiffness modulates as a function of muscle activation, fiber contraction, and interacting tendon dynamics. Furthermore, it describes how this differs from currently available stiffness definitions, including quasi-stiffness and short-range stiffness. This work offers a theoretical and computational basis for describing and investigating the neuromuscular mechanisms underlying human locomotion. Copyright © 2015 the American Physiological Society.

Wind Farm LES Simulations Using an Overset Methodology

NASA Astrophysics Data System (ADS)

Ananthan, Shreyas; Yellapantula, Shashank

2017-11-01

Accurate simulation of wind farm wakes under realistic atmospheric inflow conditions and complex terrain requires modeling a wide range of length and time scales. The computational domain can span several kilometers while requiring mesh resolutions in O(10-6) to adequately resolve the boundary layer on the blade surface. Overset mesh methodology offers an attractive option to address the disparate range of length scales; it allows embedding body-confirming meshes around turbine geomtries within nested wake capturing meshes of varying resolutions necessary to accurately model the inflow turbulence and the resulting wake structures. Dynamic overset hole-cutting algorithms permit relative mesh motion that allow this nested mesh structure to track unsteady inflow direction changes, turbine control changes (yaw and pitch), and wake propagation. An LES model with overset mesh for localized mesh refinement is used to analyze wind farm wakes and performance and compared with local mesh refinements using non-conformal (hanging node) unstructured meshes. Turbine structures will be modeled using both actuator line approaches and fully-resolved structures to test the efficacy of overset methods for wind farm applications. Exascale Computing Project (ECP), Project Number: 17-SC-20-SC, a collaborative effort of two DOE organizations - the Office of Science and the National Nuclear Security Administration.

Prediction of a Densely Loaded Particle-Laden Jet using a Euler-Lagrange Dense Spray Model

NASA Astrophysics Data System (ADS)

Pakseresht, Pedram; Apte, Sourabh V.

2017-11-01

Modeling of a dense spray regime using an Euler-Lagrange discrete-element approach is challenging because of local high volume loading. A subgrid cluster of droplets can lead to locally high void fractions for the disperse phase. Under these conditions, spatio-temporal changes in the carrier phase volume fractions, which are commonly neglected in spray simulations in an Euler-Lagrange two-way coupling model, could become important. Accounting for the carrier phase volume fraction variations, leads to zero-Mach number, variable density governing equations. Using pressure-based solvers, this gives rise to a source term in the pressure Poisson equation and a non-divergence free velocity field. To test the validity and predictive capability of such an approach, a round jet laden with solid particles is investigated using Direct Numerical Simulation and compared with available experimental data for different loadings. Various volume fractions spanning from dilute to dense regimes are investigated with and without taking into account the volume displacement effects. The predictions of the two approaches are compared and analyzed to investigate the effectiveness of the dense spray model. Financial support was provided by National Aeronautics and Space Administration (NASA).

Modeling and simulating the neuromuscular mechanisms regulating ankle and knee joint stiffness during human locomotion

PubMed Central

Maculan, Marco; Pizzolato, Claudio; Reggiani, Monica; Farina, Dario

2015-01-01

This work presents an electrophysiologically and dynamically consistent musculoskeletal model to predict stiffness in the human ankle and knee joints as derived from the joints constituent biological tissues (i.e., the spanning musculotendon units). The modeling method we propose uses electromyography (EMG) recordings from 13 muscle groups to drive forward dynamic simulations of the human leg in five healthy subjects during overground walking and running. The EMG-driven musculoskeletal model estimates musculotendon and resulting joint stiffness that is consistent with experimental EMG data as well as with the experimental joint moments. This provides a framework that allows for the first time observing 1) the elastic interplay between the knee and ankle joints, 2) the individual muscle contribution to joint stiffness, and 3) the underlying co-contraction strategies. It provides a theoretical description of how stiffness modulates as a function of muscle activation, fiber contraction, and interacting tendon dynamics. Furthermore, it describes how this differs from currently available stiffness definitions, including quasi-stiffness and short-range stiffness. This work offers a theoretical and computational basis for describing and investigating the neuromuscular mechanisms underlying human locomotion. PMID:26245321

Implementation of a Marauding Insect Module (MIM, version 1.0) in the Integrated BIosphere Simulator (IBIS, version 2.6b4) dynamic vegetation-land surface model

NASA Astrophysics Data System (ADS)

Landry, Jean-Sébastien; Price, David T.; Ramankutty, Navin; Parrott, Lael; Damon Matthews, H.

2016-04-01

Insects defoliate and kill plants in many ecosystems worldwide. The consequences of these natural processes on terrestrial ecology and nutrient cycling are well established, and their potential climatic effects resulting from modified land-atmosphere exchanges of carbon, energy, and water are increasingly being recognized. We developed a Marauding Insect Module (MIM) to quantify, in the Integrated BIosphere Simulator (IBIS), the consequences of insect activity on biogeochemical and biogeophysical fluxes, also accounting for the effects of altered vegetation dynamics. MIM can simulate damage from three different insect functional types: (1) defoliators on broadleaf deciduous trees, (2) defoliators on needleleaf evergreen trees, and (3) bark beetles on needleleaf evergreen trees, with the resulting impacts being estimated by IBIS based on the new, insect-modified state of the vegetation. MIM further accounts for the physical presence and gradual fall of insect-killed dead standing trees. The design of MIM should facilitate the addition of other insect types besides the ones already included and could guide the development of similar modules for other process-based vegetation models. After describing IBIS-MIM, we illustrate the usefulness of the model by presenting results spanning daily to centennial timescales for vegetation dynamics and cycling of carbon, energy, and water in a simplified setting and for bark beetles only. More precisely, we simulated 100 % mortality events from the mountain pine beetle for three locations in western Canada. We then show that these simulated impacts agree with many previous studies based on field measurements, satellite data, or modelling. MIM and similar tools should therefore be of great value in assessing the wide array of impacts resulting from insect-induced plant damage in the Earth system.

The importance of wind-flux feedbacks during the November CINDY-DYNAMO MJO event

NASA Astrophysics Data System (ADS)

Riley Dellaripa, Emily; Maloney, Eric; van den Heever, Susan

2015-04-01

High-resolution, large-domain cloud resolving model (CRM) simulations probing the importance of wind-flux feedbacks to Madden-Julian Oscillation (MJO) convection are performed for the November 2011 CINDY-DYNAMO MJO event. The work is motivated by observational analysis from RAMA buoys in the Indian Ocean and TRMM precipitation retrievals that show a positive correlation between MJO precipitation and wind-induced surface fluxes, especially latent heat fluxes, during and beyond the CINDY-DYNAMO time period. Simulations are done using Colorado State University's Regional Atmospheric Modeling System (RAMS). The domain setup is oceanic and spans 1000 km x 1000 km with 1.5 km horizontal resolution and 65 stretched vertical levels centered on the location of Gan Island - one of the major CINDY-DYNAMO observation points. The model is initialized with ECMWF reanalysis and Aqua MODIS sea surface temperatures. Nudging from ECMWF reanalysis is applied at the domain periphery to encourage realistic evolution of MJO convection. The control experiment is run for the entire month of November so both suppressed and active, as well as, transitional phases of the MJO are modeled. In the control experiment, wind-induced surface fluxes are activated through the surface bulk aerodynamic formula and allowed to evolve organically. Sensitivity experiments are done by restarting the control run one week into the simulation and controlling the wind-induced flux feedbacks. In one sensitivity experiment, wind-induced surface flux feedbacks are completely denied, while in another experiment the winds are kept constant at the control simulations mean surface wind speed. The evolution of convection, especially on the mesoscale, is compared between the control and sensitivity simulations.

Insight Center | Computational Science | NREL

Science.gov Websites

effectively convey information and illustrate research findings to stakeholders and visitors. The -turbine array simulations. Observational data span from the nanostructures of biomass pretreatments to the

Atomic Detail Visualization of Photosynthetic Membranes with GPU-Accelerated Ray Tracing

PubMed Central

Vandivort, Kirby L.; Barragan, Angela; Singharoy, Abhishek; Teo, Ivan; Ribeiro, João V.; Isralewitz, Barry; Liu, Bo; Goh, Boon Chong; Phillips, James C.; MacGregor-Chatwin, Craig; Johnson, Matthew P.; Kourkoutis, Lena F.; Hunter, C. Neil

2016-01-01

The cellular process responsible for providing energy for most life on Earth, namely photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. We present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. We describe the techniques that were used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers. PMID:27274603

Simulation of bio-locomotion by a momentum redistribution technique for self-propulsion

NASA Astrophysics Data System (ADS)

Curet, Oscar; Shirgaonkar, Anup; Patankar, Neelesh; Maciver, Malcolm

2007-11-01

We have developed a general purpose computational approach for self-propulsion based on a momentum redistribution concept. In this poster, our primary goal is to show that the technique can simulate swimming of various organisms without using reduced order models for fluid dynamics. The approach fully resolves the motion of the organism and the surrounding fluid. Thus, it is an effective tool to obtain forces, flow fields, as well as the swimming velocity when the deformation kinematics of the organism are available from observational data. We will present images of computational flow fields for several examples including the aquatic locomotion of sperm, jellyfish, eel, and blackghost knifefish. These examples span a range of body configurations, swimming gaits, and Reynolds numbers in their natural environments. Peculiarities of various modes of swimming will be highlighted.

Antimisting kerosene atomization and flammability

NASA Technical Reports Server (NTRS)

Fleeter, R.; Petersen, R. A.; Toaz, R. D.; Jakub, A.; Sarohia, V.

1982-01-01

Various parameters found to affect the flammability of antimisting kerosene (Jet A + polymer additive) are investigated. Digital image processing was integrated into a technique for measurement of fuel spray characteristics. This technique was developed to avoid many of the error sources inherent to other spray assessment techniques and was applied to the study of engine fuel nozzle atomization performance with Jet A and antimisting fuel. Aircraft accident fuel spill and ignition dynamics were modeled in a steady state simulator allowing flammability to be measured as a function of airspeed, fuel flow rate, fuel jet Reynolds number and polymer concentration. The digital imaging technique was employed to measure spray characteristics in this simulation and these results were related to flammability test results. Scaling relationships were investigated through correlation of experimental results with characteristic dimensions spanning more than two orders of magnitude.

An Experimental Study of the Aerodynamics of a Swept and Unswept Semispan Wing with a Simulated Glaze Ice Accretion

NASA Technical Reports Server (NTRS)

Bragg, Michael B.

1994-01-01

Two semispan wings, one with a rectangular planform and one with 30 degrees of leading edge sweep were tested. Both had a NACA 0012 airfoil section, and both were tested clean and with simulated glaze ice shapes on their leading edges. Several surface roughness were tested. Each model geometry is documented and each surface roughness is explained. Aerodynamic performance of the wing in the form of sectional lift and integrated three-dimensional lift is documented through pressure measurements obtained from rows of surface pressure taps placed at five span locations on the wing. For the rectangular wing, sectional drag near the midspan is obtained from wake total pressure profiles. The data is presented in tabular and graphical form and is also available on computer disk.

Mapping and uncertainty analysis of energy and pitch angle phase space in the DIII-D fast ion loss detector.

PubMed

Pace, D C; Pipes, R; Fisher, R K; Van Zeeland, M A

2014-11-01

New phase space mapping and uncertainty analysis of energetic ion loss data in the DIII-D tokamak provides experimental results that serve as valuable constraints in first-principles simulations of energetic ion transport. Beam ion losses are measured by the fast ion loss detector (FILD) diagnostic system consisting of two magnetic spectrometers placed independently along the outer wall. Monte Carlo simulations of mono-energetic and single-pitch ions reaching the FILDs are used to determine the expected uncertainty in the measurements. Modeling shows that the variation in gyrophase of 80 keV beam ions at the FILD aperture can produce an apparent measured energy signature spanning across 50-140 keV. These calculations compare favorably with experiments in which neutral beam prompt loss provides a well known energy and pitch distribution.

Communication: Polymer entanglement dynamics: Role of attractive interactions

DOE PAGES

Grest, Gary S.

2016-10-10

The coupled dynamics of entangled polymers, which span broad time and length scales, govern their unique viscoelastic properties. To follow chain mobility by numerical simulations from the intermediate Rouse and reptation regimes to the late time diffusive regime, highly coarse grained models with purely repulsive interactions between monomers are widely used since they are computationally the most efficient. In this paper, using large scale molecular dynamics simulations, the effect of including the attractive interaction between monomers on the dynamics of entangled polymer melts is explored for the first time over a wide temperature range. Attractive interactions have little effect onmore » the local packing for all temperatures T and on the chain mobility for T higher than about twice the glass transition T g. Finally, these results, across a broad range of molecular weight, show that to study the dynamics of entangled polymer melts, the interactions can be treated as pure repulsive, confirming a posteriori the validity of previous studies and opening the way to new large scale numerical simulations.« less

Model-based functional neuroimaging using dynamic neural fields: An integrative cognitive neuroscience approach

PubMed Central

Wijeakumar, Sobanawartiny; Ambrose, Joseph P.; Spencer, John P.; Curtu, Rodica

2017-01-01

A fundamental challenge in cognitive neuroscience is to develop theoretical frameworks that effectively span the gap between brain and behavior, between neuroscience and psychology. Here, we attempt to bridge this divide by formalizing an integrative cognitive neuroscience approach using dynamic field theory (DFT). We begin by providing an overview of how DFT seeks to understand the neural population dynamics that underlie cognitive processes through previous applications and comparisons to other modeling approaches. We then use previously published behavioral and neural data from a response selection Go/Nogo task as a case study for model simulations. Results from this study served as the ‘standard’ for comparisons with a model-based fMRI approach using dynamic neural fields (DNF). The tutorial explains the rationale and hypotheses involved in the process of creating the DNF architecture and fitting model parameters. Two DNF models, with similar structure and parameter sets, are then compared. Both models effectively simulated reaction times from the task as we varied the number of stimulus-response mappings and the proportion of Go trials. Next, we directly simulated hemodynamic predictions from the neural activation patterns from each model. These predictions were tested using general linear models (GLMs). Results showed that the DNF model that was created by tuning parameters to capture simultaneously trends in neural activation and behavioral data quantitatively outperformed a Standard GLM analysis of the same dataset. Further, by using the GLM results to assign functional roles to particular clusters in the brain, we illustrate how DNF models shed new light on the neural populations’ dynamics within particular brain regions. Thus, the present study illustrates how an interactive cognitive neuroscience model can be used in practice to bridge the gap between brain and behavior. PMID:29118459

Representation of Precipitation in a Decade-long Continental-Scale Convection-Resolving Climate Simulation

NASA Astrophysics Data System (ADS)

Leutwyler, D.; Fuhrer, O.; Ban, N.; Lapillonne, X.; Lüthi, D.; Schar, C.

2017-12-01

The representation of moist convection in climate models represents a major challenge, due to the small scales involved. Regional climate simulations using horizontal resolutions of O(1km) allow to explicitly resolve deep convection leading to an improved representation of the water cycle. However, due to their extremely demanding computational requirements, they have so far been limited to short simulations and/or small computational domains. A new version of the Consortium for Small-Scale Modeling weather and climate model (COSMO) is capable of exploiting new supercomputer architectures employing GPU accelerators, and allows convection-resolving climate simulations on computational domains spanning continents and time periods up to one decade. We present results from a decade-long, convection-resolving climate simulation on a European-scale computational domain. The simulation has a grid spacing of 2.2 km, 1536x1536x60 grid points, covers the period 1999-2008, and is driven by the ERA-Interim reanalysis. Specifically we present an evaluation of hourly rainfall using a wide range of data sets, including several rain-gauge networks and a remotely-sensed lightning data set. Substantial improvements are found in terms of the diurnal cycles of precipitation amount, wet-hour frequency and all-hour 99th percentile. However the results also reveal substantial differences between regions with and without strong orographic forcing. Furthermore we present an index for deep-convective activity based on the statistics of vertical motion. Comparison of the index with lightning data shows that the convection-resolving climate simulations are able to reproduce important features of the annual cycle of deep convection in Europe. Leutwyler D., D. Lüthi, N. Ban, O. Fuhrer, and C. Schär (2017): Evaluation of the Convection-Resolving Climate Modeling Approach on Continental Scales , J. Geophys. Res. Atmos., 122, doi:10.1002/2016JD026013.

Including Finite Surface Span Effects in Empirical Jet-Surface Interaction Noise Models

NASA Technical Reports Server (NTRS)

Brown, Clifford A.

2016-01-01

The effect of finite span on the jet-surface interaction noise source and the jet mixing noise shielding and reflection effects is considered using recently acquired experimental data. First, the experimental setup and resulting data are presented with particular attention to the role of surface span on far-field noise. These effects are then included in existing empirical models that have previously assumed that all surfaces are semi-infinite. This extended abstract briefly describes the experimental setup and data leaving the empirical modeling aspects for the final paper.

Coseismic slip on the southern Cascadia megathrust implied by tsunami deposits in an Oregon lake and earthquake-triggered marine turbidites

NASA Astrophysics Data System (ADS)

Witter, Robert C.; Zhang, Yinglong; Wang, Kelin; Goldfinger, Chris; Priest, George R.; Allan, Jonathan C.

2012-10-01

We test hypothetical tsunami scenarios against a 4,600-year record of sandy deposits in a southern Oregon coastal lake that offer minimum inundation limits for prehistoric Cascadia tsunamis. Tsunami simulations constrain coseismic slip estimates for the southern Cascadia megathrust and contrast with slip deficits implied by earthquake recurrence intervals from turbidite paleoseismology. We model the tsunamigenic seafloor deformation using a three-dimensional elastic dislocation model and test three Cascadia earthquake rupture scenarios: slip partitioned to a splay fault; slip distributed symmetrically on the megathrust; and slip skewed seaward. Numerical tsunami simulations use the hydrodynamic finite element model, SELFE, that solves nonlinear shallow-water wave equations on unstructured grids. Our simulations of the 1700 Cascadia tsunami require >12-13 m of peak slip on the southern Cascadia megathrust offshore southern Oregon. The simulations account for tidal and shoreline variability and must crest the ˜6-m-high lake outlet to satisfy geological evidence of inundation. Accumulating this slip deficit requires ≥360-400 years at the plate convergence rate, exceeding the 330-year span of two earthquake cycles preceding 1700. Predecessors of the 1700 earthquake likely involved >8-9 m of coseismic slip accrued over >260 years. Simple slip budgets constrained by tsunami simulations allow an average of 5.2 m of slip per event for 11 additional earthquakes inferred from the southern Cascadia turbidite record. By comparison, slip deficits inferred from time intervals separating earthquake-triggered turbidites are poor predictors of coseismic slip because they meet geological constraints for only 4 out of 12 (˜33%) Cascadia tsunamis.

Carbon storage, timber production, and biodiversity: comparing ecosystem services with multi-criteria decision analysis

USGS Publications Warehouse

Schwenk, W. Scott; Donovan, Therese; Keeton, William S.; Nunery, Jared S.

2012-01-01

Increasingly, land managers seek ways to manage forests for multiple ecosystem services and functions, yet considerable challenges exist in comparing disparate services and balancing trade-offs among them. We applied multi-criteria decision analysis (MCDA) and forest simulation models to simultaneously consider three objectives: (1) storing carbon, (2) producing timber and wood products, and (3) sustaining biodiversity. We used the Forest Vegetation Simulator (FVS) applied to 42 northern hardwood sites to simulate forest development over 100 years and to estimate carbon storage and timber production. We estimated biodiversity implications with occupancy models for 51 terrestrial bird species that were linked to FVS outputs. We simulated four alternative management prescriptions that spanned a range of harvesting intensities and forest structure retention. We found that silvicultural approaches emphasizing less frequent harvesting and greater structural retention could be expected to achieve the greatest net carbon storage but also produce less timber. More intensive prescriptions would enhance biodiversity because positive responses of early successional species exceeded negative responses of late successional species within the heavily forested study area. The combinations of weights assigned to objectives had a large influence on which prescriptions were scored as optimal. Overall, we found that a diversity of silvicultural approaches is likely to be preferable to any single approach, emphasizing the need for landscape-scale management to provide a full range of ecosystem goods and services. Our analytical framework that combined MCDA with forest simulation modeling was a powerful tool in understanding trade-offs among management objectives and how they can be simultaneously accommodated.

Scalar mixing in LES/PDF of a high-Ka premixed turbulent jet flame

NASA Astrophysics Data System (ADS)

You, Jiaping; Yang, Yue

2016-11-01

We report a large-eddy simulation (LES)/probability density function (PDF) study of a high-Ka premixed turbulent flame in the Lund University Piloted Jet (LUPJ) flame series, which has been investigated using direct numerical simulation (DNS) and experiments. The target flame, featuring broadened preheat and reaction zones, is categorized into the broken reaction zone regime. In the present study, three widely used mixing modes, namely the Interaction by Exchange with the Mean (IEM), Modified Curl (MC), and Euclidean Minimum Spanning Tree (EMST) models are applied to assess their performance through detailed a posteriori comparisons with DNS. A dynamic model for the time scale of scalar mixing is formulated to describe the turbulent mixing of scalars at small scales. Better quantitative agreement for the mean temperature and mean mass fractions of major and minor species are obtained with the MC and EMST models than with the IEM model. The multi-scalar mixing in composition space with the three models are analyzed to assess the modeling of the conditional molecular diffusion term. In addition, we demonstrate that the product of OH and CH2O concentrations can be a good surrogate of the local heat release rate in this flame. This work is supported by the National Natural Science Foundation of China (Grant Nos. 11521091 and 91541204).

Materials integrity in microsystems: a framework for a petascale predictive-science-based multiscale modeling and simulation system

NASA Astrophysics Data System (ADS)

To, Albert C.; Liu, Wing Kam; Olson, Gregory B.; Belytschko, Ted; Chen, Wei; Shephard, Mark S.; Chung, Yip-Wah; Ghanem, Roger; Voorhees, Peter W.; Seidman, David N.; Wolverton, Chris; Chen, J. S.; Moran, Brian; Freeman, Arthur J.; Tian, Rong; Luo, Xiaojuan; Lautenschlager, Eric; Challoner, A. Dorian

2008-09-01

Microsystems have become an integral part of our lives and can be found in homeland security, medical science, aerospace applications and beyond. Many critical microsystem applications are in harsh environments, in which long-term reliability needs to be guaranteed and repair is not feasible. For example, gyroscope microsystems on satellites need to function for over 20 years under severe radiation, thermal cycling, and shock loading. Hence a predictive-science-based, verified and validated computational models and algorithms to predict the performance and materials integrity of microsystems in these situations is needed. Confidence in these predictions is improved by quantifying uncertainties and approximation errors. With no full system testing and limited sub-system testings, petascale computing is certainly necessary to span both time and space scales and to reduce the uncertainty in the prediction of long-term reliability. This paper presents the necessary steps to develop predictive-science-based multiscale modeling and simulation system. The development of this system will be focused on the prediction of the long-term performance of a gyroscope microsystem. The environmental effects to be considered include radiation, thermo-mechanical cycling and shock. Since there will be many material performance issues, attention is restricted to creep resulting from thermal aging and radiation-enhanced mass diffusion, material instability due to radiation and thermo-mechanical cycling and damage and fracture due to shock. To meet these challenges, we aim to develop an integrated multiscale software analysis system that spans the length scales from the atomistic scale to the scale of the device. The proposed software system will include molecular mechanics, phase field evolution, micromechanics and continuum mechanics software, and the state-of-the-art model identification strategies where atomistic properties are calibrated by quantum calculations. We aim to predict the long-term (in excess of 20 years) integrity of the resonator, electrode base, multilayer metallic bonding pads, and vacuum seals in a prescribed mission. Although multiscale simulations are efficient in the sense that they focus the most computationally intensive models and methods on only the portions of the space time domain needed, the execution of the multiscale simulations associated with evaluating materials and device integrity for aerospace microsystems will require the application of petascale computing. A component-based software strategy will be used in the development of our massively parallel multiscale simulation system. This approach will allow us to take full advantage of existing single scale modeling components. An extensive, pervasive thrust in the software system development is verification, validation, and uncertainty quantification (UQ). Each component and the integrated software system need to be carefully verified. An UQ methodology that determines the quality of predictive information available from experimental measurements and packages the information in a form suitable for UQ at various scales needs to be developed. Experiments to validate the model at the nanoscale, microscale, and macroscale are proposed. The development of a petascale predictive-science-based multiscale modeling and simulation system will advance the field of predictive multiscale science so that it can be used to reliably analyze problems of unprecedented complexity, where limited testing resources can be adequately replaced by petascale computational power, advanced verification, validation, and UQ methodologies.

On a sparse pressure-flow rate condensation of rigid circulation models

PubMed Central

Schiavazzi, D. E.; Hsia, T. Y.; Marsden, A. L.

2015-01-01

Cardiovascular simulation has shown potential value in clinical decision-making, providing a framework to assess changes in hemodynamics produced by physiological and surgical alterations. State-of-the-art predictions are provided by deterministic multiscale numerical approaches coupling 3D finite element Navier Stokes simulations to lumped parameter circulation models governed by ODEs. Development of next-generation stochastic multiscale models whose parameters can be learned from available clinical data under uncertainty constitutes a research challenge made more difficult by the high computational cost typically associated with the solution of these models. We present a methodology for constructing reduced representations that condense the behavior of 3D anatomical models using outlet pressure-flow polynomial surrogates, based on multiscale model solutions spanning several heart cycles. Relevance vector machine regression is compared with maximum likelihood estimation, showing that sparse pressure/flow rate approximations offer superior performance in producing working surrogate models to be included in lumped circulation networks. Sensitivities of outlets flow rates are also quantified through a Sobol’ decomposition of their total variance encoded in the orthogonal polynomial expansion. Finally, we show that augmented lumped parameter models including the proposed surrogates accurately reproduce the response of multiscale models they were derived from. In particular, results are presented for models of the coronary circulation with closed loop boundary conditions and the abdominal aorta with open loop boundary conditions. PMID:26671219

Modelling ecosystem service flows under uncertainty with stochiastic SPAN

USGS Publications Warehouse

Johnson, Gary W.; Snapp, Robert R.; Villa, Ferdinando; Bagstad, Kenneth J.

2012-01-01

Ecosystem service models are increasingly in demand for decision making. However, the data required to run these models are often patchy, missing, outdated, or untrustworthy. Further, communication of data and model uncertainty to decision makers is often either absent or unintuitive. In this work, we introduce a systematic approach to addressing both the data gap and the difficulty in communicating uncertainty through a stochastic adaptation of the Service Path Attribution Networks (SPAN) framework. The SPAN formalism assesses ecosystem services through a set of up to 16 maps, which characterize the services in a study area in terms of flow pathways between ecosystems and human beneficiaries. Although the SPAN algorithms were originally defined deterministically, we present them here in a stochastic framework which combines probabilistic input data with a stochastic transport model in order to generate probabilistic spatial outputs. This enables a novel feature among ecosystem service models: the ability to spatially visualize uncertainty in the model results. The stochastic SPAN model can analyze areas where data limitations are prohibitive for deterministic models. Greater uncertainty in the model inputs (including missing data) should lead to greater uncertainty expressed in the model’s output distributions. By using Bayesian belief networks to fill data gaps and expert-provided trust assignments to augment untrustworthy or outdated information, we can account for uncertainty in input data, producing a model that is still able to run and provide information where strictly deterministic models could not. Taken together, these attributes enable more robust and intuitive modelling of ecosystem services under uncertainty.

Improved spectral comparisons of paleoclimate models and observations via proxy system modeling: Implications for multi-decadal variability

NASA Astrophysics Data System (ADS)

Dee, S. G.; Parsons, L. A.; Loope, G. R.; Overpeck, J. T.; Ault, T. R.; Emile-Geay, J.

2017-10-01

The spectral characteristics of paleoclimate observations spanning the last millennium suggest the presence of significant low-frequency (multi-decadal to centennial scale) variability in the climate system. Since this low-frequency climate variability is critical for climate predictions on societally-relevant scales, it is essential to establish whether General Circulation models (GCMs) are able to simulate it faithfully. Recent studies find large discrepancies between models and paleoclimate data at low frequencies, prompting concerns surrounding the ability of GCMs to predict long-term, high-magnitude variability under greenhouse forcing (Laepple and Huybers, 2014a, 2014b). However, efforts to ground climate model simulations directly in paleoclimate observations are impeded by fundamental differences between models and the proxy data: proxy systems often record a multivariate and/or nonlinear response to climate, precluding a direct comparison to GCM output. In this paper we bridge this gap via a forward proxy modeling approach, coupled to an isotope-enabled GCM. This allows us to disentangle the various contributions to signals embedded in ice cores, speleothem calcite, coral aragonite, tree-ring width, and tree-ring cellulose. The paper addresses the following questions: (1) do forward-modeled ;pseudoproxies; exhibit variability comparable to proxy data? (2) if not, which processes alter the shape of the spectrum of simulated climate variability, and are these processes broadly distinguishable from climate? We apply our method to representative case studies, and broaden these insights with an analysis of the PAGES2k database (PAGES2K Consortium, 2013). We find that current proxy system models (PSMs) can help resolve model-data discrepancies on interannual to decadal timescales, but cannot account for the mismatch in variance on multi-decadal to centennial timescales. We conclude that, specific to this set of PSMs and isotope-enabled model, the paleoclimate record may exhibit larger low-frequency variability than GCMs currently simulate, indicative of incomplete physics and/or forcings.

In Search of the Physics: NASA's Approach to Airframe Noise

NASA Technical Reports Server (NTRS)

Macaraeg, Michele G.; Lockard, David P.; Streett, Craig L.

1999-01-01

An extensive numerical and experimental study of airframe noise mechanisms associated with a subsonic high-lift system has been performed at NASA Langley Research Center (LaRC). Investigations involving both steady and unsteady computations and experiments on small-scale models with part-span flaps and full-span flaps are presented. Both surface (steady and unsteady pressure measurements, hot films, oil flows, pressure sensitive paint) and off-surface (5 holeprobe, particle-imaged velocimetry, laser velocimetry, laser light sheet measurements) were taken in the LaRC Quiet Flow Facility (QFF) and several hard-wall tunnels. Experiments in the Low Turbulence Pressure Tunnel (LTPT) included Reynolds number variations up to flight conditions. Successful microphone array measurements were also taken providing both acoustic source maps on the model, and quantitative spectra. Critical directivity measurements were obtained in the QFF. NASA Langley unstructured and structured Reynolds-Averaged Navier-Stokes codes modeled the steady aspects of the flows. Excellent comparisons with surface and off-surface experimental data were obtained. Subsequently, these meanflow calculations were utilized in both linear stability and direct numerical simulations of the flow fields to calculate unsteady surface pressures and farfield acoustic spectra. Accurate calculations were critical in obtaining not only noise source characteristics, but shear layer correction data as well. Techniques utilized in these investigations as well as brief overviews of the results are given.

Effects of injection pressure variation on mixing in a cold supersonic combustor with kerosene fuel

NASA Astrophysics Data System (ADS)

Liu, Wei-Lai; Zhu, Lin; Qi, Yin-Yin; Ge, Jia-Ru; Luo, Feng; Zou, Hao-Ran; Wei, Min; Jen, Tien-Chien

2017-10-01

Spray jet in cold kerosene-fueled supersonic flow has been characterized under different injection pressures to assess the effects of the pressure variation on the mixing between incident shock wave and transverse cavity injection. Based on the real scramjet combustor, a detailed computational fluid dynamics model is developed. The injection pressures are specified as 0.5, 1.0, 2.0, 3.0 and 4.0 MPa, respectively, with the other constant operation parameters (such as the injection diameter, angle and velocity). A three dimensional Couple Level Set & Volume of Fluids approach incorporating an improved Kelvin-Helmholtz & Rayleigh-Taylor model is used to investigate the interaction between kerosene and supersonic air. The numerical simulations primarily concentrate on penetration depth, span expansion area, angle of shock wave and sauter mean diameter distribution of the kerosene droplets with/without evaporation. Validation has been implemented by comparing the calculated against the measured in literature with good qualitative agreement. Results show that the penetration depth, span-wise angle and expansion area of the transverse cavity jet are all increased with the injection pressure. However, when the injection pressure is further increased, the value in either penetration depth or expansion area increases appreciably. This study demonstrates the feasibility and effectiveness of the combination of Couple Level Set & Volume of Fluids approach and an improved Kelvin-Helmholtz & Rayleigh-Taylor model, in turn providing insights into scramjet design improvement.

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

NASA Astrophysics Data System (ADS)

Haarsma, Reindert J.; Roberts, Malcolm J.; Vidale, Pier Luigi; Senior, Catherine A.; Bellucci, Alessio; Bao, Qing; Chang, Ping; Corti, Susanna; Fučkar, Neven S.; Guemas, Virginie; von Hardenberg, Jost; Hazeleger, Wilco; Kodama, Chihiro; Koenigk, Torben; Leung, L. Ruby; Lu, Jian; Luo, Jing-Jia; Mao, Jiafu; Mizielinski, Matthew S.; Mizuta, Ryo; Nobre, Paulo; Satoh, Masaki; Scoccimarro, Enrico; Semmler, Tido; Small, Justin; von Storch, Jin-Song

2016-11-01

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950-2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, "what are the origins and consequences of systematic model biases?", but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

Investigation of index finger triggering force using a cadaver experiment: Effects of trigger grip span, contact location, and internal tendon force.

PubMed

Chang, Joonho; Freivalds, Andris; Sharkey, Neil A; Kong, Yong-Ku; Mike Kim, H; Sung, Kiseok; Kim, Dae-Min; Jung, Kihyo

2017-11-01

A cadaver study was conducted to investigate the effects of triggering conditions (trigger grip span, contact location, and internal tendon force) on index finger triggering force and the force efficiency of involved tendons. Eight right human cadaveric hands were employed, and a motion simulator was built to secure and control the specimens. Index finger triggering forces were investigated as a function of different internal tendon forces (flexor digitorum profundus + flexor digitorum superficialis = 40, 70, and 100 N), trigger grip spans (40, 50, and 60 mm), and contact locations between the index finger and a trigger. Triggering forces significantly increased when internal tendon forces increased from 40 to 100 N. Also, trigger grip spans and contact locations had significant effects on triggering forces; maximum triggering forces were found at a 50 mm span and the most proximal contact location. The results revealed that only 10-30% of internal tendon forces were converted to their external triggering forces. Copyright © 2017 Elsevier Ltd. All rights reserved.

Stratospheric O3 changes during 2001-2010: the small role of solar flux variations in a chemical transport model

NASA Astrophysics Data System (ADS)

Dhomse, S. S.; Chipperfield, M. P.; Feng, W.; Ball, W. T.; Unruh, Y. C.; Haigh, J. D.; Krivova, N. A.; Solanki, S. K.; Smith, A. K.

2013-10-01

Solar spectral fluxes (or irradiance) measured by the SOlar Radiation and Climate Experiment (SORCE) show different variability at ultraviolet (UV) wavelengths compared to other irradiance measurements and models (e.g. NRL-SSI, SATIRE-S). Some modelling studies have suggested that stratospheric/lower mesospheric O3 changes during solar cycle 23 (1996-2008) can only be reproduced if SORCE solar fluxes are used. We have used a 3-D chemical transport model (CTM), forced by meteorology from the European Centre for Medium-Range Weather Forecasts (ECMWF), to simulate middle atmospheric O3 using three different solar flux data sets (SORCE, NRL-SSI and SATIRE-S). Simulated O3 changes are compared with Microwave Limb Sounder (MLS) and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite data. Modelled O3 anomalies from all solar flux data sets show good agreement with the observations, despite the different flux variations. The off-line CTM reproduces these changes through dynamical information contained in the analyses. A notable feature during this period is a robust positive solar signal in the tropical middle stratosphere, which is due to realistic dynamical changes in our simulations. Ozone changes in the lower mesosphere cannot be used to discriminate between solar flux data sets due to large uncertainties and the short time span of the observations. Overall this study suggests that, in a CTM, the UV variations detected by SORCE are not necessary to reproduce observed stratospheric O3 changes during 2001-2010.

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

Dierickx, Marion I. P.; Loeb, Abraham, E-mail: mdierickx@cfa.harvard.edu, E-mail: aloeb@cfa.harvard.edu

The extensive span of the Sagittarius (Sgr) stream makes it a promising tool for studying the gravitational potential of the Milky Way (MW). Characterizing its stellar kinematics can constrain halo properties and provide a benchmark for the paradigm of galaxy formation from cold dark matter. Accurate models of the disruption dynamics of the Sgr progenitor are necessary to employ this tool. Using a combination of analytic modeling and N -body simulations, we build a new model of the Sgr orbit and resulting stellar stream. In contrast to previous models, we simulate the full infall trajectory of the Sgr progenitor frommore » the time it first crossed the MW virial radius 8 Gyr ago. An exploration of the parameter space of initial phase-space conditions yields tight constraints on the angular momentum of the Sgr progenitor. Our best-fit model is the first to accurately reproduce existing data on the 3D positions and radial velocities of the debris detected 100 kpc away in the MW halo. In addition to replicating the mapped stream, the simulation also predicts the existence of several arms of the Sgr stream extending to hundreds of kiloparsecs. The two most distant stars known in the MW halo coincide with the predicted structure. Additional stars in the newly predicted arms can be found with future data from the Large Synoptic Survey Telescope. Detecting a statistical sample of stars in the most distant Sgr arms would provide an opportunity to constrain the MW potential out to unprecedented Galactocentric radii.« less

A nonparametric statistical technique for combining global precipitation datasets: development and hydrological evaluation over the Iberian Peninsula

NASA Astrophysics Data System (ADS)

Abul Ehsan Bhuiyan, Md; Nikolopoulos, Efthymios I.; Anagnostou, Emmanouil N.; Quintana-Seguí, Pere; Barella-Ortiz, Anaïs

2018-02-01

This study investigates the use of a nonparametric, tree-based model, quantile regression forests (QRF), for combining multiple global precipitation datasets and characterizing the uncertainty of the combined product. We used the Iberian Peninsula as the study area, with a study period spanning 11 years (2000-2010). Inputs to the QRF model included three satellite precipitation products, CMORPH, PERSIANN, and 3B42 (V7); an atmospheric reanalysis precipitation and air temperature dataset; satellite-derived near-surface daily soil moisture data; and a terrain elevation dataset. We calibrated the QRF model for two seasons and two terrain elevation categories and used it to generate ensemble for these conditions. Evaluation of the combined product was based on a high-resolution, ground-reference precipitation dataset (SAFRAN) available at 5 km 1 h-1 resolution. Furthermore, to evaluate relative improvements and the overall impact of the combined product in hydrological response, we used the generated ensemble to force a distributed hydrological model (the SURFEX land surface model and the RAPID river routing scheme) and compared its streamflow simulation results with the corresponding simulations from the individual global precipitation and reference datasets. We concluded that the proposed technique could generate realizations that successfully encapsulate the reference precipitation and provide significant improvement in streamflow simulations, with reduction in systematic and random error on the order of 20-99 and 44-88 %, respectively, when considering the ensemble mean.

Recent Improvements in Semi-Span Testing at the National Transonic Facility (Invited)

NASA Technical Reports Server (NTRS)

Gatlin, G. M.; Tomek, W. G.; Payne, F. M.; Griffiths, R. C.

2006-01-01

Three wind tunnel investigations of a commercial transport, high-lift, semi-span configuration have recently been conducted in the National Transonic Facility at the NASA Langley Research Center. Throughout the course of these investigations multiple improvements have been developed in the facility semi-span test capability. The primary purpose of the investigations was to assess Reynolds number scale effects on a modern commercial transport configuration up to full-scale flight test conditions (Reynolds numbers on the order of 27 million). The tests included longitudinal aerodynamic studies at subsonic takeoff and landing conditions across a range of Reynolds numbers from that available in conventional wind tunnels up to flight conditions. The purpose of this paper is to discuss lessons learned and improvements incorporated into the semi-span testing process. Topics addressed include enhanced thermal stabilization and moisture reduction procedures, assessments and improvements in model sealing techniques, compensation of model reference dimensions due to test temperature, significantly improved semi-span model access capability, and assessments of data repeatability.

Model simulations of flood and debris flow timing in steep catchments after wildfire

NASA Astrophysics Data System (ADS)

Rengers, F. K.; McGuire, L. A.; Kean, J. W.; Staley, D. M.; Hobley, D. E. J.

2016-08-01

Debris flows are a typical hazard on steep slopes after wildfire, but unlike debris flows that mobilize from landslides, most postwildfire debris flows are generated from water runoff. The majority of existing debris flow modeling has focused on landslide-triggered debris flows. In this study we explore the potential for using process-based rainfall-runoff models to simulate the timing of water flow and runoff-generated debris flows in recently burned areas. Two different spatially distributed hydrologic models with differing levels of complexity were used: the full shallow water equations and the kinematic wave approximation. Model parameter values were calibrated in two different watersheds, spanning two orders of magnitude in drainage area. These watersheds were affected by the 2009 Station Fire in the San Gabriel Mountains, CA, USA. Input data for the numerical models were constrained by time series of soil moisture, flow stage, and rainfall collected at field sites, as well as high-resolution lidar-derived digital elevation models. The calibrated parameters were used to model a third watershed in the burn area, and the results show a good match with observed timing of flow peaks. The calibrated roughness parameter (Manning's n) was generally higher when using the kinematic wave approximation relative to the shallow water equations, and decreased with increasing spatial scale. The calibrated effective watershed hydraulic conductivity was low for both models, even for storms occurring several months after the fire, suggesting that wildfire-induced changes to soil-water infiltration were retained throughout that time. Overall, the two model simulations were quite similar suggesting that a kinematic wave model, which is simpler and more computationally efficient, is a suitable approach for predicting flood and debris flow timing in steep, burned watersheds.

Model simulations of flood and debris flow timing in steep catchments after wildfire

USGS Publications Warehouse

Rengers, Francis K.; McGuire, Luke; Kean, Jason W.; Staley, Dennis M.; Hobley, D.E.J

2016-01-01

Debris flows are a typical hazard on steep slopes after wildfire, but unlike debris flows that mobilize from landslides, most post-wildfire debris flows are generated from water runoff. The majority of existing debris-flow modeling has focused on landslide-triggered debris flows. In this study we explore the potential for using process-based rainfall-runoff models to simulate the timing of water flow and runoff-generated debris flows in recently burned areas. Two different spatially distributed hydrologic models with differing levels of complexity were used: the full shallow water equations and the kinematic wave approximation. Model parameter values were calibrated in two different watersheds, spanning two orders of magnitude in drainage area. These watersheds were affected by the 2009 Station Fire in the San Gabriel Mountains, CA, USA. Input data for the numerical models were constrained by time series of soil moisture, flow stage, and rainfall collected at field sites, as well as high-resolution lidar-derived digital elevation models. The calibrated parameters were used to model a third watershed in the burn area, and the results show a good match with observed timing of flow peaks. The calibrated roughness parameter (Manning's $n$) was generally higher when using the kinematic wave approximation relative to the shallow water equations, and decreased with increasing spatial scale. The calibrated effective watershed hydraulic conductivity was low for both models, even for storms occurring several months after the fire, suggesting that wildfire-induced changes to soil-water infiltration were retained throughout that time. Overall the two model simulations were quite similar suggesting that a kinematic wave model, which is simpler and more computationally efficient, is a suitable approach for predicting flood and debris flow timing in steep, burned watersheds.

Numerical simulation of earth fissures caused by overly aquifer exploitation at Guangming Village, China

NASA Astrophysics Data System (ADS)

Ye, S.; Franceschini, A.; Zhang, Y.; Janna, C.; Gong, X.; Yu, J.; Teatini, P.

2017-12-01

Earth fissures accompanying anthropogenic land subsidence due to overly aquifer exploitation create significant geohazards in China. In the framework of an efficient and safe management of groundwater, numerical models represent a unique scientific approach to predict the generation and development of earth fissures. However, the common geomechanical simulators fail to reproduce fissure development because, due to compatibility conditions, they cannot be effectively applied in discontinuous mechanics. We present an innovative modelling approach for the simulation of fissure development. Firstly, a regional 3D groundwater model is calibrated on available piezometric records; secondly, the regional outcome is used to define the boundary conditions of a local 3D groundwater model developed at the fissure scale and implementing a refined discretization of the local hydrogeologic setting; finally, the pressure change are used as forcing factor in a local 3D geomechanical model, which combines Finite Elements and Interface Elements to simulate the deformation of the continuous aquifer system and the generation and sliding/opening of earth fissures The approach has been applied to simulate the earth fissure at Guangming Village in Wuxi, China with land subsidence of more than 1 m caused by the overexploitation of the second confined aquifer. The first earth fissure was observed in 1998. It developed fast from 1998 to 2007. The domain addressed by the local simulations is 2 km wide and 5 km long. The thickness of the aquifer system ranges from 0 m, in the proximity of a mountain ridge southward, to 210 m northward and includes a phreatic aquifer, the first and second confined aquifers, and four aquitards. The simulations spanned the period from 1980, i.e. before the inception of large groundwater withdrawals, to 2015. The modelling results highlight that the earth fissures at Guangming Village have been caused by tension and shear, which developed from the land surface downward. The main factors contributing to fissure generation include the shallow bedrock with a sharp ridge, the asymmetric shape of ridge, the uneven thickness of the sedimentary deposits at the west and east side of the buried ridge together with a different pressure change that have been responsible for differential land subsidence.

A Nonlinear Dynamical Systems based Model for Stochastic Simulation of Streamflow

NASA Astrophysics Data System (ADS)

Erkyihun, S. T.; Rajagopalan, B.; Zagona, E. A.

2014-12-01

Traditional time series methods model the evolution of the underlying process as a linear or nonlinear function of the autocorrelation. These methods capture the distributional statistics but are incapable of providing insights into the dynamics of the process, the potential regimes, and predictability. This work develops a nonlinear dynamical model for stochastic simulation of streamflows. In this, first a wavelet spectral analysis is employed on the flow series to isolate dominant orthogonal quasi periodic timeseries components. The periodic bands are added denoting the 'signal' component of the time series and the residual being the 'noise' component. Next, the underlying nonlinear dynamics of this combined band time series is recovered. For this the univariate time series is embedded in a d-dimensional space with an appropriate lag T to recover the state space in which the dynamics unfolds. Predictability is assessed by quantifying the divergence of trajectories in the state space with time, as Lyapunov exponents. The nonlinear dynamics in conjunction with a K-nearest neighbor time resampling is used to simulate the combined band, to which the noise component is added to simulate the timeseries. We demonstrate this method by applying it to the data at Lees Ferry that comprises of both the paleo reconstructed and naturalized historic annual flow spanning 1490-2010. We identify interesting dynamics of the signal in the flow series and epochal behavior of predictability. These will be of immense use for water resources planning and management.

Numerical Investigation of Compressor Non-Synchronous Vibration with Full Annulus Rotor-Stator Interaction

NASA Astrophysics Data System (ADS)

Espinal, Daniel

The objective of this research is to investigate and confirm the periodicity of the Non-Synchronous Vibration (NSV) mechanism of a GE axial compressor with a full-annulus simulation. A second objective is to develop a high fidelity single-passage tool with time-accurate unsteady capabilities able to capture rotor-stator interactions and NSV excitation response. A high fidelity methodology for axial turbomachinery simulation is developed using the low diffusion shock-capturing Riemann solver with high order schemes, the Spalart-Allmaras turbulence closure model, the fully conservative unsteady sliding BC for rotor-stator interaction with extension to full-annulus and single-passage configurations, and the phase lag boundary conditions applied to rotor-stator interface and circumferential BC. A URANS solver is used and captures the NSV flow excitation frequency of 2439 Hz, which agrees reasonably well with the measured NSV frequency of 2600 Hz from strain gage test data. It is observed that the circumferentially traveling vortex formed in the vicinity of the rotor tip propagates at the speed of a non-engine order frequency and causes the NSV. The vortex travels along the suction surface of the blade and crosses the passage outlet near blade trailing edge. Such a vortex motion trajectory repeats in each blade passage and generates two low pressure regions due to the vortex core positions, one at the leading edge and one at the trailing edge, both are oscillating due to the vortex coming and leaving. These two low pressure regions create a pair of coupling forces that generates a torsion moment causing NSV. The full-annulus simulation shows that the circumferentially traveling vortex has fairly periodical behavior and is a full annulus structure. Also, frequencies below the NSV excitation frequency of 2439 Hz with large amplitudes in response to flow-separation related phenomena are present. This behavior is consistent with experimental measurements. For circumferentially averaged parameters like total pressure ratio, NSV is observed to have an effect, particularly at radial locations above 70% span. Therefore, to achieve similar or better total pressure ratio a design with a smaller loading of the upper blade span and a higher loading of the mid blade spans should be considered. A fully-conservative sliding interface boundary condition (BC) is implemented with phase-lag capabilities using the Direct Store method for single-passage simulations. Also Direct Store phase-lag was applied to the circumferential BCs to enforce longer disturbance wavelengths. The unsteady simulation using single-blade-passage with periodic BC for an inlet guide vane (IGV)-rotor configuration captures a 2291 Hz NSV excitation frequency and an IGV-rotor-stator configuration predicts a 2365 Hz NSV excitation frequency with a significantly higher amplitude above 90% span. This correlates closely to the predicted NSV excitation frequency of 2439 Hz for the full-annulus configuration. The two-blade-row configuration exhibits the same vortex structures captured in the full-annulus study. The three-blade-row configuration only captures a tip vortex shedding at the leading edge, which can be attributed to the reflective nature of the BCs causing IGV-rotor-stator interactions to be augmented, becoming dominant and shifting NSV excitation response to engine order regime. Phase-lag simulations with a Nodal Diameter (ND) of 5 is enforced for the circumferential BCs for the three-blade-row configuration, and the results exactly matched the frequency response and flow structures of the periodic simulation, illustrating the small effect that phase-lag has on strongly periodic flow disturbances. A ND of 7 is enforced at the sliding interface, however the NSV excitation completely disappears and only the wake propagation from IGV-Rotor-Stator interactions are captured. Rotor blade passage exhibits a circumferentially travelling vortex similar to those observed in the full-annulus and two-blade-row simulations. This can occur when the rotating instability responsible for the NSV no longer maintains a pressure variation with a characteristic frequency signature as it rotates relative to the rotor rotation, and now has become the beginning of a spike-type stall cell. In this scenario the travelling vortex has become evidence of part-stall of the upper spans of the rotor blade, but stalling is contained maintaining stable operation. In conclusion, an efficient method of capturing NSV excitation has been proposed in a high-fidelity manner, where only 2% of the computational resources used in a full-annulus simulation are required for an accurate single-blade-passage multi-stage simulation.

Constraints on High Northern Photosynthesis Increase Using Earth System Models and a Set of Independent Observations

NASA Astrophysics Data System (ADS)

Winkler, A. J.; Brovkin, V.; Myneni, R.; Alexandrov, G.

2017-12-01

Plant growth in the northern high latitudes benefits from increasing temperature (radiative effect) and CO2 fertilization as a consequence of rising atmospheric CO2 concentration. This enhanced gross primary production (GPP) is evident in large scale increase in summer time greening over the 36-year record of satellite observations. In this time period also various global ecosystem models simulate a greening trend in terms of increasing leaf area index (LAI). We also found a persistent greening trend analyzing historical simulations of Earth system models (ESM) participating in Phase 5 of the Coupled Model Intercomparison Project (CMIP5). However, these models span a large range in strength of the LAI trend, expressed as sensitivity to both key environmental factors, temperature and CO2 concentration. There is also a wide spread in magnitude of the associated increase of terrestrial GPP among the ESMs, which contributes to pronounced uncertainties in projections of future climate change. Here we demonstrate that there is a linear relationship across the CMIP5 model ensemble between projected GPP changes and historical LAI sensitivity, which allows using the observed LAI sensitivity as an "emerging constraint" on GPP estimation at future CO2 concentration. This constrained estimate of future GPP is substantially higher than the traditional multi-model mean suggesting that the majority of current ESMs may be significantly underestimating carbon fixation by vegetation in NHL. We provide three independent lines of evidence in analyzing observed and simulated CO2 amplitude as well as atmospheric CO2 inversion products to arrive at the same conclusion.

Control of information in working memory: Encoding and removal of distractors in the complex-span paradigm.

PubMed

Oberauer, Klaus; Lewandowsky, Stephan

2016-11-01

The article reports four experiments with complex-span tasks in which encoding of memory items alternates with processing of distractors. The experiments test two assumptions of a computational model of complex span, SOB-CS: (1) distractor processing impairs memory because distractors are encoded into working memory, thereby interfering with memoranda; and (2) free time following distractors is used to remove them from working memory by unbinding their representations from list context. Experiment 1 shows that distractors are erroneously chosen for recall more often than not-presented stimuli, demonstrating that distractors are encoded into memory. Distractor intrusions declined with longer free time, as predicted by distractor removal. Experiment 2 shows these effects even when distractors precede the memory list, ruling out an account based on selective rehearsal of memoranda during free time. Experiments 3 and 4 test the notion that distractors decay over time. Both experiments show that, contrary to the notion of distractor decay, the chance of a distractor intruding at test does not decline with increasing time since encoding of that distractor. Experiment 4 provides additional evidence against the prediction from distractor decay that distractor intrusions decline over an unfilled retention interval. Taken together, the results support SOB-CS and rule out alternative explanations. Data and simulation code are available on Open Science Framework: osf.io/3ewh7. Copyright © 2016 Elsevier B.V. All rights reserved.

A Novel Quasi-3D Method for Cascade Flow Considering Axial Velocity Density Ratio

NASA Astrophysics Data System (ADS)

Chen, Zhiqiang; Zhou, Ming; Xu, Quanyong; Huang, Xudong

2018-03-01

A novel quasi-3D Computational Fluid Dynamics (CFD) method of mid-span flow simulation for compressor cascades is proposed. Two dimension (2D) Reynolds-Averaged Navier-Stokes (RANS) method is shown facing challenge in predicting mid-span flow with a unity Axial Velocity Density Ratio (AVDR). Three dimension (3D) RANS solution also shows distinct discrepancies if the AVDR is not predicted correctly. In this paper, 2D and 3D CFD results discrepancies are analyzed and a novel quasi-3D CFD method is proposed. The new quasi-3D model is derived by reducing 3D RANS Finite Volume Method (FVM) discretization over a one-spanwise-layer structured mesh cell. The sidewall effect is considered by two parts. The first part is explicit interface fluxes of mass, momentum and energy as well as turbulence. The second part is a cell boundary scaling factor representing sidewall boundary layer contraction. The performance of the novel quasi-3D method is validated on mid-span pressure distribution, pressure loss and shock prediction of two typical cascades. The results show good agreement with the experiment data on cascade SJ301-20 and cascade AC6-10 at all test condition. The proposed quasi-3D method shows superior accuracy over traditional 2D RANS method and 3D RANS method in performance prediction of compressor cascade.

The relations between network-operation and topological-property in a scale-free and small-world network with community structure

NASA Astrophysics Data System (ADS)

Ma, Fei; Yao, Bing

2017-10-01

It is always an open, demanding and difficult task for generating available model to simulate dynamical functions and reveal inner principles from complex systems and networks. In this article, due to lots of real-life and artificial networks are built from series of simple and small groups (components), we discuss some interesting and helpful network-operation to generate more realistic network models. In view of community structure (modular topology), we present a class of sparse network models N(t , m) . At the moment, we capture the fact the N(t , 4) has not only scale-free feature, which means that the probability that a randomly selected vertex with degree k decays as a power-law, following P(k) ∼k-γ, where γ is the degree exponent, but also small-world property, which indicates that the typical distance between two uniform randomly chosen vertices grows proportionally to logarithm of the order of N(t , 4) , namely, relatively shorter diameter and lower average path length, simultaneously displays higher clustering coefficient. Next, as a new topological parameter correlating to reliability, synchronization capability and diffusion properties of networks, the number of spanning trees over a network is studied in more detail, an exact analytical solution for the number of spanning trees of the N(t , 4) is obtained. Based on the network-operation, part hub-vertex linking with each other will be helpful for structuring various network models and investigating the rules related with real-life networks.

Analysis of cavern and well stability at the West Hackberry SPR site using a full-dome model.

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

Sobolik, Steven R.

2015-08-01

This report presents computational analyses that simulate the structural response of caverns at the Strategic Petroleum Reserve (SPR) West Hackberry site. The cavern field comprises 22 caverns. Five caverns (6, 7, 8, 9, 11) were acquired from industry and have unusual shapes and a history dating back to 1946. The other 17 caverns (101-117) were leached according to SPR standards in the mid-1980s and have tall cylindrical shapes. The history of the caverns and their shapes are simulated in a three-dimensional geomechanics model of the site that predicts deformations, strains, and stresses. Future leaching scenarios corresponding to oil drawdowns usingmore » fresh water are also simulated by increasing the volume of the caverns. Cavern pressures are varied in the model to capture operational practices in the field. The results of the finite element model are interpreted to provide information on the current and future status of subsidence, well integrity, and cavern stability. The most significant results in this report are relevant to Cavern 6. The cavern is shaped like a bowl with a large ceiling span and is in close proximity to Cavern 9. The analyses predict tensile stresses at the edge of the ceiling during repressurization of Cavern 6 following workover conditions. During a workover the cavern is at low pressure to service a well. The wellhead pressures are atmospheric. When the workover is complete, the cavern is repressurized. The resulting elastic stresses are sufficient to cause tension around the edge of the large ceiling span. With time, these stresses relax to a compressive state because of salt creep. However, the potential for salt fracture and propagation exists, particularly towards Cavern 9. With only 200 feet of salt between the caverns, the operational consequences must be examined if the two caverns become connected. A critical time may be during a workover of Cavern 9 in part because of the operational vulnerabilities, but also because dilatant damage is predicted under the ledge that forms the lower lobe in the cavern. The remaining caverns have no significant issues regarding cavern stability and may be safely enlarged during subsequent oil drawdowns. Predicted well strains and subsidence are significant and consequently future remedial actions may be necessary. These predicted well strains certainly suggest appropriate monitoring through a well-logging program. Subsidence is currently being monitored.« less

A model linking immediate serial recall, the Hebb repetition effect and the learning of phonological word forms

PubMed Central

Page, M. P. A.; Norris, D.

2009-01-01

We briefly review the considerable evidence for a common ordering mechanism underlying both immediate serial recall (ISR) tasks (e.g. digit span, non-word repetition) and the learning of phonological word forms. In addition, we discuss how recent work on the Hebb repetition effect is consistent with the idea that learning in this task is itself a laboratory analogue of the sequence-learning component of phonological word-form learning. In this light, we present a unifying modelling framework that seeks to account for ISR and Hebb repetition effects, while being extensible to word-form learning. Because word-form learning is performed in the service of later word recognition, our modelling framework also subsumes a mechanism for word recognition from continuous speech. Simulations of a computational implementation of the modelling framework are presented and are shown to be in accordance with data from the Hebb repetition paradigm. PMID:19933143

Flight Investigation at High Speeds of the Drag of Three Airfoils and a Circular Cylinder Representing Full-Scale Propeller Shanks

NASA Technical Reports Server (NTRS)

Barlow, William H

1946-01-01

Tests have been made at high speeds to determine the drag of models, simulating propeller shanks, in the form of a circular cylinder and three airfoils, the NACA 16-025, the NACA 16-040, and the NACA 16-040 with the rear 25 percent chord cut off. All the models had a maximum thickness of 4 1/2 inches to conform with average propeller-shank dimensions and a span of 20 1/4 inches. For the tests the models were supported perpendicular to the lower surface of the wing of an XP-51 airplane. A wake-survey rake mounted below the wing directly behind the models was used to determine profile drag of Mach numbers of 0.3 to 0.8 over a small range of angle of attack. The drag of the cylinder was also determined from pressure-distribution and force measurements.

Scientific Visualization to Study Flux Transfer Events at the Community Coordinated Modeling Center

NASA Technical Reports Server (NTRS)

Rastatter, Lutz; Kuznetsova, Maria M.; Sibeck, David G.; Berrios, David H.

2011-01-01

In this paper we present results of modeling of reconnection at the dayside magnetopause with subsequent development of flux transfer event signatures. The tools used include new methods that have been added to the suite of visualization methods that are used at the Community Coordinated Modeling Center (CCMC). Flux transfer events result from localized reconnection that connect magnetosheath magnetic field and plasma with magnetospheric fields and plasma and results in flux rope structures that span the dayside magnetopause. The onset of flux rope formation and the three-dimensional structure of flux ropes are studied as they have been modeled by high-resolution magnetohydrodynamic simulations of the dayside magnetosphere of the Earth. We show that flux transfer events are complex three-dimensional structures that require modern visualization and analysis techniques. Two suites of visualization methods are presented and we demonstrate the usefulness of those methods through the CCMC web site to the general science user.

The kinetics of heterogeneous nucleation and growth: an approach based on a grain explicit model

NASA Astrophysics Data System (ADS)

Rouet-Leduc, B.; Maillet, J.-B.; Denoual, C.

2014-04-01

A model for phase transitions initiated on grain boundaries is proposed and tested against numerical simulations: this approach, based on a grain explicit model, allows us to consider the granular structure, resulting in accurate predictions for a wide span of nucleation processes. Comparisons are made with classical models of homogeneous (JMAK: Johnson and Mehl 1939 Trans. Am. Inst. Min. Eng. 135 416; Avrami 1939 J. Chem. Phys. 7 1103; Kolmogorov 1937 Bull. Acad. Sci. USSR, Mat. Ser. 1 335) as well as heterogeneous (Cahn 1996 Thermodynamics and Kinetics of Phase Transformations Im et al (Pittsburgh: Materials Research Society)) nucleation. A transition scale based on material properties is proposed, allowing us to discriminate between random and site-saturated regimes. Finally, we discuss the relationship between an Avrami-type exponent and the transition regime, establishing conditions for its extraction from experiments.

Synthetic River Valleys

NASA Astrophysics Data System (ADS)

Brown, R.; Pasternack, G. B.

2011-12-01

The description of fluvial form has evolved from anecdotal descriptions to artistic renderings to 2D plots of cross section or longitudinal profiles and more recently 3D digital models. Synthetic river valleys, artificial 3D topographic models of river topography, have a plethora of potential applications in fluvial geomorphology, and the earth sciences in general, as well as in computer science and ecology. Synthetic river channels have existed implicitly since approximately the 1970s and can be simulated from a variety of approaches spanning the artistic and numerical. An objective method of synthesizing 3D stream topography based on reach scale attributes would be valuable for sizing 3D flumes in the physical and numerical realms, as initial input topography for morphodynamic models, stream restoration design, historical reconstruction, and mechanistic testing of interactions of channel geometric elements. Quite simply - simulation of synthetic channel geometry of prescribed conditions can allow systematic evaluation of the dominant relationships between river flow and geometry. A new model, the control curve method, is presented that uses hierarchically scaled parametric curves in over-lapping 2D planes to create synthetic river valleys. The approach is able to simulate 3D stream geometry from paired 2D descriptions and can allow experimental insight into form-process relationships in addition to visualizing past measurements of channel form that are limited to two dimension descriptions. Results are presented that illustrate the models ability to simulate fluvial topography representative of real world rivers as well as how channel geometric elements can be adjusted. The testing of synthetic river valleys would open up a wealth of knowledge as to why some 3D attributes of river channels are more prevalent than others as well as bridging the gap between the 2D descriptions that have dominated fluvial geomorphology the past century and modern, more complete, 3D treatments.

Percolation of binary disk systems: Modeling and theory

DOE PAGES

Meeks, Kelsey; Tencer, John; Pantoya, Michelle L.

2017-01-12

The dispersion and connectivity of particles with a high degree of polydispersity is relevant to problems involving composite material properties and reaction decomposition prediction and has been the subject of much study in the literature. This paper utilizes Monte Carlo models to predict percolation thresholds for a two-dimensional systems containing disks of two different radii. Monte Carlo simulations and spanning probability are used to extend prior models into regions of higher polydispersity than those previously considered. A correlation to predict the percolation threshold for binary disk systems is proposed based on the extended dataset presented in this work and comparedmore » to previously published correlations. Finally, a set of boundary conditions necessary for a good fit is presented, and a condition for maximizing percolation threshold for binary disk systems is suggested.« less

Calibration Results and Anticipated Science Ops for the Parker Solar Probe's SWEAP/SPAN-Electron Analyzers

NASA Astrophysics Data System (ADS)

Whittlesey, P. L.; Larson, D. E.; Livi, R.; Abiad, R.; Parker, C.; Halekas, J. S.; Kasper, J. C.; Korreck, K. E.

2017-12-01

We present the SPAN-E calibration results and science operation plans this instrument on the Parker Solar Probe mission. SPAN-E is a pair of highly configurable ESA sensors, one on the RAM side of the spacecraft (SPAN-Ae) and one on anti-RAM (SPAN-B). Together, SPAN-E will jointly measure the full 3D thermal and suprathermal electron distribution function at cadences as fast as 4.58Hz. Joined with the SPAN-Ai and SPC instruments that are part of the Solar Wind Electrons, Alphas, and Protons (SWEAP) suite, SPAN-E will measure the solar coronal plasma across a range of energies and densities with a FOV over >90% of the sky, returning data over a 7 year long PSP mission lifetime. The SPAN-E instruments have passed environmental testing at the instrument level, and the final instrument calibrations are complete. This presentation details the final instrument calibration results as performed at UCB/SSL after environmental testing, and details the planned configurations for PSP's first orbit. In addition, the PSP spacecraft's magnetic fields are expected to distort the measured electron VDFs at low energies, thus we present a novel computer vision method of measuring and modeling the spacecraft magnetic fields as seen during an observatory-level "swing" test. Ultimately, the model will feed into an algorithm for ground corrections to electron VDFs distorted by these stray spacecraft magnetic fields.

Simulation of the Impact of Climate Variability on Malaria Transmission in the Sahel

NASA Astrophysics Data System (ADS)

Bomblies, A.; Eltahir, E.; Duchemin, J.

2007-12-01

A coupled hydrology and entomology model for simulation of malaria transmission and malaria transmitting mosquito population dynamics is presented. Model development and validation is done using field data and observations collected at Banizoumbou and Zindarou, Niger spanning three wet seasons, from 2005 through 2007. The primary model objective is the accurate determination of climate variability effects on village scale malaria transmission. Malaria transmission dependence on climate variables is highly nonlinear and complex. Temperature and humidity affect mosquito longevity, temperature controls parasite development rates in the mosquito as well as subadult mosquito development rates, and precipitation determines the formation and persistence of adequate breeding pools. Moreover, unsaturated zone hydrology influences overland flow, and climate controlled evapotranspiration rates and root zone uptake therefore also influence breeding pool formation. High resolution distributed hydrologic simulation allows representation of the small-scale ephemeral pools that constitute the primary habitat of Anopheles gambiae mosquitoes, the dominant malaria vectors in the Niger Sahel. Remotely sensed soil type, vegetation type, and microtopography rasters are used to assign the distributed parameter fields for simulation of the land surface hydrologic response to precipitation and runoff generation. Predicted runoff from each cell flows overland and into topographic depressions, with explicit representation of infiltration and evapotranspiration. The model's entomology component interacts with simulated pools. Subadult (aquatic stage) mosquito breeding is simulated in the pools, and water temperature dependent stage advancement rates regulate adult mosquito emergence into the model domain. Once emerged, adult mosquitoes are tracked as independent individual agents that interact with their immediate environment. Attributes relevant to malaria transmission such as gonotrophic state, infected and infectious states, age, and location relative to human population are tracked for each individual. The model operates at a resolution consistent with the characteristic scale of relevant ecological processes. Microhabitat exploitation and spatial structure of the mosquito population surrounding villages is reproduced in this manner. The resulting coupled model predicts not only malaria transmission's response to interannual climate variability, but can also evaluate land use change effects on malaria transmission. The late Professor Andrew Spielman of the Harvard School of Public Health provided medical entomology expertise and was a part of this effort.

Evaluating the skills of isotope-enabled general circulation models against in situ atmospheric water vapor isotope observations

NASA Astrophysics Data System (ADS)

Steen-Larsen, H. C.; Risi, C.; Werner, M.; Yoshimura, K.; Masson-Delmotte, V.

2017-01-01

The skills of isotope-enabled general circulation models are evaluated against atmospheric water vapor isotopes. We have combined in situ observations of surface water vapor isotopes spanning multiple field seasons (2010, 2011, and 2012) from the top of the Greenland Ice Sheet (NEEM site: 77.45°N, 51.05°W, 2484 m above sea level) with observations from the marine boundary layer of the North Atlantic and Arctic Ocean (Bermuda Islands 32.26°N, 64.88°W, year: 2012; south coast of Iceland 63.83°N, 21.47°W, year: 2012; South Greenland 61.21°N, 47.17°W, year: 2012; Svalbard 78.92°N, 11.92°E, year: 2014). This allows us to benchmark the ability to simulate the daily water vapor isotope variations from five different simulations using isotope-enabled general circulation models. Our model-data comparison documents clear isotope biases both on top of the Greenland Ice Sheet (1-11‰ for δ18O and 4-19‰ for d-excess depending on model and season) and in the marine boundary layer (maximum differences for the following: Bermuda δ18O = 1‰, d-excess = 3‰; South coast of Iceland δ18O = 2‰, d-excess = 5‰; South Greenland δ18O = 4‰, d-excess = 7‰; Svalbard δ18O = 2‰, d-excess = 7‰). We find that the simulated isotope biases are not just explained by simulated biases in temperature and humidity. Instead, we argue that these isotope biases are related to a poor simulation of the spatial structure of the marine boundary layer water vapor isotopic composition. Furthermore, we specifically show that the marine boundary layer water vapor isotopes of the Baffin Bay region show strong influence on the water vapor isotopes at the NEEM deep ice core-drilling site in northwest Greenland. Our evaluation of the simulations using isotope-enabled general circulation models also documents wide intermodel spatial variability in the Arctic. This stresses the importance of a coordinated water vapor isotope-monitoring network in order to discriminate amongst these model behaviors.

Quasi-Dynamic Versus Fully-Dynamic Simulations of Slip Accumulation on Faults with Enhanced Dynamic Weakening

NASA Astrophysics Data System (ADS)

Lapusta, N.; Thomas, M.; Noda, H.; Avouac, J.

2012-12-01

Long-term simulations that incorporate both seismic events and aseismic slip are quite important for studies of earthquake physics but challenging computationally. To study long deformation histories, most simulation methods do not incorporate full inertial effects (wave propagation) during simulated earthquakes, using quasi-dynamic approximations instead. Here we compare the results of quasi-dynamic simulations to the fully dynamic ones for a range of problems to determine the applicability of the quasi-dynamic approach. Intuitively, the quasi-dynamic approach should do relatively well in problems where wave-mediated effects are relatively simple but should have substantially different (and hence wrong) response when the wave-mediated stress transfers dominate the character of the seismic events. This is exactly what we observe in our simulations. We consider a 2D model of a rate-and-state fault with a seismogenic (steady-state velocity-weakening) zone surrounded by creeping (steady-state velocity-strengthening) areas. If the seismogenic zone is described by the standard Dieterich-Ruina rate-and-state friction, the resulting earthquake sequences consist of relatively simple crack-like ruptures, and the inclusion of true wave-propagation effects mostly serves to concentrate stress more efficiently at the rupture front. Hence, in such models, rupture speeds and slip rates are significantly (several times) lower in the quasi-dynamic simulations compared to the fully dynamic ones, but the total slip, the crack-like nature of seismic events, and the overall pattern of earthquake sequences is comparable, consistently with prior studies. Such behavior can be classified as qualitatively similar but quantitatively different, and it motivates the popularity of the quasi-dynamic methods in simulations. However, the comparison changes dramatically once we consider a model with enhanced dynamic weakening in the seismogenic zone in the form of flash heating. In this case, the fully dynamic simulations produce seismic ruptures in the form of short-duration slip pulses, where the pulses form due to a combination of enhanced weakening and wave effects. The quasi-dynamic simulations in the same model produce completely different results, with large crack-like ruptures, different total slips, different rupture patterns, and different prestress state before large, model-spanning events. Such qualitative differences between the quasi-dynamic and fully-dynamic simulation should result in any model where inertial effects lead to qualitative differences, such as cases with supershear transition or fault with different materials on the two sides. We will present results on our current work on how the quasi-dynamic and fully dynamic simulations compare for the cases with heterogeneous fault properties.

Prediction of shallow landslide occurrence: Validation of a physically-based approach through a real case study.

PubMed

Schilirò, Luca; Montrasio, Lorella; Scarascia Mugnozza, Gabriele

2016-11-01

In recent years, physically-based numerical models have frequently been used in the framework of early-warning systems devoted to rainfall-induced landslide hazard monitoring and mitigation. For this reason, in this work we describe the potential of SLIP (Shallow Landslides Instability Prediction), a simplified physically-based model for the analysis of shallow landslide occurrence. In order to test the reliability of this model, a back analysis of recent landslide events occurred in the study area (located SW of Messina, northeastern Sicily, Italy) on October 1st, 2009 was performed. The simulation results have been compared with those obtained for the same event by using TRIGRS, another well-established model for shallow landslide prediction. Afterwards, a simulation over a 2-year span period has been performed for the same area, with the aim of evaluating the performance of SLIP as early warning tool. The results confirm the good predictive capability of the model, both in terms of spatial and temporal prediction of the instability phenomena. For this reason, we recommend an operating procedure for the real-time definition of shallow landslide triggering scenarios at the catchment scale, which is based on the use of SLIP calibrated through a specific multi-methodological approach. Copyright © 2016 Elsevier B.V. All rights reserved.

Hypersonic separated flows about "tick" configurations with sensitivity to model design

NASA Astrophysics Data System (ADS)

Moss, J. N.; O'Byrne, S.; Gai, S. L.

2014-12-01

This paper presents computational results obtained by applying the direct simulation Monte Carlo (DSMC) method for hypersonic nonequilibrium flow about "tick-shaped" model configurations. These test models produces a complex flow where the nonequilibrium and rarefied aspects of the flow are initially enhanced as the flow passes over an expansion surface, and then the flow encounters a compression surface that can induce flow separation. The resulting flow is such that meaningful numerical simulations must have the capability to account for a significant range of rarefaction effects; hence the application of the DSMC method in the current study as the flow spans several flow regimes, including transitional, slip, and continuum. The current focus is to examine the sensitivity of both the model surface response (heating, friction and pressure) and flowfield structure to assumptions regarding surface boundary conditions and more extensively the impact of model design as influenced by leading edge configuration as well as the geometrical features of the expansion and compression surfaces. Numerical results indicate a strong sensitivity to both the extent of the leading edge sharpness and the magnitude of the leading edge bevel angle. Also, the length of the expansion surface for a fixed compression surface has a significant impact on the extent of separated flow.

Emulation and Sensitivity Analysis of the Community Multiscale Air Quality Model for a UK Ozone Pollution Episode.

PubMed

Beddows, Andrew V; Kitwiroon, Nutthida; Williams, Martin L; Beevers, Sean D

2017-06-06

Gaussian process emulation techniques have been used with the Community Multiscale Air Quality model, simulating the effects of input uncertainties on ozone and NO 2 output, to allow robust global sensitivity analysis (SA). A screening process ranked the effect of perturbations in 223 inputs, isolating the 30 most influential from emissions, boundary conditions (BCs), and reaction rates. Community Multiscale Air Quality (CMAQ) simulations of a July 2006 ozone pollution episode in the UK were made with input values for these variables plus ozone dry deposition velocity chosen according to a 576 point Latin hypercube design. Emulators trained on the output of these runs were used in variance-based SA of the model output to input uncertainties. Performing these analyses for every hour of a 21 day period spanning the episode and several days on either side allowed the results to be presented as a time series of sensitivity coefficients, showing how the influence of different input uncertainties changed during the episode. This is one of the most complex models to which these methods have been applied, and here, they reveal detailed spatiotemporal patterns of model sensitivities, with NO and isoprene emissions, NO 2 photolysis, ozone BCs, and deposition velocity being among the most influential input uncertainties.

Hypersonic Separated Flows About "Tick" Configurations With Sensitivity to Model Design

NASA Technical Reports Server (NTRS)

Moss, J. N.; O'Byrne, S.; Gai, S. L.

2014-01-01

This paper presents computational results obtained by applying the direct simulation Monte Carlo (DSMC) method for hypersonic nonequilibrium flow about "tick-shaped" model configurations. These test models produces a complex flow where the nonequilibrium and rarefied aspects of the flow are initially enhanced as the flow passes over an expansion surface, and then the flow encounters a compression surface that can induce flow separation. The resulting flow is such that meaningful numerical simulations must have the capability to account for a significant range of rarefaction effects; hence the application of the DSMC method in the current study as the flow spans several flow regimes, including transitional, slip, and continuum. The current focus is to examine the sensitivity of both the model surface response (heating, friction and pressure) and flowfield structure to assumptions regarding surface boundary conditions and more extensively the impact of model design as influenced by leading edge configuration as well as the geometrical features of the expansion and compression surfaces. Numerical results indicate a strong sensitivity to both the extent of the leading edge sharpness and the magnitude of the leading edge bevel angle. Also, the length of the expansion surface for a fixed compression surface has a significant impact on the extent of separated flow.

Jet-boundary and Plan-form Corrections for Partial-Span Models with Reflection-Plane, End-Plate, or No End-Plate in a Closed Circular Wind Tunnel

NASA Technical Reports Server (NTRS)

Sivells, James C; Deters, Owen J

1946-01-01

A method is presented for determining the jet-boundary and plan-form corrections necessary for application to test data for a partial-span model with a reflection plane, an end plate, or no end plate in a closed circular wind tunnel. Examples are worked out for a partial-span model with each of the three end conditions in the Langley 19-foot pressure tunnel and the corrections are applied to measured values of lift, drag, pitching-moment, rolling-moment, and yawing-moment coefficients.

Do dynamic global vegetation models capture the seasonality of carbon fluxes in the Amazon basin? A data-model intercomparison

DOE PAGES

Restrepo-Coupe, Natalia; Levine, Naomi M.; Christoffersen, Bradley O.; ...

2016-08-29

To predict forest response to long-term climate change with high confidence requires that dynamic global vegetation models (DGVMs) be successfully tested against ecosystem response to short-term variations in environmental drivers, including regular seasonal patterns. Here, we used an integrated dataset from four forests in the Brasil flux network, spanning a range of dry-season intensities and lengths, to determine how well four state-of-the-art models (IBIS, ED2, JULES, and CLM3.5) simulated the seasonality of carbon exchanges in Amazonian tropical forests. We found that most DGVMs poorly represented the annual cycle of gross primary productivity (GPP), of photosynthetic capacity (Pc), and of othermore » fluxes and pools. Models simulated consistent dry-season declines in GPP in the equatorial Amazon (Manaus K34, Santarem K67, and Caxiuanã CAX); a contrast to observed GPP increases. Model simulated dry-season GPP reductions were driven by an external environmental factor, ‘soil water stress’ and consequently by a constant or decreasing photosynthetic infrastructure (Pc), while observed dry-season GPP resulted from a combination of internal biological (leaf-flush and abscission and increased Pc) and environmental (incoming radiation) causes. Moreover, we found models generally overestimated observed seasonal net ecosystem exchange (NEE) and respiration (Re) at equatorial locations. In contrast, a southern Amazon forest (Jarú RJA) exhibited dry-season declines in GPP and Re consistent with most DGVMs simulations. While water limitation was represented in models and the primary driver of seasonal photosynthesis in southern Amazonia, changes in internal biophysical processes, light-harvesting adaptations (e.g., variations in leaf area index (LAI) and increasing leaf-level assimilation rate related to leaf demography), and allocation lags between leaf and wood, dominated equatorial Amazon carbon flux dynamics and were deficient or absent from current model formulations. In conclusion, correctly simulating flux seasonality at tropical forests requires a greater understanding and the incorporation of internal biophysical mechanisms in future model developments.« less

Do dynamic global vegetation models capture the seasonality of carbon fluxes in the Amazon basin? A data-model intercomparison

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

Restrepo-Coupe, Natalia; Levine, Naomi M.; Christoffersen, Bradley O.

To predict forest response to long-term climate change with high confidence requires that dynamic global vegetation models (DGVMs) be successfully tested against ecosystem response to short-term variations in environmental drivers, including regular seasonal patterns. Here, we used an integrated dataset from four forests in the Brasil flux network, spanning a range of dry-season intensities and lengths, to determine how well four state-of-the-art models (IBIS, ED2, JULES, and CLM3.5) simulated the seasonality of carbon exchanges in Amazonian tropical forests. We found that most DGVMs poorly represented the annual cycle of gross primary productivity (GPP), of photosynthetic capacity (Pc), and of othermore » fluxes and pools. Models simulated consistent dry-season declines in GPP in the equatorial Amazon (Manaus K34, Santarem K67, and Caxiuanã CAX); a contrast to observed GPP increases. Model simulated dry-season GPP reductions were driven by an external environmental factor, ‘soil water stress’ and consequently by a constant or decreasing photosynthetic infrastructure (Pc), while observed dry-season GPP resulted from a combination of internal biological (leaf-flush and abscission and increased Pc) and environmental (incoming radiation) causes. Moreover, we found models generally overestimated observed seasonal net ecosystem exchange (NEE) and respiration (Re) at equatorial locations. In contrast, a southern Amazon forest (Jarú RJA) exhibited dry-season declines in GPP and Re consistent with most DGVMs simulations. While water limitation was represented in models and the primary driver of seasonal photosynthesis in southern Amazonia, changes in internal biophysical processes, light-harvesting adaptations (e.g., variations in leaf area index (LAI) and increasing leaf-level assimilation rate related to leaf demography), and allocation lags between leaf and wood, dominated equatorial Amazon carbon flux dynamics and were deficient or absent from current model formulations. In conclusion, correctly simulating flux seasonality at tropical forests requires a greater understanding and the incorporation of internal biophysical mechanisms in future model developments.« less

Finite Elements Analysis of a Composite Semi-Span Test Article With and Without Discrete Damage

NASA Technical Reports Server (NTRS)

Lovejoy, Andrew E.; Jegley, Dawn C. (Technical Monitor)

2000-01-01

AS&M Inc. performed finite element analysis, with and without discrete damage, of a composite semi-span test article that represents the Boeing 220-passenger transport aircraft composite semi-span test article. A NASTRAN bulk data file and drawings of the test mount fixtures and semi-span components were utilized to generate the baseline finite element model. In this model, the stringer blades are represented by shell elements, and the stringer flanges are combined with the skin. Numerous modeling modifications and discrete source damage scenarios were applied to the test article model throughout the course of the study. This report details the analysis method and results obtained from the composite semi-span study. Analyses were carried out for three load cases: Braked Roll, LOG Down-Bending and 2.5G Up-Bending. These analyses included linear and nonlinear static response, as well as linear and nonlinear buckling response. Results are presented in the form of stress and strain plots. factors of safety for failed elements, buckling loads and modes, deflection prediction tables and plots, and strainage prediction tables and plots. The collected results are presented within this report for comparison to test results.

A Quantitative Study on the Correlation between Grade Span Configuration of Sixth Grade Students in Private Florida Schools and Academic Achievement on Standardized Achievement Scores

ERIC Educational Resources Information Center

Rantin, Deborah

2017-01-01

The applied dissertation was designed to investigate the three models of grade span configurations of sixth grade and the effects grade span configuration has on results of the standardized achievement scores of sixth grade students in private, Florida schools. Studies that have been conducted on sixth grade students and grade span configuration…

Study of Semi-Span Model Testing Techniques

NASA Technical Reports Server (NTRS)

Gatlin, Gregory M.; McGhee, Robert J.

1996-01-01

An investigation has been conducted in the NASA Langley 14- by 22-Foot Subsonic Tunnel in order to further the development of semi-span testing capabilities. A twin engine, energy efficient transport (EET) model with a four-element wing in a takeoff configuration was used for this investigation. Initially a full span configuration was tested and force and moment data, wing and fuselage surface pressure data, and fuselage boundary layer measurements were obtained as a baseline data set. The semi-span configurations were then mounted on the wind tunnel floor, and the effects of fuselage standoff height and shape as well as the effects of the tunnel floor boundary layer height were investigated. The effectiveness of tangential blowing at the standoff/floor juncture as an active boundary-layer control technique was also studied. Results indicate that the semi-span configuration was more sensitive to variations in standoff height than to variations in floor boundary layer height. A standoff height equivalent to 30 percent of the fuselage radius resulted in better correlation with full span data than no standoff or the larger standoff configurations investigated. Undercut standoff leading edges or the use of tangential blowing in the standoff/ floor juncture improved correlation of semi-span data with full span data in the region of maximum lift coefficient.

Investigation of Seal-to-Floor Effects on Semi-Span Transonic Models

NASA Technical Reports Server (NTRS)

Sleppy, Mark A.; Engel, Eric A.; Watson, Kevin T.; Atler, Douglas M.

2009-01-01

In an effort to achieve the maximum possible Reynolds number (Re) when conducting production testing for flight loads aerodynamic databases, it has been the preferred practice of The Boeing Company / Commercial Airplanes (BCA) -- Loads and Dynamics Group since the early 1990's to test large scale semi-span models in the 11- By 11-Foot Transonic Wind Tunnel (TWT) leg of the Unitary Plan Wind Tunnel (UPWT) at the NASA Ames Research Center (ARC). There are many problems related to testing large scale semi-span models of high aspect ratio flexible transport wings, such as; floor boundary layer effects, wing spanwise wall effects, solid blockage buoyancy effects, floor mechanical interference effects, airflow under the model effects, or tunnel flow gradient effects. For most of these issues, BCA has developed and implemented either standard testing methods or numerical correction schemes and these will not be discussed in this document. Other researchers have reported on semi-span transonic testing correction issues, however most of the reported research has been for low Mach testing. Some of the reports for low Mach testing address the difficult problem of preventing undesirable airflow under a semi-span model while ensuring unrestricted main balance functionality, however, for transonic models this issue has gone unresolved. BCA has been cognizant for sometime that there are marked differences in wing pressure distributions from semi-span transonic model testing than from full model or flight testing. It has been suspected that these differences are at least in part due to airflow under the model. Previous efforts by BCA to address this issue have proven to be ineffective or inconclusive and in one situation resulted in broken hardware. This paper reports on a Boeing-NASA collaborative investigation based on a series of small tests conducted between June 2006 and November 2007 in the 11 by 11 foot Transonic Wind Tunnel at NASA Ames on three large commercial jet transport configurations to assess the effects of sealing a semi-span model to the floor and to investigate efficient sealing and testing techniques. This document will show how sealing the model to the floor has a small but remarkably far reaching spanwise effect on wing pressures, wing local section forces and wing force summations.

Finite element model updating of multi-span steel-arch-steel-girder bridges based on ambient vibrations

NASA Astrophysics Data System (ADS)

Hou, Tsung-Chin; Gao, Wei-Yuan; Chang, Chia-Sheng; Zhu, Guan-Rong; Su, Yu-Min

2017-04-01

The three-span steel-arch-steel-girder Jiaxian Bridge was newly constructed in 2010 to replace the former one that has been destroyed by Typhoon Sinlaku (2008, Taiwan). It was designed and built to continue the domestic service requirement, as well as to improve the tourism business of the Kaohsiung city government, Taiwan. This study aimed at establishing the baseline model of Jiaxian Bridge for hazardous scenario simulation such as typhoons, floods and earthquakes. Necessities of these precaution works were attributed to the inherent vulnerability of the sites: near fault and river cross. The uncalibrated baseline bridge model was built with structural finite element in accordance with the blueprints. Ambient vibration measurements were performed repeatedly to acquire the elastic dynamic characteristics of the bridge structure. Two frequency domain system identification algorithms were employed to extract the measured operational modal parameters. Modal shapes, frequencies, and modal assurance criteria (MAC) were configured as the fitting targets so as to calibrate/update the structural parameters of the baseline model. It has been recognized that different types of structural parameters contribute distinguishably to the fitting targets, as this study has similarly explored. For steel-arch-steel-girder bridges in particular this case, joint rigidity of the steel components was found to be dominant while material properties and section geometries relatively minor. The updated model was capable of providing more rational elastic responses of the bridge superstructure under normal service conditions as well as hazardous scenarios, and can be used for manage the health conditions of the bridge structure.

Coarse-grained stochastic processes and kinetic Monte Carlo simulators for the diffusion of interacting particles

NASA Astrophysics Data System (ADS)

Katsoulakis, Markos A.; Vlachos, Dionisios G.

2003-11-01

We derive a hierarchy of successively coarse-grained stochastic processes and associated coarse-grained Monte Carlo (CGMC) algorithms directly from the microscopic processes as approximations in larger length scales for the case of diffusion of interacting particles on a lattice. This hierarchy of models spans length scales between microscopic and mesoscopic, satisfies a detailed balance, and gives self-consistent fluctuation mechanisms whose noise is asymptotically identical to the microscopic MC. Rigorous, detailed asymptotics justify and clarify these connections. Gradient continuous time microscopic MC and CGMC simulations are compared under far from equilibrium conditions to illustrate the validity of our theory and delineate the errors obtained by rigorous asymptotics. Information theory estimates are employed for the first time to provide rigorous error estimates between the solutions of microscopic MC and CGMC, describing the loss of information during the coarse-graining process. Simulations under periodic boundary conditions are used to verify the information theory error estimates. It is shown that coarse-graining in space leads also to coarse-graining in time by q2, where q is the level of coarse-graining, and overcomes in part the hydrodynamic slowdown. Operation counting and CGMC simulations demonstrate significant CPU savings in continuous time MC simulations that vary from q3 for short potentials to q4 for long potentials. Finally, connections of the new coarse-grained stochastic processes to stochastic mesoscopic and Cahn-Hilliard-Cook models are made.

Simulating the elimination of sleeping sickness with an agent-based model.

PubMed

Grébaut, Pascal; Girardin, Killian; Fédérico, Valentine; Bousquet, François

2016-01-01

Although Human African Trypanosomiasis is largely considered to be in the process of extinction today, the persistence of human and animal reservoirs, as well as the vector, necessitates a laborious elimination process. In this context, modeling could be an effective tool to evaluate the ability of different public health interventions to control the disease. Using the Cormas ® system, we developed HATSim, an agent-based model capable of simulating the possible endemic evolutions of sleeping sickness and the ability of National Control Programs to eliminate the disease. This model takes into account the analysis of epidemiological, entomological, and ecological data from field studies conducted during the last decade, making it possible to predict the evolution of the disease within this area over a 5-year span. In this article, we first present HATSim according to the Overview, Design concepts, and Details (ODD) protocol that is classically used to describe agent-based models, then, in a second part, we present predictive results concerning the evolution of Human African Trypanosomiasis in the village of Lambi (Cameroon), in order to illustrate the interest of such a tool. Our results are consistent with what was observed in the field by the Cameroonian National Control Program (CNCP). Our simulations also revealed that regular screening can be sufficient, although vector control applied to all areas with human activities could be significantly more efficient. Our results indicate that the current model can already help decision-makers in planning the elimination of the disease in foci. © P. Grébaut et al., published by EDP Sciences, 2016.

Influence of Tooth Spacing Error on Gears With and Without Profile Modifications

NASA Technical Reports Server (NTRS)

Padmasolala, Giri; Lin, Hsiang H.; Oswald, Fred B.

2000-01-01

A computer simulation was conducted to investigate the effectiveness of profile modification for reducing dynamic loads in gears with different tooth spacing errors. The simulation examined varying amplitudes of spacing error and differences in the span of teeth over which the error occurs. The modification considered included both linear and parabolic tip relief. The analysis considered spacing error that varies around most of the gear circumference (similar to a typical sinusoidal error pattern) as well as a shorter span of spacing errors that occurs on only a few teeth. The dynamic analysis was performed using a revised version of a NASA gear dynamics code, modified to add tooth spacing errors to the analysis. Results obtained from the investigation show that linear tip relief is more effective in reducing dynamic loads on gears with small spacing errors but parabolic tip relief becomes more effective as the amplitude of spacing error increases. In addition, the parabolic modification is more effective for the more severe error case where the error is spread over a longer span of teeth. The findings of this study can be used to design robust tooth profile modification for improving dynamic performance of gear sets with different tooth spacing errors.

Precipitation forecast verification over Brazilian watersheds on present and future climate

NASA Astrophysics Data System (ADS)

Xavier, L.; Bruyere, C. L.; Rotunno, O.

2016-12-01

Evaluating the quality of precipitation forecast is an essential step for hydrological studies, among other applications, which is particularly relevant when taking into account climate change and the consequent likely modification of precipitation patterns. In this study we analyzed daily precipitation forecasts given by the global model CESM and the regional model WRF on present and future climate. For present runs, CESM data have been considered from 1980 to 2005, and WRF data from 1990 to 2000. CESM future runs were available for 3 RCP scenarios (4.5, 6.0 and 8.5), over 2005-2100 period; for WRF, future runs spanned 4 different 11-year periods (2020-2030, 2030-2040, 2050-2060 and 2080-2090). WRF simulations had been driven by bias-corrected forcings, and had been done on present climate for a 24 members ensemble created by varying the adopted parameterization schemes. On WRF future climate simulations, data from 3 members out of the original ensemble were available. Precipitation data have been spatially averaged over some large Brazilian watersheds (Amazon and subbasins, Tocantins, Sao Francisco, 4 of Parana`s subbasins) and have been evaluated for present climate against a gauge gridded dataset and ERA Interim data both spanning the 1980-2013 period. The evaluation was focused on the analysis of precipitation forecasts probabilities distribution. Taking into account daily and monthly mean precipitation aggregated on 3-month periods (DJF,MAM,JJA,SON), we adopted some skill measures, amongst them, the Perkins Skill Score (PSS). From the results we verified that on present climate WRF ensemble mean led to clearly better results when compared with CESM data for Amazon, Tocantins and Sao Francisco, but model was not as skillful to the other basins, which could be also been observed for future climate. PSS results from future runs showed that few changes would be observed over the different periods for the considered basins.

Mechanistic investigation into sunlight-facilitated photodegradation of pyrene in seawater with oil dispersants.

PubMed

Fu, Jie; Gong, Yanyan; Cai, Zhengqing; O'Reilly, S E; Zhao, Dongye

2017-01-30

This study investigated the effects of 3 model oil dispersants (Corexit EC9500A, Corexit EC9527A and SPC 1000) on photodegradation of pyrene under simulated sunlight. Both Corexit dispersants enhanced photodegradation of pyrene, while SPC1000 slightly inhibited the reaction. Span 80 and Tween 85 were the key ingredients causing the effects, though the underlying mechanisms differed. Span 80 enriches pyrene in the upper layer of water column, whereas Tween 85 induces a photosensitization process. Two reactive oxygen species, 1 O 2 and O 2 - , were found responsible for pyrene photodegradation, though the presence of EC9500A suppressed the 1 O 2 pathway. In terms of photodegradation products, EC9500A enhanced generation of polyaromatic intermediates, i.e., phenaleno[1,9-cd][1,2]dioxine, 1-hydroxypyrene, and 1,8-pyrenequinone, but did not alter the classical photodegradation pathway. The Corexit dispersants were more prone to photochemical decomposition, with multiple by-products detected. The information aids in our understanding of the effects of dispersants on photochemical weathering of oil compositions. Copyright © 2016 Elsevier Ltd. All rights reserved.

Information rates of probabilistically shaped coded modulation for a multi-span fiber-optic communication system with 64QAM

NASA Astrophysics Data System (ADS)

Fehenberger, Tobias

2018-02-01

This paper studies probabilistic shaping in a multi-span wavelength-division multiplexing optical fiber system with 64-ary quadrature amplitude modulation (QAM) input. In split-step fiber simulations and via an enhanced Gaussian noise model, three figures of merit are investigated, which are signal-to-noise ratio (SNR), achievable information rate (AIR) for capacity-achieving forward error correction (FEC) with bit-metric decoding, and the information rate achieved with low-density parity-check (LDPC) FEC. For the considered system parameters and different shaped input distributions, shaping is found to decrease the SNR by 0.3 dB yet simultaneously increases the AIR by up to 0.4 bit per 4D-symbol. The information rates of LDPC-coded modulation with shaped 64QAM input are improved by up to 0.74 bit per 4D-symbol, which is larger than the shaping gain when considering AIRs. This increase is attributed to the reduced coding gap of the higher-rate code that is used for decoding the nonuniform QAM input.

A Comparative Study of Simulated and Measured Gear-Flap Flow Interaction

NASA Technical Reports Server (NTRS)

Khorrami, Mehdi R.; Mineck, Raymond E.; Yao, Chungsheng; Jenkins, Luther N.; Fares, Ehab

2015-01-01

The ability of two CFD solvers to accurately characterize the transient, complex, interacting flowfield asso-ciated with a realistic gear-flap configuration is assessed via comparison of simulated flow with experimental measurements. The simulated results, obtained with NASA's FUN3D and Exa's PowerFLOW® for a high-fidelity, 18% scale semi-span model of a Gulfstream aircraft in landing configuration (39 deg flap deflection, main landing gear on and off) are compared to two-dimensional and stereo particle image velocimetry measurements taken within the gear-flap flow interaction region during wind tunnel tests of the model. As part of the bench-marking process, direct comparisons of the mean and fluctuating velocity fields are presented in the form of planar contour plots and extracted line profiles at measurement planes in various orientations stationed in the main gear wake. The measurement planes in the vicinity of the flap side edge and downstream of the flap trailing edge are used to highlight the effects of gear presence on tip vortex development and the ability of the computational tools to accurately capture such effects. The present study indicates that both computed datasets contain enough detail to construct a relatively accurate depiction of gear-flap flow interaction. Such a finding increases confidence in using the simulated volumetric flow solutions to examine the behavior of pertinent aer-odynamic mechanisms within the gear-flap interaction zone.

Folding and insertion thermodynamics of the transmembrane WALP peptide

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

Bereau, Tristan, E-mail: bereau@mpip-mainz.mpg.de; Bennett, W. F. Drew; Pfaendtner, Jim

The anchor of most integral membrane proteins consists of one or several helices spanning the lipid bilayer. The WALP peptide, GWW(LA){sub n} (L)WWA, is a common model helix to study the fundamentals of protein insertion and folding, as well as helix-helix association in the membrane. Its structural properties have been illuminated in a large number of experimental and simulation studies. In this combined coarse-grained and atomistic simulation study, we probe the thermodynamics of a single WALP peptide, focusing on both the insertion across the water-membrane interface, as well as folding in both water and a membrane. The potential of meanmore » force characterizing the peptide’s insertion into the membrane shows qualitatively similar behavior across peptides and three force fields. However, the Martini force field exhibits a pronounced secondary minimum for an adsorbed interfacial state, which may even become the global minimum—in contrast to both atomistic simulations and the alternative PLUM force field. Even though the two coarse-grained models reproduce the free energy of insertion of individual amino acids side chains, they both underestimate its corresponding value for the full peptide (as compared with atomistic simulations), hinting at cooperative physics beyond the residue level. Folding of WALP in the two environments indicates the helix as the most stable structure, though with different relative stabilities and chain-length dependence.« less

Folding and insertion thermodynamics of the transmembrane WALP peptide

NASA Astrophysics Data System (ADS)

Bereau, Tristan; Bennett, W. F. Drew; Pfaendtner, Jim; Deserno, Markus; Karttunen, Mikko

2015-12-01

The anchor of most integral membrane proteins consists of one or several helices spanning the lipid bilayer. The WALP peptide, GWW(LA)n (L)WWA, is a common model helix to study the fundamentals of protein insertion and folding, as well as helix-helix association in the membrane. Its structural properties have been illuminated in a large number of experimental and simulation studies. In this combined coarse-grained and atomistic simulation study, we probe the thermodynamics of a single WALP peptide, focusing on both the insertion across the water-membrane interface, as well as folding in both water and a membrane. The potential of mean force characterizing the peptide's insertion into the membrane shows qualitatively similar behavior across peptides and three force fields. However, the Martini force field exhibits a pronounced secondary minimum for an adsorbed interfacial state, which may even become the global minimum—in contrast to both atomistic simulations and the alternative PLUM force field. Even though the two coarse-grained models reproduce the free energy of insertion of individual amino acids side chains, they both underestimate its corresponding value for the full peptide (as compared with atomistic simulations), hinting at cooperative physics beyond the residue level. Folding of WALP in the two environments indicates the helix as the most stable structure, though with different relative stabilities and chain-length dependence.

Simulation of streamflow and sediment transport in two surface-coal-mined basins in Fayette County, Pennsylvania

USGS Publications Warehouse

Sams, J. I.; Witt, E. C.

1995-01-01

The Hydrological Simulation Program - Fortran (HSPF) was used to simulate streamflow and sediment transport in two surface-mined basins of Fayette County, Pa. Hydrologic data from the Stony Fork Basin (0.93 square miles) was used to calibrate HSPF parameters. The calibrated parameters were applied to an HSPF model of the Poplar Run Basin (8.83 square miles) to evaluate the transfer value of model parameters. The results of this investigation provide information to the Pennsylvania Department of Environmental Resources, Bureau of Mining and Reclamation, regarding the value of the simulated hydrologic data for use in cumulative hydrologic-impact assessments of surface-mined basins. The calibration period was October 1, 1985, through September 30, 1988 (water years 1986-88). The simulated data were representative of the observed data from the Stony Fork Basin. Mean simulated streamflow was 1.64 cubic feet per second compared to measured streamflow of 1.58 cubic feet per second for the 3-year period. The difference between the observed and simulated peak stormflow ranged from 4.0 to 59.7 percent for 12 storms. The simulated sediment load for the 1987 water year was 127.14 tons (0.21 ton per acre), which compares to a measured sediment load of 147.09 tons (0.25 ton per acre). The total simulated suspended-sediment load for the 3-year period was 538.2 tons (0.30 ton per acre per year), which compares to a measured sediment load of 467.61 tons (0.26 ton per acre per year). The model was verified by comparing observed and simulated data from October 1, 1988, through September 30, 1989. The results obtained were comparable to those from the calibration period. The simulated mean daily discharge was representative of the range of data observed from the basin and of the frequency with which specific discharges were equalled or exceeded. The calibrated and verified parameters from the Stony Fork model were applied to an HSPF model of the Poplar Run Basin. The two basins are in a similar physical setting. Data from October 1, 1987, through September 30, 1989, were used to evaluate the Poplar Run model. In general, the results from the Poplar Run model were comparable to those obtained from the Stony Fork model. The difference between observed and simulated total streamflow was 1.1 percent for the 2-year period. The mean annual streamflow simulated by the Poplar Run model was 18.3 cubic feet per second. This compares to an observed streamflow of 18.15 cubic feet per second. For the 2-year period, the simulated sediment load was 2,754 tons (0.24 ton per acre per year), which compares to a measured sediment load of 3,051.2 tons (0.27 ton per acre per year) for the Poplar Run Basin. Cumulative frequency-distribution curves of the observed and simulated streamflow compared well. The comparison between observed and simulated data improved as the time span increased. Simulated annual means and totals were more representative of the observed data than hourly data used in comparing storm events. The structure and organization of the HSPF model facilitated the simulation of a wide range of hydrologic processes. The simulation results from this investigation indicate that model parameters may be transferred to ungaged basins to generate representative hydrologic data through modeling techniques.

Simulating large atmospheric phase screens using a woofer-tweeter algorithm.

PubMed

Buscher, David F

2016-10-03

We describe an algorithm for simulating atmospheric wavefront perturbations over ranges of spatial and temporal scales spanning more than 4 orders of magnitude. An open-source implementation of the algorithm written in Python can simulate the evolution of the perturbations more than an order-of-magnitude faster than real time. Testing of the implementation using metrics appropriate to adaptive optics systems and long-baseline interferometers show accuracies at the few percent level or better.

Further investigations of the aeroelastic behavior of the AFW wind-tunnel model using transonic small disturbance theory

NASA Technical Reports Server (NTRS)

Silva, Walter A.; Bennett, Robert M.

1992-01-01

The Computational Aeroelasticity Program-Transonic Small Disturbance (CAP-TSD) code, developed at LaRC, is applied to the active flexible wing wind-tunnel model for prediction of transonic aeroelastic behavior. A semi-span computational model is used for evaluation of symmetric motions, and a full-span model is used for evaluation of antisymmetric motions, and a full-span model is used for evaluation of antisymmetric motions. Static aeroelastic solutions using CAP-TSD are computed. Dynamic deformations are presented as flutter boundaries in terms of Mach number and dynamic pressure. Flutter boundaries that take into account modal refinements, vorticity and entropy corrections, antisymmetric motion, and sensitivity to the modeling of the wing tip ballast stores are also presented with experimental flutter results.

Understanding the Kinetics of Protein-Nanoparticle Corona Formation.

PubMed

Vilanova, Oriol; Mittag, Judith J; Kelly, Philip M; Milani, Silvia; Dawson, Kenneth A; Rädler, Joachim O; Franzese, Giancarlo

2016-12-27

When a pristine nanoparticle (NP) encounters a biological fluid, biomolecules spontaneously form adsorption layers around the NP, called "protein corona". The corona composition depends on the time-dependent environmental conditions and determines the NP's fate within living organisms. Understanding how the corona evolves is fundamental in nanotoxicology as well as medical applications. However, the process of corona formation is challenging due to the large number of molecules involved and to the large span of relevant time scales ranging from 100 μs, hard to probe in experiments, to hours, out of reach of all-atoms simulations. Here we combine experiments, simulations, and theory to study (i) the corona kinetics (over 10 -3 -10 3 s) and (ii) its final composition for silica NPs in a model plasma made of three blood proteins (human serum albumin, transferrin, and fibrinogen). When computer simulations are calibrated by experimental protein-NP binding affinities measured in single-protein solutions, the theoretical model correctly reproduces competitive protein replacement as proven by independent experiments. When we change the order of administration of the three proteins, we observe a memory effect in the final corona composition that we can explain within our model. Our combined experimental and computational approach is a step toward the development of systematic prediction and control of protein-NP corona composition based on a hierarchy of equilibrium protein binding constants.

Relations between Preschool Attention Span-Persistence and Age 25 Educational Outcomes

ERIC Educational Resources Information Center

McClelland, Megan M.; Acock, Alan C.; Piccinin, Andrea; Rhea, Sally Ann; Stallings, Michael C.

2013-01-01

This study examined relations between children's attention span-persistence in preschool and later school achievement and college completion. Children were drawn from the Colorado Adoption Project using adopted and non-adopted children (N = 430). Results of structural equation modeling indicated that children's age 4 attention span-persistence…

A Motivational Theory of Life-Span Development

ERIC Educational Resources Information Center

Heckhausen, Jutta; Wrosch, Carsten; Schulz, Richard

2010-01-01

This article had four goals. First, the authors identified a set of general challenges and questions that a life-span theory of development should address. Second, they presented a comprehensive account of their Motivational Theory of Life-Span Development. They integrated the model of optimization in primary and secondary control and the…

Life Span as the Measure of Performance and Learning in a Business Gaming Simulation

ERIC Educational Resources Information Center

Thavikulwat, Precha

2012-01-01

This study applies the learning curve method of measuring learning to participants of a computer-assisted business gaming simulation that includes a multiple-life-cycle feature. The study involved 249 participants. It verified the workability of the feature and estimated the participants' rate of learning at 17.4% for every doubling of experience.…

Distributed-observer-based cooperative control for synchronization of linear discrete-time multi-agent systems.

PubMed

Liang, Hongjing; Zhang, Huaguang; Wang, Zhanshan

2015-11-01

This paper considers output synchronization of discrete-time multi-agent systems with directed communication topologies. The directed communication graph contains a spanning tree and the exosystem as its root. Distributed observer-based consensus protocols are proposed, based on the relative outputs of neighboring agents. A multi-step algorithm is presented to construct the observer-based protocols. In light of the discrete-time algebraic Riccati equation and internal model principle, synchronization problem is completed. At last, numerical simulation is provided to verify the effectiveness of the theoretical results. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Atomic detail visualization of photosynthetic membranes with GPU-accelerated ray tracing

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

Stone, John E.; Sener, Melih; Vandivort, Kirby L.

The cellular process responsible for providing energy for most life on Earth, namely, photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. In this paper, we present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. Finally, we describemore » the techniques that were used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers.« less

Mapping and uncertainty analysis of energy and pitch angle phase space in the DIII-D fast ion loss detector

DOE PAGES

Pace, D. C.; Pipes, R.; Fisher, R. K.; ...

2014-08-05

New phase space mapping and uncertainty analysis of energetic ion loss data in the DIII-D tokamak provides experimental results that serve as valuable constraints in first-principles simulations of energetic ion transport. Beam ion losses are measured by the fast ion loss detector (FILD) diagnostic system consisting of two magnetic spectrometers placed independently along the outer wall. Monte Carlo simulations of mono-energetic and single-pitch ions reaching the FILDs are used to determine the expected uncertainty in the measurements. Modeling shows that the variation in gyrophase of 80 keV beam ions at the FILD aperture can produce an apparent measured energy signaturemore » spanning across 50-140 keV. As a result, these calculations compare favorably with experiments in which neutral beam prompt loss provides a well known energy and pitch distribution.« less

Atomic detail visualization of photosynthetic membranes with GPU-accelerated ray tracing

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

Stone, John E.; Sener, Melih; Vandivort, Kirby L.

The cellular process responsible for providing energy for most life on Earth, namely, photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. We present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. We describe the techniques that weremore » used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers.« less

Atomic detail visualization of photosynthetic membranes with GPU-accelerated ray tracing

DOE PAGES

Stone, John E.; Sener, Melih; Vandivort, Kirby L.; ...

2015-12-12

The cellular process responsible for providing energy for most life on Earth, namely, photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. In this paper, we present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. Finally, we describemore » the techniques that were used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers.« less

Impact of tropical cyclones on the intensity and phase propagation of fall Wyrtki jets

NASA Astrophysics Data System (ADS)

Sreenivas, P.; Chowdary, J. S.; Gnanaseelan, C.

2012-11-01

Observations and model simulations are used to study the impact of tropical cyclones (TC) on the fall Wyrtki jets (WJ). These strong narrow equatorial currents peak during November and play a vital role in the energy and mass transport in the tropical Indian Ocean (TIO). Maximum number of TCs is observed over TIO during November with longer than normal life span (8-15 days). These TCs enhance equatorial westerly winds (surface) and amplify monthly mean WJs (both at surface and subsurface) by 0.4 ms-1 (anomalies exceed 0.7 ms-1 during TC), which is about half of the climatological amplitude. Intensified WJs increase the heat content of eastern TIO and modulate air-sea interaction. It is also shown that movement of TCs is mainly responsible for the westward phase propagation of WJs, a previously unexplored mechanism. These features are evident in ECCO2 simulations as well.

Docking 90Sr radionuclide in cement: An atomistic modeling study

NASA Astrophysics Data System (ADS)

Youssef, Mostafa; Pellenq, Roland J.-M.; Yildiz, Bilge

Cementitious materials are considered to be a waste form for the ultimate disposal of radioactive materials in geological repositories. We investigated by means of atomistic simulations the encapsulation of strontium-90, an important radionuclide, in calcium-silicate-hydrate (C-S-H) and its crystalline analog, the 9 Å-tobermorite. C-S-H is the major binding phase of cement. Strontium was shown to energetically favor substituting calcium in the interlayer sites in C-S-H and 9 Å-tobermorite with the trend more pronounced in the latter. The integrity of the silicate chains in both cementitious waste forms were not affected by strontium substitution within the time span of molecular dynamics simulation. Finally, we observed a limited degradation of the mechanical properties in the strontium-containing cementitious waste form with the increasing strontium concentration. These results suggest the cement hydrate as a good candidate for immobilizing radioactive strontium.

Simulation of hydraulic characteristics in the white sturgeon spawning habitat of the Kootenai River near Bonners Ferry, Idaho

USGS Publications Warehouse

Berenbrock, Charles

2005-01-01

Hydraulic characterization of the Kootenai River, especially in the white sturgeon spawning habitat reach, is needed by the Kootenai River White Sturgeon Recovery Team to promote hydraulic conditions that improve spawning conditions for the white sturgeon (Acipenser transmontanus) in the Kootenai River. The decreasing population and spawning failure of white sturgeon has led to much concern. Few wild juvenile sturgeons are found in the river today. Determining the location of the transition between backwater and free-flowing water in the Kootenai River is a primary focus for biologists who believe that hydraulic changes at the transition affect the location where the sturgeon choose to spawn. The Kootenai River begins in British Columbia, Canada, and flows through Montana, Idaho, and back into British Columbia. The 65.6-mile reach of the Kootenai River in Idaho was studied. The study area encompasses the white sturgeon spawning reach that has been designated as a critical habitat. A one-dimensional hydraulic-flow model of the study reach was developed, calibrated, and used to develop relations between hydraulic characteristics and water-surface elevation, discharge, velocity, and backwater extent. The model used 164 cross sections, most of which came from a previous river survey conducted in 2002-03. The model was calibrated to water-surface elevations at specific discharges at five gaging stations. Calibrated water-surface elevations ranged from about 1,743 to about 1,759 feet, and discharges used in calibration ranged from 5,000 to 47,500 cubic feet per second. Model calibration was considered acceptable when the difference between measured and simulated water-surface elevations was ?0.15 foot or less. Measured and simulated average velocities also were compared. These comparisons indicated agreement between measured and simulated values. The location of the transition between backwater and free-flowing water was determined using the calibrated model. The model was used to simulate hydraulic characteristics for a range of water-surface elevations from 1,741 to 1,762 feet and discharges from 4,000 to 75,000 cubic feet per second. These simulated hydraulic characteristics were used to develop a three-parameter relation-discharge in the study reach, water-surface elevation at Kootenai River at Porthill gaging station (12322000), and the location of the transition between backwater and free-flowing water. Simulated hydraulic characteristics produced backwater locations ranging from river mile (RM) 105.6 (Porthill) to RM 158 (near Crossport), a span of about 52 miles. However, backwater locations from measured data ranged primarily from RM 152 to RM 157, a 5-mile span. The average backwater location from measured data was at about RM 154. Three-parameter relations also were developed for determining the amount of discharge in the Shorty Island side channel and average velocity at selected cross sections in the study reach. Simulated discharge for the side channel relative to measured data ranged from 0 to about 5,500 cubic feet per second, and simulated average velocity relative to measured data ranged from 0 to about 3.5 feet per second. Relations using other hydraulic, sediment/incipient motion, ecological, and biological characteristics also could be developed. The relations also can be used in real time by accessing data from the Web. Discharge and stage data for two gaging stations, Tribal Hatchery (12310100) and Porthill (12322500), are available from the Idaho U.S. Geological Survey web page (URL: http://waterdata.usgs.gov/id/nwis/current/?type=flow). Because the coordinate axes of the three-parameter relations use discharge from the Tribal Hatchery gaging station and water-surface elevation from the Porthill gaging station, the location of the transition between backwater and free-flowing water can be determined for current conditions using the real-time data. Similarly, discharge in the Shorty Island side channel and (or) average velocity at selected cross sections also can be determined for current conditions.

Implementation of local grid refinement (LGR) for the Lake Michigan Basin regional groundwater-flow model

USGS Publications Warehouse

Hoard, C.J.

2010-01-01

The U.S. Geological Survey is evaluating water availability and use within the Great Lakes Basin. This is a pilot effort to develop new techniques and methods to aid in the assessment of water availability. As part of the pilot program, a regional groundwater-flow model for the Lake Michigan Basin was developed using SEAWAT-2000. The regional model was used as a framework for assessing local-scale water availability through grid-refinement techniques. Two grid-refinement techniques, telescopic mesh refinement and local grid refinement, were used to illustrate the capability of the regional model to evaluate local-scale problems. An intermediate model was developed in central Michigan spanning an area of 454 square miles (mi2) using telescopic mesh refinement. Within the intermediate model, a smaller local model covering an area of 21.7 mi2 was developed and simulated using local grid refinement. Recharge was distributed in space and time using a daily output from a modified Thornthwaite-Mather soil-water-balance method. The soil-water-balance method derived recharge estimates from temperature and precipitation data output from an atmosphere-ocean coupled general-circulation model. The particular atmosphere-ocean coupled general-circulation model used, simulated climate change caused by high global greenhouse-gas emissions to the atmosphere. The surface-water network simulated in the regional model was refined and simulated using a streamflow-routing package for MODFLOW. The refined models were used to demonstrate streamflow depletion and potential climate change using five scenarios. The streamflow-depletion scenarios include (1) natural conditions (no pumping), (2) a pumping well near a stream; the well is screened in surficial glacial deposits, (3) a pumping well near a stream; the well is screened in deeper glacial deposits, and (4) a pumping well near a stream; the well is open to a deep bedrock aquifer. Results indicated that a range of 59 to 50 percent of the water pumped originated from the stream for the shallow glacial and deep bedrock pumping scenarios, respectively. The difference in streamflow reduction between the shallow and deep pumping scenarios was compensated for in the deep well by deriving more water from regional sources. The climate-change scenario only simulated natural conditions from 1991-2044, so there was no pumping stress simulated. Streamflows were calculated for the simulated period and indicated that recharge over the period generally increased from the start of the simulation until approximately 2017, and decreased from then to the end of the simulation. Streamflow was highly correlated with recharge so that the lowest streamflows occurred in the later stress periods of the model when recharge was lowest.

An Investigation of the Wing and the Wing-Fuselage Combination of a Full-Scale Model of the Republic XP-91 Airplane in the Ames 40-by 80-Foot Wing Tunnel

NASA Technical Reports Server (NTRS)

Hunton, Lynn W.; Dew, Joseph K.

1948-01-01

Wind-tunnel tests of a full-scale model of the Republic XP-91 airplane were conducted to determine the longitudinal and lateral characteristics of the wing alone and the wing-fuselage combination, the characteristics of the aileron, and the damping in roll af the wing alone. Various high-lift devices were investigated including trailing-edge split flaps and partial- and full-span leading-edge slats and Krueger-type nose flaps. Results of this investigation showed that a very significant gain in maximum lift could be achieved through use of the proper leading-edge device, The maximum lift coefficient of the model with split flaps and the original partial-span straight slats was only 1.2; whereas a value of approximately 1.8 was obtained by drooping the slat and extending it full span, Improvement in maximum lift of approximately the same amount resulted when a full-span nose flap was substituted for the original partial-span slat.

A COMPARISON OF STATIC AND DYNAMIC OPTIMIZATION MUSCLE FORCE PREDICTIONS DURING WHEELCHAIR PROPULSION

PubMed Central

Morrow, Melissa M.; Rankin, Jeffery W.; Neptune, Richard R.; Kaufman, Kenton R.

2014-01-01

The primary purpose of this study was to compare static and dynamic optimization muscle force and work predictions during the push phase of wheelchair propulsion. A secondary purpose was to compare the differences in predicted shoulder and elbow kinetics and kinematics and handrim forces. The forward dynamics simulation minimized differences between simulated and experimental data (obtained from 10 manual wheelchair users) and muscle co-contraction. For direct comparison between models, the shoulder and elbow muscle moment arms and net joint moments from the dynamic optimization were used as inputs into the static optimization routine. RMS errors between model predictions were calculated to quantify model agreement. There was a wide range of individual muscle force agreement that spanned from poor (26.4 % Fmax error in the middle deltoid) to good (6.4 % Fmax error in the anterior deltoid) in the prime movers of the shoulder. The predicted muscle forces from the static optimization were sufficient to create the appropriate motion and joint moments at the shoulder for the push phase of wheelchair propulsion, but showed deviations in the elbow moment, pronation-supination motion and hand rim forces. These results suggest the static approach does not produce results similar enough to be a replacement for forward dynamics simulations, and care should be taken in choosing the appropriate method for a specific task and set of constraints. Dynamic optimization modeling approaches may be required for motions that are greatly influenced by muscle activation dynamics or that require significant co-contraction. PMID:25282075

How Much Dust Does Enceladus eject?

NASA Astrophysics Data System (ADS)

Kempf, S.; Southworth, B.; Srama, R.; Schmidt, J.; Postberg, F.

2016-12-01

There is an ongoing argument how much dust per second the ice volcanoes on Saturn's ice moon eject. By adjusting their plume model to the dust flux measured by the Cassini dust detector during the close Enceladus flyby in 2005, Schmidt et al. (2008) obtained a total dust production rate in the plumes of about 􏱱5 kg/s. On the other hand, Ingersoll and Ewald (2005) derived a dust production rate of 51 kg/s from the total plume brightness. Knowledge of the production rate is essential for estimating the dust to gas mass ratio, which in turn is an important constraint for finding the plume source mechanism. Here we report on measurements of the plume dust density during the last close Cassini flyby at Enceladus in October 2015. The data match our numerical model for the Enceladus plume. The model is based on a large number of dynamical simulations including gravity and Lorentz force to investigate the earliest phase of the ring particle life span. The evolution of the electrostatic charge carried by the initially uncharged grains is treated self-consistently. Our numerical simulations reproduce all Enceladus data sets obtained by Cassini's Cosmic Dust Analyzer (CDA). Our model calculations together with the new density data constrain the Enceladus dust source rate to < 5 kg/s. Based on our simulation results we are able to draw conclusions about the emission of plume particles along the fractures in the south polar terrain.

Characterisation of minimal-span plane Couette turbulence with pressure gradients

NASA Astrophysics Data System (ADS)

Sekimoto, Atsushi; Atkinson, Callum; Soria, Julio

2018-04-01

The turbulence statistics and dynamics in the spanwise-minimal plane Couette flow with pressure gradients, so-called, Couette-Poiseuille (C-P) flow, are investigated using direct numerical simulation. The large-scale motion is limited in the spanwise box dimension as in the minimal-span channel turbulence of Flores & Jiménez (Phys. Fluids, vol. 22, 2010, 071704). The effect of the top wall, where normal pressure-driven Poiseuille flow is realised, is distinguished from the events on the bottom wall, where the pressure gradient results in mild or almost-zero wall-shear stress. A proper scaling of turbulence statistics in minimal-span C-P flows is presented. Also the ‘shear-less’ wall-bounded turbulence, where the Corrsin shear parameter is very weak compared to normal wall-bounded turbulence, represents local separation, which is also observed as spanwise streaks of reversed flow in full-size plane C-P turbulence. The local separation is a multi-scale event, which grows up to the order of the channel height even in the minimal-span geometry.

Several Families of Sequences with Low Correlation and Large Linear Span

NASA Astrophysics Data System (ADS)

Zeng, Fanxin; Zhang, Zhenyu

In DS-CDMA systems and DS-UWB radios, low correlation of spreading sequences can greatly help to minimize multiple access interference (MAI) and large linear span of spreading sequences can reduce their predictability. In this letter, new sequence sets with low correlation and large linear span are proposed. Based on the construction Trm1[Trnm(αbt+γiαdt)]r for generating p-ary sequences of period pn-1, where n=2m, d=upm±v, b=u±v, γi∈GF(pn), and p is an arbitrary prime number, several methods to choose the parameter d are provided. The obtained sequences with family size pn are of four-valued, five-valued, six-valued or seven-valued correlation and the maximum nontrivial correlation value is (u+v-1)pm-1. The simulation by a computer shows that the linear span of the new sequences is larger than that of the sequences with Niho-type and Welch-type decimations, and similar to that of [10].

Simulations of the future precipitation climate of the Central Andes using a coupled regional climate model

NASA Astrophysics Data System (ADS)

Nicholls, S.; Mohr, K. I.

2014-12-01

The meridional extent and complex orography of the South American continent contributes to a wide diversity of climate regimes ranging from hyper-arid deserts to tropical rainforests to sub-polar highland regions. Global climate models, although capable of resolving synoptic-scale South American climate features, are inadequate for fully-resolving the strong gradients between climate regimes and the complex orography which define the Tropical Andes given their low spatial and temporal resolution. Recent computational advances now make practical regional climate modeling with prognostic mesoscale atmosphere-ocean coupled models, such as the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system, to climate research. Previous work has shown COAWST to reasonably simulate the both the entire 2003-2004 wet season (Dec-Feb) as validated against both satellite and model analysis data. More recently, COAWST simulations have also been shown to sensibly reproduce the entire annual cycle of rainfall (Oct 2003 - Oct 2004) with historical climate model input. Using future global climate model input for COAWST, the present work involves year-long cycle spanning October to October for the years 2031, 2059, and 2087 assuming the most likely regional climate pathway (RCP): RCP 6.0. COAWST output is used to investigate how global climate change impacts the spatial distribution, precipitation rates, and diurnal cycle of precipitation patterns in the Central Andes vary in these yearly "snapshots". Initial results show little change to precipitation coverage or its diurnal cycle, however precipitation amounts did tend drier over the Brazilian Plateau and wetter over the Western Amazon and Central Andes. These results suggest potential adjustments to large-scale climate features (such as the Bolivian High).

Individual differences in working memory and reasoning-remembering relationships in solving class-inclusion problems.

PubMed

Howe, M L; Rabinowitz, F M; Powell, T L

1998-09-01

In the present experiment, we evaluated the effects of individual differences in reading span and variation in memory demands on class-inclusion performance. One hundred twenty college students whose reading spans ranged from low to medium to high (as indexed by a computerized version of the Daneman and Carpenter [1980] reading-span task) solved 48 class-inclusion problems. Half of the subjects had the solution information available when the problems were presented; the other half performed a detection task between solution information and problem presentation. The results from both standard statistical analyses and from a mathematical model indicated that differences in reading span and memory load had predictable, similar effects. Specifically, the sophistication of reasoning strategies declined when memory demands increased or when reading spans decreased. Surprisingly, these effects were primarily additive. The results were interpreted in terms of global resource models and findings from the developmental literature.

Multiple loop conformations of peptides predicted by molecular dynamics simulations are compatible with nuclear magnetic resonance.

PubMed

Carstens, Heiko; Renner, Christian; Milbradt, Alexander G; Moroder, Luis; Tavan, Paul

2005-03-29

The affinity and selectivity of protein-protein interactions can be fine-tuned by varying the size, flexibility, and amino acid composition of involved surface loops. As a model for such surface loops, we study the conformational landscape of an octapeptide, whose flexibility is chemically steered by a covalent ring closure integrating an azobenzene dye into and by a disulfide bridge additionally constraining the peptide backbone. Because the covalently integrated azobenzene dyes can be switched by light between a bent cis state and an elongated trans state, six cyclic peptide models of strongly different flexibilities are obtained. The conformational states of these peptide models are sampled by NMR and by unconstrained molecular dynamics (MD) simulations. Prototypical conformations and the free-energy landscapes in the high-dimensional space spanned by the phi/psi angles at the peptide backbone are obtained by clustering techniques from the MD trajectories. Multiple open-loop conformations are shown to be predicted by MD particularly in the very flexible cases and are shown to comply with the NMR data despite the fact that such open-loop conformations are missing in the refined NMR structures.

Eurasian Heat Waves: Mechanisms and Predictability

NASA Technical Reports Server (NTRS)

Wang, Hailan; Schubert, Siegfried; Koster, Randal; Suarez, Max

2012-01-01

This study uses the NASA MERRA reanalysis and GEOS 5 model simulations to examine the causes of Eurasian heat waves including the recent extreme events that occurred in Europe during 2003 and in Russia during 2010. The focus is on the warm season and the part of the Eurasian continent that extends north of about 40oN, or roughly to the north of the mean upper tropospheric jet stream. The results show that such extreme events are an amplification of natural patterns of atmospheric variability that develop over the Eurasian continent as a result of internal atmospheric forcing. The amplification occurs when the wave occasionally becomes locked in place for several weeks to months resulting in extreme heat and drying with the location depending on the phase of the upper atmospheric wave. An ensemble of very long GEOS-S model simulations (spanning the 20th century) forced with observed SST and greenhouse gases show that the model is capable of generating very similar heat waves, and that they have become more intense in the last thirty years as a result of the overall warming of the Asian continent. Sensitivity experiments with perturbed initial conditions indicate that these events have limited predictability.

Magnetic cluster expansion model for random and ordered magnetic face-centered cubic Fe-Ni-Cr alloys

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

Lavrentiev, M. Yu., E-mail: Mikhail.Lavrentiev@ukaea.uk; Nguyen-Manh, D.; Dudarev, S. L.

A Magnetic Cluster Expansion model for ternary face-centered cubic Fe-Ni-Cr alloys has been developed, using DFT data spanning binary and ternary alloy configurations. Using this Magnetic Cluster Expansion model Hamiltonian, we perform Monte Carlo simulations and explore magnetic structures of alloys over the entire range of compositions, considering both random and ordered alloy structures. In random alloys, the removal of magnetic collinearity constraint reduces the total magnetic moment but does not affect the predicted range of compositions where the alloys adopt low-temperature ferromagnetic configurations. During alloying of ordered fcc Fe-Ni compounds with Cr, chromium atoms tend to replace nickel rathermore » than iron atoms. Replacement of Ni by Cr in ordered alloys with high iron content increases the Curie temperature of the alloys. This can be explained by strong antiferromagnetic Fe-Cr coupling, similar to that found in bcc Fe-Cr solutions, where the Curie temperature increase, predicted by simulations as a function of Cr concentration, is confirmed by experimental observations. In random alloys, both magnetization and the Curie temperature decrease abruptly with increasing chromium content, in agreement with experiment.« less

Nanomechanics of slip avalanches in amorphous plasticity

NASA Astrophysics Data System (ADS)

Cao, Penghui; Dahmen, Karin A.; Kushima, Akihiro; Wright, Wendelin J.; Park, Harold S.; Short, Michael P.; Yip, Sidney

2018-05-01

Discrete stress relaxations (slip avalanches) in a model metallic glass under uniaxial compression are studied using a metadynamics algorithm for molecular simulation at experimental strain rates. The onset of yielding is observed at the first major stress drop, accompanied, upon analysis, by the formation of a single localized shear band region spanning the entire system. During the elastic response prior to yielding, low concentrations of shear transformation deformation events appear intermittently and spatially uncorrelated. During serrated flow following yielding, small stress drops occur interspersed between large drops. The simulation results point to a threshold value of stress dissipation as a characteristic feature separating major and minor avalanches consistent with mean-field modeling analysis and mechanical testing experiments. We further interpret this behavior to be a consequence of a nonlinear interplay of two prevailing mechanisms of amorphous plasticity, thermally activated atomic diffusion and stress-induced shear transformations, originally proposed by Spaepen and Argon, respectively. Probing the atomistic processes at widely separate strain rates gives insight to different modes of shear band formation: percolation of shear transformations versus crack-like propagation. Additionally a focus on crossover avalanche size has implications for nanomechanical modeling of spatially and temporally heterogeneous dynamics.

Time Constraints and Resource Sharing in Adults' Working Memory Spans

ERIC Educational Resources Information Center

Barrouillet, Pierre; Bernardin, Sophie; Camos, Valerie

2004-01-01

This article presents a new model that accounts for working memory spans in adults, the time-based resource-sharing model. The model assumes that both components (i.e., processing and maintenance) of the main working memory tasks require attention and that memory traces decay as soon as attention is switched away. Because memory retrievals are…

Modeling and simulation of protein-surface interactions: achievements and challenges.

PubMed

Ozboyaci, Musa; Kokh, Daria B; Corni, Stefano; Wade, Rebecca C

2016-01-01

Understanding protein-inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein-surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein-surface interactions and discuss the successes and limitations of current approaches.

A Tale of Two Forcings: Present-Day Coupled Antarctic Ice-sheet/Southern Ocean dynamics using the POPSICLES model.

NASA Astrophysics Data System (ADS)

Martin, Daniel; Asay-Davis, Xylar; Cornford, Stephen; Price, Stephen; Ng, Esmond; Collins, William

2015-04-01

We present POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period 1990 to 2010 resulting from two different choices of climate forcing: a 'normal-year' climatology and the CORE v. 2 interannual forcing data (Large and Yeager 2008). Simulations are performed at 0.1o (~5 km) ocean resolution and adaptive ice sheet resolution as fine as 500 m. We compare time-averaged melt rates below a number of major ice shelves with those reported by Rignot et al. (2013) as well as other recent studies. We also present seasonal variability and decadal melting trends from several Antarctic regions, along with the response of the ice shelves and consequent dynamics of the grounded ice sheet. POPSICLES couples the POP2x ocean model, a modified version of the Parallel Ocean Program (Smith and Gent, 2002), and the BISICLES ice-sheet model (Cornford et al., 2012). POP2x includes sub-ice-shelf circulation using partial top cells (Losch, 2008) and boundary layer physics following Holland and Jenkins (1999), Jenkins (2001), and Jenkins et al. (2010). Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP; Losch, 2008) and other continental-scale simulations and melt-rate observations (Kimura et al., 2013; Rignot et al., 2013). BISICLES makes use of adaptive mesh refinement and a 1st-order accurate momentum balance similar to the L1L2 model of Schoof and Hindmarsh (2009) to accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests (MISMIP-3d; Pattyn et al., 2013) and realistic configurations (Favier et al. 2014).

Study on Wind-induced Vibration and Fatigue Life of Cable-stayed Flexible Antenna

NASA Astrophysics Data System (ADS)

He, Kongde; He, Xuehui; Fang, Zifan; Zheng, Xiaowei; Yu, Hongchang

2018-03-01

The cable-stayed flexible antenna is a large-span space structure composed of flexible multibody, with low frequency of vibration, vortex-induced resonance can occur under the action of Stochastic wind, and a larger amplitude is generated when resonance occurs. To solve this problem, based on the theory of vortex-induced vibration, this paper analyzes the vortex-induced vibration of a cable-stayed flexible antenna under the action of Wind. Based on the sinusoidal force model and Autoregressive Model (AR) method, the vortex-induced force is simulated, then the fatigue analysis of the structure is based on the linear fatigue cumulative damage principle and the rain-flow method. The minimum fatigue life of the structure is calculated to verify the vibration fatigue performance of the structure.

Design-Based Comparison of Spine Surgery Simulators: Optimizing Educational Features of Surgical Simulators.

PubMed

Ryu, Won Hyung A; Mostafa, Ahmed E; Dharampal, Navjit; Sharlin, Ehud; Kopp, Gail; Jacobs, W Bradley; Hurlbert, R John; Chan, Sonny; Sutherland, Garnette R

2017-10-01

Simulation-based education has made its entry into surgical residency training, particularly as an adjunct to hands-on clinical experience. However, one of the ongoing challenges to wide adoption is the capacity of simulators to incorporate educational features required for effective learning. The aim of this study was to identify strengths and limitations of spine simulators to characterize design elements that are essential in enhancing resident education. We performed a mixed qualitative and quantitative cohort study with a focused survey and interviews of stakeholders in spine surgery pertaining to their experiences on 3 spine simulators. Ten participants were recruited spanning all levels of training and expertise until qualitative analysis reached saturation of themes. Participants were asked to perform lumbar pedicle screw insertion on 3 simulators. Afterward, a 10-item survey was administrated and a focused interview was conducted to explore topics pertaining to the design features of the simulators. Overall impressions of the simulators were positive with regards to their educational benefit, but our qualitative analysis revealed differing strengths and limitations. Main design strengths of the computer-based simulators were incorporation of procedural guidance and provision of performance feedback. The synthetic model excelled in achieving more realistic haptic feedback and incorporating use of actual surgical tools. Stakeholders from trainees to experts acknowledge the growing role of simulation-based education in spine surgery. However, different simulation modalities have varying design elements that augment learning in distinct ways. Characterization of these design characteristics will allow for standardization of simulation curricula in spinal surgery, optimizing educational benefit. Copyright © 2017 Elsevier Inc. All rights reserved.

Combining super-ensembles and statistical emulation to improve a regional climate and vegetation model

NASA Astrophysics Data System (ADS)

Hawkins, L. R.; Rupp, D. E.; Li, S.; Sarah, S.; McNeall, D. J.; Mote, P.; Betts, R. A.; Wallom, D.

2017-12-01

Changing regional patterns of surface temperature, precipitation, and humidity may cause ecosystem-scale changes in vegetation, altering the distribution of trees, shrubs, and grasses. A changing vegetation distribution, in turn, alters the albedo, latent heat flux, and carbon exchanged with the atmosphere with resulting feedbacks onto the regional climate. However, a wide range of earth-system processes that affect the carbon, energy, and hydrologic cycles occur at sub grid scales in climate models and must be parameterized. The appropriate parameter values in such parameterizations are often poorly constrained, leading to uncertainty in predictions of how the ecosystem will respond to changes in forcing. To better understand the sensitivity of regional climate to parameter selection and to improve regional climate and vegetation simulations, we used a large perturbed physics ensemble and a suite of statistical emulators. We dynamically downscaled a super-ensemble (multiple parameter sets and multiple initial conditions) of global climate simulations using a 25-km resolution regional climate model HadRM3p with the land-surface scheme MOSES2 and dynamic vegetation module TRIFFID. We simultaneously perturbed land surface parameters relating to the exchange of carbon, water, and energy between the land surface and atmosphere in a large super-ensemble of regional climate simulations over the western US. Statistical emulation was used as a computationally cost-effective tool to explore uncertainties in interactions. Regions of parameter space that did not satisfy observational constraints were eliminated and an ensemble of parameter sets that reduce regional biases and span a range of plausible interactions among earth system processes were selected. This study demonstrated that by combining super-ensemble simulations with statistical emulation, simulations of regional climate could be improved while simultaneously accounting for a range of plausible land-atmosphere feedback strengths.

Modeling possible spreadings of a buoyant surface plume with lagrangian and eulerian approaches at different resolutions using flow syntheses from 1992-2007 - a Gulf of Mexico study

NASA Astrophysics Data System (ADS)

Tulloch, R.; Hill, C. N.; Jahn, O.

2010-12-01

We present results from an ensemble of BP oil spill simulations. The oil spill slick is modeled as a buoyant surface plume that is transported by ocean currents modulated, in some experiments, by surface winds. Ocean currents are taken from ECCO2 project (see http://ecco2.org ) observationally constrained state estimates spanning 1992-2007. In this work we (i) explore the role of increased resolution of ocean eddies, (ii) compare inferences from particle based, lagrangian, approaches with eulerian, field based, approaches and (ii) examine the impact of differential response of oil particles and water to normal and extreme, hurricane derived, wind stress. We focus on three main questions. Is the simulated response to an oil spill markedly different for different years, depending on ocean circulation and wind forcing? Does the simulated response depend heavily on resolution and are lagrangian and eulerian estimates comparable? We start from two regional configurations of the MIT General Circulation Model (MITgcm - see http://mitgcm.org ) at 16km and 4km resolutions respectively, both covering the Gulf of Mexico and western North Atlantic regions. The simulations are driven at open boundaries with momentum and hydrographic fields from ECCO2 observationally constrained global circulation estimates. The time dependent surface flow fields from these simulations are used to transport a dye that can optionally decay over time (approximating biological breakdown) and to transport lagrangian particles. Using these experiments we examine the robustness of conclusions regarding the fate of a buoyant slick, injected at a single point. In conclusion we discuss how future drilling operations could use similar approaches to better anticipate outcomes of accidents both in this region and elsewhere.

Large-scale comparison of protein essential dynamics from molecular dynamics simulations and coarse-grained normal mode analyses.

PubMed

Ahmed, Aqeel; Villinger, Saskia; Gohlke, Holger

2010-12-01

A large-scale comparison of essential dynamics (ED) modes from molecular dynamic simulations and normal modes from coarse-grained normal mode methods (CGNM) was performed on a dataset of 335 proteins. As CGNM methods, the elastic network model (ENM) and the rigid cluster normal mode analysis (RCNMA) were used. Low-frequency normal modes from ENM correlate very well with ED modes in terms of directions of motions and relative amplitudes of motions. Notably, a similar performance was found if normal modes from RCNMA were used, despite a higher level of coarse graining. On average, the space spanned by the first quarter of ENM modes describes 84% of the space spanned by the five ED modes. Furthermore, no prominent differences for ED and CGNM modes among different protein structure classes (CATH classification) were found. This demonstrates the general potential of CGNM approaches for describing intrinsic motions of proteins with little computational cost. For selected cases, CGNM modes were found to be more robust among proteins that have the same topology or are of the same homologous superfamily than ED modes. In view of recent evidence regarding evolutionary conservation of vibrational dynamics, this suggests that ED modes, in some cases, might not be representative of the underlying dynamics that are characteristic of a whole family, probably due to insufficient sampling of some of the family members by MD. Copyright © 2010 Wiley-Liss, Inc.

Development of Semi-Span Model Test Techniques

NASA Technical Reports Server (NTRS)

Milholen, William E., II; Chokani, Ndaona; McGhee, Robert J.

1996-01-01

A computational investigation was performed to support the development of a semispan model test capability in the NASA Langley Research Center's National Transonic Facility. This capability is desirable for the testing of advanced subsonic transport aircraft at full-scale Reynolds numbers. A state-of-the-art three-dimensional Navier-Stokes solver was used to examine methods to improve the flow over a semi-span configuration. First, a parametric study is conducted to examine the influence of the stand-off height on the flow over the semispan model. It is found that decreasing the stand-off height, below the maximum fuselage radius, improves the aerodynamic characteristics of the semi-span model. Next, active sidewall boundary layer control techniques are examined. Juncture region blowing jets, upstream tangential blowing, and sidewall suction are found to improve the flow over the aft portion of the semispan model. Both upstream blowing and suction are found to reduce the sidewall boundary layer separation. The resulting near surface streamline patterns are improved, and found to be quite similar to the full-span results. Both techniques however adversely affect the pitching moment coefficient.

Multiphase High-Frequency Isolated DC-DC Converter for Industrial Applications

NASA Astrophysics Data System (ADS)

Maurya, Rakesh; Srivastava, S. P.; Agarwal, Pramod

2014-01-01

Industrial applications such as welding, plasma cutting, and surface hardening require a large DC current at low voltage. In such applications, the rating of power supply varies from few kilowatts to hundreds of kilowatts. The power supply employs in such applications particularly in arc welding process is expected to operate from open-circuit (no-load) to short-circuit (when the electrode sticks to the workpiece for a short span of time) quickly. In this paper, high-frequency isolated multiphase DC-DC converter is proposed which is well suited for aforementioned applications. Based on mathematical analysis, a simulation study with 5 kW, 5 V/1,000 A proposed model is carried out using Simulink block set and Sim Power System tool box and its performances are evaluated under symmetrical control methods. To verify the simulation results, scaled prototype model of rating 1.5 V/100 A is developed and tested with aforementioned control method under different operating conditions. In comparison with conventional welding power supply employed in many industries, the performance of proposed converter is improved significantly in terms of size and weight, efficiency and dynamic response.

Prolongation of SMAP to Spatio-temporally Seamless Coverage of Continental US Using a Deep Learning Neural Network

NASA Astrophysics Data System (ADS)

Fang, K.; Shen, C.; Kifer, D.; Yang, X.

2017-12-01

The Soil Moisture Active Passive (SMAP) mission has delivered high-quality and valuable sensing of surface soil moisture since 2015. However, its short time span, coarse resolution, and irregular revisit schedule have limited its use. Utilizing a state-of-the-art deep-in-time neural network, Long Short-Term Memory (LSTM), we created a system that predicts SMAP level-3 soil moisture data using climate forcing, model-simulated moisture, and static physical attributes as inputs. The system removes most of the bias with model simulations and also improves predicted moisture climatology, achieving a testing accuracy of 0.025 to 0.03 in most parts of Continental United States (CONUS). As the first application of LSTM in hydrology, we show that it is more robust than simpler methods in either temporal or spatial extrapolation tests. We also discuss roles of different predictors, the effectiveness of regularization algorithms and impacts of training strategies. With high fidelity to SMAP products, our data can aid various applications including data assimilation, weather forecasting, and soil moisture hindcasting.

Introduction and application of the multiscale coefficient of variation analysis.

PubMed

Abney, Drew H; Kello, Christopher T; Balasubramaniam, Ramesh

2017-10-01

Quantifying how patterns of behavior relate across multiple levels of measurement typically requires long time series for reliable parameter estimation. We describe a novel analysis that estimates patterns of variability across multiple scales of analysis suitable for time series of short duration. The multiscale coefficient of variation (MSCV) measures the distance between local coefficient of variation estimates within particular time windows and the overall coefficient of variation across all time samples. We first describe the MSCV analysis and provide an example analytical protocol with corresponding MATLAB implementation and code. Next, we present a simulation study testing the new analysis using time series generated by ARFIMA models that span white noise, short-term and long-term correlations. The MSCV analysis was observed to be sensitive to specific parameters of ARFIMA models varying in the type of temporal structure and time series length. We then apply the MSCV analysis to short time series of speech phrases and musical themes to show commonalities in multiscale structure. The simulation and application studies provide evidence that the MSCV analysis can discriminate between time series varying in multiscale structure and length.

Integrated Path Differential Absorption Lidar Optimizations Based on Pre-Analyzed Atmospheric Data for ASCENDS Mission Applications

NASA Technical Reports Server (NTRS)

Pliutau, Denis; Prasad, Narasimha S.

2012-01-01

In this paper a modeling method based on data reductions is investigated which includes pre analyzed MERRA atmospheric fields for quantitative estimates of uncertainties introduced in the integrated path differential absorption methods for the sensing of various molecules including CO2. This approach represents the extension of our existing lidar modeling framework previously developed and allows effective on- and offline wavelength optimizations and weighting function analysis to minimize the interference effects such as those due to temperature sensitivity and water vapor absorption. The new simulation methodology is different from the previous implementation in that it allows analysis of atmospheric effects over annual spans and the entire Earth coverage which was achieved due to the data reduction methods employed. The effectiveness of the proposed simulation approach is demonstrated with application to the mixing ratio retrievals for the future ASCENDS mission. Independent analysis of multiple accuracy limiting factors including the temperature, water vapor interferences, and selected system parameters is further used to identify favorable spectral regions as well as wavelength combinations facilitating the reduction in total errors in the retrieved XCO2 values.

Dynamics of stellar black holes in young star clusters with different metallicities - II. Black hole-black hole binaries

NASA Astrophysics Data System (ADS)

Ziosi, Brunetto Marco; Mapelli, Michela; Branchesi, Marica; Tormen, Giuseppe

2014-07-01

In this paper, we study the formation and dynamical evolution of black hole-black hole (BH-BH) binaries in young star clusters (YSCs), by means of N-body simulations. The simulations include metallicity-dependent recipes for stellar evolution and stellar winds, and have been run for three different metallicities (Z = 0.01, 0.1 and 1 Z⊙). Following recent theoretical models of wind mass-loss and core-collapse supernovae, we assume that the mass of the stellar remnants depends on the metallicity of the progenitor stars. We find that BH-BH binaries form efficiently because of dynamical exchanges: in our simulations, we find about 10 times more BH-BH binaries than double neutron star binaries. The simulated BH-BH binaries form earlier in metal-poor YSCs, which host more massive black holes (BHs) than in metal-rich YSCs. The simulated BH-BH binaries have very large chirp masses (up to 80 M⊙), because the BH mass is assumed to depend on metallicity, and because BHs can grow in mass due to the merger with stars. The simulated BH-BH binaries span a wide range of orbital periods (10-3-107 yr), and only a small fraction of them (0.3 per cent) is expected to merge within a Hubble time. We discuss the estimated merger rate from our simulations and the implications for Advanced VIRGO and LIGO.

Transport regimes spanning magnetization-coupling phase space

NASA Astrophysics Data System (ADS)

Baalrud, Scott D.; Daligault, Jérôme

2017-10-01

The manner in which transport properties vary over the entire parameter-space of coupling and magnetization strength is explored. Four regimes are identified based on the relative size of the gyroradius compared to other fundamental length scales: the collision mean free path, Debye length, distance of closest approach, and interparticle spacing. Molecular dynamics simulations of self-diffusion and temperature anisotropy relaxation spanning the parameter space are found to agree well with the predicted boundaries. Comparison with existing theories reveals regimes where they succeed, where they fail, and where no theory has yet been developed.

Shear strength and microstructure of polydisperse packings: The effect of size span and shape of particle size distribution.

PubMed

Azéma, Emilien; Linero, Sandra; Estrada, Nicolas; Lizcano, Arcesio

2017-08-01

By means of extensive contact dynamics simulations, we analyzed the effect of particle size distribution (PSD) on the strength and microstructure of sheared granular materials composed of frictional disks. The PSDs are built by means of a normalized β function, which allows the systematic investigation of the effects of both, the size span (from almost monodisperse to highly polydisperse) and the shape of the PSD (from linear to pronouncedly curved). We show that the shear strength is independent of the size span, which substantiates previous results obtained for uniform distributions by packing fraction. Notably, the shear strength is also independent of the shape of the PSD, as shown previously for systems composed of frictionless disks. In contrast, the packing fraction increases with the size span, but decreases with more pronounced PSD curvature. At the microscale, we analyzed the connectivity and anisotropies of the contacts and forces networks. We show that the invariance of the shear strength with the PSD is due to a compensation mechanism which involves both geometrical sources of anisotropy. In particular, contact orientation anisotropy decreases with the size span and increases with PSD curvature, while the branch length anisotropy behaves inversely.

The MICE grand challenge lightcone simulation - I. Dark matter clustering

NASA Astrophysics Data System (ADS)

Fosalba, P.; Crocce, M.; Gaztañaga, E.; Castander, F. J.

2015-04-01

We present a new N-body simulation from the Marenostrum Institut de Ciències de l'Espai (MICE) collaboration, the MICE Grand Challenge (MICE-GC), containing about 70 billion dark matter particles in a (3 Gpc h-1)3 comoving volume. Given its large volume and fine spatial resolution, spanning over five orders of magnitude in dynamic range, it allows an accurate modelling of the growth of structure in the universe from the linear through the highly non-linear regime of gravitational clustering. We validate the dark matter simulation outputs using 3D and 2D clustering statistics, and discuss mass-resolution effects in the non-linear regime by comparing to previous simulations and the latest numerical fits. We show that the MICE-GC run allows for a measurement of the BAO feature with per cent level accuracy and compare it to state-of-the-art theoretical models. We also use sub-arcmin resolution pixelized 2D maps of the dark matter counts in the lightcone to make tomographic analyses in real and redshift space. Our analysis shows the simulation reproduces the Kaiser effect on large scales, whereas we find a significant suppression of power on non-linear scales relative to the real space clustering. We complete our validation by presenting an analysis of the three-point correlation function in this and previous MICE simulations, finding further evidence for mass-resolution effects. This is the first of a series of three papers in which we present the MICE-GC simulation, along with a wide and deep mock galaxy catalogue built from it. This mock is made publicly available through a dedicated web portal, http://cosmohub.pic.es.

Evaluation and projected changes of precipitation statistics in convection-permitting WRF climate simulations over Central Europe

NASA Astrophysics Data System (ADS)

Knist, Sebastian; Goergen, Klaus; Simmer, Clemens

2018-02-01

We perform simulations with the WRF regional climate model at 12 and 3 km grid resolution for the current and future climates over Central Europe and evaluate their added value with a focus on the daily cycle and frequency distribution of rainfall and the relation between extreme precipitation and air temperature. First, a 9 year period of ERA-Interim driven simulations is evaluated against observations; then global climate model runs (MPI-ESM-LR RCP4.5 scenario) are downscaled and analyzed for three 12-year periods: a control, a mid-of-century and an end-of-century projection. The higher resolution simulations reproduce both the diurnal cycle and the hourly intensity distribution of precipitation more realistically compared to the 12 km simulation. Moreover, the observed increase of the temperature-extreme precipitation scaling from the Clausius-Clapeyron (C-C) scaling rate of 7% K-1 to a super-adiabatic scaling rate for temperatures above 11 °C is reproduced only by the 3 km simulation. The drop of the scaling rates at high temperatures under moisture limited conditions differs between sub-regions. For both future scenario time spans both simulations suggest a slight decrease in mean summer precipitation and an increase in hourly heavy and extreme precipitation. This increase is stronger in the 3 km runs. Temperature-extreme precipitation scaling curves in the future climate are projected to shift along the 7% K-1 trajectory to higher peak extreme precipitation values at higher temperatures. The curves keep their typical shape of C-C scaling followed by super-adiabatic scaling and a drop-off at higher temperatures due to moisture limitation.

Comparisons of Simulated and Observed Sub-Auroral Polarization Stream (SAPS) during the 17 March 2013 Storm

NASA Astrophysics Data System (ADS)

Chen, M.; Lemon, C.; Sazykin, S. Y.; Wolf, R.; Anderson, P. C.

2016-12-01

Sub-Auroral Polarization Streams (SAPS), characterized by large subauroral E x B velocities that span from dusk to the early morning sector for high magnetic activity, result from strong magnetosphere-ionosphere coupling. We investigate how electron and ion precipitation and the ionospheric conductance affect the simulated development of the SAPS electric field for the 17 March 2013 storm. Our approach is to use the magnetically and electrically self-consistent Rice Convection Model - Equilibrium (RCM-E) of the inner magnetosphere to simulate the SAPS. We use parameterized rates of whistler-generated electron pitch-angle scattering from Orlova and Shprits [JGR, 2014] that depend on equatorial radial distance, magnetic activity (Kp), and magnetic local time (MLT) outside the simulated plasmasphere. Inside the plasmasphere, parameterized scattering rates due to hiss [Orlova et al., GRL, 2014] are used. Ions are scattered at a fraction of strong pitch-angle scattering where the fraction is scaled by epsilon, the ratio of the gyroradius to the field-line radius of curvature, when epsilon is greater than 0.1. The electron and proton contributions to the auroral conductance in the RCM-E are calculated using the empirical Robinson et al. [JGR, 1987] and Galand and Richmond [JGR, 2001] equations, respectively. The "background" ionospheric conductance is based on parameters from the International Reference Ionosphere [Bilitza and Reinisch, JASR, 2008] but modified to include the effect of specified ionospheric troughs. Parameterized simulations will aid in understanding the underlying physical process. We compare simulated precipitating particle energy flux and E x B velocities with DMSP observations where SAPS are observed during the 17 March 2013 storm. Analysis of discerpancies between the simulation results and data will aid us in assessing needed improvements in the model.

Active Magnetic Regenerative Liquefier

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

Barclay, John A.; Oseen-Send, Kathryn; Ferguson, Luke

2016-01-12

This final report for the DOE Project entitled Active Magnetic Regenerative Liquefier (AMRL) funded under Grant DE-FG36-08GO18064 to Heracles Energy Corporation d.b.a. Prometheus Energy (Heracles/Prometheus) describes an active magnetic regenerative refrigerator (AMRR) prototype designed and built during the period from July 2008 through May 2011. The primary goal of this project was to make significant technical advances toward highly efficient liquefaction of hydrogen. Conventional hydrogen liquefiers at any scale have a maximum FOM of ~0.35 due primarily to the intrinsic difficulty of rapid, efficient compression of either hydrogen or helium working gases. Numerical simulation modeling of high performance AMRL designsmore » indicates certain designs have promise to increase thermodynamic efficiency from a FOM of ~0.35 toward ~0.5 to ~0.6. The technical approach was the use of solid magnetic working refrigerants cycled in and out of high magnetic fields to build an efficient active regenerative magnetic refrigeration module providing cooling power for AMRL. A single-stage reciprocating AMRR with a design temperature span from ~290 K to ~120 K was built and tested with dual magnetic regenerators moving in and out of the conductively-cooled superconducting magnet subsystem. The heat transfer fluid (helium) was coupled to the process stream (refrigeration/liquefaction load) via high performance heat exchangers. In order to maximize AMRR efficiency a helium bypass loop with adjustable flow was incorporated in the design because the thermal mass of magnetic refrigerants is higher in low magnetic field than in high magnetic field. Heracles/Prometheus designed experiments to measure AMRR performance under a variety of different operational parameters such as cycle frequency, magnetic field strength, heat transfer fluid flow rate, amount of bypass flow of the heat transfer fluid while measuring work input, temperature span, cooling capability as a function of cold temperature as a function of the amount of bypass flow of the heat transfer fluid. The operational AMRR prototype can be used to answer key questions such as the best recipe for multiple layers of different magnetic refrigerants in one or more integrated regenerators with varying amounts of bypass flow of the heat transfer fluid. Layered regenerators are necessary to span the AMRR range from 290 K to 120K. Our AMRR performance simulation model predicts that ~10-15 % of bypass flow should significantly improve the thermodynamic performance. Initial results obtained with regenerators made of gadolinium spheres were very encouraging; a temperature span of ~ 50 K (between 295K and 245 K) across both regenerators was achieved with zero bypass flow of the heat transfer fluid and with the magnetic field strength of ~4 T.« less

Exploring the dimensionality of digit span.

PubMed

Bowden, Stephen C; Petrauskas, Vilija M; Bardenhagen, Fiona J; Meade, Catherine E; Simpson, Leonie C

2013-04-01

The Digit Span subtest from the Wechsler Scales is used to measure Freedom from Distractibility or Working Memory. Some published research suggests that Digit Span forward should be interpreted differently from Digit Span backward. The present study explored the dimensionality of the Wechsler Memory Scale-III Digit Span (forward and backward) items in a sample of heterogeneous neuroscience patients (n = 267) using confirmatory factor analysis (CFA) for dichotomous items. Results suggested that four correlated factors underlie Digit Span, reflecting easy and hard items in both forward and backward presentation orders. The model for Digit Span was then cross-validated in a seizure disorders sample (n = 223) by replication of the CFA and by examination of measurement invariance. Measurement invariance tests of the precise numerical generalization of trait estimation across groups. Results supported measurement invariance and it was concluded that forward and backward digit span scores should be interpreted as measures of the same cognitive ability.

A Holistic Model for Wellness and Prevention over the Life Span.

ERIC Educational Resources Information Center

Witmer, J. Melvin; Sweeney, Thomas J.

1992-01-01

Presents integrated paradigm for wellness and prevention over the life span for purpose of theory building, research, clinical application, education, advocacy, and consciousness raising. Model described includes 11 characteristics desirable for optimal health and functioning. Notes characteristics are expressed through five life tasks of…

Approximate models for the ion-kinetic regime in inertial-confinement-fusion capsule implosions

DOE PAGES

Hoffman, Nelson M.; Zimmerman, George B.; Molvig, Kim; ...

2015-05-19

“Reduced” (i.e., simplified or approximate) ion-kinetic (RIK) models in radiation-hydrodynamic simulations permit a useful description of inertial-confinement-fusion (ICF) implosions where kinetic deviations from hydrodynamic behavior are important. For implosions in or near the kinetic regime (i.e., when ion mean free paths are comparable to the capsule size), simulations using a RIK model give a detailed picture of the time- and space-dependent structure of imploding capsules, allow an assessment of the relative importance of various kinetic processes during the implosion, enable explanations of past and current observations, and permit predictions of the results of future experiments. The RIK simulation method describedmore » here uses moment-based reduced kinetic models for transport of mass, momentum, and energy by long-mean-free-path ions, a model for the decrease of fusion reactivity owing to the associated modification of the ion distribution function, and a model of hydrodynamic turbulent mixing. Transport models are based on local gradient-diffusion approximations for the transport of moments of the ion distribution functions, with coefficients to impose flux limiting or account for transport modification. After calibration against a reference set of ICF implosions spanning the hydrodynamic-to-kinetic transition, the method has useful, quantifiable predictive ability over a broad range of capsule parameter space. Calibrated RIK simulations show that an important contributor to ion species separation in ICF capsule implosions is the preferential flux of longer-mean-free-path species out of the fuel and into the shell, leaving the fuel relatively enriched in species with shorter mean free paths. Also, the transport of ion thermal energy is enhanced in the kinetic regime, causing the fuel region to have a more uniform, lower ion temperature, extending over a larger volume, than implied by clean simulations. Furthermore, we expect that the success of our simple approach will motivate continued theoretical research into the development of first-principles-based, comprehensive, self-consistent, yet useable models of kinetic multispecies ion behavior in ICF plasmas.« less

The impact of anti-congestion policies on fuel consumption, carbon dioxide emissions and urban sprawl: Application of RELU-TRAN2, a CGE model

NASA Astrophysics Data System (ADS)

Hiramatsu, Tomoru

RELU-TRAN (Regional Economy and Land Use and Transportation) is a numerically solvable general equilibrium model (Anas and Liu, 2007), which treats in a unified manner the regional economy, urban land use and urban personal transportation sectors. In this dissertation, the model is extended by adding the consumer-workers' choice of private vehicle type according to the vehicle's fuel economy, by treating congestion on local roads as well as on major roads and by introducing car fuel consumption as a function of congested vehicle speed. By making the extensions, the model becomes more suitable to analyze the fuel consumption and CO2 emission consequences of urban development. The model is calibrated and simulated for the Chicago metropolitan area. By adjusting the model to the longer time span gradually, the shortand long-run price elasticities of fuel consumption are examined. As the time span becomes longer, fuel consumption becomes more elastic with respect to gasoline price, but when technological improvements in car fuel economy over comparable time spans are introduced exogenously, then the elasticity of fuel with respect to gasoline price becomes similar to that estimated in the econometric literature. Comparative statics exercises show that, if travel by auto becomes relatively more attractive in terms of travel time or travel cost than travel by public transit, then the Chicago MSA becomes more sprawled in total developed land area, whereas if public transit travel becomes relatively more attractive, then the Chicago MSA becomes more centralized. To mitigate fuel consumption and CO2 emissions, relative effectiveness of quasi-Pigouvian congestion tolls, a fuel tax on gasoline, a cordon toll around the downtown and a downtown parking fee are tested. All of these policies successfully reduce the aggregate fuel consumption and CO2. The urban growth boundary (UGB) is an alternative policy tested by the model. The UGB directly makes the Chicago MSA more centralized by prohibiting the development into urban use of a part of the vacant land in the suburban areas. The UGB also reduces aggregate fuel and CO2 emissions, but the impact is much smaller than the quasi-Pigouvian toll. Although Chicago MSA is centralized by both the UGB and the quasi-Pigouvian toll, the auto travel is directly discouraged by quasi-Pigouvian toll and but not by the UGB.

Computer Modeling of the Earliest Cellular Structures and Functions

NASA Technical Reports Server (NTRS)

Pohorille, Andrew; Chipot, Christophe; Schweighofer, Karl

2000-01-01

In the absence of extinct or extant record of protocells (the earliest ancestors of contemporary cells). the most direct way to test our understanding of the origin of cellular life is to construct laboratory models of protocells. Such efforts are currently underway in the NASA Astrobiology Program. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs for molecules that perform proto-cellular functions. Many of these functions, such as import of nutrients, capture and storage of energy. and response to changes in the environment are carried out by proteins bound to membrane< We will discuss a series of large-scale, molecular-level computer simulations which demonstrate (a) how small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes, (b) how these peptides aggregate to form membrane-spanning structures (eg. channels), and (c) by what mechanisms such aggregates perform essential proto-cellular functions, such as proton transport of protons across cell walls, a key step in cellular bioenergetics. The simulations were performed using the molecular dynamics method, in which Newton's equations of motion for each item in the system are solved iteratively. The problems of interest required simulations on multi-nanosecond time scales, which corresponded to 10(exp 6)-10(exp 8) time steps.

Towards understanding the variability in biospheric CO2 fluxes: using FTIR spectrometry and a chemical transport model to investigate the sources and sinks of carbonyl sulfide and its link to CO2

NASA Astrophysics Data System (ADS)

Wang, Yuting; Deutscher, Nicholas M.; Palm, Mathias; Warneke, Thorsten; Notholt, Justus; Baker, Ian; Berry, Joe; Suntharalingam, Parvadha; Jones, Nicholas; Mahieu, Emmanuel; Lejeune, Bernard; Hannigan, James; Conway, Stephanie; Mendonca, Joseph; Strong, Kimberly; Campbell, J. Elliott; Wolf, Adam; Kremser, Stefanie

2016-02-01

Understanding carbon dioxide (CO2) biospheric processes is of great importance because the terrestrial exchange drives the seasonal and interannual variability of CO2 in the atmosphere. Atmospheric inversions based on CO2 concentration measurements alone can only determine net biosphere fluxes, but not differentiate between photosynthesis (uptake) and respiration (production). Carbonyl sulfide (OCS) could provide an important additional constraint: it is also taken up by plants during photosynthesis but not emitted during respiration, and therefore is a potential means to differentiate between these processes. Solar absorption Fourier Transform InfraRed (FTIR) spectrometry allows for the retrievals of the atmospheric concentrations of both CO2 and OCS from measured solar absorption spectra. Here, we investigate co-located and quasi-simultaneous FTIR measurements of OCS and CO2 performed at five selected sites located in the Northern Hemisphere. These measurements are compared to simulations of OCS and CO2 using a chemical transport model (GEOS-Chem). The coupled biospheric fluxes of OCS and CO2 from the simple biosphere model (SiB) are used in the study. The CO2 simulation with SiB fluxes agrees with the measurements well, while the OCS simulation reproduced a weaker drawdown than FTIR measurements at selected sites, and a smaller latitudinal gradient in the Northern Hemisphere during growing season when comparing with HIPPO (HIAPER Pole-to-Pole Observations) data spanning both hemispheres. An offset in the timing of the seasonal cycle minimum between SiB simulation and measurements is also seen. Using OCS as a photosynthesis proxy can help to understand how the biospheric processes are reproduced in models and to further understand the carbon cycle in the real world.

Holocene constraints on simulated tropical Pacific climate

NASA Astrophysics Data System (ADS)

Emile-Geay, J.; Cobb, K. M.; Carre, M.; Braconnot, P.; Leloup, J.; Zhou, Y.; Harrison, S. P.; Correge, T.; Mcgregor, H. V.; Collins, M.; Driscoll, R.; Elliot, M.; Schneider, B.; Tudhope, A. W.

2015-12-01

The El Niño-Southern Oscillation (ENSO) influences climate and weather worldwide, so uncertainties in its response to external forcings contribute to the spread in global climate projections. Theoretical and modeling studies have argued that such forcings may affect ENSO either via the seasonal cycle, the mean state, or extratropical influences, but these mechanisms are poorly constrained by the short instrumental record. Here we synthesize a pan-Pacific network of high-resolution marine biocarbonates spanning discrete snapshots of the Holocene (past 10, 000 years of Earth's history), which we use to constrain a set of global climate model (GCM) simulations via a forward model and a consistent treatment of uncertainty. Observations suggest important reductions in ENSO variability throughout the interval, most consistently during 3-5 kyBP, when approximately 2/3 reductions are inferred. The magnitude and timing of these ENSO variance reductions bear little resemblance to those sim- ulated by GCMs, or to equatorial insolation. The central Pacific witnessed a mid-Holocene increase in seasonality, at odds with the reductions simulated by GCMs. Finally, while GCM aggregate behavior shows a clear inverse relationship between seasonal amplitude and ENSO-band variance in sea-surface temperature, in agreement with many previous studies, such a relationship is not borne out by these observations. Our synthesis suggests that tropical Pacific climate is highly variable, but exhibited millennia-long periods of reduced ENSO variability whose origins, whether forced or unforced, contradict existing explanations. It also points to deficiencies in the ability of current GCMs to simulate forced changes in the tropical Pacific seasonal cycle and its interaction with ENSO, highlighting a key area of growth for future modeling efforts.

Dark matter self-interactions and small scale structure

NASA Astrophysics Data System (ADS)

Tulin, Sean; Yu, Hai-Bo

2018-02-01

We review theories of dark matter (DM) beyond the collisionless paradigm, known as self-interacting dark matter (SIDM), and their observable implications for astrophysical structure in the Universe. Self-interactions are motivated, in part, due to the potential to explain long-standing (and more recent) small scale structure observations that are in tension with collisionless cold DM (CDM) predictions. Simple particle physics models for SIDM can provide a universal explanation for these observations across a wide range of mass scales spanning dwarf galaxies, low and high surface brightness spiral galaxies, and clusters of galaxies. At the same time, SIDM leaves intact the success of ΛCDM cosmology on large scales. This report covers the following topics: (1) small scale structure issues, including the core-cusp problem, the diversity problem for rotation curves, the missing satellites problem, and the too-big-to-fail problem, as well as recent progress in hydrodynamical simulations of galaxy formation; (2) N-body simulations for SIDM, including implications for density profiles, halo shapes, substructure, and the interplay between baryons and self-interactions; (3) semi-analytic Jeans-based methods that provide a complementary approach for connecting particle models with observations; (4) merging systems, such as cluster mergers (e.g., the Bullet Cluster) and minor infalls, along with recent simulation results for mergers; (5) particle physics models, including light mediator models and composite DM models; and (6) complementary probes for SIDM, including indirect and direct detection experiments, particle collider searches, and cosmological observations. We provide a summary and critical look for all current constraints on DM self-interactions and an outline for future directions.