Longitudinal train dynamics model for a rail transit simulation system

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

Wang, Jinghui; Rakha, Hesham A.

The paper develops a longitudinal train dynamics model in support of microscopic railway transportation simulation. The model can be calibrated without any mechanical data making it ideal for implementation in transportation simulators. The calibration and validation work is based on data collected from the Portland light rail train fleet. The calibration procedure is mathematically formulated as a constrained non-linear optimization problem. The validity of the model is assessed by comparing instantaneous model predictions against field observations, and also evaluated in the domains of acceleration/deceleration versus speed and acceleration/deceleration versus distance. A test is conducted to investigate the adequacy of themore » model in simulation implementation. The results demonstrate that the proposed model can adequately capture instantaneous train dynamics, and provides good performance in the simulation test. Thus, the model provides a simple theoretical foundation for microscopic simulators and will significantly support the planning, management and control of railway transportation systems.« less

Face, content, and construct validity of human placenta as a haptic training tool in neurointerventional surgery.

PubMed

Ribeiro de Oliveira, Marcelo Magaldi; Nicolato, Arthur; Santos, Marcilea; Godinho, Joao Victor; Brito, Rafael; Alvarenga, Alexandre; Martins, Ana Luiza Valle; Prosdocimi, André; Trivelato, Felipe Padovani; Sabbagh, Abdulrahman J; Reis, Augusto Barbosa; Maestro, Rolando Del

2016-05-01

OBJECT The development of neurointerventional treatments of central nervous system disorders has resulted in the need for adequate training environments for novice interventionalists. Virtual simulators offer anatomical definition but lack adequate tactile feedback. Animal models, which provide more lifelike training, require an appropriate infrastructure base. The authors describe a training model for neurointerventional procedures using the human placenta (HP), which affords haptic training with significantly fewer resource requirements, and discuss its validation. METHODS Twelve HPs were prepared for simulated endovascular procedures. Training exercises performed by interventional neuroradiologists and novice fellows were placental angiography, stent placement, aneurysm coiling, and intravascular liquid embolic agent injection. RESULTS The endovascular training exercises proposed can be easily reproduced in the HP. Face, content, and construct validity were assessed by 6 neurointerventional radiologists and 6 novice fellows in interventional radiology. CONCLUSIONS The use of HP provides an inexpensive training model for the training of neurointerventionalists. Preliminary validation results show that this simulation model has face and content validity and has demonstrated construct validity for the interventions assessed in this study.

Simulating traffic for incident management and ITS investment decisions

DOT National Transportation Integrated Search

1998-08-01

UTPS-type models were designed to adequately support planning activities typical of the 1960s and 1970s. However, these packages were not designed to model intelligent transportation systems (ITS) and support incident management planning. To ov...

"The Effect of Alternative Representations of Lake ...

EPA Pesticide Factsheets

Lakes can play a significant role in regional climate, modulating inland extremes in temperature and enhancing precipitation. Representing these effects becomes more important as regional climate modeling (RCM) efforts focus on simulating smaller scales. When using the Weather Research and Forecasting (WRF) model to downscale future global climate model (GCM) projections into RCM simulations, model users typically must rely on the GCM to represent temperatures at all water points. However, GCMs have insufficient resolution to adequately represent even large inland lakes, such as the Great Lakes. Some interpolation methods, such as setting lake surface temperatures (LSTs) equal to the nearest water point, can result in inland lake temperatures being set from sea surface temperatures (SSTs) that are hundreds of km away. In other cases, a single point is tasked with representing multiple large, heterogeneous lakes. Similar consequences can result from interpolating ice from GCMs to inland lake points, resulting in lakes as large as Lake Superior freezing completely in the space of a single timestep. The use of a computationally-efficient inland lake model can improve RCM simulations where the input data is too coarse to adequately represent inland lake temperatures and ice (Gula and Peltier 2012). This study examines three scenarios under which ice and LSTs can be set within the WRF model when applied as an RCM to produce 2-year simulations at 12 km gri

Knowledge representation requirements for model sharing between model-based reasoning and simulation in process flow domains

NASA Technical Reports Server (NTRS)

Throop, David R.

1992-01-01

The paper examines the requirements for the reuse of computational models employed in model-based reasoning (MBR) to support automated inference about mechanisms. Areas in which the theory of MBR is not yet completely adequate for using the information that simulations can yield are identified, and recent work in these areas is reviewed. It is argued that using MBR along with simulations forces the use of specific fault models. Fault models are used so that a particular fault can be instantiated into the model and run. This in turn implies that the component specification language needs to be capable of encoding any fault that might need to be sensed or diagnosed. It also means that the simulation code must anticipate all these faults at the component level.

Three Dimensional CFD Analysis of the GTX Combustor

NASA Technical Reports Server (NTRS)

Steffen, C. J., Jr.; Bond, R. B.; Edwards, J. R.

2002-01-01

The annular combustor geometry of a combined-cycle engine has been analyzed with three-dimensional computational fluid dynamics. Both subsonic combustion and supersonic combustion flowfields have been simulated. The subsonic combustion analysis was executed in conjunction with a direct-connect test rig. Two cold-flow and one hot-flow results are presented. The simulations compare favorably with the test data for the two cold flow calculations; the hot-flow data was not yet available. The hot-flow simulation indicates that the conventional ejector-ramjet cycle would not provide adequate mixing at the conditions tested. The supersonic combustion ramjet flowfield was simulated with frozen chemistry model. A five-parameter test matrix was specified, according to statistical design-of-experiments theory. Twenty-seven separate simulations were used to assemble surrogate models for combustor mixing efficiency and total pressure recovery. ScramJet injector design parameters (injector angle, location, and fuel split) as well as mission variables (total fuel massflow and freestream Mach number) were included in the analysis. A promising injector design has been identified that provides good mixing characteristics with low total pressure losses. The surrogate models can be used to develop performance maps of different injector designs. Several complex three-way variable interactions appear within the dataset that are not adequately resolved with the current statistical analysis.

Near-wall k-epsilon turbulence modeling

NASA Technical Reports Server (NTRS)

Mansour, N. N.; Kim, J.; Moin, P.

1987-01-01

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Transonic flow about a thick circular-arc airfoil

NASA Technical Reports Server (NTRS)

Mcdevitt, J. B.; Levy, L. L., Jr.; Deiwert, G. S.

1975-01-01

An experimental and theoretical study of transonic flow over a thick airfoil, prompted by a need for adequately documented experiments that could provide rigorous verification of viscous flow simulation computer codes, is reported. Special attention is given to the shock-induced separation phenomenon in the turbulent regime. Measurements presented include surface pressures, streamline and flow separation patterns, and shadowgraphs. For a limited range of free-stream Mach numbers the airfoil flow field is found to be unsteady. Dynamic pressure measurements and high-speed shadowgraph movies were taken to investigate this phenomenon. Comparisons of experimentally determined and numerically simulated steady flows using a new viscous-turbulent code are also included. The comparisons show the importance of including an accurate turbulence model. When the shock-boundary layer interaction is weak the turbulence model employed appears adequate, but when the interaction is strong, and extensive regions of separation are present, the model is inadequate and needs further development.

Application of Probabilistic Analysis to Aircraft Impact Dynamics

NASA Technical Reports Server (NTRS)

Lyle, Karen H.; Padula, Sharon L.; Stockwell, Alan E.

2003-01-01

Full-scale aircraft crash simulations performed with nonlinear, transient dynamic, finite element codes can incorporate structural complexities such as: geometrically accurate models; human occupant models; and advanced material models to include nonlinear stressstrain behaviors, laminated composites, and material failure. Validation of these crash simulations is difficult due to a lack of sufficient information to adequately determine the uncertainty in the experimental data and the appropriateness of modeling assumptions. This paper evaluates probabilistic approaches to quantify the uncertainty in the simulated responses. Several criteria are used to determine that a response surface method is the most appropriate probabilistic approach. The work is extended to compare optimization results with and without probabilistic constraints.

Population pharmacokinetics of rifampin in the treatment of Mycobacterium tuberculosis in Asian elephants.

PubMed

Egelund, E F; Isaza, R; Brock, A P; Alsultan, A; An, G; Peloquin, C A

2015-04-01

The objective of this study was to develop a population pharmacokinetic model for rifampin in elephants. Rifampin concentration data from three sources were pooled to provide a total of 233 oral concentrations from 37 Asian elephants. The population pharmacokinetic models were created using Monolix (version 4.2). Simulations were conducted using ModelRisk. We examined the influence of age, food, sex, and weight as model covariates. We further optimized the dosing of rifampin based upon simulations using the population pharmacokinetic model. Rifampin pharmacokinetics were best described by a one-compartment open model including first-order absorption with a lag time and first-order elimination. Body weight was a significant covariate for volume of distribution, and food intake was a significant covariate for lag time. The median Cmax of 6.07 μg/mL was below the target range of 8-24 μg/mL. Monte Carlo simulations predicted the highest treatable MIC of 0.25 μg/mL with the current initial dosing recommendation of 10 mg/kg, based upon a previously published target AUC0-24/MIC > 271 (fAUC > 41). Simulations from the population model indicate that the current dose of 10 mg/kg may be adequate for MICs up to 0.25 μg/mL. While the targeted AUC/MIC may be adequate for most MICs, the median Cmax for all elephants is below the human and elephant targeted ranges. © 2014 John Wiley & Sons Ltd.

A Study of Fan Stage/Casing Interaction Models

NASA Technical Reports Server (NTRS)

Lawrence, Charles; Carney, Kelly; Gallardo, Vicente

2003-01-01

The purpose of the present study is to investigate the performance of several existing and new, blade-case interactions modeling capabilities that are compatible with the large system simulations used to capture structural response during blade-out events. Three contact models are examined for simulating the interactions between a rotor bladed disk and a case: a radial and linear gap element and a new element based on a hydrodynamic formulation. The first two models are currently available in commercial finite element codes such as NASTRAN and have been showed to perform adequately for simulating rotor-case interactions. The hydrodynamic model, although not readily available in commercial codes, may prove to be better able to characterize rotor-case interactions.

Scale-Resolving simulations (SRS): How much resolution do we really need?

NASA Astrophysics Data System (ADS)

Pereira, Filipe M. S.; Girimaji, Sharath

2017-11-01

Scale-resolving simulations (SRS) are emerging as the computational approach of choice for many engineering flows with coherent structures. The SRS methods seek to resolve only the most important features of the coherent structures and model the remainder of the flow field with canonical closures. With reference to a typical Large-Eddy Simulation (LES), practical SRS methods aim to resolve a considerably narrower range of scales (reduced physical resolution) to achieve an adequate degree of accuracy at reasonable computational effort. While the objective of SRS is well-founded, the criteria for establishing the optimal degree of resolution required to achieve an acceptable level of accuracy are not clear. This study considers the canonical case of the flow around a circular cylinder to address the issue of `optimal' resolution. Two important criteria are developed. The first condition addresses the issue of adequate resolution of the flow field. The second guideline provides an assessment of whether the modeled field is canonical (stochastic) turbulence amenable to closure-based computations.

Empirical models of wind conditions on Upper Klamath Lake, Oregon

USGS Publications Warehouse

Buccola, Norman L.; Wood, Tamara M.

2010-01-01

Upper Klamath Lake is a large (230 square kilometers), shallow (mean depth 2.8 meters at full pool) lake in southern Oregon. Lake circulation patterns are driven largely by wind, and the resulting currents affect the water quality and ecology of the lake. To support hydrodynamic modeling of the lake and statistical investigations of the relation between wind and lake water-quality measurements, the U.S. Geological Survey has monitored wind conditions along the lakeshore and at floating raft sites in the middle of the lake since 2005. In order to make the existing wind archive more useful, this report summarizes the development of empirical wind models that serve two purposes: (1) to fill short (on the order of hours or days) wind data gaps at raft sites in the middle of the lake, and (2) to reconstruct, on a daily basis, over periods of months to years, historical wind conditions at U.S. Geological Survey sites prior to 2005. Empirical wind models based on Artificial Neural Network (ANN) and Multivariate-Adaptive Regressive Splines (MARS) algorithms were compared. ANNs were better suited to simulating the 10-minute wind data that are the dependent variables of the gap-filling models, but the simpler MARS algorithm may be adequate to accurately simulate the daily wind data that are the dependent variables of the historical wind models. To further test the accuracy of the gap-filling models, the resulting simulated winds were used to force the hydrodynamic model of the lake, and the resulting simulated currents were compared to measurements from an acoustic Doppler current profiler. The error statistics indicated that the simulation of currents was degraded as compared to when the model was forced with observed winds, but probably is adequate for short gaps in the data of a few days or less. Transport seems to be less affected by the use of the simulated winds in place of observed winds. The simulated tracer concentration was similar between model results when simulated winds were used to force the model, and when observed winds were used to force the model, and differences between the two results did not accumulate over time.

Assessment of the Ability of Contemporary Climate Models to Assess Adequately the Risk of Possible Regional Anomalies and Trends

NASA Astrophysics Data System (ADS)

Mokhov, I. I.

2018-04-01

The results describing the ability of contemporary global and regional climate models not only to assess the risk of general trends of changes but also to predict qualitatively new regional effects are presented. In particular, model simulations predicted spatially inhomogeneous changes in the wind and wave conditions in the Arctic basins, which have been confirmed in recent years. According to satellite and reanalysis data, a qualitative transition to the regime predicted by model simulations occurred about a decade ago.

Modelling and simulation techniques for membrane biology.

PubMed

Burrage, Kevin; Hancock, John; Leier, André; Nicolau, Dan V

2007-07-01

One of the most important aspects of Computational Cell Biology is the understanding of the complicated dynamical processes that take place on plasma membranes. These processes are often so complicated that purely temporal models cannot always adequately capture the dynamics. On the other hand, spatial models can have large computational overheads. In this article, we review some of these issues with respect to chemistry, membrane microdomains and anomalous diffusion and discuss how to select appropriate modelling and simulation paradigms based on some or all the following aspects: discrete, continuous, stochastic, delayed and complex spatial processes.

Ocean-Atmosphere Coupled Model Simulations of Precipitation in the Central Andes

NASA Technical Reports Server (NTRS)

Nicholls, Stephen D.; Mohr, Karen I.

2015-01-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. In addition, South American meteorology and climate are also made further complicated by ENSO, a powerful coupled ocean-atmosphere phenomenon. Modelling studies in this region have typically resorted to either atmospheric mesoscale or atmosphere-ocean coupled global climate models. The latter offers full physics and high spatial resolution, but it is computationally inefficient typically lack an interactive ocean, whereas the former offers high computational efficiency and ocean-atmosphere coupling, but it lacks adequate spatial and temporal resolution to adequate resolve the complex orography and explicitly simulate precipitation. Explicit simulation of precipitation is vital in the Central Andes where rainfall rates are light (0.5-5 mm hr-1), there is strong seasonality, and most precipitation is associated with weak mesoscale-organized convection. Recent increases in both computational power and model development have led to the advent of coupled ocean-atmosphere mesoscale models for both weather and climate study applications. These modelling systems, while computationally expensive, include two-way ocean-atmosphere coupling, high resolution, and explicit simulation of precipitation. In this study, we use the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST), a fully-coupled mesoscale atmosphere-ocean modeling system. Previous work has shown COAWST to reasonably simulate the entire 2003-2004 wet season (Dec-Feb) as validated against both satellite and model analysis data when ECMWF interim analysis data were used for boundary conditions on a 27-9-km grid configuration (Outer grid extent: 60.4S to 17.7N and 118.6W to 17.4W).

Application of Probability Methods to Assess Crash Modeling Uncertainty

NASA Technical Reports Server (NTRS)

Lyle, Karen H.; Stockwell, Alan E.; Hardy, Robin C.

2003-01-01

Full-scale aircraft crash simulations performed with nonlinear, transient dynamic, finite element codes can incorporate structural complexities such as: geometrically accurate models; human occupant models; and advanced material models to include nonlinear stress-strain behaviors, and material failure. Validation of these crash simulations is difficult due to a lack of sufficient information to adequately determine the uncertainty in the experimental data and the appropriateness of modeling assumptions. This paper evaluates probabilistic approaches to quantify the effects of finite element modeling assumptions on the predicted responses. The vertical drop test of a Fokker F28 fuselage section will be the focus of this paper. The results of a probabilistic analysis using finite element simulations will be compared with experimental data.

Application of Probability Methods to Assess Crash Modeling Uncertainty

NASA Technical Reports Server (NTRS)

Lyle, Karen H.; Stockwell, Alan E.; Hardy, Robin C.

2007-01-01

Full-scale aircraft crash simulations performed with nonlinear, transient dynamic, finite element codes can incorporate structural complexities such as: geometrically accurate models; human occupant models; and advanced material models to include nonlinear stress-strain behaviors, and material failure. Validation of these crash simulations is difficult due to a lack of sufficient information to adequately determine the uncertainty in the experimental data and the appropriateness of modeling assumptions. This paper evaluates probabilistic approaches to quantify the effects of finite element modeling assumptions on the predicted responses. The vertical drop test of a Fokker F28 fuselage section will be the focus of this paper. The results of a probabilistic analysis using finite element simulations will be compared with experimental data.

Modeling tree-level fuel connectivity to evaluate the effectiveness of thinning treatments for reducing crown fire potential

Treesearch

Marco A. Contreras; Russell A. Parsons; Woodam Chung

2012-01-01

Land managers have been using fire behavior and simulation models to assist in several fire management tasks. These widely-used models use average attributes to make stand-level predictions without considering spatial variability of fuels within a stand. Consequently, as the existing models have limitations in adequately modeling crown fire initiation and propagation,...

Using Transom Jack in the Human Engineering Analysis of the Materials Science Research Rack-1 and Quench Module Insert

NASA Technical Reports Server (NTRS)

Dunn, Mariea C.; Alves, Jeffrey R.; Hutchinson, Sonya L.

1999-01-01

This paper describes the human engineering analysis performed on the Materials Science Research Rack-1 and Quench Module Insert (MSRR-1/QMI) using Transom Jack (Jack) software. The Jack software was used to model a virtual environment consisting of the MSRR-1/QMI hardware configuration and human figures representing the 95th percentile male and 5th percentile female. The purpose of the simulation was to assess the human interfaces in the design for their ability to meet the requirements of the Pressurized Payloads Interface Requirements Document - International Space Program, Revision C (SSP 57000). Jack was used in the evaluation because of its ability to correctly model anthropometric body measurements and the physical behavior of astronauts working in microgravity, which is referred to as the neutral body posture. The Jack model allows evaluation of crewmember interaction with hardware through task simulation including but not limited to collision avoidance behaviors, hand/eye coordination, reach path planning, and automatic grasping to part contours. Specifically, this virtual simulation depicts the human figures performing the QMI installation and check-out, sample cartridge insertion and removal, and gas bottle drawer removal. These tasks were evaluated in terms of adequate clearance in reach envelopes, adequate accessibility in work envelopes, appropriate line of sight in visual envelopes, and accommodation of full size range for male and female stature maneuverability. The results of the human engineering analysis virtual simulation indicate that most of the associated requirements of SSP 57000 were met. However, some hardware design considerations and crew procedures modifications are recommended to improve accessibility, provide an adequate work envelope, reduce awkward body posture, and eliminate permanent protrusions.

Assessment of Alternative Conceptual Models Using Reactive Transport Modeling with Monitoring Data

NASA Astrophysics Data System (ADS)

Dai, Z.; Price, V.; Heffner, D.; Hodges, R.; Temples, T.; Nicholson, T.

2005-12-01

Monitoring data proved very useful in evaluating alternative conceptual models, simulating contaminant transport behavior, and reducing uncertainty. A graded approach using three alternative conceptual site models was formulated to simulate a field case of tetrachloroethene (PCE) transport and biodegradation. These models ranged from simple to complex in their representation of subsurface heterogeneities. The simplest model was a single-layer homogeneous aquifer that employed an analytical reactive transport code, BIOCHLOR (Aziz et al., 1999). Due to over-simplification of the aquifer structure, this simulation could not reproduce the monitoring data. The second model consisted of a multi-layer conceptual model, in combination with numerical modules, MODFLOW and RT3D within GMS, to simulate flow and reactive transport. Although the simulation results from the second model were comparatively better than those from the simple model, they still did not adequately reproduce the monitoring well concentrations because the geological structures were still inadequately defined. Finally, a more realistic conceptual model was formulated that incorporated heterogeneities and geologic structures identified from well logs and seismic survey data using the Petra and PetraSeis software. This conceptual model included both a major channel and a younger channel that were detected in the PCE source area. In this model, these channels control the local ground-water flow direction and provide a preferential chemical transport pathway. Simulation results using this conceptual site model proved compatible with the monitoring concentration data. This study demonstrates that the bias and uncertainty from inadequate conceptual models are much larger than those introduced from an inadequate choice of model parameter values (Neuman and Wierenga, 2003; Meyer et al., 2004; Ye et al., 2004). This case study integrated conceptual and numerical models, based on interpreted local hydrogeologic and geochemical data, with detailed monitoring plume data. It provided key insights for confirming alternative conceptual site models and assessing the performance of monitoring networks. A monitoring strategy based on this graded approach for assessing alternative conceptual models can provide the technical bases for identifying critical monitoring locations, adequate monitoring frequency, and performance indicator parameters for performance monitoring involving ground-water levels and PCE concentrations.

Behavior of the gypsy moth life system model and development of synoptic model formulations

Treesearch

J. J. Colbert; Xu Rumei

1991-01-01

Aims of the research: The gypsy moth life system model (GMLSM) is a complex model which incorporates numerous components (both biotic and abiotic) and ecological processes. It is a detailed simulation model which has much biological reality. However, it has not yet been tested with life system data. For such complex models, evaluation and testing cannot be adequately...

Simulating potato gas exchange as influenced by CO2 and irrigation

USDA-ARS?s Scientific Manuscript database

Recent research suggests that an energy balance approach is required for crop models to adequately respond to current and future climatic conditions associated with elevated CO2, higher temperatures, and water scarcity. More realistic models are needed in order to understand the impact of, and deve...

Techniques of Australian forest planning

Treesearch

Australian Forestry Council

1978-01-01

Computer modeling has been extensively adopted for Australian forest planning over the last ten years. It has been confined almost entirely to the plantations of fast-growing species for which adequate inventory, growth, and experimental data are available. Stand simulation models have replaced conventional yield tables and enabled a wide range of alternative...

Simulation tools for robotics research and assessment

NASA Astrophysics Data System (ADS)

Fields, MaryAnne; Brewer, Ralph; Edge, Harris L.; Pusey, Jason L.; Weller, Ed; Patel, Dilip G.; DiBerardino, Charles A.

2016-05-01

The Robotics Collaborative Technology Alliance (RCTA) program focuses on four overlapping technology areas: Perception, Intelligence, Human-Robot Interaction (HRI), and Dexterous Manipulation and Unique Mobility (DMUM). In addition, the RCTA program has a requirement to assess progress of this research in standalone as well as integrated form. Since the research is evolving and the robotic platforms with unique mobility and dexterous manipulation are in the early development stage and very expensive, an alternate approach is needed for efficient assessment. Simulation of robotic systems, platforms, sensors, and algorithms, is an attractive alternative to expensive field-based testing. Simulation can provide insight during development and debugging unavailable by many other means. This paper explores the maturity of robotic simulation systems for applications to real-world problems in robotic systems research. Open source (such as Gazebo and Moby), commercial (Simulink, Actin, LMS), government (ANVEL/VANE), and the RCTA-developed RIVET simulation environments are examined with respect to their application in the robotic research domains of Perception, Intelligence, HRI, and DMUM. Tradeoffs for applications to representative problems from each domain are presented, along with known deficiencies and disadvantages. In particular, no single robotic simulation environment adequately covers the needs of the robotic researcher in all of the domains. Simulation for DMUM poses unique constraints on the development of physics-based computational models of the robot, the environment and objects within the environment, and the interactions between them. Most current robot simulations focus on quasi-static systems, but dynamic robotic motion places an increased emphasis on the accuracy of the computational models. In order to understand the interaction of dynamic multi-body systems, such as limbed robots, with the environment, it may be necessary to build component-level computational models to provide the necessary simulation fidelity for accuracy. However, the Perception domain remains the most problematic for adequate simulation performance due to the often cartoon nature of computer rendering and the inability to model realistic electromagnetic radiation effects, such as multiple reflections, in real-time.

Multiresponse modeling of variably saturated flow and isotope tracer transport for a hillslope experiment at the Landscape Evolution Observatory

NASA Astrophysics Data System (ADS)

Scudeler, Carlotta; Pangle, Luke; Pasetto, Damiano; Niu, Guo-Yue; Volkmann, Till; Paniconi, Claudio; Putti, Mario; Troch, Peter

2016-10-01

This paper explores the challenges of model parameterization and process representation when simulating multiple hydrologic responses from a highly controlled unsaturated flow and transport experiment with a physically based model. The experiment, conducted at the Landscape Evolution Observatory (LEO), involved alternate injections of water and deuterium-enriched water into an initially very dry hillslope. The multivariate observations included point measures of water content and tracer concentration in the soil, total storage within the hillslope, and integrated fluxes of water and tracer through the seepage face. The simulations were performed with a three-dimensional finite element model that solves the Richards and advection-dispersion equations. Integrated flow, integrated transport, distributed flow, and distributed transport responses were successively analyzed, with parameterization choices at each step supported by standard model performance metrics. In the first steps of our analysis, where seepage face flow, water storage, and average concentration at the seepage face were the target responses, an adequate match between measured and simulated variables was obtained using a simple parameterization consistent with that from a prior flow-only experiment at LEO. When passing to the distributed responses, it was necessary to introduce complexity to additional soil hydraulic parameters to obtain an adequate match for the point-scale flow response. This also improved the match against point measures of tracer concentration, although model performance here was considerably poorer. This suggests that still greater complexity is needed in the model parameterization, or that there may be gaps in process representation for simulating solute transport phenomena in very dry soils.

Simulations using patient data to evaluate systematic errors that may occur in 4D treatment planning: a proof of concept study.

PubMed

St James, Sara; Seco, Joao; Mishra, Pankaj; Lewis, John H

2013-09-01

The purpose of this work is to present a framework to evaluate the accuracy of four-dimensional treatment planning in external beam radiation therapy using measured patient data and digital phantoms. To accomplish this, 4D digital phantoms of two model patients were created using measured patient lung tumor positions. These phantoms were used to simulate a four-dimensional computed tomography image set, which in turn was used to create a 4D Monte Carlo (4DMC) treatment plan. The 4DMC plan was evaluated by simulating the delivery of the treatment plan over approximately 5 min of tumor motion measured from the same patient on a different day. Unique phantoms accounting for the patient position (tumor position and thorax position) at 2 s intervals were used to represent the model patients on the day of treatment delivery and the delivered dose to the tumor was determined using Monte Carlo simulations. For Patient 1, the tumor was adequately covered with 95.2% of the tumor receiving the prescribed dose. For Patient 2, the tumor was not adequately covered and only 74.3% of the tumor received the prescribed dose. This study presents a framework to evaluate 4D treatment planning methods and demonstrates a potential limitation of 4D treatment planning methods. When systematic errors are present, including when the imaging study used for treatment planning does not represent all potential tumor locations during therapy, the treatment planning methods may not adequately predict the dose to the tumor. This is the first example of a simulation study based on patient tumor trajectories where systematic errors that occur due to an inaccurate estimate of tumor motion are evaluated.

An Improved MUSIC Model for Gibbsite Surfaces

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

Mitchell, Scott C.; Bickmore, Barry R.; Tadanier, Christopher J.

2004-06-01

Here we use gibbsite as a model system with which to test a recently published, bond-valence method for predicting intrinsic pKa values for surface functional groups on oxides. At issue is whether the method is adequate when valence parameters for the functional groups are derived from ab initio structure optimization of surfaces terminated by vacuum. If not, ab initio molecular dynamics (AIMD) simulations of solvated surfaces (which are much more computationally expensive) will have to be used. To do this, we had to evaluate extant gibbsite potentiometric titration data that where some estimate of edge and basal surface area wasmore » available. Applying BET and recently developed atomic force microscopy methods, we found that most of these data sets were flawed, in that their surface area estimates were probably wrong. Similarly, there may have been problems with many of the titration procedures. However, one data set was adequate on both counts, and we applied our method of surface pKa int prediction to fitting a MUSIC model to this data with considerable success—several features of the titration data were predicted well. However, the model fit was certainly not perfect, and we experienced some difficulties optimizing highly charged, vacuum-terminated surfaces. Therefore, we conclude that we probably need to do AIMD simulations of solvated surfaces to adequately predict intrinsic pKa values for surface functional groups.« less

Quadrupedal locomotor simulation: producing more realistic gaits using dual-objective optimization

PubMed Central

Hirasaki, Eishi

2018-01-01

In evolutionary biomechanics it is often considered that gaits should evolve to minimize the energetic cost of travelling a given distance. In gait simulation this goal often leads to convincing gait generation. However, as the musculoskeletal models used get increasingly sophisticated, it becomes apparent that such a single goal can lead to extremely unrealistic gait patterns. In this paper, we explore the effects of requiring adequate lateral stability and show how this increases both energetic cost and the realism of the generated walking gait in a high biofidelity chimpanzee musculoskeletal model. We also explore the effects of changing the footfall sequences in the simulation so it mimics both the diagonal sequence walking gaits that primates typically use and also the lateral sequence walking gaits that are much more widespread among mammals. It is apparent that adding a lateral stability criterion has an important effect on the footfall phase relationship, suggesting that lateral stability may be one of the key drivers behind the observed footfall sequences in quadrupedal gaits. The observation that single optimization goals are no longer adequate for generating gait in current models has important implications for the use of biomimetic virtual robots to predict the locomotor patterns in fossil animals. PMID:29657790

Experiences in integrating auto-translated state-chart designs for model checking

NASA Technical Reports Server (NTRS)

Pingree, P. J.; Benowitz, E. G.

2003-01-01

In the complex environment of JPL's flight missions with increasing dependency on advanced software designs, traditional software validation methods of simulation and testing are being stretched to adequately cover the needs of software development.

Hydrological simulation of Sperchios River basin in Central Greece using the MIKE SHE model and geographic information systems

NASA Astrophysics Data System (ADS)

Paparrizos, Spyridon; Maris, Fotios

2017-05-01

The MIKE SHE model is able to simulate the entire stream flow which includes direct and basic flow. Many models either do not simulate or use simplistic methods to determine the basic flow. The MIKE SHE model takes into account many hydrological data. Since this study was directed towards the simulation of surface runoff and infiltration into saturated and unsaturated zone, the MIKE SHE is an appropriate model for reliable conclusions. In the current research, the MIKE SHE model was used to simulate runoff in the area of Sperchios River basin. Meteorological data from eight rainfall stations within the Sperchios River basin were used as inputs. Vegetation as well as geological data was used to perform the calibration and validation of the physical processes of the model. Additionally, ArcGIS program was used. The results indicated that the model was able to simulate the surface runoff satisfactorily, representing all the hydrological data adequately. Some minor differentiations appeared which can be eliminated with the appropriate adjustments that can be decided by the researcher's experience.

Process Modeling of Composite Materials for Wind-Turbine Rotor Blades: Experiments and Numerical Modeling

PubMed Central

Wieland, Birgit; Ropte, Sven

2017-01-01

The production of rotor blades for wind turbines is still a predominantly manual process. Process simulation is an adequate way of improving blade quality without a significant increase in production costs. This paper introduces a module for tolerance simulation for rotor-blade production processes. The investigation focuses on the simulation of temperature distribution for one-sided, self-heated tooling and thick laminates. Experimental data from rotor-blade production and down-scaled laboratory tests are presented. Based on influencing factors that are identified, a physical model is created and implemented as a simulation. This provides an opportunity to simulate temperature and cure-degree distribution for two-dimensional cross sections. The aim of this simulation is to support production processes. Hence, it is modelled as an in situ simulation with direct input of temperature data and real-time capability. A monolithic part of the rotor blade, the main girder, is used as an example for presenting the results. PMID:28981458

Process Modeling of Composite Materials for Wind-Turbine Rotor Blades: Experiments and Numerical Modeling.

PubMed

Wieland, Birgit; Ropte, Sven

2017-10-05

The production of rotor blades for wind turbines is still a predominantly manual process. Process simulation is an adequate way of improving blade quality without a significant increase in production costs. This paper introduces a module for tolerance simulation for rotor-blade production processes. The investigation focuses on the simulation of temperature distribution for one-sided, self-heated tooling and thick laminates. Experimental data from rotor-blade production and down-scaled laboratory tests are presented. Based on influencing factors that are identified, a physical model is created and implemented as a simulation. This provides an opportunity to simulate temperature and cure-degree distribution for two-dimensional cross sections. The aim of this simulation is to support production processes. Hence, it is modelled as an in situ simulation with direct input of temperature data and real-time capability. A monolithic part of the rotor blade, the main girder, is used as an example for presenting the results.

Simulation of advective flow under steady-state and transient recharge conditions, Camp Edwards, Massachusetts Military Reservation, Cape Cod, Massachusetts

USGS Publications Warehouse

Walter, Donald A.; Masterson, John P.

2003-01-01

The U.S. Geological Survey has developed several ground-water models in support of an investigation of ground-water contamination being conducted by the Army National Guard Bureau at Camp Edwards, Massachusetts Military Reservation on western Cape Cod, Massachusetts. Regional and subregional steady-state models and regional transient models were used to (1) improve understanding of the hydrologic system, (2) simulate advective transport of contaminants, (3) delineate recharge areas to municipal wells, and (4) evaluate how model discretization and time-varying recharge affect simulation results. A water-table mound dominates ground-water-flow patterns. Near the top of the mound, which is within Camp Edwards, hydraulic gradients are nearly vertically downward and horizontal gradients are small. In downgradient areas that are further from the top of the water-table mound, the ratio of horizontal to vertical gradients is larger and horizontal flow predominates. The steady-state regional model adequately simulates advective transport in some areas of the aquifer; however, simulation of ground-water flow in areas with local hydrologic boundaries, such as ponds, requires more finely discretized subregional models. Subregional models also are needed to delineate recharge areas to municipal wells that are inadequately represented in the regional model or are near other pumped wells. Long-term changes in recharge rates affect hydraulic heads in the aquifer and shift the position of the top of the water-table mound. Hydraulic-gradient directions do not change over time in downgradient areas, whereas they do change substantially with temporal changes in recharge near the top of the water-table mound. The assumption of steady-state hydraulic conditions is valid in downgradient area, where advective transport paths change little over time. In areas closer to the top of the water-table mound, advective transport paths change as a function of time, transient and steady-state paths do not coincide, and the assumption of steady-state conditions is not valid. The simulation results indicate that several modeling tools are needed to adequately simulate ground-water flow at the site and that the utility of a model varies according to hydrologic conditions in the specific areas of interest.

Simulation of Triple Oxidation Ditch Wastewater Treatment Process

NASA Astrophysics Data System (ADS)

Yang, Yue; Zhang, Jinsong; Liu, Lixiang; Hu, Yongfeng; Xu, Ziming

2010-11-01

This paper presented the modeling mechanism and method of a sewage treatment system. A triple oxidation ditch process of a WWTP was simulated based on activated sludge model ASM2D with GPS-X software. In order to identify the adequate model structure to be implemented into the GPS-X environment, the oxidation ditch was divided into several completely stirred tank reactors depended on the distribution of aeration devices and dissolved oxygen concentration. The removal efficiency of COD, ammonia nitrogen, total nitrogen, total phosphorus and SS were simulated by GPS-X software with influent quality data of this WWTP from June to August 2009, to investigate the differences between the simulated results and the actual results. The results showed that, the simulated values could well reflect the actual condition of the triple oxidation ditch process. Mathematical modeling method was appropriate in effluent quality predicting and process optimizing.

Recent progress towards predicting aircraft ground handling performance

NASA Technical Reports Server (NTRS)

Yager, T. J.; White, E. J.

1981-01-01

Capability implemented in simulating aircraft ground handling performance is reviewed and areas for further expansion and improvement are identified. Problems associated with providing necessary simulator input data for adequate modeling of aircraft tire/runway friction behavior are discussed and efforts to improve tire/runway friction definition, and simulator fidelity are described. Aircraft braking performance data obtained on several wet runway surfaces are compared to ground vehicle friction measurements. Research to improve methods of predicting tire friction performance are discussed.

Modeling laser velocimeter signals as triply stochastic Poisson processes

NASA Technical Reports Server (NTRS)

Mayo, W. T., Jr.

1976-01-01

Previous models of laser Doppler velocimeter (LDV) systems have not adequately described dual-scatter signals in a manner useful for analysis and simulation of low-level photon-limited signals. At low photon rates, an LDV signal at the output of a photomultiplier tube is a compound nonhomogeneous filtered Poisson process, whose intensity function is another (slower) Poisson process with the nonstationary rate and frequency parameters controlled by a random flow (slowest) process. In the present paper, generalized Poisson shot noise models are developed for low-level LDV signals. Theoretical results useful in detection error analysis and simulation are presented, along with measurements of burst amplitude statistics. Computer generated simulations illustrate the difference between Gaussian and Poisson models of low-level signals.

The Co-simulation of Humanoid Robot Based on Solidworks, ADAMS and Simulink

NASA Astrophysics Data System (ADS)

Song, Dalei; Zheng, Lidan; Wang, Li; Qi, Weiwei; Li, Yanli

A simulation method of adaptive controller is proposed for the humanoid robot system based on co-simulation of Solidworks, ADAMS and Simulink. A complex mathematical modeling process is avoided by this method, and the real time dynamic simulating function of Simulink would be exerted adequately. This method could be generalized to other complicated control system. This method is adopted to build and analyse the model of humanoid robot. The trajectory tracking and adaptive controller design also proceed based on it. The effect of trajectory tracking is evaluated by fitting-curve theory of least squares method. The anti-interference capability of the robot is improved a lot through comparative analysis.

From soilscapes to landscapes: A landscape-oriented approach to simulate soil organic carbon dynamics in intensively managed landscapes

USDA-ARS?s Scientific Manuscript database

Most available biogeochemical models focus within a soil profile and cannot adequately resolve contributions of the lighter size fractions of organic rich soils for Enrichment Ratio (ER) estimates, thereby causing unintended errors in Soil Organic Carbon (SOC) storage predictions. These models set E...

A visual-environment simulator with variable contrast

NASA Astrophysics Data System (ADS)

Gusarova, N. F.; Demin, A. V.; Polshchikov, G. V.

1987-01-01

A visual-environment simulator is proposed in which the image contrast can be varied continuously up to the reversal of the image. Contrast variability can be achieved by using two independently adjustable light sources to simultaneously illuminate the carrier of visual information (e.g., a slide or a cinematographic film). It is shown that such a scheme makes it possible to adequately model a complex visual environment.

Mind the Noise When Identifying Computational Models of Cognition from Brain Activity.

PubMed

Kolossa, Antonio; Kopp, Bruno

2016-01-01

The aim of this study was to analyze how measurement error affects the validity of modeling studies in computational neuroscience. A synthetic validity test was created using simulated P300 event-related potentials as an example. The model space comprised four computational models of single-trial P300 amplitude fluctuations which differed in terms of complexity and dependency. The single-trial fluctuation of simulated P300 amplitudes was computed on the basis of one of the models, at various levels of measurement error and at various numbers of data points. Bayesian model selection was performed based on exceedance probabilities. At very low numbers of data points, the least complex model generally outperformed the data-generating model. Invalid model identification also occurred at low levels of data quality and under low numbers of data points if the winning model's predictors were closely correlated with the predictors from the data-generating model. Given sufficient data quality and numbers of data points, the data-generating model could be correctly identified, even against models which were very similar to the data-generating model. Thus, a number of variables affects the validity of computational modeling studies, and data quality and numbers of data points are among the main factors relevant to the issue. Further, the nature of the model space (i.e., model complexity, model dependency) should not be neglected. This study provided quantitative results which show the importance of ensuring the validity of computational modeling via adequately prepared studies. The accomplishment of synthetic validity tests is recommended for future applications. Beyond that, we propose to render the demonstration of sufficient validity via adequate simulations mandatory to computational modeling studies.

The strengths and weaknesses of inverted pendulum models of human walking.

PubMed

McGrath, Michael; Howard, David; Baker, Richard

2015-02-01

An investigation into the kinematic and kinetic predictions of two "inverted pendulum" (IP) models of gait was undertaken. The first model consisted of a single leg, with anthropometrically correct mass and moment of inertia, and a point mass at the hip representing the rest of the body. A second model incorporating the physiological extension of a head-arms-trunk (HAT) segment, held upright by an actuated hip moment, was developed for comparison. Simulations were performed, using both models, and quantitatively compared with empirical gait data. There was little difference between the two models' predictions of kinematics and ground reaction force (GRF). The models agreed well with empirical data through mid-stance (20-40% of the gait cycle) suggesting that IP models adequately simulate this phase (mean error less than one standard deviation). IP models are not cyclic, however, and cannot adequately simulate double support and step-to-step transition. This is because the forces under both legs augment each other during double support to increase the vertical GRF. The incorporation of an actuated hip joint was the most novel change and added a new dimension to the classic IP model. The hip moment curve produced was similar to those measured during experimental walking trials. As a result, it was interpreted that the primary role of the hip musculature in stance is to keep the HAT upright. Careful consideration of the differences between the models throws light on what the different terms within the GRF equation truly represent. Copyright © 2014 Elsevier B.V. All rights reserved.

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

Taylor, Paul A.; Cooper, Candice Frances; Burnett, Damon J.

Light body armor development for the warfighter is based on trial-and-error testing of prototype designs against ballistic projectiles. Torso armor testing against blast is virtually nonexistent but necessary to ensure adequate protection against injury to the heart and lungs. In this report, we discuss the development of a high-fidelity human torso model, it's merging with the existing Sandia Human Head-Neck Model, and development of the modeling & simulation (M&S) capabilities necessary to simulate wound injury scenarios. Using the new Sandia Human Torso Model, we demonstrate the advantage of virtual simulation in the investigation of wound injury as it relates tomore » the warfighter experience. We present the results of virtual simulations of blast loading and ballistic projectile impact to the tors o with and without notional protective armor. In this manner, we demonstrate the ad vantages of applying a modeling and simulation approach to the investigation of wound injury and relative merit assessments of protective body armor without the need for trial-and-error testing.« less

Pharmacokinetic/Pharmacodynamic Modeling and Simulation of Cefiderocol, a Parenteral Siderophore Cephalosporin, for Dose Adjustment Based on Renal Function.

PubMed

Katsube, Takayuki; Wajima, Toshihiro; Ishibashi, Toru; Arjona Ferreira, Juan Camilo; Echols, Roger

2017-01-01

Cefiderocol, a novel parenteral siderophore cephalosporin, exhibits potent efficacy against most Gram-negative bacteria, including carbapenem-resistant strains. Since cefiderocol is excreted primarily via the kidneys, this study was conducted to develop a population pharmacokinetics (PK) model to determine dose adjustment based on renal function. Population PK models were developed based on data for cefiderocol concentrations in plasma, urine, and dialysate with a nonlinear mixed-effects model approach. Monte-Carlo simulations were conducted to calculate the probability of target attainment (PTA) of fraction of time during the dosing interval where the free drug concentration in plasma exceeds the MIC (T f >MIC ) for an MIC range of 0.25 to 16 μg/ml. For the simulations, dose regimens were selected to compare cefiderocol exposure among groups with different levels of renal function. The developed models well described the PK of cefiderocol for each renal function group. A dose of 2 g every 8 h with 3-h infusions provided >90% PTA for 75% T f >MIC for an MIC of ≤4 μg/ml for patients with normal renal function, while a more frequent dose (every 6 h) could be used for patients with augmented renal function. A reduced dose and/or extended dosing interval was selected for patients with impaired renal function. A supplemental dose immediately after intermittent hemodialysis was proposed for patients requiring intermittent hemodialysis. The PK of cefiderocol could be adequately modeled, and the modeling-and-simulation approach suggested dose regimens based on renal function, ensuring drug exposure with adequate bactericidal effect. Copyright © 2016 American Society for Microbiology.

Adapting the Water Erosion Prediction Project (WEPP) model for forest applications

Treesearch

Shuhui Dun; Joan Q. Wu; William J. Elliot; Peter R. Robichaud; Dennis C. Flanagan; James R. Frankenberger; Robert E. Brown; Arthur C. Xu

2009-01-01

There has been an increasing public concern over forest stream pollution by excessive sedimentation due to natural or human disturbances. Adequate erosion simulation tools are needed for sound management of forest resources. The Water Erosion Prediction Project (WEPP) watershed model has proved useful in forest applications where Hortonian flow is the major form of...

Computer simulations of liquid crystals: Defects, deformations and dynamics

NASA Astrophysics Data System (ADS)

Billeter, Jeffrey Lee

1999-11-01

Computer simulations play an increasingly important role in investigating fundamental issues in the physics of liquid crystals. Presented here are the results of three projects which utilize the unique power of simulations to probe questions which neither theory nor experiment can adequately answer. Throughout, we use the (generalized) Gay-Berne model, a widely-used phenomenological potential which captures the essential features of the anisotropic mesogen shapes and interactions. First, we used a Molecular Dynamics simulation with 65536 Gay-Berne particles to study the behaviors of topological defects in a quench from the isotropic to the nematic phase. Twist disclination loops were the dominant defects, and we saw evidence for dynamical scaling. We observed the loops separating, combining and collapsing, and we also observed numerous non-singular type-1 lines which appeared to be intimately involved with many of the loop processes. Second, we used a Molecular Dynamics simulation of a sphere embedded in a system of 2048 Gay-Berne particles to study the effects of radial anchoring of the molecules at the sphere's surface. A saturn ring defect configuration was observed, and the ring caused a driven sphere (modelling the falling ball experiment) to experience an increased resistance as it moved through the nematic. Deviations from a linear relationship between the driving force and the terminal speed are attributed to distortions of the saturn ring which we observed. The existence of the saturn ring confirms theoretical predictions for small spheres. Finally, we constructed a model for wedge-shaped molecules and used a linear response approach in a Monte Carlo simulation to investigate the flexoelectric behavior of a system of 256 such wedges. Novel potential models as well as novel analytical and visualization techniques were developed for these projects. Once again, the emphasis throughout was to investigate questions which simulations alone can adequately answer.

Local rules simulation of the kinetics of virus capsid self-assembly.

PubMed

Schwartz, R; Shor, P W; Prevelige, P E; Berger, B

1998-12-01

A computer model is described for studying the kinetics of the self-assembly of icosahedral viral capsids. Solution of this problem is crucial to an understanding of the viral life cycle, which currently cannot be adequately addressed through laboratory techniques. The abstract simulation model employed to address this is based on the local rules theory of. Proc. Natl. Acad. Sci. USA. 91:7732-7736). It is shown that the principle of local rules, generalized with a model of kinetics and other extensions, can be used to simulate complicated problems in self-assembly. This approach allows for a computationally tractable molecular dynamics-like simulation of coat protein interactions while retaining many relevant features of capsid self-assembly. Three simple simulation experiments are presented to illustrate the use of this model. These show the dependence of growth and malformation rates on the energetics of binding interactions, the tolerance of errors in binding positions, and the concentration of subunits in the examples. These experiments demonstrate a tradeoff within the model between growth rate and fidelity of assembly for the three parameters. A detailed discussion of the computational model is also provided.

Flight Testing an Iced Business Jet for Flight Simulation Model Validation

NASA Technical Reports Server (NTRS)

Ratvasky, Thomas P.; Barnhart, Billy P.; Lee, Sam; Cooper, Jon

2007-01-01

A flight test of a business jet aircraft with various ice accretions was performed to obtain data to validate flight simulation models developed through wind tunnel tests. Three types of ice accretions were tested: pre-activation roughness, runback shapes that form downstream of the thermal wing ice protection system, and a wing ice protection system failure shape. The high fidelity flight simulation models of this business jet aircraft were validated using a software tool called "Overdrive." Through comparisons of flight-extracted aerodynamic forces and moments to simulation-predicted forces and moments, the simulation models were successfully validated. Only minor adjustments in the simulation database were required to obtain adequate match, signifying the process used to develop the simulation models was successful. The simulation models were implemented in the NASA Ice Contamination Effects Flight Training Device (ICEFTD) to enable company pilots to evaluate flight characteristics of the simulation models. By and large, the pilots confirmed good similarities in the flight characteristics when compared to the real airplane. However, pilots noted pitch up tendencies at stall with the flaps extended that were not representative of the airplane and identified some differences in pilot forces. The elevator hinge moment model and implementation of the control forces on the ICEFTD were identified as a driver in the pitch ups and control force issues, and will be an area for future work.

Spherical harmonic analysis of a model-generated climatology

NASA Technical Reports Server (NTRS)

Christidis, Z. D.; Spar, J.

1981-01-01

Monthly mean fields of 850 mb temperature (T850), 500 mb geopotential height (G500) and sea level pressure (SLP) were generated in the course of a five-year climate simulation run with a global general circulation model. Both the model-generated climatology and an observed climatology were subjected to spherical harmonic analysis, with separate analyses of the globe and the Northern Hemisphere. Comparison of the dominant harmonics of the two climatologies indicates that more than 95% of the area-weighted spatial variance of G500 and more than 90% of that of T850 are explained by fewer than three components, and that the model adequately simulates these large-scale characteristics. On the other hand, as many as 25 harmonics are needed to explain 95% of the observed variance of SLP, and the model simulation of this field is much less satisfactory. The model climatology is also evaluated in terms of the annual cycles of the dominant harmonics.

A simple electric circuit model for proton exchange membrane fuel cells

NASA Astrophysics Data System (ADS)

Lazarou, Stavros; Pyrgioti, Eleftheria; Alexandridis, Antonio T.

A simple and novel dynamic circuit model for a proton exchange membrane (PEM) fuel cell suitable for the analysis and design of power systems is presented. The model takes into account phenomena like activation polarization, ohmic polarization, and mass transport effect present in a PEM fuel cell. The proposed circuit model includes three resistors to approach adequately these phenomena; however, since for the PEM dynamic performance connection or disconnection of an additional load is of crucial importance, the proposed model uses two saturable inductors accompanied by an ideal transformer to simulate the double layer charging effect during load step changes. To evaluate the effectiveness of the proposed model its dynamic performance under load step changes is simulated. Experimental results coming from a commercial PEM fuel cell module that uses hydrogen from a pressurized cylinder at the anode and atmospheric oxygen at the cathode, clearly verify the simulation results.

Rainfall runoff modelling of the Upper Ganga and Brahmaputra basins using PERSiST.

PubMed

Futter, M N; Whitehead, P G; Sarkar, S; Rodda, H; Crossman, J

2015-06-01

There are ongoing discussions about the appropriate level of complexity and sources of uncertainty in rainfall runoff models. Simulations for operational hydrology, flood forecasting or nutrient transport all warrant different levels of complexity in the modelling approach. More complex model structures are appropriate for simulations of land-cover dependent nutrient transport while more parsimonious model structures may be adequate for runoff simulation. The appropriate level of complexity is also dependent on data availability. Here, we use PERSiST; a simple, semi-distributed dynamic rainfall-runoff modelling toolkit to simulate flows in the Upper Ganges and Brahmaputra rivers. We present two sets of simulations driven by single time series of daily precipitation and temperature using simple (A) and complex (B) model structures based on uniform and hydrochemically relevant land covers respectively. Models were compared based on ensembles of Bayesian Information Criterion (BIC) statistics. Equifinality was observed for parameters but not for model structures. Model performance was better for the more complex (B) structural representations than for parsimonious model structures. The results show that structural uncertainty is more important than parameter uncertainty. The ensembles of BIC statistics suggested that neither structural representation was preferable in a statistical sense. Simulations presented here confirm that relatively simple models with limited data requirements can be used to credibly simulate flows and water balance components needed for nutrient flux modelling in large, data-poor basins.

A molecular dynamics description of the conformational flexibility of the L-iduronate ring in glycosaminoglycans.

PubMed

Angulo, Jesús; Nieto, Pedro M; Martín-Lomas, Manuel

2003-07-07

For a synthetic hexasaccharide model it is shown that the conformational flexibility of the L-iduronate ring in glycosaminoglycans can be adequately described by using the PME methodology together with simulation protocols suitable for highly charged systems.

Numerical Simulation of Slag Eye Formation and Slag Entrapment in a Bottom-Blown Argon-Stirred Ladle

NASA Astrophysics Data System (ADS)

Liu, Wei; Tang, Haiyan; Yang, Shufeng; Wang, Minghui; Li, Jingshe; Liu, Qing; Liu, Jianhui

2018-06-01

A transient mathematical model is developed for simulating the bubble-steel-slag-top gas four-phase flow in a bottom-blown argon-stirred ladle with a 70-ton capacity. The Lagrangian discrete phase model (DPM) is used for describing the moving behavior of bubbles in the steel and slag. To observe the formation process of slag eye, the volume of fluid (VOF) model is used to track the interfaces between three incompressible phases: metal/slag, metal/gas, and slag/gas. The complex multiphase turbulent flow induced by bubble-liquid interactions is solved by a large eddy simulation (LES) model. Slag eye area and slag droplet dispersion are investigated under different gas flow rates. The results show that the movement of bubbles, formation and collapse of slag eye, volatility of steel/slag interface and behavior of slag entrapment can be properly predicted in the current model. When the gas flow rate is 300 L/min, the circulation driven by the bubble plume will stir the entire ladle adequately and form a slag eye of the right size. At the same time, it will not cause strong erosion to the ladle wall, and the fluctuation of the interface is of adequate intensity, which will be helpful for improving the desulfurization efficiency; the slag entrapment behavior can also be decreased. Interestingly, with the motion of liquid steel circulation, the collision and coalescence of dispersed slag droplets occur during the floating process in the vicinity of the wall.

Application of an Integrated HPC Reliability Prediction Framework to HMMWV Suspension System

DTIC Science & Technology

2010-09-13

model number M966 (TOW Missle Carrier, Basic Armor without weapons), since they were available. Tires used for all simulations were the bias-type...vehicle fleet, including consideration of all kinds of uncertainty, especially including model uncertainty. The end result will be a tool to use...building an adequate vehicle reliability prediction framework for military vehicles is the accurate modeling of the integration of various types of

Analysis of the impact of simulation model simplifications on the quality of low-energy buildings simulation results

NASA Astrophysics Data System (ADS)

Klimczak, Marcin; Bojarski, Jacek; Ziembicki, Piotr; Kęskiewicz, Piotr

2017-11-01

The requirements concerning energy performance of buildings and their internal installations, particularly HVAC systems, have been growing continuously in Poland and all over the world. The existing, traditional calculation methods following from the static heat exchange model are frequently not sufficient for a reasonable heating design of a building. Both in Poland and elsewhere in the world, methods and software are employed which allow a detailed simulation of the heating and moisture conditions in a building, and also an analysis of the performance of HVAC systems within a building. However, these systems are usually difficult in use and complex. In addition, the development of a simulation model that is sufficiently adequate to the real building requires considerable time involvement of a designer, is time-consuming and laborious. A simplification of the simulation model of a building renders it possible to reduce the costs of computer simulations. The paper analyses in detail the effect of introducing a number of different variants of the simulation model developed in Design Builder on the quality of final results obtained. The objective of this analysis is to find simplifications which allow obtaining simulation results which have an acceptable level of deviations from the detailed model, thus facilitating a quick energy performance analysis of a given building.

Numerical simulation of small-scale thermal convection in the atmosphere

NASA Technical Reports Server (NTRS)

Somerville, R. C. J.

1973-01-01

A Boussinesq system is integrated numerically in three dimensions and time in a study of nonhydrostatic convection in the atmosphere. Simulation of cloud convection is achieved by the inclusion of parametrized effects of latent heat and small-scale turbulence. The results are compared with the cell structure observed in Rayleigh-Benard laboratory conversion experiments in air. At a Rayleigh number of 4000, the numerical model adequately simulates the experimentally observed evolution, including some prominent transients of a flow from a randomly perturbed initial conductive state into the final state of steady large-amplitude two-dimensional rolls. At Rayleigh number 9000, the model reproduces the experimentally observed unsteady equilibrium of vertically coherent oscillatory waves superimposed on rolls.

Nucleic acid polymeric properties and electrostatics: Directly comparing theory and simulation with experiment.

PubMed

Sim, Adelene Y L

2016-06-01

Nucleic acids are biopolymers that carry genetic information and are also involved in various gene regulation functions such as gene silencing and protein translation. Because of their negatively charged backbones, nucleic acids are polyelectrolytes. To adequately understand nucleic acid folding and function, we need to properly describe its i) polymer/polyelectrolyte properties and ii) associating ion atmosphere. While various theories and simulation models have been developed to describe nucleic acids and the ions around them, many of these theories/simulations have not been well evaluated due to complexities in comparison with experiment. In this review, I discuss some recent experiments that have been strategically designed for straightforward comparison with theories and simulation models. Such data serve as excellent benchmarks to identify limitations in prevailing theories and simulation parameters. Copyright © 2015 Elsevier B.V. All rights reserved.

Nursing simulation: a community experience.

PubMed

Gunowa, Neesha Oozageer; Elliott, Karen; McBride, Michelle

2018-04-02

The education sector faces major challenges in providing learning experiences so that newly qualified nurses feel adequately prepared to work in a community setting. With this in mind, higher education institutions need to develop more innovative ways to deliver the community-nurse experience to student nurses. This paper presents and explores how simulation provides an opportunity for educators to support and evaluate student performance in an environment that models a complete patient encounter in the community. Following the simulation, evaluative data were collated and the answers analysed to identify key recommendations.

Microscopic transport model animation visualisation on KML base

NASA Astrophysics Data System (ADS)

Yatskiv, I.; Savrasovs, M.

2012-10-01

By reading classical literature devoted to the simulation theory it could be found that one of the greatest possibilities of simulation is the ability to present processes inside the system by animation. This gives to the simulation model additional value during presentation of simulation results for the public and authorities who are not familiar enough with simulation. That is why most of universal and specialised simulation tools have the ability to construct 2D and 3D representation of the model. Usually the development of such representation could take much time and there must be put a lot forces into creating an adequate 3D representation of the model. For long years such well-known microscopic traffic flow simulation software tools as VISSIM, AIMSUN and PARAMICS have had a possibility to produce 2D and 3D animation. But creation of realistic 3D model of the place where traffic flows are simulated, even in these professional software tools it is a hard and time consuming action. The goal of this paper is to describe the concepts of use the existing on-line geographical information systems for visualisation of animation produced by simulation software. For demonstration purposes the following technologies and tools have been used: PTV VISION VISSIM, KML and Google Earth.

Simulation of Chronic Liver Injury Due to Environmental Chemicals

EPA Science Inventory

US EPA Virtual Liver (v-Liver) is a cellular systems model of hepatic tissues to predict the effects of chronic exposure to chemicals. Tens of thousands of chemicals are currently in commerce and hundreds more are introduced every year. Few of these chemicals have been adequate...

A subsurface drip irrigation system for weighing lysimetry

USDA-ARS?s Scientific Manuscript database

Large, precision weighing lysimeters can have accuracies as good as 0.04 mm equivalent depth of water, adequate for hourly and even half-hourly determinations of evapotranspiration (ET) rate from crops. Such data are important for testing and improving simulation models of the complex interactions o...

[The history of development of evolutionary methods in St. Petersburg school of computer simulation in biology].

PubMed

Menshutkin, V V; Kazanskiĭ, A B; Levchenko, V F

2010-01-01

The history of rise and development of evolutionary methods in Saint Petersburg school of biological modelling is traced and analyzed. Some pioneering works in simulation of ecological and evolutionary processes, performed in St.-Petersburg school became an exemplary ones for many followers in Russia and abroad. The individual-based approach became the crucial point in the history of the school as an adequate instrument for construction of models of biological evolution. This approach is natural for simulation of the evolution of life-history parameters and adaptive processes in populations and communities. In some cases simulated evolutionary process was used for solving a reverse problem, i. e., for estimation of uncertain life-history parameters of population. Evolutionary computations is one more aspect of this approach application in great many fields. The problems and vistas of ecological and evolutionary modelling in general are discussed.

The accuracy of climate models' simulated season lengths and the effectiveness of grid scale correction factors

DOE PAGES

Winterhalter, Wade E.

2011-09-01

Global climate change is expected to impact biological populations through a variety of mechanisms including increases in the length of their growing season. Climate models are useful tools for predicting how season length might change in the future. However, the accuracy of these models tends to be rather low at regional geographic scales. Here, I determined the ability of several atmosphere and ocean general circulating models (AOGCMs) to accurately simulate historical season lengths for a temperate ectotherm across the continental United States. I also evaluated the effectiveness of regional-scale correction factors to improve the accuracy of these models. I foundmore » that both the accuracy of simulated season lengths and the effectiveness of the correction factors to improve the model's accuracy varied geographically and across models. These results suggest that regional specific correction factors do not always adequately remove potential discrepancies between simulated and historically observed environmental parameters. As such, an explicit evaluation of the correction factors' effectiveness should be included in future studies of global climate change's impact on biological populations.« less

Lunar gravitational field estimation and the effects of mismodeling upon lunar satellite orbit prediction. M.S. Thesis

NASA Technical Reports Server (NTRS)

Davis, John H.

1993-01-01

Lunar spherical harmonic gravity coefficients are estimated from simulated observations of a near-circular low altitude polar orbiter disturbed by lunar mascons. Lunar gravity sensing missions using earth-based nearside observations with and without satellite-based far-side observations are simulated and least squares maximum likelihood estimates are developed for spherical harmonic expansion fit models. Simulations and parameter estimations are performed by a modified version of the Smithsonian Astrophysical Observatory's Planetary Ephemeris Program. Two different lunar spacecraft mission phases are simulated to evaluate the estimated fit models. Results for predicting state covariances one orbit ahead are presented along with the state errors resulting from the mismodeled gravity field. The position errors from planning a lunar landing maneuver with a mismodeled gravity field are also presented. These simulations clearly demonstrate the need to include observations of satellite motion over the far side in estimating the lunar gravity field. The simulations also illustrate that the eighth degree and order expansions used in the simulated fits were unable to adequately model lunar mascons.

STS/DBS power subsystem end-to-end stability margin

NASA Astrophysics Data System (ADS)

Devaux, R. N.; Vattimo, R. J.; Peck, S. R.; Baker, W. E.

Attention is given to a full-up end-to-end subsystem stability test which was performed with a flight solar array providing power to a fully operational spacecraft. The solar array simulator is described, and a comparison is made between test results obtained with the simulator and those obtained with the actual array. It is concluded that stability testing with a fully integrated spacecraft is necessary to ensure that all elements have been adequately modeled.

Predictor characteristics necessary for building a clinically useful risk prediction model: a simulation study.

PubMed

Schummers, Laura; Himes, Katherine P; Bodnar, Lisa M; Hutcheon, Jennifer A

2016-09-21

Compelled by the intuitive appeal of predicting each individual patient's risk of an outcome, there is a growing interest in risk prediction models. While the statistical methods used to build prediction models are increasingly well understood, the literature offers little insight to researchers seeking to gauge a priori whether a prediction model is likely to perform well for their particular research question. The objective of this study was to inform the development of new risk prediction models by evaluating model performance under a wide range of predictor characteristics. Data from all births to overweight or obese women in British Columbia, Canada from 2004 to 2012 (n = 75,225) were used to build a risk prediction model for preeclampsia. The data were then augmented with simulated predictors of the outcome with pre-set prevalence values and univariable odds ratios. We built 120 risk prediction models that included known demographic and clinical predictors, and one, three, or five of the simulated variables. Finally, we evaluated standard model performance criteria (discrimination, risk stratification capacity, calibration, and Nagelkerke's r 2 ) for each model. Findings from our models built with simulated predictors demonstrated the predictor characteristics required for a risk prediction model to adequately discriminate cases from non-cases and to adequately classify patients into clinically distinct risk groups. Several predictor characteristics can yield well performing risk prediction models; however, these characteristics are not typical of predictor-outcome relationships in many population-based or clinical data sets. Novel predictors must be both strongly associated with the outcome and prevalent in the population to be useful for clinical prediction modeling (e.g., one predictor with prevalence ≥20 % and odds ratio ≥8, or 3 predictors with prevalence ≥10 % and odds ratios ≥4). Area under the receiver operating characteristic curve values of >0.8 were necessary to achieve reasonable risk stratification capacity. Our findings provide a guide for researchers to estimate the expected performance of a prediction model before a model has been built based on the characteristics of available predictors.

Molecular-dynamics simulation of mutual diffusion in nonideal liquid mixtures

NASA Astrophysics Data System (ADS)

Rowley, R. L.; Stoker, J. M.; Giles, N. F.

1991-05-01

The mutual-diffusion coefficients, D 12, of n-hexane, n-heptane, and n-octane in chloroform were modeled using equilibrium molecular-dynamics (MD) simulations of simple Lennard-Jones (LJ) fluids. Pure-component LJ parameters were obtained by comparison of simulations to experimental self-diffusion coefficients. While values of “effective” LJ parameters are not expected to simulate accurately diverse thermophysical properties over a wide range of conditions, it was recently shown that effective parameters obtained from pure self-diffusion coefficients can accurately model mutual diffusion in ideal, liquid mixtures. In this work, similar simulations are used to model diffusion in nonideal mixtures. The same combining rules used in the previous study for the cross-interaction parameters were found to be adequate to represent the composition dependence of D 12. The effect of alkane chain length on D 12 is also correctly predicted by the simulations. A commonly used assumption in empirical correlations of D 12, that its kinetic portion is a simple, compositional average of the intradiffusion coefficients, is inconsistent with the simulation results. In fact, the value of the kinetic portion of D 12 was often outside the range of values bracketed by the two intradiffusion coefficients for the nonideal system modeled here.

The use of discrete-event simulation modeling to compare handwritten and electronic prescribing systems.

PubMed

Ghany, Ahmad; Vassanji, Karim; Kuziemsky, Craig; Keshavjee, Karim

2013-01-01

Electronic prescribing (e-prescribing) is expected to bring many benefits to Canadian healthcare, such as a reduction in errors and adverse drug reactions. As there currently is no functioning e-prescribing system in Canada that is completely electronic, we are unable to evaluate the performance of a live system. An alternative approach is to use simulation modeling for evaluation. We developed two discrete-event simulation models, one of the current handwritten prescribing system and one of a proposed e-prescribing system, to compare the performance of these two systems. We were able to compare the number of processes in each model, workflow efficiency, and the distribution of patients or prescriptions. Although we were able to compare these models to each other, using discrete-event simulation software was challenging. We were limited in the number of variables we could measure. We discovered non-linear processes and feedback loops in both models that could not be adequately represented using discrete-event simulation software. Finally, interactions between entities in both models could not be modeled using this type of software. We have come to the conclusion that a more appropriate approach to modeling both the handwritten and electronic prescribing systems would be to use a complex adaptive systems approach using agent-based modeling or systems-based modeling.

Toolbox for Urban Mobility Simulation: High Resolution Population Dynamics for Global Cities

NASA Astrophysics Data System (ADS)

Bhaduri, B. L.; Lu, W.; Liu, C.; Thakur, G.; Karthik, R.

2015-12-01

In this rapidly urbanizing world, unprecedented rate of population growth is not only mirrored by increasing demand for energy, food, water, and other natural resources, but has detrimental impacts on environmental and human security. Transportation simulations are frequently used for mobility assessment in urban planning, traffic operation, and emergency management. Previous research, involving purely analytical techniques to simulations capturing behavior, has investigated questions and scenarios regarding the relationships among energy, emissions, air quality, and transportation. Primary limitations of past attempts have been availability of input data, useful "energy and behavior focused" models, validation data, and adequate computational capability that allows adequate understanding of the interdependencies of our transportation system. With increasing availability and quality of traditional and crowdsourced data, we have utilized the OpenStreetMap roads network, and has integrated high resolution population data with traffic simulation to create a Toolbox for Urban Mobility Simulations (TUMS) at global scale. TUMS consists of three major components: data processing, traffic simulation models, and Internet-based visualizations. It integrates OpenStreetMap, LandScanTM population, and other open data (Census Transportation Planning Products, National household Travel Survey, etc.) to generate both normal traffic operation and emergency evacuation scenarios. TUMS integrates TRANSIMS and MITSIM as traffic simulation engines, which are open-source and widely-accepted for scalable traffic simulations. Consistent data and simulation platform allows quick adaption to various geographic areas that has been demonstrated for multiple cities across the world. We are combining the strengths of geospatial data sciences, high performance simulations, transportation planning, and emissions, vehicle and energy technology development to design and develop a simulation framework to assist decision makers at all levels - local, state, regional, and federal. Using Cleveland, Tennessee as an example, in this presentation, we illustrate how emerging cities could easily assess future land use scenario driven impacts on energy and environment utilizing such a capability.

Density and white light brightness in looplike coronal mass ejections - Temporal evolution

NASA Technical Reports Server (NTRS)

Steinolfson, R. S.; Hundhausen, A. J.

1988-01-01

Three ambient coronal models suitable for studies of time-dependent phenomena were used to investigate the propagation of coronal mass ejections initiated in each atmosphere by an identical energy source. These models included those of a static corona with a dipole magnetic field, developed by Dryer et al. (1979); a steady polytropic corona with an equatorial coronal streamer, developed by Steinolfson et al. (1982); and Steinolfson's (1988) model of heated corona with an equatorial coronal streamer. The results indicated that the first model does not adequately represent the general characteristics of observed looplike mass ejections, and the second model simulated only some of the observed features. Only the third model, which included a heating term and a streamer, was found to yield accurate simulation of the mess ejection observations.

Defining metrics of the Quasi-Biennial Oscillation in global climate models

NASA Astrophysics Data System (ADS)

Schenzinger, Verena; Osprey, Scott; Gray, Lesley; Butchart, Neal

2017-06-01

As the dominant mode of variability in the tropical stratosphere, the Quasi-Biennial Oscillation (QBO) has been subject to extensive research. Though there is a well-developed theory of this phenomenon being forced by wave-mean flow interaction, simulating the QBO adequately in global climate models still remains difficult. This paper presents a set of metrics to characterize the morphology of the QBO using a number of different reanalysis datasets and the FU Berlin radiosonde observation dataset. The same metrics are then calculated from Coupled Model Intercomparison Project 5 and Chemistry-Climate Model Validation Activity 2 simulations which included a representation of QBO-like behaviour to evaluate which aspects of the QBO are well captured by the models and which ones remain a challenge for future model development.

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

Zhang, Fan; Parker, Jack C.; Watson, David B

This study investigates uranium and technetium sorption onto aluminum and iron hydroxides during titration of acidic groundwater. The contaminated groundwater exhibits oxic conditions with high concentrations of NO{sub 3}{sup -}, SO{sub 4}{sup 2-}, U, Tc, and various metal cations. More than 90% of U and Tc was removed from the aqueous phase as Al and Fe precipitated above pH 5.5, but was partially resolublized at higher pH values. An equilibrium hydrolysis and precipitation reaction model adequately described variations in aqueous concentrations of metal cations. An anion exchange reaction model was incorporated to simulate sulfate, U and Tc sorption onto variablymore » charged (pH-dependent) Al and Fe hydroxides. Modeling results indicate that competitive sorption/desorption on mixed mineral phases needs to be considered to adequately predict U and Tc mobility. The model could be useful for future studies of the speciation of U, Tc and co-existing ions during pre- and post-groundwater treatment practices.« less

Water injection into vapor- and liquid-dominated reservoirs: Modeling of heat transfer and mass transport

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

Pruess, K.; Oldenburg, C.; Moridis, G.

1997-12-31

This paper summarizes recent advances in methods for simulating water and tracer injection, and presents illustrative applications to liquid- and vapor-dominated geothermal reservoirs. High-resolution simulations of water injection into heterogeneous, vertical fractures in superheated vapor zones were performed. Injected water was found to move in dendritic patterns, and to experience stronger lateral flow effects than predicted from homogeneous medium models. Higher-order differencing methods were applied to modeling water and tracer injection into liquid-dominated systems. Conventional upstream weighting techniques were shown to be adequate for predicting the migration of thermal fronts, while higher-order methods give far better accuracy for tracer transport.more » A new fluid property module for the TOUGH2 simulator is described which allows a more accurate description of geofluids, and includes mineral dissolution and precipitation effects with associated porosity and permeability change. Comparisons between numerical simulation predictions and data for laboratory and field injection experiments are summarized. Enhanced simulation capabilities include a new linear solver package for TOUGH2, and inverse modeling techniques for automatic history matching and optimization.« less

Controls/CFD Interdisciplinary Research Software Generates Low-Order Linear Models for Control Design From Steady-State CFD Results

NASA Technical Reports Server (NTRS)

Melcher, Kevin J.

1997-01-01

The NASA Lewis Research Center is developing analytical methods and software tools to create a bridge between the controls and computational fluid dynamics (CFD) disciplines. Traditionally, control design engineers have used coarse nonlinear simulations to generate information for the design of new propulsion system controls. However, such traditional methods are not adequate for modeling the propulsion systems of complex, high-speed vehicles like the High Speed Civil Transport. To properly model the relevant flow physics of high-speed propulsion systems, one must use simulations based on CFD methods. Such CFD simulations have become useful tools for engineers that are designing propulsion system components. The analysis techniques and software being developed as part of this effort are an attempt to evolve CFD into a useful tool for control design as well. One major aspect of this research is the generation of linear models from steady-state CFD results. CFD simulations, often used during the design of high-speed inlets, yield high resolution operating point data. Under a NASA grant, the University of Akron has developed analytical techniques and software tools that use these data to generate linear models for control design. The resulting linear models have the same number of states as the original CFD simulation, so they are still very large and computationally cumbersome. Model reduction techniques have been successfully applied to reduce these large linear models by several orders of magnitude without significantly changing the dynamic response. The result is an accurate, easy to use, low-order linear model that takes less time to generate than those generated by traditional means. The development of methods for generating low-order linear models from steady-state CFD is most complete at the one-dimensional level, where software is available to generate models with different kinds of input and output variables. One-dimensional methods have been extended somewhat so that linear models can also be generated from two- and three-dimensional steady-state results. Standard techniques are adequate for reducing the order of one-dimensional CFD-based linear models. However, reduction of linear models based on two- and three-dimensional CFD results is complicated by very sparse, ill-conditioned matrices. Some novel approaches are being investigated to solve this problem.

Development of the CCP-200 mathematical model for Syzran CHPP using the Thermolib software package

NASA Astrophysics Data System (ADS)

Usov, S. V.; Kudinov, A. A.

2016-04-01

Simplified cycle diagram of the CCP-200 power generating unit of Syzran CHPP containing two gas turbines PG6111FA with generators, two steam recovery boilers KUP-110/15-8.0/0.7-540/200, and one steam turbine Siemens SST-600 (one-cylinder with two variable heat extraction units of 60/75 MW in heatextraction and condensing modes, accordingly) with S-GEN5-100 generators was presented. Results of experimental guarantee tests of the CCP-200 steam-gas unit are given. Brief description of the Thermolib application for the MatLab Simulink software package is given. Basic equations used in Thermolib for modeling thermo-technical processes are given. Mathematical models of gas-turbine plant, heat-recovery steam generator, steam turbine and integrated plant for power generating unit CCP-200 of Syzran CHPP were developed with the help of MatLab Simulink and Thermolib. The simulation technique at different ambient temperature values was used in order to get characteristics of the developed mathematical model. Graphic comparison of some characteristics of the CCP-200 simulation model (gas temperature behind gas turbine, gas turbine and combined cycle plant capacity, high and low pressure steam consumption and feed water consumption for high and low pressure economizers) with actual characteristics of the steam-gas unit received at experimental (field) guarantee tests at different ambient temperature are shown. It is shown that the chosen degrees of complexity, characteristics of the CCP-200 simulation model, developed by Thermolib, adequately correspond to the actual characteristics of the steam-gas unit received at experimental (field) guarantee tests; this allows considering the developed mathematical model as adequate and acceptable it for further work.

Robot, computer problem solving system

NASA Technical Reports Server (NTRS)

Becker, J. D.

1972-01-01

The development of a computer problem solving system is reported that considers physical problems faced by an artificial robot moving around in a complex environment. Fundamental interaction constraints with a real environment are simulated for the robot by visual scan and creation of an internal environmental model. The programming system used in constructing the problem solving system for the simulated robot and its simulated world environment is outlined together with the task that the system is capable of performing. A very general framework for understanding the relationship between an observed behavior and an adequate description of that behavior is included.

Real-time simulation for intra-operative navigation in robotic surgery. Using a mass spring system for a basic study of organ deformation.

PubMed

Kawamura, Kazuya; Kobayashi, Yo; Fujie, Masakatsu G

2007-01-01

Medical technology has advanced with the introduction of robot technology, making previous medical treatments that were very difficult far more possible. However, operation of a surgical robot demands substantial training and continual practice on the part of the surgeon because it requires difficult techniques that are different from those of traditional surgical procedures. We focused on a simulation technology based on the physical characteristics of organs. In this research, we proposed the development of surgical simulation, based on a physical model, for intra-operative navigation by a surgeon. In this paper, we describe the design of our system, in particular our organ deformation calculator. The proposed simulation system consists of an organ deformation calculator and virtual slave manipulators. We obtained adequate experimental results of a target node at a nearby point of interaction, because this point ensures better accuracy for our simulation model. The next research step would be to focus on a surgical environment in which internal organ models would be integrated into a slave simulation system.

Optimisation and evaluation of pre-design models for offshore wind turbines with jacket support structures and their influence on integrated load simulations

NASA Astrophysics Data System (ADS)

Schafhirt, S.; Kaufer, D.; Cheng, P. W.

2014-12-01

In recent years many advanced load simulation tools, allowing an aero-servo-hydroelastic analyses of an entire offshore wind turbine, have been developed and verified. Nowadays, even an offshore wind turbine with a complex support structure such as a jacket can be analysed. However, the computational effort rises significantly with an increasing level of details. This counts especially for offshore wind turbines with lattice support structures, since those models do naturally have a higher number of nodes and elements than simpler monopile structures. During the design process multiple load simulations are demanded to obtain an optimal solution. In the view of pre-design tasks it is crucial to apply load simulations which keep the simulation quality and the computational effort in balance. The paper will introduce a reference wind turbine model consisting of the REpower5M wind turbine and a jacket support structure with a high level of detail. In total twelve variations of this reference model are derived and presented. Main focus is to simplify the models of the support structure and the foundation. The reference model and the simplified models are simulated with the coupled simulation tool Flex5-Poseidon and analysed regarding frequencies, fatigue loads, and ultimate loads. A model has been found which reaches an adequate increase of simulation speed while holding the results in an acceptable range compared to the reference results.

Conditioning geostatistical simulations of a heterogeneous paleo-fluvial bedrock aquifer using lithologs and pumping tests

NASA Astrophysics Data System (ADS)

Niazi, A.; Bentley, L. R.; Hayashi, M.

2016-12-01

Geostatistical simulations are used to construct heterogeneous aquifer models. Optimally, such simulations should be conditioned with both lithologic and hydraulic data. We introduce an approach to condition lithologic geostatistical simulations of a paleo-fluvial bedrock aquifer consisting of relatively high permeable sandstone channels embedded in relatively low permeable mudstone using hydraulic data. The hydraulic data consist of two-hour single well pumping tests extracted from the public water well database for a 250-km2 watershed in Alberta, Canada. First, lithologic models of the entire watershed are simulated and conditioned with hard lithological data using transition probability - Markov chain geostatistics (TPROGS). Then, a segment of the simulation around a pumping well is used to populate a flow model (FEFLOW) with either sand or mudstone. The values of the hydraulic conductivity and specific storage of sand and mudstone are then adjusted to minimize the difference between simulated and actual pumping test data using the parameter estimation program PEST. If the simulated pumping test data do not adequately match the measured data, the lithologic model is updated by locally deforming the lithology distribution using the probability perturbation method and the model parameters are again updated with PEST. This procedure is repeated until the simulated and measured data agree within a pre-determined tolerance. The procedure is repeated for each well that has pumping test data. The method creates a local groundwater model that honors both the lithologic model and pumping test data and provides estimates of hydraulic conductivity and specific storage. Eventually, the simulations will be integrated into a watershed-scale groundwater model.

Using simulators to teach pediatric airway procedures in an international setting.

PubMed

Schwartz, Marissa A; Kavanagh, Katherine R; Frampton, Steven J; Bruce, Iain A; Valdez, Tulio A

2018-01-01

There has been a growing shift towards endoscopic management of laryngeal procedures in pediatric otolaryngology. There still appears to be a shortage of pediatric otolaryngology programs and children's hospitals worldwide where physicians can learn and practice these skills. Laryngeal simulation models have the potential to be part of the educational training of physicians who lack exposure to relatively uncommon pediatric otolaryngologic pathology. The objective of this study was to assess the utility of pediatric laryngeal models to teach laryngeal pathology to physicians at an international meeting. Pediatric laryngeal models were assessed by participants at an international pediatric otolaryngology meeting. Participants provided demographic information and previous experience with pediatric airways. Participants then performed simulated surgery on these models and evaluated them using both a previously validated Tissue Likeness Scale and a pre-simulation to post-simulation confidence scale. Participants reported significant subjective improvement in confidence level after use of the simulation models (p < 0.05). Participants reported realistic representations of human anatomy and pathology. The models' tissue mechanics were adequate to practice operative technique including the ability to incise, suture, and suspend models. The pediatric laryngeal models demonstrate high quality anatomy, which is easy manipulated with surgical instruments. These models allow both trainees and surgeons to practice time-sensitive airway surgeries in a safe and controlled environment. Copyright © 2017 Elsevier B.V. All rights reserved.

Small-scale multi-axial hybrid simulation of a shear-critical reinforced concrete frame

NASA Astrophysics Data System (ADS)

Sadeghian, Vahid; Kwon, Oh-Sung; Vecchio, Frank

2017-10-01

This study presents a numerical multi-scale simulation framework which is extended to accommodate hybrid simulation (numerical-experimental integration). The framework is enhanced with a standardized data exchange format and connected to a generalized controller interface program which facilitates communication with various types of laboratory equipment and testing configurations. A small-scale experimental program was conducted using a six degree-of-freedom hydraulic testing equipment to verify the proposed framework and provide additional data for small-scale testing of shearcritical reinforced concrete structures. The specimens were tested in a multi-axial hybrid simulation manner under a reversed cyclic loading condition simulating earthquake forces. The physical models were 1/3.23-scale representations of a beam and two columns. A mixed-type modelling technique was employed to analyze the remainder of the structures. The hybrid simulation results were compared against those obtained from a large-scale test and finite element analyses. The study found that if precautions are taken in preparing model materials and if the shear-related mechanisms are accurately considered in the numerical model, small-scale hybrid simulations can adequately simulate the behaviour of shear-critical structures. Although the findings of the study are promising, to draw general conclusions additional test data are required.

Discrete-time pilot model. [human dynamics and digital simulation

NASA Technical Reports Server (NTRS)

Cavalli, D.

1978-01-01

Pilot behavior is considered as a discrete-time process where the decision making has a sequential nature. This model differs from both the quasilinear model which follows from classical control theory and from the optimal control model which considers the human operator as a Kalman estimator-predictor. An additional factor considered is that the pilot's objective may not be adequately formulated as a quadratic cost functional to be minimized, but rather as a more fuzzy measure of the closeness with which the aircraft follows a reference trajectory. All model parameters, in the digital program simulating the pilot's behavior, were successfully compared in terms of standard-deviation and performance with those of professional pilots in IFR configuration. The first practical application of the model was in the study of its performance degradation when the aircraft model static margin decreases.

78 FR 72612 - Criteria for the Certification and Recertification of the Waste Isolation Pilot Plant's...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-12-03

... WIPP PA process culminates in a series of computer simulations that model the physical attributes of... and Processes LWA Land Withdrawal Act MSHA Mine Safety and Health Administration NMED New Mexico... Agency's technical review process was to determine whether, with the new design, the WIPP adequately...

Simulation of κ-Carbide Precipitation Kinetics in Aged Low-Density Fe-Mn-Al-C Steels and Its Effects on Strengthening

NASA Astrophysics Data System (ADS)

Lee, Jaeeun; Park, Siwook; Kim, Hwangsun; Park, Seong-Jun; Lee, Keunho; Kim, Mi-Young; Madakashira, Phaniraj P.; Han, Heung Nam

2018-03-01

Fe-Al-Mn-C alloy systems are low-density austenite-based steels that show excellent mechanical properties. After aging such steels at adequate temperatures for adequate time, nano-scale precipitates such as κ-carbide form, which have profound effects on the mechanical properties. Therefore, it is important to predict the amount and size of the generated κ-carbide precipitates in order to control the mechanical properties of low-density steels. In this study, the microstructure and mechanical properties of aged low-density austenitic steel were characterized. Thermo-kinetic simulations of the aging process were used to predict the size and phase fraction of κ-carbide after different aging periods, and these results were validated by comparison with experimental data derived from dark-field transmission electron microscopy images. Based on these results, models for precipitation strengthening based on different mechanisms were assessed. The measured increase in the strength of aged specimens was compared with that calculated from the models to determine the exact precipitation strengthening mechanism.

Modelling of slaughterhouse solid waste anaerobic digestion: determination of parameters and continuous reactor simulation.

PubMed

López, Iván; Borzacconi, Liliana

2010-10-01

A model based on the work of Angelidaki et al. (1993) was applied to simulate the anaerobic biodegradation of ruminal contents. In this study, two fractions of solids with different biodegradation rates were considered. A first-order kinetic was used for the easily biodegradable fraction and a kinetic expression that is function of the extracellular enzyme concentration was used for the slowly biodegradable fraction. Batch experiments were performed to obtain an accumulated methane curve that was then used to obtain the model parameters. For this determination, a methodology derived from the "multiple-shooting" method was successfully used. Monte Carlo simulations allowed a confidence range to be obtained for each parameter. Simulations of a continuous reactor were performed using the optimal set of model parameters. The final steady-states were determined as functions of the operational conditions (solids load and residence time). The simulations showed that methane flow peaked at a flow rate of 0.5-0.8 Nm(3)/d/m(reactor)(3) at a residence time of 10-20 days. Simulations allow the adequate selection of operating conditions of a continuous reactor. (c) 2010 Elsevier Ltd. All rights reserved.

Knowledge representation and qualitative simulation of salmon redd functioning. Part II: qualitative model of redds.

PubMed

Guerrin, F; Dumas, J

2001-02-01

This paper describes a qualitative model of the functioning of salmon redds (spawning areas of salmon) and its impact on mortality rates of early stages. For this, we use Qsim, a qualitative simulator, which appeared adequate for representing available qualitative knowledge of freshwater ecology experts (see Part I of this paper). Since the number of relevant variables was relatively large, it appeared necessary to decompose the model into two parts, corresponding to processes occurring at separate time-scales. A qualitative clock allows us to submit the simulation of salmon developmental stages to the calculation of accumulated daily temperatures (degree-days), according to the clock ticks and a water temperature regime set by the user. Therefore, this introduces some way of real-time dating and duration in a purely qualitative model. Simulating both sub-models, either separately or by means of alternate transitions, allows us to generate the evolutions of variables of interest, such as the mortality rates according to two factors (flow of oxygenated water and plugging of gravel interstices near the bed surface), under various scenarios.

A formal language for the specification and verification of synchronous and asynchronous circuits

NASA Technical Reports Server (NTRS)

Russinoff, David M.

1993-01-01

A formal hardware description language for the intended application of verifiable asynchronous communication is described. The language is developed within the logical framework of the Nqthm system of Boyer and Moore and is based on the event-driven behavioral model of VHDL, including the basic VHDL signal propagation mechanisms, the notion of simulation deltas, and the VHDL simulation cycle. A core subset of the language corresponds closely with a subset of VHDL and is adequate for the realistic gate-level modeling of both combinational and sequential circuits. Various extensions to this subset provide means for convenient expression of behavioral circuit specifications.

Energy decay in a granular gas collapse

NASA Astrophysics Data System (ADS)

Almazán, Lidia; Serero, Dan; Salueña, Clara; Pöschel, Thorsten

2017-01-01

An inelastic hard ball bouncing repeatedly off the ground comes to rest in finite time by performing an infinite number of collisions. Similarly, a granular gas under the influence of external gravity, condenses at the bottom of the confinement due to inelastic collisions. By means of hydrodynamical simulations, we find that the condensation process of a granular gas reveals a similar dynamics as the bouncing ball. Our result is in agreement with both experiments and particle simulations, but disagrees with earlier simplified hydrodynamical description. Analyzing the result in detail, we find that the adequate modeling of pressure plays a key role in continuum modeling of granular matter.

Simulation loop between cad systems, GEANT-4 and GeoModel: Implementation and results

NASA Astrophysics Data System (ADS)

Sharmazanashvili, A.; Tsutskiridze, Niko

2016-09-01

Compare analysis of simulation and as-built geometry descriptions of detector is important field of study for data_vs_Monte-Carlo discrepancies. Shapes consistency and detalization is not important while adequateness of volumes and weights of detector components are essential for tracking. There are 2 main reasons of faults of geometry descriptions in simulation: (1) Difference between simulated and as-built geometry descriptions; (2) Internal inaccuracies of geometry transformations added by simulation software infrastructure itself. Georgian Engineering team developed hub on the base of CATIA platform and several tools enabling to read in CATIA different descriptions used by simulation packages, like XML->CATIA; VP1->CATIA; Geo-Model->CATIA; Geant4->CATIA. As a result it becomes possible to compare different descriptions with each other using the full power of CATIA and investigate both classes of reasons of faults of geometry descriptions. Paper represents results of case studies of ATLAS Coils and End-Cap toroid structures.

Nonequilibrium Nonideal Nanoplasma Generated by a Fast Single Ion in Condensed Matter

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

Faenov, A. Ya.; Kansai Photon Science Institut, Japan Atomic Energy Agency; Lankin, A. V.

A plasma model of relaxation of a medium in heavy ion tracks in condensed matter is proposed. The model is based on three assumptions: the Maxwell distribution of plasma electrons, localization of plasma inside the track nanochannel and constant values of the plasma electron density and temperature during the X-ray irradiation. It is demonstrated that the plasma relaxation model adequately describes the X-ray spectra observed upon interaction of a fast ion with condensed target. Preassumptions of plasma relaxation model are validated by the molecular dynamics modeling and simulation.

Moving base simulation evaluation of translational rate command systems for STOVL aircraft in hover

NASA Technical Reports Server (NTRS)

Franklin, James A.; Stortz, Michael W.

1996-01-01

Using a generalized simulation model, a moving-base simulation of a lift-fan short takeoff/vertical landing fighter aircraft has been conducted on the Vertical Motion Simulator at Ames Research Center. Objectives of the experiment were to determine the influence of system bandwidth and phase delay on flying qualities for translational rate command and vertical velocity command systems. Assessments were made for precision hover control and for landings aboard an LPH type amphibious assault ship in the presence of winds and rough seas. Results obtained define the boundaries between satisfactory and adequate flying qualities for these design features for longitudinal and lateral translational rate command and for vertical velocity command.

Hydrological and water quality processes simulation by the integrated MOHID model

NASA Astrophysics Data System (ADS)

Epelde, Ane; Antiguedad, Iñaki; Brito, David; Eduardo, Jauch; Neves, Ramiro; Sauvage, Sabine; Sánchez-Pérez, José Miguel

2016-04-01

Different modelling approaches have been used in recent decades to study the water quality degradation caused by non-point source pollution. In this study, the MOHID fully distributed and physics-based model has been employed to simulate hydrological processes and nitrogen dynamics in a nitrate vulnerable zone: the Alegria River watershed (Basque Country, Northern Spain). The results of this study indicate that the MOHID code is suitable for hydrological processes simulation at the watershed scale, as the model shows satisfactory performance at simulating the discharge (with NSE: 0.74 and 0.76 during calibration and validation periods, respectively). The agronomical component of the code, allowed the simulation of agricultural practices, which lead to adequate crop yield simulation in the model. Furthermore, the nitrogen exportation also shows satisfactory performance (with NSE: 0.64 and 0.69 during calibration and validation periods, respectively). While the lack of field measurements do not allow to evaluate the nutrient cycling processes in depth, it has been observed that the MOHID model simulates the annual denitrification according to general ranges established for agricultural watersheds (in this study, 9 kg N ha-1 year-1). In addition, the model has simulated coherently the spatial distribution of the denitrification process, which is directly linked to the simulated hydrological conditions. Thus, the model has localized the highest rates nearby the discharge zone of the aquifer and also where the aquifer thickness is low. These results evidence the strength of this model to simulate watershed scale hydrological processes as well as the crop production and the agricultural activity derived water quality degradation (considering both nutrient exportation and nutrient cycling processes).

Large-Eddy Simulation of Waked Turbines in a Scaled Wind Farm Facility

NASA Astrophysics Data System (ADS)

Wang, J.; McLean, D.; Campagnolo, F.; Yu, T.; Bottasso, C. L.

2017-05-01

The aim of this paper is to present the numerical simulation of waked scaled wind turbines operating in a boundary layer wind tunnel. The simulation uses a LES-lifting-line numerical model. An immersed boundary method in conjunction with an adequate wall model is used to represent the effects of both the wind turbine nacelle and tower, which are shown to have a considerable effect on the wake behavior. Multi-airfoil data calibrated at different Reynolds numbers are used to account for the lift and drag characteristics at the low and varying Reynolds conditions encountered in the experiments. The present study focuses on low turbulence inflow conditions and inflow non-uniformity due to wind tunnel characteristics, while higher turbulence conditions are considered in a separate study. The numerical model is validated by using experimental data obtained during test campaigns conducted with the scaled wind farm facility. The simulation and experimental results are compared in terms of power capture, rotor thrust, downstream velocity profiles and turbulence intensity.

A sophisticated simulation for the fracture behavior of concrete material using XFEM

NASA Astrophysics Data System (ADS)

Zhai, Changhai; Wang, Xiaomin; Kong, Jingchang; Li, Shuang; Xie, Lili

2017-10-01

The development of a powerful numerical model to simulate the fracture behavior of concrete material has long been one of the dominant research areas in earthquake engineering. A reliable model should be able to adequately represent the discontinuous characteristics of cracks and simulate various failure behaviors under complicated loading conditions. In this paper, a numerical formulation, which incorporates a sophisticated rigid-plastic interface constitutive model coupling cohesion softening, contact, friction and shear dilatation into the XFEM, is proposed to describe various crack behaviors of concrete material. An effective numerical integration scheme for accurately assembling the contribution to the weak form on both sides of the discontinuity is introduced. The effectiveness of the proposed method has been assessed by simulating several well-known experimental tests. It is concluded that the numerical method can successfully capture the crack paths and accurately predict the fracture behavior of concrete structures. The influence of mode-II parameters on the mixed-mode fracture behavior is further investigated to better determine these parameters.

Recent Progress Towards Predicting Aircraft Ground Handling Performance

NASA Technical Reports Server (NTRS)

Yager, T. J.; White, E. J.

1981-01-01

The significant progress which has been achieved in development of aircraft ground handling simulation capability is reviewed and additional improvements in software modeling identified. The problem associated with providing necessary simulator input data for adequate modeling of aircraft tire/runway friction behavior is discussed and efforts to improve this complex model, and hence simulator fidelity, are described. Aircraft braking performance data obtained on several wet runway surfaces is compared to ground vehicle friction measurements and, by use of empirically derived methods, good agreement between actual and estimated aircraft braking friction from ground vehilce data is shown. The performance of a relatively new friction measuring device, the friction tester, showed great promise in providing data applicable to aircraft friction performance. Additional research efforts to improve methods of predicting tire friction performance are discussed including use of an instrumented tire test vehicle to expand the tire friction data bank and a study of surface texture measurement techniques.

Proceedings of the Augmented VIsual Display (AVID) Research Workshop

NASA Technical Reports Server (NTRS)

Kaiser, Mary K. (Editor); Sweet, Barbara T. (Editor)

1993-01-01

The papers, abstracts, and presentations were presented at a three day workshop focused on sensor modeling and simulation, and image enhancement, processing, and fusion. The technical sessions emphasized how sensor technology can be used to create visual imagery adequate for aircraft control and operations. Participants from industry, government, and academic laboratories contributed to panels on Sensor Systems, Sensor Modeling, Sensor Fusion, Image Processing (Computer and Human Vision), and Image Evaluation and Metrics.

Aurally-adequate time-frequency analysis for scattered sound in auditoria

NASA Astrophysics Data System (ADS)

Norris, Molly K.; Xiang, Ning; Kleiner, Mendel

2005-04-01

The goal of this work was to apply an aurally-adequate time-frequency analysis technique to the analysis of sound scattering effects in auditoria. Time-frequency representations were developed as a motivated effort that takes into account binaural hearing, with a specific implementation of interaural cross-correlation process. A model of the human auditory system was implemented in the MATLAB platform based on two previous models [A. Härmä and K. Palomäki, HUTear, Espoo, Finland; and M. A. Akeroyd, A. Binaural Cross-correlogram Toolbox for MATLAB (2001), University of Sussex, Brighton]. These stages include proper frequency selectivity, the conversion of the mechanical motion of the basilar membrane to neural impulses, and binaural hearing effects. The model was then used in the analysis of room impulse responses with varying scattering characteristics. This paper discusses the analysis results using simulated and measured room impulse responses. [Work supported by the Frank H. and Eva B. Buck Foundation.

About one counterexample of applying method of splitting in modeling of plating processes

NASA Astrophysics Data System (ADS)

Solovjev, D. S.; Solovjeva, I. A.; Litovka, Yu V.; Korobova, I. L.

2018-05-01

The paper presents the main factors that affect the uniformity of the thickness distribution of plating on the surface of the product. The experimental search for the optimal values of these factors is expensive and time-consuming. The problem of adequate simulation of coating processes is very relevant. The finite-difference approximation using seven-point and five-point templates in combination with the splitting method is considered as solution methods for the equations of the model. To study the correctness of the solution of equations of the mathematical model by these methods, the experiments were conducted on plating with a flat anode and cathode, which relative position was not changed in the bath. The studies have shown that the solution using the splitting method was up to 1.5 times faster, but it did not give adequate results due to the geometric features of the task under the given boundary conditions.

Different modelling approaches to evaluate nitrogen transport and turnover at the watershed scale

NASA Astrophysics Data System (ADS)

Epelde, Ane Miren; Antiguedad, Iñaki; Brito, David; Jauch, Eduardo; Neves, Ramiro; Garneau, Cyril; Sauvage, Sabine; Sánchez-Pérez, José Miguel

2016-08-01

This study presents the simulation of hydrological processes and nutrient transport and turnover processes using two integrated numerical models: Soil and Water Assessment Tool (SWAT) (Arnold et al., 1998), an empirical and semi-distributed numerical model; and Modelo Hidrodinâmico (MOHID) (Neves, 1985), a physics-based and fully distributed numerical model. This work shows that both models reproduce satisfactorily water and nitrate exportation at the watershed scale at annual and daily basis, MOHID providing slightly better results. At the watershed scale, both SWAT and MOHID simulated similarly and satisfactorily the denitrification amount. However, as MOHID numerical model was the only one able to reproduce adequately the spatial variation of the soil hydrological conditions and water table level fluctuation, it proved to be the only model able of reproducing the spatial variation of the nutrient cycling processes that are dependent to the soil hydrological conditions such as the denitrification process. This evidences the strength of the fully distributed and physics-based models to simulate the spatial variability of nutrient cycling processes that are dependent to the hydrological conditions of the soils.

Shape-based approach for the estimation of individual facial mimics in craniofacial surgery planning

NASA Astrophysics Data System (ADS)

Gladilin, Evgeny; Zachow, Stefan; Deuflhard, Peter; Hege, Hans-Christian

2002-05-01

Besides the static soft tissue prediction, the estimation of basic facial emotion expressions is another important criterion for the evaluation of craniofacial surgery planning. For a realistic simulation of facial mimics, an adequate biomechanical model of soft tissue including the mimic musculature is needed. In this work, we present an approach for the modeling of arbitrarily shaped muscles and the estimation of basic individual facial mimics, which is based on the geometrical model derived from the individual tomographic data and the general finite element modeling of soft tissue biomechanics.

Design of Bi-Directional Hydrofoils for Tidal Current Turbines

NASA Astrophysics Data System (ADS)

Nedyalkov, Ivaylo; Wosnik, Martin

2015-11-01

Tidal Current Turbines operate in flows which reverse direction. Bi-directional hydrofoils have rotational symmetry and allow such turbines to operate without the need for pitch or yaw control, decreasing the initial and maintenance costs. A numerical test-bed was developed to automate the simulations of hydrofoils in OpenFOAM and was utilized to simulate the flow over eleven classes of hydrofoils comprising a total of 700 foil shapes at different angles of attack. For promising candidate foil shapes physical models of 75 mm chord and 150 mm span were fabricated and tested in the University of New Hampshire High-Speed Cavitation Tunnel (HiCaT). The experimental results were compared to the simulations for model validation. The numerical test-bed successfully generated simulations for a wide range of foil shapes, although, as expected, the k - ω - SST turbulence model employed here was not adequate for some of the foils and for large angles of attack at which separation occurred. An optimization algorithm is currently being coupled with the numerical test-bed and additional turbulence models will be implemented in the future.

SPH-based numerical simulations of flow slides in municipal solid waste landfills.

PubMed

Huang, Yu; Dai, Zili; Zhang, Weijie; Huang, Maosong

2013-03-01

Most municipal solid waste (MSW) is disposed of in landfills. Over the past few decades, catastrophic flow slides have occurred in MSW landfills around the world, causing substantial economic damage and occasionally resulting in human victims. It is therefore important to predict the run-out, velocity and depth of such slides in order to provide adequate mitigation and protection measures. To overcome the limitations of traditional numerical methods for modelling flow slides, a mesh-free particle method entitled smoothed particle hydrodynamics (SPH) is introduced in this paper. The Navier-Stokes equations were adopted as the governing equations and a Bingham model was adopted to analyse the relationship between material stress rates and particle motion velocity. The accuracy of the model is assessed using a series of verifications, and then flow slides that occurred in landfills located in Sarajevo and Bandung were simulated to extend its applications. The simulated results match the field data well and highlight the capability of the proposed SPH modelling method to simulate such complex phenomena as flow slides in MSW landfills.

Opportunities and Challenges in Supply-Side Simulation: Physician-Based Models

PubMed Central

Gresenz, Carole Roan; Auerbach, David I; Duarte, Fabian

2013-01-01

Objective To provide a conceptual framework and to assess the availability of empirical data for supply-side microsimulation modeling in the context of health care. Data Sources Multiple secondary data sources, including the American Community Survey, Health Tracking Physician Survey, and SK&A physician database. Study Design We apply our conceptual framework to one entity in the health care market—physicians—and identify, assess, and compare data available for physician-based simulation models. Principal Findings Our conceptual framework describes three broad types of data required for supply-side microsimulation modeling. Our assessment of available data for modeling physician behavior suggests broad comparability across various sources on several dimensions and highlights the need for significant integration of data across multiple sources to provide a platform adequate for modeling. A growing literature provides potential estimates for use as behavioral parameters that could serve as the models' engines. Sources of data for simulation modeling that account for the complex organizational and financial relationships among physicians and other supply-side entities are limited. Conclusions A key challenge for supply-side microsimulation modeling is optimally combining available data to harness their collective power. Several possibilities also exist for novel data collection. These have the potential to serve as catalysts for the next generation of supply-side-focused simulation models to inform health policy. PMID:23347041

Early Validation of Failure Detection, Isolation, and Recovery Design Using Heterogeneous Modelling and Simulation

NASA Astrophysics Data System (ADS)

van der Plas, Peter; Guerriero, Suzanne; Cristiano, Leorato; Rugina, Ana

2012-08-01

Modelling and simulation can support a number of use cases across the spacecraft development life-cycle. Given the increasing complexity of space missions, the observed general trend is for a more extensive usage of simulation already in the early phases. A major perceived advantage is that modelling and simulation can enable the validation of critical aspects of the spacecraft design before the actual development is started, as such reducing the risk in later phases.Failure Detection, Isolation, and Recovery (FDIR) is one of the areas with a high potential to benefit from early modelling and simulation. With the increasing level of required spacecraft autonomy, FDIR specifications can grow in such a way that the traditional document-based review process soon becomes inadequate.This paper shows that FDIR modelling and simulation in a system context can provide a powerful tool to support the FDIR verification process. It is highlighted that FDIR modelling at this early stage requires heterogeneous modelling tools and languages, in order to provide an adequate functional description of the different components (i.e. FDIR functions, environment, equipment, etc.) to be modelled.For this reason, an FDIR simulation framework is proposed in this paper. This framework is based on a number of tools already available in the Avionics Systems Laboratory at ESTEC, which are the Avionics Test Bench Functional Engineering Simulator (ATB FES), Matlab/Simulink, TASTE, and Real Time Developer Studio (RTDS).The paper then discusses the application of the proposed simulation framework to a real case-study, i.e. the FDIR modelling of a satellite in support of actual ESA mission. Challenges and benefits of the approach are described. Finally, lessons learned and the generality of the proposed approach are discussed.

Using High-Resolution Satellite Observations for Evaluation of Cloud and Precipitation Statistics from Cloud-Resolving Model Simulations. Part I: South China Sea Monsoon Experiment

NASA Astrophysics Data System (ADS)

Zhou, Y.; Hou, A.; Lau, W. K.; Shie, C.; Tao, W.; Lin, X.; Chou, M.; Olson, W. S.; Grecu, M.

2006-05-01

The cloud and precipitation statistics simulated by 3D Goddard Cumulus Ensemble (GCE) model during the South China Sea Monsoon Experiment (SCSMEX) is compared with Tropical Rainfall Measuring Mission (TRMM) TMI and PR rainfall measurements and the Earth's Radiant Energy System (CERES) single scanner footprint (SSF) radiation and cloud retrievals. It is found that GCE is capable of simulating major convective system development and reproducing total surface rainfall amount as compared with rainfall estimated from the soundings. Mesoscale organization is adequately simulated except when environmental wind shear is very weak. The partitions between convective and stratiform rain are also close to TMI and PR classification. However, the model simulated rain spectrum is quite different from either TMI or PR measurements. The model produces more heavy rains and light rains (less than 0.1 mm/hr) than the observations. The model also produces heavier vertical hydrometer profiles of rain, graupel when compared with TMI retrievals and PR radar reflectivity. Comparing GCE simulated OLR and cloud properties with CERES measurements found that the model has much larger domain averaged OLR due to smaller total cloud fraction and a much skewed distribution of OLR and cloud top than CERES observations, indicating that the model's cloud field is not wide spread, consistent with the model's precipitation activity. These results will be used as guidance for improving the model's microphysics.

A Lagrangian stochastic model for aerial spray transport above an oak forest

USGS Publications Warehouse

Wang, Yansen; Miller, David R.; Anderson, Dean E.; McManus, Michael L.

1995-01-01

An aerial spray droplets' transport model has been developed by applying recent advances in Lagrangian stochastic simulation of heavy particles. A two-dimensional Lagrangian stochastic model was adopted to simulate the spray droplet dispersion in atmospheric turbulence by adjusting the Lagrangian integral time scale along the drop trajectory. The other major physical processes affecting the transport of spray droplets above a forest canopy, the aircraft wingtip vortices and the droplet evaporation, were also included in each time step of the droplets' transport.The model was evaluated using data from an aerial spray field experiment. In generally neutral stability conditions, the accuracy of the model predictions varied from run-to-run as expected. The average root-mean-square error was 24.61 IU cm−2, and the average relative error was 15%. The model prediction was adequate in two-dimensional steady wind conditions, but was less accurate in variable wind condition. The results indicated that the model can simulate successfully the ensemble; average transport of aerial spray droplets under neutral, steady atmospheric wind conditions.

Using High Resolution Simulations with WRF/SSiB Regional Climate Model Constrained by In Situ Observations to Assess the Impacts of Dust in Snow in the Upper Colorado River Basin

NASA Astrophysics Data System (ADS)

Oaida, C. M.; Skiles, M.; Painter, T. H.; Xue, Y.

2015-12-01

The mountain snowpack is an essential resource for both the environment as well as society. Observational and energy balance modeling work have shown that dust on snow (DOS) in western U.S. (WUS) is a major contributor to snow processes, including snowmelt timing and runoff amount in regions like the Upper Colorado River Basin (UCRB). In order to accurately estimate the impact of DOS to the hydrologic cycle and water resources, now and under a changing climate, we need to be able to (1) adequately simulate the snowpack (accumulation), and (2) realistically represent DOS processes in models. Energy balance models do not capture the impact on a broader local or regional scale, nor the land-atmosphere feedbacks, while GCM studies cannot resolve orographic-related precipitation processes, and therefore snowpack accumulation, owing to coarse spatial resolution and smoother terrain. All this implies the impacts of dust on snow on the mountain snowpack and other hydrologic processes are likely not well captured in current modeling studies. Recent increase in computing power allows for RCMs to be used at higher spatial resolutions, while recent in situ observations of dust in snow properties can help constrain modeling simulations. Therefore, in the work presented here, we take advantage of these latest resources to address the some of the challenges outlined above. We employ the newly enhanced WRF/SSiB regional climate model at 4 km horizontal resolution. This scale has been shown by others to be adequate in capturing orographic processes over WUS. We also constrain the magnitude of dust deposition provided by a global chemistry and transport model, with in situ measurements taken at sites in the UCRB. Furthermore, we adjust the dust absorptive properties based on observed values at these sites, as opposed to generic global ones. This study aims to improve simulation of the impact of dust in snow on the hydrologic cycle and related water resources.

Effects of magnetospheric lobe cell convection on dayside upper thermospheric winds at high latitudes

NASA Astrophysics Data System (ADS)

Zhang, B.; Wang, W.; Wu, Q.; Knipp, D.; Kilcommons, L.; Brambles, O. J.; Liu, J.; Wiltberger, M.; Lyon, J. G.; Häggström, I.

2016-08-01

This paper investigates a possible physical mechanism of the observed dayside high-latitude upper thermospheric wind using numerical simulations from the coupled magnetosphere-ionosphere-thermosphere (CMIT) model. Results show that the CMIT model is capable of reproducing the unexpected afternoon equatorward winds in the upper thermosphere observed by the High altitude Interferometer WIND observation (HIWIND) balloon. Models that lack adequate coupling produce poleward winds. The modeling study suggests that ion drag driven by magnetospheric lobe cell convection is another possible mechanism for turning the climatologically expected dayside poleward winds to the observed equatorward direction. The simulation results are validated by HIWIND, European Incoherent Scatter, and Defense Meteorological Satellite Program. The results suggest a strong momentum coupling between high-latitude ionospheric plasma circulation and thermospheric neutral winds in the summer hemisphere during positive IMF Bz periods, through the formation of magnetospheric lobe cell convection driven by persistent positive IMF By. The CMIT simulation adds important insight into the role of dayside coupling during intervals of otherwise quiet geomagnetic activity

Vibronic coupling simulations for linear and nonlinear optical processes: Simulation results

NASA Astrophysics Data System (ADS)

Silverstein, Daniel W.; Jensen, Lasse

2012-02-01

A vibronic coupling model based on time-dependent wavepacket approach is applied to simulate linear optical processes, such as one-photon absorbance and resonance Raman scattering, and nonlinear optical processes, such as two-photon absorbance and resonance hyper-Raman scattering, on a series of small molecules. Simulations employing both the long-range corrected approach in density functional theory and coupled cluster are compared and also examined based on available experimental data. Although many of the small molecules are prone to anharmonicity in their potential energy surfaces, the harmonic approach performs adequately. A detailed discussion of the non-Condon effects is illustrated by the molecules presented in this work. Linear and nonlinear Raman scattering simulations allow for the quantification of interference between the Franck-Condon and Herzberg-Teller terms for different molecules.

Multislice spiral CT simulator for dynamic cardiopulmonary studies

NASA Astrophysics Data System (ADS)

De Francesco, Silvia; Ferreira da Silva, Augusto M.

2002-04-01

We've developed a Multi-slice Spiral CT Simulator modeling the acquisition process of a real tomograph over a 4-dimensional phantom (4D MCAT) of the human thorax. The simulator allows us to visually characterize artifacts due to insufficient temporal sampling and a priori evaluate the quality of the images obtained in cardio-pulmonary studies (both with single-/multi-slice and ECG gated acquisition processes). The simulating environment allows both for conventional and spiral scanning modes and includes a model of noise in the acquisition process. In case of spiral scanning, reconstruction facilities include longitudinal interpolation methods (360LI and 180LI both for single and multi-slice). Then, the reconstruction of the section is performed through FBP. The reconstructed images/volumes are affected by distortion due to insufficient temporal sampling of the moving object. The developed simulating environment allows us to investigate the nature of the distortion characterizing it qualitatively and quantitatively (using, for example, Herman's measures). Much of our work is focused on the determination of adequate temporal sampling and sinogram regularization techniques. At the moment, the simulator model is limited to the case of multi-slice tomograph, being planned as a next step of development the extension to cone beam or area detectors.

Simulation tools for particle-based reaction-diffusion dynamics in continuous space

PubMed Central

2014-01-01

Particle-based reaction-diffusion algorithms facilitate the modeling of the diffusional motion of individual molecules and the reactions between them in cellular environments. A physically realistic model, depending on the system at hand and the questions asked, would require different levels of modeling detail such as particle diffusion, geometrical confinement, particle volume exclusion or particle-particle interaction potentials. Higher levels of detail usually correspond to increased number of parameters and higher computational cost. Certain systems however, require these investments to be modeled adequately. Here we present a review on the current field of particle-based reaction-diffusion software packages operating on continuous space. Four nested levels of modeling detail are identified that capture incrementing amount of detail. Their applicability to different biological questions is discussed, arching from straight diffusion simulations to sophisticated and expensive models that bridge towards coarse grained molecular dynamics. PMID:25737778

A summary of computational experience at GE Aircraft Engines for complex turbulent flows in gas turbines

NASA Astrophysics Data System (ADS)

Zerkle, Ronald D.; Prakash, Chander

1995-03-01

This viewgraph presentation summarizes some CFD experience at GE Aircraft Engines for flows in the primary gaspath of a gas turbine engine and in turbine blade cooling passages. It is concluded that application of the standard k-epsilon turbulence model with wall functions is not adequate for accurate CFD simulation of aerodynamic performance and heat transfer in the primary gas path of a gas turbine engine. New models are required in the near-wall region which include more physics than wall functions. The two-layer modeling approach appears attractive because of its computational complexity. In addition, improved CFD simulation of film cooling and turbine blade internal cooling passages will require anisotropic turbulence models. New turbulence models must be practical in order to have a significant impact on the engine design process. A coordinated turbulence modeling effort between NASA centers would be beneficial to the gas turbine industry.

A summary of computational experience at GE Aircraft Engines for complex turbulent flows in gas turbines

NASA Technical Reports Server (NTRS)

Zerkle, Ronald D.; Prakash, Chander

1995-01-01

This viewgraph presentation summarizes some CFD experience at GE Aircraft Engines for flows in the primary gaspath of a gas turbine engine and in turbine blade cooling passages. It is concluded that application of the standard k-epsilon turbulence model with wall functions is not adequate for accurate CFD simulation of aerodynamic performance and heat transfer in the primary gas path of a gas turbine engine. New models are required in the near-wall region which include more physics than wall functions. The two-layer modeling approach appears attractive because of its computational complexity. In addition, improved CFD simulation of film cooling and turbine blade internal cooling passages will require anisotropic turbulence models. New turbulence models must be practical in order to have a significant impact on the engine design process. A coordinated turbulence modeling effort between NASA centers would be beneficial to the gas turbine industry.

Opportunities and pitfalls in clinical proof-of-concept: principles and examples.

PubMed

Chen, Chao

2018-04-01

Clinical proof-of-concept trials crucially inform major resource deployment decisions. This paper discusses several mechanisms for enhancing their rigour and efficiency. The importance of careful consideration when using a surrogate endpoint is illustrated; situational effectiveness of run-in patient enrichment is explored; a versatile tool is introduced to ensure a strong pharmacological underpinning; the benefits of dose-titration are revealed by simulation; and the importance of adequately scheduled observations is shown. The general process of model-based trial design and analysis is described and several examples demonstrate the value in historical data, simulation-guided design, model-based analysis and trial adaptation informed by interim analysis. Copyright © 2018 Elsevier Ltd. All rights reserved.

Relaxing the rule of ten events per variable in logistic and Cox regression.

PubMed

Vittinghoff, Eric; McCulloch, Charles E

2007-03-15

The rule of thumb that logistic and Cox models should be used with a minimum of 10 outcome events per predictor variable (EPV), based on two simulation studies, may be too conservative. The authors conducted a large simulation study of other influences on confidence interval coverage, type I error, relative bias, and other model performance measures. They found a range of circumstances in which coverage and bias were within acceptable levels despite less than 10 EPV, as well as other factors that were as influential as or more influential than EPV. They conclude that this rule can be relaxed, in particular for sensitivity analyses undertaken to demonstrate adequate control of confounding.

Effect of Turbulence Modeling on an Excited Jet

NASA Technical Reports Server (NTRS)

Brown, Clifford A.; Hixon, Ray

2010-01-01

The flow dynamics in a high-speed jet are dominated by unsteady turbulent flow structures in the plume. Jet excitation seeks to control these flow structures through the natural instabilities present in the initial shear layer of the jet. Understanding and optimizing the excitation input, for jet noise reduction or plume mixing enhancement, requires many trials that may be done experimentally or computationally at a significant cost savings. Numerical simulations, which model various parts of the unsteady dynamics to reduce the computational expense of the simulation, must adequately capture the unsteady flow dynamics in the excited jet for the results are to be used. Four CFD methods are considered for use in an excited jet problem, including two turbulence models with an Unsteady Reynolds Averaged Navier-Stokes (URANS) solver, one Large Eddy Simulation (LES) solver, and one URANS/LES hybrid method. Each method is used to simulate a simplified excited jet and the results are evaluated based on the flow data, computation time, and numerical stability. The knowledge gained about the effect of turbulence modeling and CFD methods from these basic simulations will guide and assist future three-dimensional (3-D) simulations that will be used to understand and optimize a realistic excited jet for a particular application.

Numerical simulation of double‐diffusive finger convection

USGS Publications Warehouse

Hughes, Joseph D.; Sanford, Ward E.; Vacher, H. Leonard

2005-01-01

A hybrid finite element, integrated finite difference numerical model is developed for the simulation of double‐diffusive and multicomponent flow in two and three dimensions. The model is based on a multidimensional, density‐dependent, saturated‐unsaturated transport model (SUTRA), which uses one governing equation for fluid flow and another for solute transport. The solute‐transport equation is applied sequentially to each simulated species. Density coupling of the flow and solute‐transport equations is accounted for and handled using a sequential implicit Picard iterative scheme. High‐resolution data from a double‐diffusive Hele‐Shaw experiment, initially in a density‐stable configuration, is used to verify the numerical model. The temporal and spatial evolution of simulated double‐diffusive convection is in good agreement with experimental results. Numerical results are very sensitive to discretization and correspond closest to experimental results when element sizes adequately define the spatial resolution of observed fingering. Numerical results also indicate that differences in the molecular diffusivity of sodium chloride and the dye used to visualize experimental sodium chloride concentrations are significant and cause inaccurate mapping of sodium chloride concentrations by the dye, especially at late times. As a result of reduced diffusion, simulated dye fingers are better defined than simulated sodium chloride fingers and exhibit more vertical mass transfer.

Aspects of intelligent electronic device based switchgear control training model application

NASA Astrophysics Data System (ADS)

Bogdanov, Dimitar; Popov, Ivaylo

2018-02-01

The design of the protection and control equipment for electrical power sector application was object of extensive advance in the last several decades. The modern technologies offer a wide range of multifunctional flexible applications, making the protection and control of facilities more sophisticated. In the same time, the advance of technology imposes the necessity of simulators, training models and tutorial laboratory equipment to be used for adequate training of students and field specialists

Computer simulation of surface and film processes

NASA Technical Reports Server (NTRS)

Tiller, W. A.; Halicioglu, M. T.

1983-01-01

Adequate computer methods, based on interactions between discrete particles, provide information leading to an atomic level understanding of various physical processes. The success of these simulation methods, however, is related to the accuracy of the potential energy function representing the interactions among the particles. The development of a potential energy function for crystalline SiO2 forms that can be employed in lengthy computer modelling procedures was investigated. In many of the simulation methods which deal with discrete particles, semiempirical two body potentials were employed to analyze energy and structure related properties of the system. Many body interactions are required for a proper representation of the total energy for many systems. Many body interactions for simulations based on discrete particles are discussed.

Bayesian Assessment of the Uncertainties of Estimates of a Conceptual Rainfall-Runoff Model Parameters

NASA Astrophysics Data System (ADS)

Silva, F. E. O. E.; Naghettini, M. D. C.; Fernandes, W.

2014-12-01

This paper evaluated the uncertainties associated with the estimation of the parameters of a conceptual rainfall-runoff model, through the use of Bayesian inference techniques by Monte Carlo simulation. The Pará River sub-basin, located in the upper São Francisco river basin, in southeastern Brazil, was selected for developing the studies. In this paper, we used the Rio Grande conceptual hydrologic model (EHR/UFMG, 2001) and the Markov Chain Monte Carlo simulation method named DREAM (VRUGT, 2008a). Two probabilistic models for the residues were analyzed: (i) the classic [Normal likelihood - r ≈ N (0, σ²)]; and (ii) a generalized likelihood (SCHOUPS & VRUGT, 2010), in which it is assumed that the differences between observed and simulated flows are correlated, non-stationary, and distributed as a Skew Exponential Power density. The assumptions made for both models were checked to ensure that the estimation of uncertainties in the parameters was not biased. The results showed that the Bayesian approach proved to be adequate to the proposed objectives, enabling and reinforcing the importance of assessing the uncertainties associated with hydrological modeling.

Numerical simulations of an advection fog event over Shanghai Pudong International Airport with the WRF model

NASA Astrophysics Data System (ADS)

Lin, Caiyan; Zhang, Zhongfeng; Pu, Zhaoxia; Wang, Fengyun

2017-10-01

A series of numerical simulations is conducted to understand the formation, evolution, and dissipation of an advection fog event over Shanghai Pudong International Airport (ZSPD) with the Weather Research and Forecasting (WRF) model. Using the current operational settings at the Meteorological Center of East China Air Traffic Management Bureau, the WRF model successfully predicts the fog event at ZSPD. Additional numerical experiments are performed to examine the physical processes associated with the fog event. The results indicate that prediction of this particular fog event is sensitive to microphysical schemes for the time of fog dissipation but not for the time of fog onset. The simulated timing of the arrival and dissipation of the fog, as well as the cloud distribution, is substantially sensitive to the planetary boundary layer and radiation (both longwave and shortwave) processes. Moreover, varying forecast lead times also produces different simulation results for the fog event regarding its onset and duration, suggesting a trade-off between more accurate initial conditions and a proper forecast lead time that allows model physical processes to spin up adequately during the fog simulation. The overall outcomes from this study imply that the complexity of physical processes and their interactions within the WRF model during fog evolution and dissipation is a key area of future research.

Modeling scintillator and WLS fiber signals for fast Monte Carlo simulations

NASA Astrophysics Data System (ADS)

Sánchez, F. A.; Medina-Tanco, G.

2010-08-01

In this work we present a fast, robust and flexible procedure to simulate electronic signals of scintillator units: plastic scintillator material embedded with a wavelength shifter optical fiber coupled to a photo-multiplier tube which, in turn, is plugged to a front-end electronic board. The simple rationale behind the simulation chain allows to adapt the procedure to a broad range of detectors based on that kind of units. We show that, in order to produce realistic results, the simulation parameters can be properly calibrated against laboratory measurements and used thereafter as input of the simulations. Simulated signals of atmospheric background cosmic ray muons are presented and their main features analyzed and validated using actual measured data. Conversely, for any given practical application, the present simulation scheme can be used to find an adequate combination of photo-multiplier tube and optical fiber at the prototyping stage.

Optimization of droplets for UV-NIL using coarse-grain simulation of resist flow

NASA Astrophysics Data System (ADS)

Sirotkin, Vadim; Svintsov, Alexander; Zaitsev, Sergey

2009-03-01

A mathematical model and numerical method are described, which make it possible to simulate ultraviolet ("step and flash") nanoimprint lithography (UV-NIL) process adequately even using standard Personal Computers. The model is derived from 3D Navier-Stokes equations with the understanding that the resist motion is largely directed along the substrate surface and characterized by ultra-low values of the Reynolds number. By the numerical approximation of the model, a special finite difference method is applied (a coarse-grain method). A coarse-grain modeling tool for detailed analysis of resist spreading in UV-NIL at the structure-scale level is tested. The obtained results demonstrate the high ability of the tool to calculate optimal dispensing for given stamp design and process parameters. This dispensing provides uniform filled areas and a homogeneous residual layer thickness in UV-NIL.

Snowmelt-runoff Model Utilizing Remotely-sensed Data

NASA Technical Reports Server (NTRS)

Rango, A.

1985-01-01

Remotely sensed snow cover information is the critical data input for the Snowmelt-Runoff Model (SRM), which was developed to simulatke discharge from mountain basins where snowmelt is an important component of runoff. Of simple structure, the model requires only input of temperature, precipitation, and snow covered area. SRM was run successfully on two widely separated basins. The simulations on the Kings River basin are significant because of the large basin area (4000 sq km) and the adequate performance in the most extreme drought year of record (1976). The performance of SRM on the Okutadami River basin was important because it was accomplished with minimum snow cover data available. Tables show: optimum and minimum conditions for model application; basin sizes and elevations where SRM was applied; and SRM strengths and weaknesses. Graphs show results of discharge simulation.

Large Eddy Simulation of Turbulent Combustion

DTIC Science & Technology

2005-10-01

a new method to automatically generate skeletal kinetic mechanisms for surrogate fuels, using the directed relation graph method with error...propagation, was developed. These mechanisms are guaranteed to match results obtained using detailed chemistry within a user- defined accuracy for any...specified target. They can be combined together to produce adequate chemical models for surrogate fuels. A library containing skeletal mechanisms of various

A simulation model for the infiltration of heterogeneous sediment into a stream bed

Treesearch

Tim Lauck; Roland Lamberson; Thomas E. Lisle

1993-01-01

Abstract - Salmonid embryos depend on the adequate flow of oxygenated water to survive and interstitial passageways to emerge from the gravel bed. Spawning gravels are initially cleaned by the spawning female, but sediment transported during subsequent high-runoff events can nfiltrate the porous substrate. In many gravel-bed channels used for spawning, most of the...

Revision of the Rawls et al. (1982) pedotransfer functions for their applicability to US croplands

USDA-ARS?s Scientific Manuscript database

Large scale environmental impact studies typically involve the use of simulation models and require a variety of inputs, some of which may need to be estimated in absence of adequate measured data. As an example, soil water retention needs to be estimated for a large number of soils that are to be u...

Population pharmacokinetic-pharmacodynamic modeling and model-based prediction of docetaxel-induced neutropenia in Japanese patients with non-small cell lung cancer.

PubMed

Fukae, Masato; Shiraishi, Yoshimasa; Hirota, Takeshi; Sasaki, Yuka; Yamahashi, Mika; Takayama, Koichi; Nakanishi, Yoichi; Ieiri, Ichiro

2016-11-01

Docetaxel is used to treat many cancers, and neutropenia is the dose-limiting factor for its clinical use. A population pharmacokinetic-pharmacodynamic (PK-PD) model was introduced to predict the development of docetaxel-induced neutropenia in Japanese patients with non-small cell lung cancer (NSCLC). Forty-seven advanced or recurrent Japanese patients with NSCLC were enrolled. Patients received 50 or 60 mg/m 2 docetaxel as monotherapy, and blood samples for a PK analysis were collected up to 24 h after its infusion. Laboratory tests including absolute neutrophil count data and demographic information were used in population PK-PD modeling. The model was built by NONMEM 7.2 with a first-order conditional estimation using an interaction method. Based on the final model, a Monte Carlo simulation was performed to assess the impact of covariates on and the predictability of neutropenia. A three-compartment model was employed to describe PK data, and the PK model adequately described the docetaxel concentrations observed. Serum albumin (ALB) was detected as a covariate of clearance (CL): CL (L/h) = 32.5 × (ALB/3.6) 0.965  × (WGHT/70) 3/4 . In population PK-PD modeling, a modified semi-mechanistic myelosuppression model was applied, and characterization of the time course of neutrophil counts was adequate. The covariate selection indicated that α1-acid glycoprotein (AAG) was a predictor of neutropenia. The model-based simulation also showed that ALB and AAG negatively correlated with the development of neutropenia and that the time course of neutrophil counts was predictable. The developed model may facilitate the prediction and care of docetaxel-induced neutropenia.

Short-stack modeling of degradation in solid oxide fuel cells. Part I. Contact degradation

NASA Astrophysics Data System (ADS)

Gazzarri, J. I.; Kesler, O.

As the first part of a two paper series, we present a two-dimensional impedance model of a working solid oxide fuel cell (SOFC) to study the effect of contact degradation on the impedance spectrum for the purpose of non-invasive diagnosis. The two dimensional modeled geometry includes the ribbed interconnect, and is adequate to represent co- and counter-flow configurations. Simulated degradation modes include: cathode delamination, interconnect oxidation, and interconnect-cathode detachment. The simulations show differences in the way each degradation mode impacts the impedance spectrum shape, suggesting that identification is possible. In Part II, we present a sensitivity analysis of the results to input parameter variability that reveals strengths and limitations of the method, as well as describing possible interactions between input parameters and concurrent degradation modes.

Electron-phonon interaction within classical molecular dynamics

DOE PAGES

Tamm, A.; Samolyuk, G.; Correa, A. A.; ...

2016-07-14

Here, we present a model for nonadiabatic classical molecular dynamics simulations that captures with high accuracy the wave-vector q dependence of the phonon lifetimes, in agreement with quantum mechanics calculations. It is based on a local view of the e-ph interaction where individual atom dynamics couples to electrons via a damping term that is obtained as the low-velocity limit of the stopping power of a moving ion in a host. The model is parameter free, as its components are derived from ab initio-type calculations, is readily extended to the case of alloys, and is adequate for large-scale molecular dynamics computermore » simulations. We also show how this model removes some oversimplifications of the traditional ionic damped dynamics commonly used to describe situations beyond the Born-Oppenheimer approximation.« less

Hand-assisted laparoscopic sigmoid colectomy skills acquisition: augmented reality simulator versus human cadaver training models.

PubMed

Leblanc, Fabien; Senagore, Anthony J; Ellis, Clyde N; Champagne, Bradley J; Augestad, Knut M; Neary, Paul C; Delaney, Conor P

2010-01-01

The aim of this study was to compare a simulator with the human cadaver model for hand-assisted laparoscopic colorectal skills acquisition training. An observational prospective comparative study was conducted to compare the laparoscopic surgery training models. The study took place during the laparoscopic colectomy training course performed at the annual scientific meeting of the American Society of Colon and Rectal Surgeons. Thirty four practicing surgeons performed hand-assisted laparoscopic sigmoid colectomy on human cadavers (n = 7) and on an augmented reality simulator (n = 27). Prior laparoscopic colorectal experience was assessed. Trainers and trainees completed independently objective structured assessment forms. Training models were compared by trainees' technical skills scores, events scores, and satisfaction. Prior laparoscopic experience was similar in both surgeon groups. Generic and specific skills scores were similar on both training models. Generic events scores were significantly better on the cadaver model. The 2 most frequent generic events occurring on the simulator were poor hand-eye coordination and inefficient use of retraction. Specific events were scored better on the simulator and reached the significance limit (p = 0.051) for trainers. The specific events occurring on the cadaver were intestinal perforation and left ureter identification difficulties. Overall satisfaction was better for the cadaver than for the simulator model (p = 0.009). With regard to skills scores, the augmented reality simulator had adequate qualities for the hand-assisted laparoscopic colectomy training. Nevertheless, events scores highlighted weaknesses of the anatomical replication on the simulator. Although improvements likely will be required to incorporate the simulator more routinely into the colorectal training, it may be useful in its current form for more junior trainees or those early on their learning curve. Copyright 2010 Association of Program Directors in Surgery. Published by Elsevier Inc. All rights reserved.

Streamflow Simulations and Percolation Estimates Using the Soil and Water Assessment Tool for Selected Basins in North-Central Nebraska, 1940-2005

USGS Publications Warehouse

Strauch, Kellan R.; Linard, Joshua I.

2009-01-01

The U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, used the Soil and Water Assessment Tool to simulate streamflow and estimate percolation in north-central Nebraska to aid development of long-term strategies for management of hydrologically connected ground and surface water. Although groundwater models adequately simulate subsurface hydrologic processes, they often are not designed to simulate the hydrologically complex processes occurring at or near the land surface. The use of watershed models such as the Soil and Water Assessment Tool, which are designed specifically to simulate surface and near-subsurface processes, can provide helpful insight into the effects of surface-water hydrology on the groundwater system. The Soil and Water Assessment Tool was calibrated for five stream basins in the Elkhorn-Loup Groundwater Model study area in north-central Nebraska to obtain spatially variable estimates of percolation. Six watershed models were calibrated to recorded streamflow in each subbasin by modifying the adjustment parameters. The calibrated parameter sets were then used to simulate a validation period; the validation period was half of the total streamflow period of record with a minimum requirement of 10 years. If the statistical and water-balance results for the validation period were similar to those for the calibration period, a model was considered satisfactory. Statistical measures of each watershed model's performance were variable. These objective measures included the Nash-Sutcliffe measure of efficiency, the ratio of the root-mean-square error to the standard deviation of the measured data, and an estimate of bias. The model met performance criteria for the bias statistic, but failed to meet statistical adequacy criteria for the other two performance measures when evaluated at a monthly time step. A primary cause of the poor model validation results was the inability of the model to reproduce the sustained base flow and streamflow response to precipitation that was observed in the Sand Hills region. The watershed models also were evaluated based on how well they conformed to the annual mass balance (precipitation equals the sum of evapotranspiration, streamflow/runoff, and deep percolation). The model was able to adequately simulate annual values of evapotranspiration, runoff, and precipitation in comparison to reported values, which indicates the model may provide reasonable estimates of annual percolation. Mean annual percolation estimated by the model as basin averages varied within the study area from a maximum of 12.9 inches in the Loup River Basin to a minimum of 1.5 inches in the Shell Creek Basin. Percolation also varied within the studied basins; basin headwaters tended to have greater percolation rates than downstream areas. This variance in percolation rates was mainly was because of the predominance of sandy, highly permeable soils in the upstream areas of the modeled basins.

Parameterizing deep convection using the assumed probability density function method

DOE PAGES

Storer, R. L.; Griffin, B. M.; Höft, J.; ...

2014-06-11

Due to their coarse horizontal resolution, present-day climate models must parameterize deep convection. This paper presents single-column simulations of deep convection using a probability density function (PDF) parameterization. The PDF parameterization predicts the PDF of subgrid variability of turbulence, clouds, and hydrometeors. That variability is interfaced to a prognostic microphysics scheme using a Monte Carlo sampling method. The PDF parameterization is used to simulate tropical deep convection, the transition from shallow to deep convection over land, and mid-latitude deep convection. These parameterized single-column simulations are compared with 3-D reference simulations. The agreement is satisfactory except when the convective forcing ismore » weak. The same PDF parameterization is also used to simulate shallow cumulus and stratocumulus layers. The PDF method is sufficiently general to adequately simulate these five deep, shallow, and stratiform cloud cases with a single equation set. This raises hopes that it may be possible in the future, with further refinements at coarse time step and grid spacing, to parameterize all cloud types in a large-scale model in a unified way.« less

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

NASA Technical Reports Server (NTRS)

Goldberg, Louis F.

1992-01-01

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Parameterizing deep convection using the assumed probability density function method

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

Storer, R. L.; Griffin, B. M.; Höft, J.

2015-01-06

Due to their coarse horizontal resolution, present-day climate models must parameterize deep convection. This paper presents single-column simulations of deep convection using a probability density function (PDF) parameterization. The PDF parameterization predicts the PDF of subgrid variability of turbulence, clouds, and hydrometeors. That variability is interfaced to a prognostic microphysics scheme using a Monte Carlo sampling method.The PDF parameterization is used to simulate tropical deep convection, the transition from shallow to deep convection over land, and midlatitude deep convection. These parameterized single-column simulations are compared with 3-D reference simulations. The agreement is satisfactory except when the convective forcing is weak.more » The same PDF parameterization is also used to simulate shallow cumulus and stratocumulus layers. The PDF method is sufficiently general to adequately simulate these five deep, shallow, and stratiform cloud cases with a single equation set. This raises hopes that it may be possible in the future, with further refinements at coarse time step and grid spacing, to parameterize all cloud types in a large-scale model in a unified way.« less

Parameterizing deep convection using the assumed probability density function method

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

Storer, R. L.; Griffin, B. M.; Hoft, Jan

2015-01-06

Due to their coarse horizontal resolution, present-day climate models must parameterize deep convection. This paper presents single-column simulations of deep convection using a probability density function (PDF) parameterization. The PDF parameterization predicts the PDF of subgrid variability of turbulence, clouds, and hydrometeors. That variability is interfaced to a prognostic microphysics scheme using a Monte Carlo sampling method.The PDF parameterization is used to simulate tropical deep convection, the transition from shallow to deep convection over land, and mid-latitude deep convection.These parameterized single-column simulations are compared with 3-D reference simulations. The agreement is satisfactory except when the convective forcing is weak. Themore » same PDF parameterization is also used to simulate shallow cumulus and stratocumulus layers. The PDF method is sufficiently general to adequately simulate these five deep, shallow, and stratiform cloud cases with a single equation set. This raises hopes that it may be possible in the future, with further refinements at coarse time step and grid spacing, to parameterize all cloud types in a large-scale model in a unified way.« less

Comparisons of CTH simulations with measured wave profiles for simple flyer plate experiments

DOE PAGES

Thomas, S. A.; Veeser, L. R.; Turley, W. D.; ...

2016-06-13

We conducted detailed 2-dimensional hydrodynamics calculations to assess the quality of simulations commonly used to design and analyze simple shock compression experiments. Such simple shock experiments also contain data where dynamic properties of materials are integrated together. We wished to assess how well the chosen computer hydrodynamic code could do at capturing both the simple parts of the experiments and the integral parts. We began with very simple shock experiments, in which we examined the effects of the equation of state and the compressional and tensile strength models. We increased complexity to include spallation in copper and iron and amore » solid-solid phase transformation in iron to assess the quality of the damage and phase transformation simulations. For experiments with a window, the response of both the sample and the window are integrated together, providing a good test of the material models. While CTH physics models are not perfect and do not reproduce all experimental details well, we find the models are useful; the simulations are adequate for understanding much of the dynamic process and for planning experiments. However, higher complexity in the simulations, such as adding in spall, led to greater differences between simulation and experiment. Lastly, this comparison of simulation to experiment may help guide future development of hydrodynamics codes so that they better capture the underlying physics.« less

On the multi-scale description of micro-structured fluids composed of aggregating rods

NASA Astrophysics Data System (ADS)

Perez, Marta; Scheuer, Adrien; Abisset-Chavanne, Emmanuelle; Ammar, Amine; Chinesta, Francisco; Keunings, Roland

2018-05-01

When addressing the flow of concentrated suspensions composed of rods, dense clusters are observed. Thus, the adequate modelling and simulation of such a flow requires addressing the kinematics of these dense clusters and their impact on the flow in which they are immersed. In a former work, we addressed a first modelling framework of these clusters, assumed so dense that they were considered rigid and their kinematics (flow-induced rotation) were totally defined by a symmetric tensor c with unit trace representing the cluster conformation. Then, the rigid nature of the clusters was relaxed, assuming them deformable, and a model giving the evolution of both the cluster shape and its microstructural orientation descriptor (the so-called shape and orientation tensors) was proposed. This paper compares the predictions coming from those models with finer-scale discrete simulations inspired from molecular dynamics modelling.

Nesting large-eddy simulations within mesoscale simulations for wind energy applications

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

Lundquist, J K; Mirocha, J D; Chow, F K

2008-09-08

With increasing demand for more accurate atmospheric simulations for wind turbine micrositing, for operational wind power forecasting, and for more reliable turbine design, simulations of atmospheric flow with resolution of tens of meters or higher are required. These time-dependent large-eddy simulations (LES), which resolve individual atmospheric eddies on length scales smaller than turbine blades and account for complex terrain, are possible with a range of commercial and open-source software, including the Weather Research and Forecasting (WRF) model. In addition to 'local' sources of turbulence within an LES domain, changing weather conditions outside the domain can also affect flow, suggesting thatmore » a mesoscale model provide boundary conditions to the large-eddy simulations. Nesting a large-eddy simulation within a mesoscale model requires nuanced representations of turbulence. Our group has improved the Weather and Research Forecasting model's (WRF) LES capability by implementing the Nonlinear Backscatter and Anisotropy (NBA) subfilter stress model following Kosovic (1997) and an explicit filtering and reconstruction technique to compute the Resolvable Subfilter-Scale (RSFS) stresses (following Chow et al, 2005). We have also implemented an immersed boundary method (IBM) in WRF to accommodate complex terrain. These new models improve WRF's LES capabilities over complex terrain and in stable atmospheric conditions. We demonstrate approaches to nesting LES within a mesoscale simulation for farms of wind turbines in hilly regions. Results are sensitive to the nesting method, indicating that care must be taken to provide appropriate boundary conditions, and to allow adequate spin-up of turbulence in the LES domain.« less

A comparison between block and smooth modeling in finite element simulations of tDCS*

PubMed Central

Indahlastari, Aprinda; Sadleir, Rosalind J.

2018-01-01

Current density distributions in five selected structures, namely, anterior superior temporal gyrus (ASTG), hippocampus (HIP), inferior frontal gyrus (IFG), occipital lobe (OCC) and pre-central gyrus (PRC) were investigated as part of a comparison between electrostatic finite element models constructed directly from MRI-resolution data (block models), and smoothed tetrahedral finite element models (smooth models). Three electrode configurations were applied, mimicking different tDCS therapies. Smooth model simulations were found to require three times longer to complete. The percentage differences between mean and median current densities of each model type in arbitrarily chosen brain structures ranged from −33.33–48.08%. No clear relationship was found between structure volumes and current density differences between the two model types. Tissue regions nearby the electrodes demonstrated the least percentage differences between block and smooth models. Therefore, block models may be adequate to predict current density values in cortical regions presumed targeted by tDCS. PMID:26737023

Efficiency of endoscopy units can be improved with use of discrete event simulation modeling.

PubMed

Sauer, Bryan G; Singh, Kanwar P; Wagner, Barry L; Vanden Hoek, Matthew S; Twilley, Katherine; Cohn, Steven M; Shami, Vanessa M; Wang, Andrew Y

2016-11-01

Background and study aims: The projected increased demand for health services obligates healthcare organizations to operate efficiently. Discrete event simulation (DES) is a modeling method that allows for optimization of systems through virtual testing of different configurations before implementation. The objective of this study was to identify strategies to improve the daily efficiencies of an endoscopy center with the use of DES. Methods: We built a DES model of a five procedure room endoscopy unit at a tertiary-care university medical center. After validating the baseline model, we tested alternate configurations to run the endoscopy suite and evaluated outcomes associated with each change. The main outcome measures included adequate number of preparation and recovery rooms, blocked inflow, delay times, blocked outflows, and patient cycle time. Results: Based on a sensitivity analysis, the adequate number of preparation rooms is eight and recovery rooms is nine for a five procedure room unit (total 3.4 preparation and recovery rooms per procedure room). Simple changes to procedure scheduling and patient arrival times led to a modest improvement in efficiency. Increasing the preparation/recovery rooms based on the sensitivity analysis led to significant improvements in efficiency. Conclusions: By applying tools such as DES, we can model changes in an environment with complex interactions and find ways to improve the medical care we provide. DES is applicable to any endoscopy unit and would be particularly valuable to those who are trying to improve on the efficiency of care and patient experience.

Efficiency of endoscopy units can be improved with use of discrete event simulation modeling

PubMed Central

Sauer, Bryan G.; Singh, Kanwar P.; Wagner, Barry L.; Vanden Hoek, Matthew S.; Twilley, Katherine; Cohn, Steven M.; Shami, Vanessa M.; Wang, Andrew Y.

2016-01-01

Background and study aims: The projected increased demand for health services obligates healthcare organizations to operate efficiently. Discrete event simulation (DES) is a modeling method that allows for optimization of systems through virtual testing of different configurations before implementation. The objective of this study was to identify strategies to improve the daily efficiencies of an endoscopy center with the use of DES. Methods: We built a DES model of a five procedure room endoscopy unit at a tertiary-care university medical center. After validating the baseline model, we tested alternate configurations to run the endoscopy suite and evaluated outcomes associated with each change. The main outcome measures included adequate number of preparation and recovery rooms, blocked inflow, delay times, blocked outflows, and patient cycle time. Results: Based on a sensitivity analysis, the adequate number of preparation rooms is eight and recovery rooms is nine for a five procedure room unit (total 3.4 preparation and recovery rooms per procedure room). Simple changes to procedure scheduling and patient arrival times led to a modest improvement in efficiency. Increasing the preparation/recovery rooms based on the sensitivity analysis led to significant improvements in efficiency. Conclusions: By applying tools such as DES, we can model changes in an environment with complex interactions and find ways to improve the medical care we provide. DES is applicable to any endoscopy unit and would be particularly valuable to those who are trying to improve on the efficiency of care and patient experience. PMID:27853739

Knowledge Based Cloud FE Simulation of Sheet Metal Forming Processes.

PubMed

Zhou, Du; Yuan, Xi; Gao, Haoxiang; Wang, Ailing; Liu, Jun; El Fakir, Omer; Politis, Denis J; Wang, Liliang; Lin, Jianguo

2016-12-13

The use of Finite Element (FE) simulation software to adequately predict the outcome of sheet metal forming processes is crucial to enhancing the efficiency and lowering the development time of such processes, whilst reducing costs involved in trial-and-error prototyping. Recent focus on the substitution of steel components with aluminum alloy alternatives in the automotive and aerospace sectors has increased the need to simulate the forming behavior of such alloys for ever more complex component geometries. However these alloys, and in particular their high strength variants, exhibit limited formability at room temperature, and high temperature manufacturing technologies have been developed to form them. Consequently, advanced constitutive models are required to reflect the associated temperature and strain rate effects. Simulating such behavior is computationally very expensive using conventional FE simulation techniques. This paper presents a novel Knowledge Based Cloud FE (KBC-FE) simulation technique that combines advanced material and friction models with conventional FE simulations in an efficient manner thus enhancing the capability of commercial simulation software packages. The application of these methods is demonstrated through two example case studies, namely: the prediction of a material's forming limit under hot stamping conditions, and the tool life prediction under multi-cycle loading conditions.

Knowledge Based Cloud FE Simulation of Sheet Metal Forming Processes

PubMed Central

Zhou, Du; Yuan, Xi; Gao, Haoxiang; Wang, Ailing; Liu, Jun; El Fakir, Omer; Politis, Denis J.; Wang, Liliang; Lin, Jianguo

2016-01-01

The use of Finite Element (FE) simulation software to adequately predict the outcome of sheet metal forming processes is crucial to enhancing the efficiency and lowering the development time of such processes, whilst reducing costs involved in trial-and-error prototyping. Recent focus on the substitution of steel components with aluminum alloy alternatives in the automotive and aerospace sectors has increased the need to simulate the forming behavior of such alloys for ever more complex component geometries. However these alloys, and in particular their high strength variants, exhibit limited formability at room temperature, and high temperature manufacturing technologies have been developed to form them. Consequently, advanced constitutive models are required to reflect the associated temperature and strain rate effects. Simulating such behavior is computationally very expensive using conventional FE simulation techniques. This paper presents a novel Knowledge Based Cloud FE (KBC-FE) simulation technique that combines advanced material and friction models with conventional FE simulations in an efficient manner thus enhancing the capability of commercial simulation software packages. The application of these methods is demonstrated through two example case studies, namely: the prediction of a material's forming limit under hot stamping conditions, and the tool life prediction under multi-cycle loading conditions. PMID:28060298

Description and validation of the Simple, Efficient, Dynamic, Global, Ecological Simulator (SEDGES v.1.0)

NASA Astrophysics Data System (ADS)

Paiewonsky, Pablo; Elison Timm, Oliver

2018-03-01

In this paper, we present a simple dynamic global vegetation model whose primary intended use is auxiliary to the land-atmosphere coupling scheme of a climate model, particularly one of intermediate complexity. The model simulates and provides important ecological-only variables but also some hydrological and surface energy variables that are typically either simulated by land surface schemes or else used as boundary data input for these schemes. The model formulations and their derivations are presented here, in detail. The model includes some realistic and useful features for its level of complexity, including a photosynthetic dependency on light, full coupling of photosynthesis and transpiration through an interactive canopy resistance, and a soil organic carbon dependence for bare-soil albedo. We evaluate the model's performance by running it as part of a simple land surface scheme that is driven by reanalysis data. The evaluation against observational data includes net primary productivity, leaf area index, surface albedo, and diagnosed variables relevant for the closure of the hydrological cycle. In this setup, we find that the model gives an adequate to good simulation of basic large-scale ecological and hydrological variables. Of the variables analyzed in this paper, gross primary productivity is particularly well simulated. The results also reveal the current limitations of the model. The most significant deficiency is the excessive simulation of evapotranspiration in mid- to high northern latitudes during their winter to spring transition. The model has a relative advantage in situations that require some combination of computational efficiency, model transparency and tractability, and the simulation of the large-scale vegetation and land surface characteristics under non-present-day conditions.

Simulation of nitrate reduction in groundwater - An upscaling approach from small catchments to the Baltic Sea basin

NASA Astrophysics Data System (ADS)

Hansen, A. L.; Donnelly, C.; Refsgaard, J. C.; Karlsson, I. B.

2018-01-01

This paper describes a modeling approach proposed to simulate the impact of local-scale, spatially targeted N-mitigation measures for the Baltic Sea Basin. Spatially targeted N-regulations aim at exploiting the considerable spatial differences in the natural N-reduction taking place in groundwater and surface water. While such measures can be simulated using local-scale physically-based catchment models, use of such detailed models for the 1.8 million km2 Baltic Sea basin is not feasible due to constraints on input data and computing power. Large-scale models that are able to simulate the Baltic Sea basin, on the other hand, do not have adequate spatial resolution to simulate some of the field-scale measures. Our methodology combines knowledge and results from two local-scale physically-based MIKE SHE catchment models, the large-scale and more conceptual E-HYPE model, and auxiliary data in order to enable E-HYPE to simulate how spatially targeted regulation of agricultural practices may affect N-loads to the Baltic Sea. We conclude that the use of E-HYPE with this upscaling methodology enables the simulation of the impact on N-loads of applying a spatially targeted regulation at the Baltic Sea basin scale to the correct order-of-magnitude. The E-HYPE model together with the upscaling methodology therefore provides a sound basis for large-scale policy analysis; however, we do not expect it to be sufficiently accurate to be useful for the detailed design of local-scale measures.

Numerical Simulation of the Detonation of Condensed Explosives

NASA Astrophysics Data System (ADS)

Wang, Cheng; Ye, Ting; Ning, Jianguo

Detonation process of a condensed explosive was simulated using a finite difference method. Euler equations were applied to describe the detonation flow field, an ignition and growth model for the chemical reaction and Jones-Wilkins-Lee (JWL) equations of state for the state of explosives and detonation products. Based on the simple mixture rule that assumes the reacting explosives to be a mixture of the reactant and product components, 1D and 2D codes were developed to simulate the detonation process of high explosive PBX9404. The numerical results are in good agreement with the experimental results, which demonstrates that the finite difference method, mixture rule and chemical reaction proposed in this paper are adequate and feasible.

Strategies for sperm chemotaxis in the siphonophores and ascidians: a numerical simulation study.

PubMed

Ishikawa, Makiko; Tsutsui, Hidekazu; Cosson, Jacky; Oka, Yoshitaka; Morisawa, Masaaki

2004-04-01

Chemotactic swimming behaviors of spermatozoa toward an egg have been reported in various species. The strategies underlying these behaviors, however, are poorly understood. We focused on two types of chemotaxis, one in the siphonophores and the second in the ascidians, and then proposed two models based on experimental data. Both models assumed that the radius of the path curvature of a swimming spermatozoon depends on [Ca(2+)](i), the intracellular calcium concentration. The chemotaxis in the siphonophores could be simulated in a model that assumes that [Ca(2+)](i) depends on the local concentration of the attractant in the vicinity of the spermatozoon and that a substantial time period is required for the clearance of transient high [Ca(2+)](i). In the case of ascidians, trajectories similar to those in experiments could be adequately simulated by a variant of this model that assumes that [Ca(2+)](i) depends on the time derivative of the attractant concentration. The properties of these strategies and future problems are discussed in relation to these models.

Pyramid algorithms as models of human cognition

NASA Astrophysics Data System (ADS)

Pizlo, Zygmunt; Li, Zheng

2003-06-01

There is growing body of experimental evidence showing that human perception and cognition involves mechanisms that can be adequately modeled by pyramid algorithms. The main aspect of those mechanisms is hierarchical clustering of information: visual images, spatial relations, and states as well as transformations of a problem. In this paper we review prior psychophysical and simulation results on visual size transformation, size discrimination, speed-accuracy tradeoff, figure-ground segregation, and the traveling salesman problem. We also present our new results on graph search and on the 15-puzzle.

Sea-ice deformation in a coupled ocean-sea-ice model and in satellite remote sensing data

NASA Astrophysics Data System (ADS)

Spreen, Gunnar; Kwok, Ron; Menemenlis, Dimitris; Nguyen, An T.

2017-07-01

A realistic representation of sea-ice deformation in models is important for accurate simulation of the sea-ice mass balance. Simulated sea-ice deformation from numerical simulations with 4.5, 9, and 18 km horizontal grid spacing and a viscous-plastic (VP) sea-ice rheology are compared with synthetic aperture radar (SAR) satellite observations (RGPS, RADARSAT Geophysical Processor System) for the time period 1996-2008. All three simulations can reproduce the large-scale ice deformation patterns, but small-scale sea-ice deformations and linear kinematic features (LKFs) are not adequately reproduced. The mean sea-ice total deformation rate is about 40 % lower in all model solutions than in the satellite observations, especially in the seasonal sea-ice zone. A decrease in model grid spacing, however, produces a higher density and more localized ice deformation features. The 4.5 km simulation produces some linear kinematic features, but not with the right frequency. The dependence on length scale and probability density functions (PDFs) of absolute divergence and shear for all three model solutions show a power-law scaling behavior similar to RGPS observations, contrary to what was found in some previous studies. Overall, the 4.5 km simulation produces the most realistic divergence, vorticity, and shear when compared with RGPS data. This study provides an evaluation of high and coarse-resolution viscous-plastic sea-ice simulations based on spatial distribution, time series, and power-law scaling metrics.

Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium

NASA Astrophysics Data System (ADS)

Zhang, Mingkan; Zhang, Ye

2015-03-01

A laboratory-generated hierarchical, fully heterogeneous aquifer model (FHM) provides a reference for developing and testing an upscaling approach that integrates large-scale connectivity mapping with flow and transport modeling. Based on the FHM, three hydrostratigraphic models (HSMs) that capture lithological (static) connectivity at different resolutions are created, each corresponding to a sedimentary hierarchy. Under increasing system lnK variances (0.1, 1.0, 4.5), flow upscaling is first conducted to calculate equivalent hydraulic conductivity for individual connectivity (or unit) of the HSMs. Given the computed flow fields, an instantaneous, conservative tracer test is simulated by all models. For the HSMs, two upscaling formulations are tested based on the advection-dispersion equation (ADE), implementing space versus time-dependent macrodispersivity. Comparing flow and transport predictions of the HSMs against those of the reference model, HSMs capturing connectivity at increasing resolutions are more accurate, although upscaling errors increase with system variance. Results suggest: (1) by explicitly modeling connectivity, an enhanced degree of freedom in representing dispersion can improve the ADE-based upscaled models by capturing non-Fickian transport of the FHM; (2) when connectivity is sufficiently resolved, the type of data conditioning used to model transport becomes less critical. Data conditioning, however, is influenced by the prediction goal; (3) when aquifer is weakly-to-moderately heterogeneous, the upscaled models adequately capture the transport simulation of the FHM, despite the existence of hierarchical heterogeneity at smaller scales. When aquifer is strongly heterogeneous, the upscaled models become less accurate because lithological connectivity cannot adequately capture preferential flows; (4) three-dimensional transport connectivities of the hierarchical aquifer differ quantitatively from those analyzed for two-dimensional systems. This article was corrected on 7 MAY 2015. See the end of the full text for details.

A mathematical model for simulating spring discharge and estimating sinkhole porosity in a karst watershed

NASA Astrophysics Data System (ADS)

Li, Guangquan; Field, Malcolm S.

2014-03-01

Documenting and understanding water balances in a karst watershed in which groundwater and surface water resources are strongly interconnected are important aspects for managing regional water resources. Assessing water balances in karst watersheds can be difficult, however, because karst watersheds are so very strongly affected by groundwater flows through solution conduits that are often connected to one or more sinkholes. In this paper we develop a mathematical model to approximate sinkhole porosity from discharge at a downstream spring. The model represents a combination of a traditional linear reservoir model with turbulent hydrodynamics in the solution conduit connecting the downstream spring with the upstream sinkhole, which allows for the simulation of spring discharges and estimation of sinkhole porosity. Noting that spring discharge is an integral of all aspects of water storage and flow, it is mainly dependent on the behavior of the karst aquifer as a whole and can be adequately simulated using the analytical model described in this paper. The model is advantageous in that it obviates the need for a sophisticated numerical model that is much more costly to calibrate and operate. The model is demonstrated using the St. Marks River Watershed in northwestern Florida.

Testing the skill of numerical hydraulic modeling to simulate spatiotemporal flooding patterns in the Logone floodplain, Cameroon

NASA Astrophysics Data System (ADS)

Fernández, Alfonso; Najafi, Mohammad Reza; Durand, Michael; Mark, Bryan G.; Moritz, Mark; Jung, Hahn Chul; Neal, Jeffrey; Shastry, Apoorva; Laborde, Sarah; Phang, Sui Chian; Hamilton, Ian M.; Xiao, Ningchuan

2016-08-01

Recent innovations in hydraulic modeling have enabled global simulation of rivers, including simulation of their coupled wetlands and floodplains. Accurate simulations of floodplains using these approaches may imply tremendous advances in global hydrologic studies and in biogeochemical cycling. One such innovation is to explicitly treat sub-grid channels within two-dimensional models, given only remotely sensed data in areas with limited data availability. However, predicting inundated area in floodplains using a sub-grid model has not been rigorously validated. In this study, we applied the LISFLOOD-FP hydraulic model using a sub-grid channel parameterization to simulate inundation dynamics on the Logone River floodplain, in northern Cameroon, from 2001 to 2007. Our goal was to determine whether floodplain dynamics could be simulated with sufficient accuracy to understand human and natural contributions to current and future inundation patterns. Model inputs in this data-sparse region include in situ river discharge, satellite-derived rainfall, and the shuttle radar topography mission (SRTM) floodplain elevation. We found that the model accurately simulated total floodplain inundation, with a Pearson correlation coefficient greater than 0.9, and RMSE less than 700 km2, compared to peak inundation greater than 6000 km2. Predicted discharge downstream of the floodplain matched measurements (Nash-Sutcliffe efficiency of 0.81), and indicated that net flow from the channel to the floodplain was modeled accurately. However, the spatial pattern of inundation was not well simulated, apparently due to uncertainties in SRTM elevations. We evaluated model results at 250, 500 and 1000-m spatial resolutions, and found that results are insensitive to spatial resolution. We also compared the model output against results from a run of LISFLOOD-FP in which the sub-grid channel parameterization was disabled, finding that the sub-grid parameterization simulated more realistic dynamics. These results suggest that analysis of global inundation is feasible using a sub-grid model, but that spatial patterns at sub-kilometer resolutions still need to be adequately predicted.

Comparison of Heat and Moisture Fluxes from a Modified Soil-plant-atmosphere Model with Observations from BOREAS. Chapter 3

NASA Technical Reports Server (NTRS)

Lee, Young-Hee; Mahrt, L.

2005-01-01

This study evaluates the prediction of heat and moisture fluxes from a new land surface scheme with eddy correlation data collected at the old aspen site during the Boreal Ecosystem-Atmosphere Study (BOREAS) in 1994. The model used in this study couples a multilayer vegetation model with a soil model. Inclusion of organic material in the upper soil layer is required to adequately simulate exchange between the soil and subcanopy air. Comparisons between the model and observations are discussed to reveal model misrepresentation of some aspects of the diurnal variation of subcanopy processes. Evapotranspiration

An analysis of airline landing flare data based on flight and training simulator measurements

NASA Technical Reports Server (NTRS)

Heffley, R. K.; Schulman, T. M.; Clement, T. M.

1982-01-01

Landings by experienced airline pilots transitioning to the DC-10, performed in flight and on a simulator, were analyzed and compared using a pilot-in-the-loop model of the landing maneuver. By solving for the effective feedback gains and pilot compensation which described landing technique, it was possible to discern fundamental differences in pilot behavior between the actual aircraft and the simulator. These differences were then used to infer simulator fidelity in terms of specific deficiencies and to quantify the effectiveness of training on the simulator as compared to training in flight. While training on the simulator, pilots exhibited larger effective lag in commanding the flare. The inability to compensate adequately for this lag was associated with hard or inconsistent landings. To some degree this deficiency was carried into flight, thus resulting in a slightly different and inferior landing technique than exhibited by pilots trained exclusively on the actual aircraft.

High fidelity studies of exploding foil initiator bridges, Part 3: ALEGRA MHD simulations

NASA Astrophysics Data System (ADS)

Neal, William; Garasi, Christopher

2017-01-01

Simulations of high voltage detonators, such as Exploding Bridgewire (EBW) and Exploding Foil Initiators (EFI), have historically been simple, often empirical, one-dimensional models capable of predicting parameters such as current, voltage, and in the case of EFIs, flyer velocity. Experimental methods have correspondingly generally been limited to the same parameters. With the advent of complex, first principles magnetohydrodynamic codes such as ALEGRA and ALE-MHD, it is now possible to simulate these components in three dimensions, and predict a much greater range of parameters than before. A significant improvement in experimental capability was therefore required to ensure these simulations could be adequately verified. In this third paper of a three part study, the experimental results presented in part 2 are compared against 3-dimensional MHD simulations. This improved experimental capability, along with advanced simulations, offer an opportunity to gain a greater understanding of the processes behind the functioning of EBW and EFI detonators.

Technology transfer of operator-in-the-loop simulation

NASA Technical Reports Server (NTRS)

Yae, K. H.; Lin, H. C.; Lin, T. C.; Frisch, H. P.

1994-01-01

The technology developed for operator-in-the-loop simulation in space teleoperation has been applied to Caterpillar's backhoe, wheel loader, and off-highway truck. On an SGI workstation, the simulation integrates computer modeling of kinematics and dynamics, real-time computational and visualization, and an interface with the operator through the operator's console. The console is interfaced with the workstation through an IBM-PC in which the operator's commands were digitized and sent through an RS-232 serial port. The simulation gave visual feedback adequate for the operator in the loop, with the camera's field of vision projected on a large screen in multiple view windows. The view control can emulate either stationary or moving cameras. This simulator created an innovative engineering design environment by integrating computer software and hardware with the human operator's interactions. The backhoe simulation has been adopted by Caterpillar in building a virtual reality tool for backhoe design.

Experimental study and simulation of space charge stimulated discharge

NASA Astrophysics Data System (ADS)

Noskov, M. D.; Malinovski, A. S.; Cooke, C. M.; Wright, K. A.; Schwab, A. J.

2002-11-01

The electrical discharge of volume distributed space charge in poly(methylmethacrylate) (PMMA) has been investigated both experimentally and by computer simulation. The experimental space charge was implanted in dielectric samples by exposure to a monoenergetic electron beam of 3 MeV. Electrical breakdown through the implanted space charge region within the sample was initiated by a local electric field enhancement applied to the sample surface. A stochastic-deterministic dynamic model for electrical discharge was developed and used in a computer simulation of these breakdowns. The model employs stochastic rules to describe the physical growth of the discharge channels, and deterministic laws to describe the electric field, the charge, and energy dynamics within the discharge channels and the dielectric. Simulated spatial-temporal and current characteristics of the expanding discharge structure during physical growth are quantitatively compared with the experimental data to confirm the discharge model. It was found that a single fixed set of physically based dielectric parameter values was adequate to simulate the complete family of experimental space charge discharges in PMMA. It is proposed that such a set of parameters also provides a useful means to quantify the breakdown properties of other dielectrics.

A model of the productivity of the northern pintail

USGS Publications Warehouse

Carlson, J.D.; Clark, W.R.; Klaas, E.E.

1993-01-01

We adapted a stochastic computer model to simulate productivity of the northern pintail (Anas acuta). Researchers at the Northern Prairie Wildlife Research Center of the U.S. Fish and Wildlife Service originally developed the model to simulate productivity of the mallard (A. platyrhynchos). We obtained data and descriptive information on the breeding biology of pintails from a literature review and from discussions with waterfowl biologists. All biological parameters in the productivity component of the mallard model (e.g, initial body weights, weight loss during laying and incubation, incubation time, clutch size, nest site selection characteristics) were compared with data on pintails and adjusted accordingly. The function in the mallard model that predicts nest initiation in response to pond conditions adequately mimicked pintail behavior and did not require adjustment.Recruitment rate was most sensitive to variations in parameters that control nest success, seasonal duckling survival rate, and yearling and adult body weight. We simulated upland and wetland habitat conditions in central North Dakota and compared simulation results with observed data. Simulated numbers were not significantly different from observed numbers of successful nests during wet, average, and dry wetland conditions. The simulated effect of predator barrier fencing in a study area in central North Dakota increased recruitment rate by an average of 18.4%. This modeling synthesized existing knowledge on the breeding biology of the northern pintail, identified necessary research, and furnished a useful tool for the examination and comparison of various management options.

An improved simulation of the 2015 El Niño event by optimally correcting the initial conditions and model parameters in an intermediate coupled model

NASA Astrophysics Data System (ADS)

Zhang, Rong-Hua; Tao, Ling-Jiang; Gao, Chuan

2017-09-01

Large uncertainties exist in real-time predictions of the 2015 El Niño event, which have systematic intensity biases that are strongly model-dependent. It is critically important to characterize those model biases so they can be reduced appropriately. In this study, the conditional nonlinear optimal perturbation (CNOP)-based approach was applied to an intermediate coupled model (ICM) equipped with a four-dimensional variational data assimilation technique. The CNOP-based approach was used to quantify prediction errors that can be attributed to initial conditions (ICs) and model parameters (MPs). Two key MPs were considered in the ICM: one represents the intensity of the thermocline effect, and the other represents the relative coupling intensity between the ocean and atmosphere. Two experiments were performed to illustrate the effects of error corrections, one with a standard simulation and another with an optimized simulation in which errors in the ICs and MPs derived from the CNOP-based approach were optimally corrected. The results indicate that simulations of the 2015 El Niño event can be effectively improved by using CNOP-derived error correcting. In particular, the El Niño intensity in late 2015 was adequately captured when simulations were started from early 2015. Quantitatively, the Niño3.4 SST index simulated in Dec. 2015 increased to 2.8 °C in the optimized simulation, compared with only 1.5 °C in the standard simulation. The feasibility and effectiveness of using the CNOP-based technique to improve ENSO simulations are demonstrated in the context of the 2015 El Niño event. The limitations and further applications are also discussed.

A mechanistic diagnosis of the simulation of soil CO2 efflux of the ACME Land Model

NASA Astrophysics Data System (ADS)

Liang, J.; Ricciuto, D. M.; Wang, G.; Gu, L.; Hanson, P. J.; Mayes, M. A.

2017-12-01

Accurate simulation of the CO2 efflux from soils (i.e., soil respiration) to the atmosphere is critical to project global biogeochemical cycles and the magnitude of climate change in Earth system models (ESMs). Currently, the simulated soil respiration by ESMs still have a large uncertainty. In this study, a mechanistic diagnosis of soil respiration in the Accelerated Climate Model for Energy (ACME) Land Model (ALM) was conducted using long-term observations at the Missouri Ozark AmeriFlux (MOFLUX) forest site in the central U.S. The results showed that the ALM default run significantly underestimated annual soil respiration and gross primary production (GPP), while incorrectly estimating soil water potential. Improved simulations of soil water potential with site-specific data significantly improved the modeled annual soil respiration, primarily because annual GPP was simultaneously improved. Therefore, accurate simulations of soil water potential must be carefully calibrated in ESMs. Despite improved annual soil respiration, the ALM continued to underestimate soil respiration during peak growing seasons, and to overestimate soil respiration during non-peak growing seasons. Simulations involving increased GPP during peak growing seasons increased soil respiration, while neither improved plant phenology nor increased temperature sensitivity affected the simulation of soil respiration during non-peak growing seasons. One potential reason for the overestimation of the soil respiration during non-peak growing seasons may be that the current model structure is substrate-limited, while microbial dormancy under stress may cause the system to become decomposer-limited. Further studies with more microbial data are required to provide adequate representation of soil respiration and to understand the underlying reasons for inaccurate model simulations.

Benchmarking of vertically-integrated CO2 flow simulations at the Sleipner Field, North Sea

NASA Astrophysics Data System (ADS)

Cowton, L. R.; Neufeld, J. A.; White, N. J.; Bickle, M. J.; Williams, G. A.; White, J. C.; Chadwick, R. A.

2018-06-01

Numerical modeling plays an essential role in both identifying and assessing sub-surface reservoirs that might be suitable for future carbon capture and storage projects. Accuracy of flow simulations is tested by benchmarking against historic observations from on-going CO2 injection sites. At the Sleipner project located in the North Sea, a suite of time-lapse seismic reflection surveys enables the three-dimensional distribution of CO2 at the top of the reservoir to be determined as a function of time. Previous attempts have used Darcy flow simulators to model CO2 migration throughout this layer, given the volume of injection with time and the location of the injection point. Due primarily to computational limitations preventing adequate exploration of model parameter space, these simulations usually fail to match the observed distribution of CO2 as a function of space and time. To circumvent these limitations, we develop a vertically-integrated fluid flow simulator that is based upon the theory of topographically controlled, porous gravity currents. This computationally efficient scheme can be used to invert for the spatial distribution of reservoir permeability required to minimize differences between the observed and calculated CO2 distributions. When a uniform reservoir permeability is assumed, inverse modeling is unable to adequately match the migration of CO2 at the top of the reservoir. If, however, the width and permeability of a mapped channel deposit are allowed to independently vary, a satisfactory match between the observed and calculated CO2 distributions is obtained. Finally, the ability of this algorithm to forecast the flow of CO2 at the top of the reservoir is assessed. By dividing the complete set of seismic reflection surveys into training and validation subsets, we find that the spatial pattern of permeability required to match the training subset can successfully predict CO2 migration for the validation subset. This ability suggests that it might be feasible to forecast migration patterns into the future with a degree of confidence. Nevertheless, our analysis highlights the difficulty in estimating reservoir parameters away from the region swept by CO2 without additional observational constraints.

Growing C4 perennial grass for bioenergy using a new Agro-BGC ecosystem model

NASA Astrophysics Data System (ADS)

di Vittorio, A. V.; Anderson, R. S.; Miller, N. L.; Running, S. W.

2009-12-01

Accurate, spatially gridded estimates of bioenergy crop yields require 1) biophysically accurate crop growth models and 2) careful parameterization of unavailable inputs to these models. To meet the first requirement we have added the capacity to simulate C4 perennial grass as a bioenergy crop to the Biome-BGC ecosystem model. This new model, hereafter referred to as Agro-BGC, includes enzyme driven C4 photosynthesis, individual live and dead leaf, stem, and root carbon/nitrogen pools, separate senescence and litter fall processes, fruit growth, optional annual seeding, flood irrigation, a growing degree day phenology with a killing frost option, and a disturbance handler that effectively simulates fertilization, harvest, fire, and incremental irrigation. There are four Agro-BGC vegetation parameters that are unavailable for Panicum virgatum (switchgrass), and to meet the second requirement we have optimized the model across multiple calibration sites to obtain representative values for these parameters. We have verified simulated switchgrass yields against observations at three non-calibration sites in IL. Agro-BGC simulates switchgrass growth and yield at harvest very well at a single site. Our results suggest that a multi-site optimization scheme would be adequate for producing regional-scale estimates of bioenergy crop yields on high spatial resolution grids.

Development and validation of a modified Hybrid-III six-year-old dummy model for simulating submarining in motor-vehicle crashes.

PubMed

Hu, Jingwen; Klinich, Kathleen D; Reed, Matthew P; Kokkolaras, Michael; Rupp, Jonathan D

2012-06-01

In motor-vehicle crashes, young school-aged children restrained by vehicle seat belt systems often suffer from abdominal injuries due to submarining. However, the current anthropomorphic test device, so-called "crash dummy", is not adequate for proper simulation of submarining. In this study, a modified Hybrid-III six-year-old dummy model capable of simulating and predicting submarining was developed using MADYMO (TNO Automotive Safety Solutions). The model incorporated improved pelvis and abdomen geometry and properties previously tested in a modified physical dummy. The model was calibrated and validated against four sled tests under two test conditions with and without submarining using a multi-objective optimization method. A sensitivity analysis using this validated child dummy model showed that dummy knee excursion, torso rotation angle, and the difference between head and knee excursions were good predictors for submarining status. It was also shown that restraint system design variables, such as lap belt angle, D-ring height, and seat coefficient of friction (COF), may have opposite effects on head and abdomen injury risks; therefore child dummies and dummy models capable of simulating submarining are crucial for future restraint system design optimization for young school-aged children. Copyright © 2011 IPEM. Published by Elsevier Ltd. All rights reserved.

Transfer of training for aerospace operations: How to measure, validate, and improve it

NASA Technical Reports Server (NTRS)

Cohen, Malcolm M.

1993-01-01

It has been a commonly accepted practice to train pilots and astronauts in expensive, extremely sophisticated, high fidelity simulators, with as much of the real-world feel and response as possible. High fidelity and high validity have often been assumed to be inextricably interwoven, although this assumption may not be warranted. The Project Mercury rate-damping task on the Naval Air Warfare Center's Human Centrifuge Dynamic Flight Simulator, the shuttle landing task on the NASA-ARC Vertical Motion Simulator, and the almost complete acceptance by the airline industry of full-up Boeing 767 flight simulators, are just a few examples of this approach. For obvious reasons, the classical models of transfer of training have never been adequately evaluated in aerospace operations, and there have been few, if any, scientifically valid replacements for the classical models. This paper reviews some of the earlier work involving transfer of training in aerospace operations, and discusses some of the methods by which appropriate criteria for assessing the validity of training may be established.

Measurements and Computations of Flow in an Urban Street System

NASA Astrophysics Data System (ADS)

Castro, Ian P.; Xie, Zheng-Tong; Fuka, V.; Robins, Alan G.; Carpentieri, M.; Hayden, P.; Hertwig, D.; Coceal, O.

2017-02-01

We present results from laboratory and computational experiments on the turbulent flow over an array of rectangular blocks modelling a typical, asymmetric urban canopy at various orientations to the approach flow. The work forms part of a larger study on dispersion within such arrays (project DIPLOS) and concentrates on the nature of the mean flow and turbulence fields within the canopy region, recognising that unless the flow field is adequately represented in computational models there is no reason to expect realistic simulations of the nature of the dispersion of pollutants emitted within the canopy. Comparisons between the experimental data and those obtained from both large-eddy simulation (LES) and direct numerical simulation (DNS) are shown and it is concluded that careful use of LES can produce generally excellent agreement with laboratory and DNS results, lending further confidence in the use of LES for such situations. Various crucial issues are discussed and advice offered to both experimentalists and those seeking to compute canopy flows with turbulence resolving models.

On kinetic modelling for solar redox thermochemical H2O and CO2 splitting over NiFe2O4 for H2, CO and syngas production.

PubMed

Dimitrakis, Dimitrios A; Syrigou, Maria; Lorentzou, Souzana; Kostoglou, Margaritis; Konstandopoulos, Athanasios G

2017-10-11

This study aims at developing a kinetic model that can adequately describe solar thermochemical water and carbon dioxide splitting with nickel ferrite powder as the active redox material. The kinetic parameters of water splitting of a previous study are revised to include transition times and new kinetic parameters for carbon dioxide splitting are developed. The computational results show a satisfactory agreement with experimental data and continuous multicycle operation under varying operating conditions is simulated. Different test cases are explored in order to improve the product yield. At first a parametric analysis is conducted, investigating the appropriate duration of the oxidation and the thermal reduction step that maximizes the hydrogen yield. Subsequently, a non-isothermal oxidation step is simulated and proven as an interesting option for increasing the hydrogen production. The kinetic model is adapted to simulate the production yields in structured solar reactor components, i.e. extruded monolithic structures, as well.

Kinetic barriers in the isomerization of substituted ureas: implications for computer-aided drug design.

PubMed

Loeffler, Johannes R; Ehmki, Emanuel S R; Fuchs, Julian E; Liedl, Klaus R

2016-05-01

Urea derivatives are ubiquitously found in many chemical disciplines. N,N'-substituted ureas may show different conformational preferences depending on their substitution pattern. The high energetic barrier for isomerization of the cis and trans state poses additional challenges on computational simulation techniques aiming at a reproduction of the biological properties of urea derivatives. Herein, we investigate energetics of urea conformations and their interconversion using a broad spectrum of methodologies ranging from data mining, via quantum chemistry to molecular dynamics simulation and free energy calculations. We find that the inversion of urea conformations is inherently slow and beyond the time scale of typical simulation protocols. Therefore, extra care needs to be taken by computational chemists to work with appropriate model systems. We find that both knowledge-driven approaches as well as physics-based methods may guide molecular modelers towards accurate starting structures for expensive calculations to ensure that conformations of urea derivatives are modeled as adequately as possible.

Analysis and Description of HOLTIN Service Provision for AECG monitoring in Complex Indoor Environments

PubMed Central

Led, Santiago; Azpilicueta, Leire; Aguirre, Erik; de Espronceda, Miguel Martínez; Serrano, Luis; Falcone, Francisco

2013-01-01

In this work, a novel ambulatory ECG monitoring device developed in-house called HOLTIN is analyzed when operating in complex indoor scenarios. The HOLTIN system is described, from the technological platform level to its functional model. In addition, by using in-house 3D ray launching simulation code, the wireless channel behavior, which enables ubiquitous operation, is performed. The effect of human body presence is taken into account by a novel simplified model embedded within the 3D Ray Launching code. Simulation as well as measurement results are presented, showing good agreement. These results may aid in the adequate deployment of this novel device to automate conventional medical processes, increasing the coverage radius and optimizing energy consumption. PMID:23584122

The Impact of Simulated Mesoscale Convective Systems on Global Precipitation: A Multiscale Modeling Study

NASA Technical Reports Server (NTRS)

Tao, Wei-Kuo; Chern, Jiun-Dar

2017-01-01

The importance of precipitating mesoscale convective systems (MCSs) has been quantified from TRMM precipitation radar and microwave imager retrievals. MCSs generate more than 50% of the rainfall in most tropical regions. MCSs usually have horizontal scales of a few hundred kilometers (km); therefore, a large domain with several hundred km is required for realistic simulations of MCSs in cloud-resolving models (CRMs). Almost all traditional global and climate models do not have adequate parameterizations to represent MCSs. Typical multi-scale modeling frameworks (MMFs) may also lack the resolution (4 km grid spacing) and domain size (128 km) to realistically simulate MCSs. In this study, the impact of MCSs on precipitation is examined by conducting model simulations using the Goddard Cumulus Ensemble (GCE) model and Goddard MMF (GMMF). The results indicate that both models can realistically simulate MCSs with more grid points (i.e., 128 and 256) and higher resolutions (1 or 2 km) compared to those simulations with fewer grid points (i.e., 32 and 64) and low resolution (4 km). The modeling results also show the strengths of the Hadley circulations, mean zonal and regional vertical velocities, surface evaporation, and amount of surface rainfall are weaker or reduced in the GMMF when using more CRM grid points and higher CRM resolution. In addition, the results indicate that large-scale surface evaporation and wind feed back are key processes for determining the surface rainfall amount in the GMMF. A sensitivity test with reduced sea surface temperatures shows both reduced surface rainfall and evaporation.

The impact of simulated mesoscale convective systems on global precipitation: A multiscale modeling study

NASA Astrophysics Data System (ADS)

Tao, Wei-Kuo; Chern, Jiun-Dar

2017-06-01

The importance of precipitating mesoscale convective systems (MCSs) has been quantified from TRMM precipitation radar and microwave imager retrievals. MCSs generate more than 50% of the rainfall in most tropical regions. MCSs usually have horizontal scales of a few hundred kilometers (km); therefore, a large domain with several hundred km is required for realistic simulations of MCSs in cloud-resolving models (CRMs). Almost all traditional global and climate models do not have adequate parameterizations to represent MCSs. Typical multiscale modeling frameworks (MMFs) may also lack the resolution (4 km grid spacing) and domain size (128 km) to realistically simulate MCSs. The impact of MCSs on precipitation is examined by conducting model simulations using the Goddard Cumulus Ensemble (GCE, a CRM) model and Goddard MMF that uses the GCEs as its embedded CRMs. Both models can realistically simulate MCSs with more grid points (i.e., 128 and 256) and higher resolutions (1 or 2 km) compared to those simulations with fewer grid points (i.e., 32 and 64) and low resolution (4 km). The modeling results also show the strengths of the Hadley circulations, mean zonal and regional vertical velocities, surface evaporation, and amount of surface rainfall are weaker or reduced in the Goddard MMF when using more CRM grid points and higher CRM resolution. In addition, the results indicate that large-scale surface evaporation and wind feedback are key processes for determining the surface rainfall amount in the GMMF. A sensitivity test with reduced sea surface temperatures shows both reduced surface rainfall and evaporation.

Alternative ways of using field-based estimates to calibrate ecosystem models and their implications for carbon cycle studies

USGS Publications Warehouse

He, Yujie; Zhuang, Qianlai; McGuire, David; Liu, Yaling; Chen, Min

2013-01-01

Model-data fusion is a process in which field observations are used to constrain model parameters. How observations are used to constrain parameters has a direct impact on the carbon cycle dynamics simulated by ecosystem models. In this study, we present an evaluation of several options for the use of observations in modeling regional carbon dynamics and explore the implications of those options. We calibrated the Terrestrial Ecosystem Model on a hierarchy of three vegetation classification levels for the Alaskan boreal forest: species level, plant-functional-type level (PFT level), and biome level, and we examined the differences in simulated carbon dynamics. Species-specific field-based estimates were directly used to parameterize the model for species-level simulations, while weighted averages based on species percent cover were used to generate estimates for PFT- and biome-level model parameterization. We found that calibrated key ecosystem process parameters differed substantially among species and overlapped for species that are categorized into different PFTs. Our analysis of parameter sets suggests that the PFT-level parameterizations primarily reflected the dominant species and that functional information of some species were lost from the PFT-level parameterizations. The biome-level parameterization was primarily representative of the needleleaf PFT and lost information on broadleaf species or PFT function. Our results indicate that PFT-level simulations may be potentially representative of the performance of species-level simulations while biome-level simulations may result in biased estimates. Improved theoretical and empirical justifications for grouping species into PFTs or biomes are needed to adequately represent the dynamics of ecosystem functioning and structure.

Investigation of different modeling approaches for computational fluid dynamics simulation of high-pressure rocket combustors

NASA Astrophysics Data System (ADS)

Ivancic, B.; Riedmann, H.; Frey, M.; Knab, O.; Karl, S.; Hannemann, K.

2016-07-01

The paper summarizes technical results and first highlights of the cooperation between DLR and Airbus Defence and Space (DS) within the work package "CFD Modeling of Combustion Chamber Processes" conducted in the frame of the Propulsion 2020 Project. Within the addressed work package, DLR Göttingen and Airbus DS Ottobrunn have identified several test cases where adequate test data are available and which can be used for proper validation of the computational fluid dynamics (CFD) tools. In this paper, the first test case, the Penn State chamber (RCM1), is discussed. Presenting the simulation results from three different tools, it is shown that the test case can be computed properly with steady-state Reynolds-averaged Navier-Stokes (RANS) approaches. The achieved simulation results reproduce the measured wall heat flux as an important validation parameter very well but also reveal some inconsistencies in the test data which are addressed in this paper.

Using travel times to simulate multi-dimensional bioreactive transport in time-periodic flows.

PubMed

Sanz-Prat, Alicia; Lu, Chuanhe; Finkel, Michael; Cirpka, Olaf A

2016-04-01

In travel-time models, the spatially explicit description of reactive transport is replaced by associating reactive-species concentrations with the travel time or groundwater age at all locations. These models have been shown adequate for reactive transport in river-bank filtration under steady-state flow conditions. Dynamic hydrological conditions, however, can lead to fluctuations of infiltration velocities, putting the validity of travel-time models into question. In transient flow, the local travel-time distributions change with time. We show that a modified version of travel-time based reactive transport models is valid if only the magnitude of the velocity fluctuates, whereas its spatial orientation remains constant. We simulate nonlinear, one-dimensional, bioreactive transport involving oxygen, nitrate, dissolved organic carbon, aerobic and denitrifying bacteria, considering periodic fluctuations of velocity. These fluctuations make the bioreactive system pulsate: The aerobic zone decreases at times of low velocity and increases at those of high velocity. For the case of diurnal fluctuations, the biomass concentrations cannot follow the hydrological fluctuations and a transition zone containing both aerobic and obligatory denitrifying bacteria is established, whereas a clear separation of the two types of bacteria prevails in the case of seasonal velocity fluctuations. We map the 1-D results to a heterogeneous, two-dimensional domain by means of the mean groundwater age for steady-state flow in both domains. The mapped results are compared to simulation results of spatially explicit, two-dimensional, advective-dispersive-bioreactive transport subject to the same relative fluctuations of velocity as in the one-dimensional model. The agreement between the mapped 1-D and the explicit 2-D results is excellent. We conclude that travel-time models of nonlinear bioreactive transport are adequate in systems of time-periodic flow if the flow direction does not change. Copyright © 2016 Elsevier B.V. All rights reserved.

Numerical simulation of volume-controlled mechanical ventilated respiratory system with 2 different lungs.

PubMed

Shi, Yan; Zhang, Bolun; Cai, Maolin; Zhang, Xiaohua Douglas

2017-09-01

Mechanical ventilation is a key therapy for patients who cannot breathe adequately by themselves, and dynamics of mechanical ventilation system is of great significance for life support of patients. Recently, models of mechanical ventilated respiratory system with 1 lung are used to simulate the respiratory system of patients. However, humans have 2 lungs. When the respiratory characteristics of 2 lungs are different, a single-lung model cannot reflect real respiratory system. In this paper, to illustrate dynamic characteristics of mechanical ventilated respiratory system with 2 different lungs, we propose a mathematical model of mechanical ventilated respiratory system with 2 different lungs and conduct experiments to verify the model. Furthermore, we study the dynamics of mechanical ventilated respiratory system with 2 different lungs. This research study can be used for improving the efficiency and safety of volume-controlled mechanical ventilation system. Copyright © 2016 John Wiley & Sons, Ltd.

Evaluating simulated functional trait patterns and quantifying modelled trait diversity effects on simulated ecosystem fluxes

NASA Astrophysics Data System (ADS)

Pavlick, R.; Schimel, D.

2014-12-01

Dynamic Global Vegetation Models (DGVMs) typically employ only a small set of Plant Functional Types (PFTs) to represent the vast diversity of observed vegetation forms and functioning. There is growing evidence, however, that this abstraction may not adequately represent the observed variation in plant functional traits, which is thought to play an important role for many ecosystem functions and for ecosystem resilience to environmental change. The geographic distribution of PFTs in these models is also often based on empirical relationships between present-day climate and vegetation patterns. Projections of future climate change, however, point toward the possibility of novel regional climates, which could lead to no-analog vegetation compositions incompatible with the PFT paradigm. Here, we present results from the Jena Diversity-DGVM (JeDi-DGVM), a novel traits-based vegetation model, which simulates a large number of hypothetical plant growth strategies constrained by functional tradeoffs, thereby allowing for a more flexible temporal and spatial representation of the terrestrial biosphere. First, we compare simulated present-day geographical patterns of functional traits with empirical trait observations (in-situ and from airborne imaging spectroscopy). The observed trait patterns are then used to improve the tradeoff parameterizations of JeDi-DGVM. Finally, focusing primarily on the simulated leaf traits, we run the model with various amounts of trait diversity. We quantify the effects of these modeled biodiversity manipulations on simulated ecosystem fluxes and stocks for both present-day conditions and transient climate change scenarios. The simulation results reveal that the coarse treatment of plant functional traits by current PFT-based vegetation models may contribute substantial uncertainty regarding carbon-climate feedbacks. Further development of trait-based models and further investment in global in-situ and spectroscopic plant trait observations are needed.

Data-driven RANS for simulations of large wind farms

NASA Astrophysics Data System (ADS)

Iungo, G. V.; Viola, F.; Ciri, U.; Rotea, M. A.; Leonardi, S.

2015-06-01

In the wind energy industry there is a growing need for real-time predictions of wind turbine wake flows in order to optimize power plant control and inhibit detrimental wake interactions. To this aim, a data-driven RANS approach is proposed in order to achieve very low computational costs and adequate accuracy through the data assimilation procedure. The RANS simulations are implemented with a classical Boussinesq hypothesis and a mixing length turbulence closure model, which is calibrated through the available data. High-fidelity LES simulations of a utility-scale wind turbine operating with different tip speed ratios are used as database. It is shown that the mixing length model for the RANS simulations can be calibrated accurately through the Reynolds stress of the axial and radial velocity components, and the gradient of the axial velocity in the radial direction. It is found that the mixing length is roughly invariant in the very near wake, then it increases linearly with the downstream distance in the diffusive region. The variation rate of the mixing length in the downstream direction is proposed as a criterion to detect the transition between near wake and transition region of a wind turbine wake. Finally, RANS simulations were performed with the calibrated mixing length model, and a good agreement with the LES simulations is observed.

Thermal modelling using discrete vasculature for thermal therapy: a review

PubMed Central

Kok, H.P.; Gellermann, J.; van den Berg, C.A.T.; Stauffer, P.R.; Hand, J.W.; Crezee, J.

2013-01-01

Reliable temperature information during clinical hyperthermia and thermal ablation is essential for adequate treatment control, but conventional temperature measurements do not provide 3D temperature information. Treatment planning is a very useful tool to improve treatment quality and substantial progress has been made over the last decade. Thermal modelling is a very important and challenging aspect of hyperthermia treatment planning. Various thermal models have been developed for this purpose, with varying complexity. Since blood perfusion is such an important factor in thermal redistribution of energy in in vivo tissue, thermal simulations are most accurately performed by modelling discrete vasculature. This review describes the progress in thermal modelling with discrete vasculature for the purpose of hyperthermia treatment planning and thermal ablation. There has been significant progress in thermal modelling with discrete vasculature. Recent developments have made real-time simulations possible, which can provide feedback during treatment for improved therapy. Future clinical application of thermal modelling with discrete vasculature in hyperthermia treatment planning is expected to further improve treatment quality. PMID:23738700

Development of a canine model to enable the preclinical assessment of pH-dependent absorption of test compounds.

PubMed

Fancher, R Marcus; Zhang, Hongjian; Sleczka, Bogdan; Derbin, George; Rockar, Richard; Marathe, Punit

2011-07-01

A preclinical canine model capable of predicting a compound's potential for pH-dependent absorption in humans was developed. This involved the surgical insertion of a gastrostomy feeding tube into the stomach of a beagle dog. The tube was sutured in position to allow frequent withdrawal of gastric fluid for pH measurement. Therefore, it was possible to measure pH in the stomach and assess the effect of gastric pH-modifying agents on the absorption of various test compounds. Fasted gastric pH in the dog showed considerable inter- and intra-animal variability. Pretreatment of pentagastrin (6 µg/kg intramuscularly) 20 min prior to test compound administration was determined to be adequate for simulating fasting stomach pH in humans. Pretreatment with famotidine [40 mg orally] 1 h prior to test compound administration was determined to be adequate for simulating human gastric pH when acid-reducing agents are coadministered. Pentagastrin and famotidine pretreatments were used to test two discovery compounds and distinct differences in their potential for pH-dependent absorption were observed. The model described herein can be used preclinically to screen out compounds, differentiate compounds, and support the assessment of various formulation- and prodrug-based strategies to mitigate the pH effect. Copyright © 2011 Wiley-Liss, Inc. and the American Pharmacists Association

A new numerical benchmark for variably saturated variable-density flow and transport in porous media

NASA Astrophysics Data System (ADS)

Guevara, Carlos; Graf, Thomas

2016-04-01

In subsurface hydrological systems, spatial and temporal variations in solute concentration and/or temperature may affect fluid density and viscosity. These variations could lead to potentially unstable situations, in which a dense fluid overlies a less dense fluid. These situations could produce instabilities that appear as dense plume fingers migrating downwards counteracted by vertical upwards flow of freshwater (Simmons et al., Transp. Porous Medium, 2002). As a result of unstable variable-density flow, solute transport rates are increased over large distances and times as compared to constant-density flow. The numerical simulation of variable-density flow in saturated and unsaturated media requires corresponding benchmark problems against which a computer model is validated (Diersch and Kolditz, Adv. Water Resour, 2002). Recorded data from a laboratory-scale experiment of variable-density flow and solute transport in saturated and unsaturated porous media (Simmons et al., Transp. Porous Medium, 2002) is used to define a new numerical benchmark. The HydroGeoSphere code (Therrien et al., 2004) coupled with PEST (www.pesthomepage.org) are used to obtain an optimized parameter set capable of adequately representing the data set by Simmons et al., (2002). Fingering in the numerical model is triggered using random hydraulic conductivity fields. Due to the inherent randomness, a large number of simulations were conducted in this study. The optimized benchmark model adequately predicts the plume behavior and the fate of solutes. This benchmark is useful for model verification of variable-density flow problems in saturated and/or unsaturated media.

Large-eddy simulations with wall models

NASA Technical Reports Server (NTRS)

Cabot, W.

1995-01-01

The near-wall viscous and buffer regions of wall-bounded flows generally require a large expenditure of computational resources to be resolved adequately, even in large-eddy simulation (LES). Often as much as 50% of the grid points in a computational domain are devoted to these regions. The dense grids that this implies also generally require small time steps for numerical stability and/or accuracy. It is commonly assumed that the inner wall layers are near equilibrium, so that the standard logarithmic law can be applied as the boundary condition for the wall stress well away from the wall, for example, in the logarithmic region, obviating the need to expend large amounts of grid points and computational time in this region. This approach is commonly employed in LES of planetary boundary layers, and it has also been used for some simple engineering flows. In order to calculate accurately a wall-bounded flow with coarse wall resolution, one requires the wall stress as a boundary condition. The goal of this work is to determine the extent to which equilibrium and boundary layer assumptions are valid in the near-wall regions, to develop models for the inner layer based on such assumptions, and to test these modeling ideas in some relatively simple flows with different pressure gradients, such as channel flow and flow over a backward-facing step. Ultimately, models that perform adequately in these situations will be applied to more complex flow configurations, such as an airfoil.

Equilibrium, kinetic, and reactive transport models for plutonium

NASA Astrophysics Data System (ADS)

Schwantes, Jon Michael

Equilibrium, kinetic, and reactive transport models for plutonium (Pu) have been developed to help meet environmental concerns posed by past war-related and present and future peacetime nuclear technologies. A thorough review of the literature identified several hurdles that needed to be overcome in order to develop capable predictive tools for Pu. These hurdles include: (1) missing or ill-defined chemical equilibrium and kinetic constants for environmentally important Pu species; (2) no adequate conceptual model describing the formation of Pu oxy/hydroxide colloids and solids; and (3) an inability of two-phase reactive transport models to adequately simulate Pu behavior in the presence of colloids. A computer program called INVRS K was developed that integrates the geochemical modeling software of PHREEQC with a nonlinear regression routine. This program provides a tool for estimating equilibrium and kinetic constants from experimental data. INVRS K was used to regress on binding constants for Pu sorbing onto various mineral and humic surfaces. These constants enhance the thermodynamic database for Pu and improve the capability of current predictive tools. Time and temperature studies of the Pu intrinsic colloid were also conducted and results of these studies were presented here. Formation constants for the fresh and aged Pu intrinsic colloid were regressed upon using INVRS K. From these results, it was possible to develop a cohesive diagenetic model that describes the formation of Pu oxy/hydroxide colloids and solids. This model provides for the first time a means of deciphering historically unexplained observations with respect to the Pu intrinsic colloid, as well as a basis for simulating the behavior within systems containing these solids. Discussion of the development and application of reactive transport models is also presented and includes: (1) the general application of a 1-D in flow, three-phase (i.e., dissolved, solid, and colloidal), reactive transport model; (2) a simulation of the effects of dissolution of PuO2 solid and radiolysis on the behavior of Pu diffusing out of a confined pore space; and (3) application of a steady-state three phase reactive transport model to groundwater at the Nevada Test Site.

User's guide for MAGIC-Meteorologic and hydrologic genscn (generate scenarios) input converter

USGS Publications Warehouse

Ortel, Terry W.; Martin, Angel

2010-01-01

Meteorologic and hydrologic data used in watershed modeling studies are collected by various agencies and organizations, and stored in various formats. Data may be in a raw, un-processed format with little or no quality control, or may be checked for validity before being made available. Flood-simulation systems require data in near real-time so that adequate flood warnings can be made. Additionally, forecasted data are needed to operate flood-control structures to potentially mitigate flood damages. Because real-time data are of a provisional nature, missing data may need to be estimated for use in floodsimulation systems. The Meteorologic and Hydrologic GenScn (Generate Scenarios) Input Converter (MAGIC) can be used to convert data from selected formats into the Hydrologic Simulation System-Fortran hourly-observations format for input to a Watershed Data Management database, for use in hydrologic modeling studies. MAGIC also can reformat the data to the Full Equations model time-series format, for use in hydraulic modeling studies. Examples of the application of MAGIC for use in the flood-simulation system for Salt Creek in northeastern Illinois are presented in this report.

EnergyPlus Run Time Analysis

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

Hong, Tianzhen; Buhl, Fred; Haves, Philip

2008-09-20

EnergyPlus is a new generation building performance simulation program offering many new modeling capabilities and more accurate performance calculations integrating building components in sub-hourly time steps. However, EnergyPlus runs much slower than the current generation simulation programs. This has become a major barrier to its widespread adoption by the industry. This paper analyzed EnergyPlus run time from comprehensive perspectives to identify key issues and challenges of speeding up EnergyPlus: studying the historical trends of EnergyPlus run time based on the advancement of computers and code improvements to EnergyPlus, comparing EnergyPlus with DOE-2 to understand and quantify the run time differences,more » identifying key simulation settings and model features that have significant impacts on run time, and performing code profiling to identify which EnergyPlus subroutines consume the most amount of run time. This paper provides recommendations to improve EnergyPlus run time from the modeler?s perspective and adequate computing platforms. Suggestions of software code and architecture changes to improve EnergyPlus run time based on the code profiling results are also discussed.« less

Sensitivity analysis of conservative and reactive stream transient storage models applied to field data from multiple-reach experiments

USGS Publications Warehouse

Gooseff, M.N.; Bencala, K.E.; Scott, D.T.; Runkel, R.L.; McKnight, Diane M.

2005-01-01

The transient storage model (TSM) has been widely used in studies of stream solute transport and fate, with an increasing emphasis on reactive solute transport. In this study we perform sensitivity analyses of a conservative TSM and two different reactive solute transport models (RSTM), one that includes first-order decay in the stream and the storage zone, and a second that considers sorption of a reactive solute on streambed sediments. Two previously analyzed data sets are examined with a focus on the reliability of these RSTMs in characterizing stream and storage zone solute reactions. Sensitivities of simulations to parameters within and among reaches, parameter coefficients of variation, and correlation coefficients are computed and analyzed. Our results indicate that (1) simulated values have the greatest sensitivity to parameters within the same reach, (2) simulated values are also sensitive to parameters in reaches immediately upstream and downstream (inter-reach sensitivity), (3) simulated values have decreasing sensitivity to parameters in reaches farther downstream, and (4) in-stream reactive solute data provide adequate data to resolve effective storage zone reaction parameters, given the model formulations. Simulations of reactive solutes are shown to be equally sensitive to transport parameters and effective reaction parameters of the model, evidence of the control of physical transport on reactive solute dynamics. Similar to conservative transport analysis, reactive solute simulations appear to be most sensitive to data collected during the rising and falling limb of the concentration breakthrough curve. ?? 2005 Elsevier Ltd. All rights reserved.

Projected strengthening of Amazonian dry season by constrained climate model simulations

NASA Astrophysics Data System (ADS)

Boisier, Juan P.; Ciais, Philippe; Ducharne, Agnès; Guimberteau, Matthieu

2015-07-01

The vulnerability of Amazonian rainforest, and the ecological services it provides, depends on an adequate supply of dry-season water, either as precipitation or stored soil moisture. How the rain-bearing South American monsoon will evolve across the twenty-first century is thus a question of major interest. Extensive savanization, with its loss of forest carbon stock and uptake capacity, is an extreme although very uncertain scenario. We show that the contrasting rainfall projections simulated for Amazonia by 36 global climate models (GCMs) can be reproduced with empirical precipitation models, calibrated with historical GCM data as functions of the large-scale circulation. A set of these simple models was therefore calibrated with observations and used to constrain the GCM simulations. In agreement with the current hydrologic trends, the resulting projection towards the end of the twenty-first century is for a strengthening of the monsoon seasonal cycle, and a dry-season lengthening in southern Amazonia. With this approach, the increase in the area subjected to lengthy--savannah-prone--dry seasons is substantially larger than the GCM-simulated one. Our results confirm the dominant picture shown by the state-of-the-art GCMs, but suggest that the `model democracy' view of these impacts can be significantly underestimated.

GCM simulations of Titan's middle and lower atmosphere and comparison to observations

NASA Astrophysics Data System (ADS)

Lora, Juan M.; Lunine, Jonathan I.; Russell, Joellen L.

2015-04-01

Simulation results are presented from a new general circulation model (GCM) of Titan, the Titan Atmospheric Model (TAM), which couples the Flexible Modeling System (FMS) spectral dynamical core to a suite of external/sub-grid-scale physics. These include a new non-gray radiative transfer module that takes advantage of recent data from Cassini-Huygens, large-scale condensation and quasi-equilibrium moist convection schemes, a surface model with "bucket" hydrology, and boundary layer turbulent diffusion. The model produces a realistic temperature structure from the surface to the lower mesosphere, including a stratopause, as well as satisfactory superrotation. The latter is shown to depend on the dynamical core's ability to build up angular momentum from surface torques. Simulated latitudinal temperature contrasts are adequate, compared to observations, and polar temperature anomalies agree with observations. In the lower atmosphere, the insolation distribution is shown to strongly impact turbulent fluxes, and surface heating is maximum at mid-latitudes. Surface liquids are unstable at mid- and low-latitudes, and quickly migrate poleward. The simulated humidity profile and distribution of surface temperatures, compared to observations, corroborate the prevalence of dry conditions at low latitudes. Polar cloud activity is well represented, though the observed mid-latitude clouds remain somewhat puzzling, and some formation alternatives are suggested.

Evaluating the performance of coupled snow-soil models in SURFEXv8 to simulate the permafrost thermal regime at a high Arctic site

NASA Astrophysics Data System (ADS)

Barrere, Mathieu; Domine, Florent; Decharme, Bertrand; Morin, Samuel; Vionnet, Vincent; Lafaysse, Matthieu

2017-09-01

Climate change projections still suffer from a limited representation of the permafrost-carbon feedback. Predicting the response of permafrost temperature to climate change requires accurate simulations of Arctic snow and soil properties. This study assesses the capacity of the coupled land surface and snow models ISBA-Crocus and ISBA-ES to simulate snow and soil properties at Bylot Island, a high Arctic site. Field measurements complemented with ERA-Interim reanalyses were used to drive the models and to evaluate simulation outputs. Snow height, density, temperature, thermal conductivity and thermal insulance are examined to determine the critical variables involved in the soil and snow thermal regime. Simulated soil properties are compared to measurements of thermal conductivity, temperature and water content. The simulated snow density profiles are unrealistic, which is most likely caused by the lack of representation in snow models of the upward water vapor fluxes generated by the strong temperature gradients within the snowpack. The resulting vertical profiles of thermal conductivity are inverted compared to observations, with high simulated values at the bottom of the snowpack. Still, ISBA-Crocus manages to successfully simulate the soil temperature in winter. Results are satisfactory in summer, but the temperature of the top soil could be better reproduced by adequately representing surface organic layers, i.e., mosses and litter, and in particular their water retention capacity. Transition periods (soil freezing and thawing) are the least well reproduced because the high basal snow thermal conductivity induces an excessively rapid heat transfer between the soil and the snow in simulations. Hence, global climate models should carefully consider Arctic snow thermal properties, and especially the thermal conductivity of the basal snow layer, to perform accurate predictions of the permafrost evolution under climate change.

A modeling and simulation approach to characterize methadone QT prolongation using pooled data from five clinical trials in MMT patients.

PubMed

Florian, J; Garnett, C E; Nallani, S C; Rappaport, B A; Throckmorton, D C

2012-04-01

Pharmacokinetic (PK)-pharmacodynamic modeling and simulation were used to establish a link between methadone dose, concentrations, and Fridericia rate-corrected QT (QTcF) interval prolongation, and to identify a dose that was associated with increased risk of developing torsade de pointes. A linear relationship between concentration and QTcF described the data from five clinical trials in patients on methadone maintenance treatment (MMT). A previously published population PK model adequately described the concentration-time data, and this model was used for simulation. QTcF was increased by a mean (90% confidence interval (CI)) of 17 (12, 22) ms per 1,000 ng/ml of methadone. Based on this model, doses >120 mg/day would increase the QTcF interval by >20 ms. The model predicts that 1-3% of patients would have ΔQTcF >60 ms, and 0.3-2.0% of patients would have QTcF >500 ms at doses of 160-200 mg/day. Our predictions are consistent with available observational data and support the need for electrocardiogram (ECG) monitoring and arrhythmia risk factor assessment in patients receiving methadone doses >120 mg/day.

A neurocomputational account of cognitive deficits in Parkinson’s disease

PubMed Central

Hélie, Sébastien; Paul, Erick J.; Ashby, F. Gregory

2014-01-01

Parkinson’s disease (PD) is caused by the accelerated death of dopamine (DA) producing neurons. Numerous studies documenting cognitive deficits of PD patients have revealed impairments in a variety of tasks related to memory, learning, visuospatial skills, and attention. While there have been several studies documenting cognitive deficits of PD patients, very few computational models have been proposed. In this article, we use the COVIS model of category learning to simulate DA depletion and show that the model suffers from cognitive symptoms similar to those of human participants affected by PD. Specifically, DA depletion in COVIS produced deficits in rule-based categorization, non-linear information-integration categorization, probabilistic classification, rule maintenance, and rule switching. These were observed by simulating results from younger controls, older controls, PD patients, and severe PD patients in five well-known tasks. Differential performance among the different age groups and clinical populations was modeled simply by changing the amount of DA available in the model. This suggests that COVIS may not only be an adequate model of the simulated tasks and phenomena but also more generally of the role of DA in these tasks and phenomena. PMID:22683450

Using micro-simulation to investigate the safety impacts of transit design alternatives at signalized intersections.

PubMed

Li, Lu; Persaud, Bhagwant; Shalaby, Amer

2017-03-01

This study investigates the use of crash prediction models and micro-simulation to develop an effective surrogate safety assessment measure at the intersection level. With the use of these tools, hypothetical scenarios can be developed and explored to evaluate the safety impacts of design alternatives in a controlled environment, in which factors not directly associated with the design alternatives can be fixed. Micro-simulation models are developed, calibrated, and validated. Traffic conflicts in the micro-simulation models are estimated and linked with observed crash frequency, which greatly alleviates the lengthy time needed to collect sufficient crash data for evaluating alternatives, due to the rare and infrequent nature of crash events. A set of generalized linear models with negative binomial error structure is developed to correlate the simulated conflicts with the observed crash frequency in Toronto, Ontario, Canada. Crash prediction models are also developed for crashes of different impact types and for transit-involved crashes. The resulting statistical significance and the goodness-of-fit of the models suggest adequate predictive ability. Based on the established correlation between simulated conflicts and observed crashes, scenarios are developed in the micro-simulation models to investigate the safety effects of individual transit line elements by making hypothetical modifications to such elements and estimating changes in crash frequency from the resulting changes in conflicts. The findings imply that the existing transit signal priority schemes can have a negative effect on safety performance, and that the existing near-side stop positioning and streetcar transit type can be safer at their current state than if they were to be replaced by their respective counterparts. Copyright © 2017 Elsevier Ltd. All rights reserved.

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

NASA Technical Reports Server (NTRS)

Hsu, Andrew T.

1991-01-01

Turbulent combustion can not be simulated adequately by conventional moment closure turbulence models. The difficulty lies in the fact that the reaction rate is in general an exponential function of the temperature, and the higher order correlations in the conventional moment closure models of the chemical source term can not be neglected, making the applications of such models impractical. The probability density function (pdf) method offers an attractive alternative: in a pdf model, the chemical source terms are closed and do not require additional models. A grid dependent Monte Carlo scheme was studied, since it is a logical alternative, wherein the number of computer operations increases only linearly with the increase of number of independent variables, as compared to the exponential increase in a conventional finite difference scheme. A new algorithm was devised that satisfies a restriction in the case of pure diffusion or uniform flow problems. Although for nonuniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

The Influences of Airmass Histories on Radical Species During POLARIS

NASA Technical Reports Server (NTRS)

Pierson, James M.; Kawa, S. R.

1998-01-01

The POLARIS mission focused on understanding the processes associated with the decrease of polar stratospheric ozone from spring to fall at high latitudes. This decrease is linked primarily to in situ photochemical destruction by reactive nitrogen species, NO and NO2, which also control other catalytic loss cycles. Steady state models have been used to test photochemistry and radical behavior but are not always adequate in simulating radical species observations. In some cases, air mass history can be important and trajectory models give an improved simulation of the radical species. Trajectory chemistry models, however, still consistently underestimate NO and NO2 abundances compared to measurements along the ER-2 flight track. The Goddard chemistry on trajectory model has been used to test updated rate constants for NO2 + OH, NO2 + O and OH + HNO3, key reactions that affect NO and NO2 abundances. We present comparisons between the modified Goddard chemistry on trajectory model, the JPL steady state model and observations from selected flights.

Evaluating WRF-Chem multi-scale model in simulating aerosol radiative properties over the tropics – A case study over India

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

Seethala, C.; Pandithurai, G.; Fast, Jerome D.

We utilized WRF-Chem multi-scale model to simulate the regional distribution of aerosols, optical properties and its effect on radiation over India for a winter month. The model is evaluated using measurements obtained from upper-air soundings, AERONET sun photometers, various satellite instruments, and pyranometers operated by the Indian Meteorological Department. The simulated downward shortwave flux was overestimated when the effect of aerosols on radiation and clouds was neglected. Downward shortwave radiation from a simulation that included aerosol-radiation interaction processes was 5 to 25 Wm{sup -2} closer to the observations, while a simulation that included aerosol-cloud interaction processes were another 1 tomore » 20 Wm{sup -2} closer to the observations. For the few observations available, the model usually underestimated particulate concentration. This is likely due to turbulent mixing, transport errors and the lack of secondary organic aerosol treatment in the model. The model efficiently captured the broad regional hotspots such as high aerosol optical depth over Indo-Gangetic basin as well as the northwestern and southern part of India. The regional distribution of aerosol optical depth compares well with AVHRR aerosol optical depth and the TOMS aerosol index. The magnitude and wavelength-dependence of simulated aerosol optical depth was also similar to the AERONET observations across India. Differences in surface shortwave radiation between simulations that included and neglected aerosol-radiation interactions were as high as -25 Wm{sup -2}, while differences in surface shortwave radiation between simulations that included and neglect aerosol-radiation-cloud interactions were as high as -30 Wm{sup -2}. The spatial variations of these differences were also compared with AVHRR observation. This study suggests that the model is able to qualitatively simulate the impact of aerosols on radiation over India; however, additional measurements of particulate mass and composition are needed to fully evaluate whether the aerosol precursor emissions are adequate when simulating radiative forcing in the region.« less

An error criterion for determining sampling rates in closed-loop control systems

NASA Technical Reports Server (NTRS)

Brecher, S. M.

1972-01-01

The determination of an error criterion which will give a sampling rate for adequate performance of linear, time-invariant closed-loop, discrete-data control systems was studied. The proper modelling of the closed-loop control system for characterization of the error behavior, and the determination of an absolute error definition for performance of the two commonly used holding devices are discussed. The definition of an adequate relative error criterion as a function of the sampling rate and the parameters characterizing the system is established along with the determination of sampling rates. The validity of the expressions for the sampling interval was confirmed by computer simulations. Their application solves the problem of making a first choice in the selection of sampling rates.

LSD (Landing System Development) Impact Simulation

NASA Astrophysics Data System (ADS)

Ullio, R.; Riva, N.; Pellegrino, P.; Deloo, P.

2012-07-01

In the frame of the Exploration Programs, a soft landing on the planet surface is foreseen. To ensure a successful final landing phase, a landing system by using leg tripod design landing legs with adequate crushable damping system was selected, capable of absorbing the residual velocities (vertical, horizontal and angular) at touch- down, insuring stability. TAS-I developed a numerical non linear dynamic methodology for the landing impact simulation of the Lander system by using a commercial explicit finite element analysis code (i.e. Altair RADIOSS). In this paper the most significant FE modeling approaches and results of the analytical simulation of landing impact are reported, especially with respect to the definition of leg dimensioning loads and the design update of selected parts (if necessary).

A spatial simulation model of hydrology and vegetation dynamics in semi-permanent prairie wetlands

USGS Publications Warehouse

Poiani, Karen A.; Johnson, W. Carter

1993-01-01

The objective of this study was to construct a spatial simulation model of the vegetation dynamics in semi-permanent prairie wetlands. A hydrologic submodel estimated water levels based on precipitation, runoff, and potential evapotranspiration. A vegetation submodel calculated the amount and distribution of emergent cover and open water using a geographic information system. The response of vegetation to water-level changes was based on seed bank composition, seedling recruitment and establishment, and plant survivorship. The model was developed and tested using data from the Cottonwood Lake study site in North Dakota. Data from semi-permanent wetland P1 were used to calibrate the model. Data from a second wetland, P4, were used to evaluate model performance. Simulation results were compared with actual water data from 1797 through 1989. Test results showed that differences between calculated and observed water levels were within 10 cm 75% of the time. Open water over the past decade ranged from 0 to 7% in wetland P4 and from 0 to 8% in submodel simulations. Several model parameters including evapotranspiration and timing of seedling germination could be improved with more complex techniques or relatively minor adjustments. Despite these differences the model adequately represented vegetation dynamics of prairie wetlands and can be used to examine wetland response to natural or human-induced climate change.

Evaluation of Aztreonam Dosing Regimens in Patients With Normal and Impaired Renal Function: A Population Pharmacokinetic Modeling and Monte Carlo Simulation Analysis.

PubMed

Xu, Hongmei; Zhou, Wangda; Zhou, Diansong; Li, Jianguo; Al-Huniti, Nidal

2017-03-01

Aztreonam is a monocyclic β-lactam antibiotic often used to treat infections caused by Enterobacteriaceae or Pseudomonas aeruginosa. Despite the long history of clinical use, population pharmacokinetic modeling of aztreonam in renally impaired patients is not yet available. The aims of this study were to assess the impact of renal impairment on aztreonam exposure and to evaluate dosing regimens for patients with renal impairment. A population model describing aztreonam pharmacokinetics following intravenous administration was developed using plasma concentrations from 42 healthy volunteers and renally impaired patients from 2 clinical studies. The final pharmacokinetic model was used to predict aztreonam plasma concentrations and evaluate the probability of pharmacodynamic target attainment (PTA) in patients with different levels of renal function. A 2-compartment model with first-order elimination adequately described aztreonam pharmacokinetics. The population mean estimates of aztreonam clearance, intercompartmental clearance, volume of distribution of the central compartment, and volume of distribution of the peripheral compartment were 4.93 L/h, 9.26 L/h, 7.43 L, and 6.44 L, respectively. Creatinine clearance and body weight were the most significant variables to explain patient variability in aztreonam clearance and volume of distribution, respectively. Simulations using the final pharmacokinetic model resulted in a clinical susceptibility break point of 4 and 8 mg/L, respectively, based on the clinical use of 1- and 2-g loading doses with the same or reduced maintenance dose every 8 hours for various renal deficiency patients. The population pharmacokinetic modeling and PTA estimation support adequate PTAs (>90% PTA) from the aztreonam label for dose adjustment of aztreonam in patients with moderate and severe renal impairment. © 2016, The American College of Clinical Pharmacology.

Low-speed aerodynamic characteristics of a 0.08-scale YF-17 airplane model at high angles of attack and sideslip

NASA Technical Reports Server (NTRS)

Petroff, D. N.; Scher, S. H.; Sutton, C. E.

1978-01-01

Data were obtained with and without the nose boom and with several strake configurations; also, data were obtained for various control surface deflections. Analysis of the results revealed that selected strake configurations adequately provided low Reynolds number simulation of the high Reynolds number characteristics. The addition of the boom in general tended to reduce the Reynolds number effects.

Guidelines for Calibration and Application of Storm.

DTIC Science & Technology

1977-12-01

combination method uses the SCS method on pervious areas and the coefficient method on impervious areas of the watershed. Storm water quality is computed...stations, it should be accomplished according to procedures outlined In Reference 7. Adequate storm water quality data are the most difficult and costly...mass discharge of pollutants is negligible. The state-of-the-art in urban storm water quality modeling precludes highly accurate simulation of

Uncertainty based modeling of rainfall-runoff: Combined differential evolution adaptive Metropolis (DREAM) and K-means clustering

NASA Astrophysics Data System (ADS)

Zahmatkesh, Zahra; Karamouz, Mohammad; Nazif, Sara

2015-09-01

Simulation of rainfall-runoff process in urban areas is of great importance considering the consequences and damages of extreme runoff events and floods. The first issue in flood hazard analysis is rainfall simulation. Large scale climate signals have been proved to be effective in rainfall simulation and prediction. In this study, an integrated scheme is developed for rainfall-runoff modeling considering different sources of uncertainty. This scheme includes three main steps of rainfall forecasting, rainfall-runoff simulation and future runoff prediction. In the first step, data driven models are developed and used to forecast rainfall using large scale climate signals as rainfall predictors. Due to high effect of different sources of uncertainty on the output of hydrologic models, in the second step uncertainty associated with input data, model parameters and model structure is incorporated in rainfall-runoff modeling and simulation. Three rainfall-runoff simulation models are developed for consideration of model conceptual (structural) uncertainty in real time runoff forecasting. To analyze the uncertainty of the model structure, streamflows generated by alternative rainfall-runoff models are combined, through developing a weighting method based on K-means clustering. Model parameters and input uncertainty are investigated using an adaptive Markov Chain Monte Carlo method. Finally, calibrated rainfall-runoff models are driven using the forecasted rainfall to predict future runoff for the watershed. The proposed scheme is employed in the case study of the Bronx River watershed, New York City. Results of uncertainty analysis of rainfall-runoff modeling reveal that simultaneous estimation of model parameters and input uncertainty significantly changes the probability distribution of the model parameters. It is also observed that by combining the outputs of the hydrological models using the proposed clustering scheme, the accuracy of runoff simulation in the watershed is remarkably improved up to 50% in comparison to the simulations by the individual models. Results indicate that the developed methodology not only provides reliable tools for rainfall and runoff modeling, but also adequate time for incorporating required mitigation measures in dealing with potentially extreme runoff events and flood hazard. Results of this study can be used in identification of the main factors affecting flood hazard analysis.

Numerical simulations of clinical focused ultrasound functional neurosurgery

NASA Astrophysics Data System (ADS)

Pulkkinen, Aki; Werner, Beat; Martin, Ernst; Hynynen, Kullervo

2014-04-01

A computational model utilizing grid and finite difference methods were developed to simulate focused ultrasound functional neurosurgery interventions. The model couples the propagation of ultrasound in fluids (soft tissues) and solids (skull) with acoustic and visco-elastic wave equations. The computational model was applied to simulate clinical focused ultrasound functional neurosurgery treatments performed in patients suffering from therapy resistant chronic neuropathic pain. Datasets of five patients were used to derive the treatment geometry. Eight sonications performed in the treatments were then simulated with the developed model. Computations were performed by driving the simulated phased array ultrasound transducer with the acoustic parameters used in the treatments. Resulting focal temperatures and size of the thermal foci were compared quantitatively, in addition to qualitative inspection of the simulated pressure and temperature fields. This study found that the computational model and the simulation parameters predicted an average of 24 ± 13% lower focal temperature elevations than observed in the treatments. The size of the simulated thermal focus was found to be 40 ± 13% smaller in the anterior-posterior direction and 22 ± 14% smaller in the inferior-superior direction than in the treatments. The location of the simulated thermal focus was off from the prescribed target by 0.3 ± 0.1 mm, while the peak focal temperature elevation observed in the measurements was off by 1.6 ± 0.6 mm. Although the results of the simulations suggest that there could be some inaccuracies in either the tissue parameters used, or in the simulation methods, the simulations were able to predict the focal spot locations and temperature elevations adequately for initial treatment planning performed to assess, for example, the feasibility of sonication. The accuracy of the simulations could be improved if more precise ultrasound tissue properties (especially of the skull bone) could be obtained.

Evaluation of rainfall simulations over West Africa in dynamically downscaled CMIP5 global circulation models

NASA Astrophysics Data System (ADS)

Akinsanola, A. A.; Ajayi, V. O.; Adejare, A. T.; Adeyeri, O. E.; Gbode, I. E.; Ogunjobi, K. O.; Nikulin, G.; Abolude, A. T.

2018-04-01

This study presents evaluation of the ability of Rossby Centre Regional Climate Model (RCA4) driven by nine global circulation models (GCMs), to skilfully reproduce the key features of rainfall climatology over West Africa for the period of 1980-2005. The seasonal climatology and annual cycle of the RCA4 simulations were assessed over three homogenous subregions of West Africa (Guinea coast, Savannah, and Sahel) and evaluated using observed precipitation data from the Global Precipitation Climatology Project (GPCP). Furthermore, the model output was evaluated using a wide range of statistical measures. The interseasonal and interannual variability of the RCA4 were further assessed over the subregions and the whole of the West Africa domain. Results indicate that the RCA4 captures the spatial and interseasonal rainfall pattern adequately but exhibits a weak performance over the Guinea coast. Findings from the interannual rainfall variability indicate that the model performance is better over the larger West Africa domain than the subregions. The largest difference across the RCA4 simulated annual rainfall was found in the Sahel. Result from the Mann-Kendall test showed no significant trend for the 1980-2005 period in annual rainfall either in GPCP observation data or in the model simulations over West Africa. In many aspects, the RCA4 simulation driven by the HadGEM2-ES perform best over the region. The use of the multimodel ensemble mean has resulted to the improved representation of rainfall characteristics over the study domain.

Testing the ability of a semidistributed hydrological model to simulate contributing area

NASA Astrophysics Data System (ADS)

Mengistu, S. G.; Spence, C.

2016-06-01

A dry climate, the prevalence of small depressions, and the lack of a well-developed drainage network are characteristics of environments with extremely variable contributing areas to runoff. These types of regions arguably present the greatest challenge to properly understanding catchment streamflow generation processes. Previous studies have shown that contributing area dynamics are important for streamflow response, but the nature of the relationship between the two is not typically understood. Furthermore, it is not often tested how well hydrological models simulate contributing area. In this study, the ability of a semidistributed hydrological model, the PDMROF configuration of Environment Canada's MESH model, was tested to determine if it could simulate contributing area. The study focused on the St. Denis Creek watershed in central Saskatchewan, Canada, which with its considerable topographic depressions, exhibits wide variation in contributing area, making it ideal for this type of investigation. MESH-PDMROF was able to replicate contributing area derived independently from satellite imagery. Daily model simulations revealed a hysteretic relationship between contributing area and streamflow not apparent from the less frequent remote sensing observations. This exercise revealed that contributing area extent can be simulated by a semi-distributed hydrological model with a scheme that assumes storage capacity distribution can be represented with a probability function. However, further investigation is needed to determine if it can adequately represent the complex relationship between streamflow and contributing area that is such a key signature of catchment behavior.

Decision method for optimal selection of warehouse material handling strategies by production companies

NASA Astrophysics Data System (ADS)

Dobos, P.; Tamás, P.; Illés, B.

2016-11-01

Adequate establishment and operation of warehouse logistics determines the companies’ competitiveness significantly because it effects greatly the quality and the selling price of the goods that the production companies produce. In order to implement and manage an adequate warehouse system, adequate warehouse position, stock management model, warehouse technology, motivated work force committed to process improvement and material handling strategy are necessary. In practical life, companies have paid small attantion to select the warehouse strategy properly. Although it has a major influence on the production in the case of material warehouse and on smooth costumer service in the case of finished goods warehouse because this can happen with a huge loss in material handling. Due to the dynamically changing production structure, frequent reorganization of warehouse activities is needed, on what the majority of the companies react basically with no reactions. This work presents a simulation test system frames for eligible warehouse material handling strategy selection and also the decision method for selection.

The impact of mesoscale convective systems on global precipitation: A modeling study

NASA Astrophysics Data System (ADS)

Tao, Wei-Kuo

2017-04-01

The importance of precipitating mesoscale convective systems (MCSs) has been quantified from TRMM precipitation radar and microwave imager retrievals. MCSs generate more than 50% of the rainfall in most tropical regions. Typical MCSs have horizontal scales of a few hundred kilometers (km); therefore, a large domain and high resolution are required for realistic simulations of MCSs in cloud-resolving models (CRMs). Almost all traditional global and climate models do not have adequate parameterizations to represent MCSs. Typical multi-scale modeling frameworks (MMFs) with 32 CRM grid points and 4 km grid spacing also might not have sufficient resolution and domain size for realistically simulating MCSs. In this study, the impact of MCSs on precipitation processes is examined by conducting numerical model simulations using the Goddard Cumulus Ensemble model (GCE) and Goddard MMF (GMMF). The results indicate that both models can realistically simulate MCSs with more grid points (i.e., 128 and 256) and higher resolutions (1 or 2 km) compared to those simulations with less grid points (i.e., 32 and 64) and low resolution (4 km). The modeling results also show that the strengths of the Hadley circulations, mean zonal and regional vertical velocities, surface evaporation, and amount of surface rainfall are either weaker or reduced in the GMMF when using more CRM grid points and higher CRM resolution. In addition, the results indicate that large-scale surface evaporation and wind feed back are key processes for determining the surface rainfall amount in the GMMF. A sensitivity test with reduced sea surface temperatures (SSTs) is conducted and results in both reduced surface rainfall and evaporation.

A one-dimensional with three-dimensional velocity space hybrid-PIC model of the discharge plasma in a Hall thruster

NASA Astrophysics Data System (ADS)

Shashkov, Andrey; Lovtsov, Alexander; Tomilin, Dmitry

2017-04-01

According to present knowledge, countless numerical simulations of the discharge plasma in Hall thrusters were conducted. However, on the one hand, adequate two-dimensional (2D) models require a lot of time to carry out numerical research of the breathing mode oscillations or the discharge structure. On the other hand, existing one-dimensional (1D) models are usually too simplistic and do not take into consideration such important phenomena as neutral-wall collisions, magnetic field induced by Hall current and double, secondary, and stepwise ionizations together. In this paper a one-dimensional with three-dimensional velocity space (1D3V) hybrid-PIC model is presented. The model is able to incorporate all the phenomena mentioned above. A new method of neutral-wall collisions simulation in described space was developed and validated. Simulation results obtained for KM-88 and KM-60 thrusters are in a good agreement with experimental data. The Bohm collision coefficient was the same for both thrusters. Neutral-wall collisions, doubly charged ions, and induced magnetic field were proved to stabilize the breathing mode oscillations in a Hall thruster under some circumstances.

Electro-osmotic flow of a model electrolyte

NASA Astrophysics Data System (ADS)

Zhu, Wei; Singer, Sherwin J.; Zheng, Zhi; Conlisk, A. T.

2005-04-01

Electro-osmotic flow is studied by nonequilibrium molecular dynamics simulations in a model system chosen to elucidate various factors affecting the velocity profile and facilitate comparison with existing continuum theories. The model system consists of spherical ions and solvent, with stationary, uniformly charged walls that make a channel with a height of 20 particle diameters. We find that hydrodynamic theory adequately describes simple pressure-driven (Poiseuille) flow in this model. However, Poisson-Boltzmann theory fails to describe the ion distribution in important situations, and therefore continuum fluid dynamics based on the Poisson-Boltzmann ion distribution disagrees with simulation results in those situations. The failure of Poisson-Boltzmann theory is traced to the exclusion of ions near the channel walls resulting from reduced solvation of the ions in that region. When a corrected ion distribution is used as input for hydrodynamic theory, agreement with numerical simulations is restored. An analytic theory is presented that demonstrates that repulsion of the ions from the channel walls increases the flow rate, and attraction to the walls has the opposite effect. A recent numerical study of electro-osmotic flow is reanalyzed in the light of our findings, and the results conform well to our conclusions for the model system.

Surface diffusion of cyclic hydrocarbons on nickel

NASA Astrophysics Data System (ADS)

Silverwood, I. P.; Armstrong, J.

2018-08-01

Surface diffusion of adsorbates is difficult to measure on realistic systems, yet it is of fundamental interest in catalysis and coating reactions. quasielastic neutron scattering (QENS) was used to investigate the diffusion of cyclohexane and benzene adsorbed on a nickel metal sponge catalyst. Molecular dynamics simulations of benzene on a model (111) nickel surface showed localised motion with diffusion by intermittent jumps. The experimental data was therefore fitted to the Singwi-Sjölander model and activation energies for diffusion of 4.0 kJ mol-1 for benzene and 4.3 kJ mol-1 for cyclohexane were calculated for the two dimensional model. Limited motion out-of plane was seen in the dynamics simulations and is discussed, although the resolution of the scattering experiment is insufficient to quantify this. Good agreement is seen between the use of a perfect crystal as a model for a disordered system over short time scales, suggesting that simple models are adequate to describe diffusion over polycrystalline metal surfaces on the timescale of QENS measurement.

Discrete particle model for cement infiltration within open-cell structures: Prevention of osteoporotic fracture.

PubMed

Ramos-Infante, Samuel Jesús; Ten-Esteve, Amadeo; Alberich-Bayarri, Angel; Pérez, María Angeles

2018-01-01

This paper proposes a discrete particle model based on the random-walk theory for simulating cement infiltration within open-cell structures to prevent osteoporotic proximal femur fractures. Model parameters consider the cement viscosity (high and low) and the desired direction of injection (vertical and diagonal). In vitro and in silico characterizations of augmented open-cell structures validated the computational model and quantified the improved mechanical properties (Young's modulus) of the augmented specimens. The cement injection pattern was successfully predicted in all the simulated cases. All the augmented specimens exhibited enhanced mechanical properties computationally and experimentally (maximum improvements of 237.95 ± 12.91% and 246.85 ± 35.57%, respectively). The open-cell structures with high porosity fraction showed a considerable increase in mechanical properties. Cement augmentation in low porosity fraction specimens resulted in a lesser increase in mechanical properties. The results suggest that the proposed discrete particle model is adequate for use as a femoroplasty planning framework.

Simulation of the inhibition of microbial sulfate reduction in a two-compartment upflow bioreactor subjected to molybdate injection.

PubMed

de Jesus, E B; de Andrade Lima, L R P

2016-08-01

Souring of oil fields during secondary oil recovery by water injection occurs mainly due to the action of sulfate-reducing bacteria (SRB) adhered to the rock surface in the vicinity of injection wells. Upflow packed-bed bioreactors have been used in petroleum microbiology because of its similarity to the oil field near the injection wells or production. However, these reactors do not realistically describe the regions near the injection wells, which are characterized by the presence of a saturated zone and a void region close to the well. In this study, the hydrodynamics of the two-compartment packing-free/packed-bed pilot bioreactor that mimics an oil reservoir was studied. The packed-free compartment was modeled using a continuous stirred tank model with mass exchange between active and stagnant zones, whereas the packed-bed compartment was modeled using a piston-dispersion-exchange model. The proposed model adequately represents the hydrodynamic of the packed-free/packed-bed bioreactor while the simulations provide important information about the characteristics of the residence time distribution (RTD) curves for different sets of model parameters. Simulations were performed to represent the control of the sulfate-reducing bacteria activity in the bioreactor with the use of molybdate in different scenarios. The simulations show that increased amounts of molybdate cause an effective inhibition of the souring sulfate-reducing bacteria activity.

Modelling and simulation of [18F]fluoromisonidazole dynamics based on histology-derived microvessel maps

NASA Astrophysics Data System (ADS)

Mönnich, David; Troost, Esther G. C.; Kaanders, Johannes H. A. M.; Oyen, Wim J. G.; Alber, Markus; Thorwarth, Daniela

2011-04-01

Hypoxia can be assessed non-invasively by positron emission tomography (PET) using radiotracers such as [18F]fluoromisonidazole (Fmiso) accumulating in poorly oxygenated cells. Typical features of dynamic Fmiso PET data are high signal variability in the first hour after tracer administration and slow formation of a consistent contrast. The purpose of this study is to investigate whether these characteristics can be explained by the current conception of the underlying microscopic processes and to identify fundamental effects. This is achieved by modelling and simulating tissue oxygenation and tracer dynamics on the microscopic scale. In simulations, vessel structures on histology-derived maps act as sources and sinks for oxygen as well as tracer molecules. Molecular distributions in the extravascular space are determined by reaction-diffusion equations, which are solved numerically using a two-dimensional finite element method. Simulated Fmiso time activity curves (TACs), though not directly comparable to PET TACs, reproduce major characteristics of clinical curves, indicating that the microscopic model and the parameter values are adequate. Evidence for dependence of the early PET signal on the vascular fraction is found. Further, possible effects leading to late contrast formation and potential implications on the quantification of Fmiso PET data are discussed.

Calibration of discrete element model parameters: soybeans

NASA Astrophysics Data System (ADS)

Ghodki, Bhupendra M.; Patel, Manish; Namdeo, Rohit; Carpenter, Gopal

2018-05-01

Discrete element method (DEM) simulations are broadly used to get an insight of flow characteristics of granular materials in complex particulate systems. DEM input parameters for a model are the critical prerequisite for an efficient simulation. Thus, the present investigation aims to determine DEM input parameters for Hertz-Mindlin model using soybeans as a granular material. To achieve this aim, widely acceptable calibration approach was used having standard box-type apparatus. Further, qualitative and quantitative findings such as particle profile, height of kernels retaining the acrylic wall, and angle of repose of experiments and numerical simulations were compared to get the parameters. The calibrated set of DEM input parameters includes the following (a) material properties: particle geometric mean diameter (6.24 mm); spherical shape; particle density (1220 kg m^{-3} ), and (b) interaction parameters such as particle-particle: coefficient of restitution (0.17); coefficient of static friction (0.26); coefficient of rolling friction (0.08), and particle-wall: coefficient of restitution (0.35); coefficient of static friction (0.30); coefficient of rolling friction (0.08). The results may adequately be used to simulate particle scale mechanics (grain commingling, flow/motion, forces, etc) of soybeans in post-harvest machinery and devices.

Estimation of pollutant loads considering dam operation in Han River Basin by BASINS/Hydrological Simulation Program-FORTRAN.

PubMed

Jung, Kwang-Wook; Yoon, Choon-G; Jang, Jae-Ho; Kong, Dong-Soo

2008-01-01

Effective watershed management often demands qualitative and quantitative predictions of the effect of future management activities as arguments for policy makers and administration. The BASINS geographic information system was developed to compute total maximum daily loads, which are helpful to establish hydrological process and water quality modeling system. In this paper the BASINS toolkit HSPF model is applied in 20,271 km(2) large watershed of the Han River Basin is used for applicability of HSPF and BMPs scenarios. For proper evaluation of watershed and stream water quality, comprehensive estimation methods are necessary to assess large amounts of point source and nonpoint-source (NPS) pollution based on the total watershed area. In this study, The Hydrological Simulation Program-FORTRAN (HSPF) was estimated to simulate watershed pollutant loads containing dam operation and applied BMPs scenarios for control NPS pollution. The 8-day monitoring data (about three years) were used in the calibration and verification processes. Model performance was in the range of "very good" and "good" based on percent difference. The water-quality simulation results were encouraging for this large sizable watershed with dam operation practice and mixed land uses; HSPF proved adequate, and its application is recommended to simulate watershed processes and BMPs evaluation. IWA Publishing 2008.

From individual to population level effects of toxicants in the tubicifid Branchiura sowerbyi using threshold effect models in a Bayesian framework.

PubMed

Ducrot, Virginie; Billoir, Elise; Péry, Alexandre R R; Garric, Jeanne; Charles, Sandrine

2010-05-01

Effects of zinc were studied in the freshwater worm Branchiura sowerbyi using partial and full life-cycle tests. Only newborn and juveniles were sensitive to zinc, displaying effects on survival, growth, and age at first brood at environmentally relevant concentrations. Threshold effect models were proposed to assess toxic effects on individuals. They were fitted to life-cycle test data using Bayesian inference and adequately described life-history trait data in exposed organisms. The daily asymptotic growth rate of theoretical populations was then simulated with a matrix population model, based upon individual-level outputs. Population-level outputs were in accordance with existing literature for controls. Working in a Bayesian framework allowed incorporating parameter uncertainty in the simulation of the population-level response to zinc exposure, thus increasing the relevance of test results in the context of ecological risk assessment.

A study of deoxyribonucleotide metabolism and its relation to DNA synthesis. Supercomputer simulation and model-system analysis.

PubMed

Heinmets, F; Leary, R H

1991-06-01

A model system (1) was established to analyze purine and pyrimidine metabolism. This system has been expanded to include macrosimulation of DNA synthesis and the study of its regulation by terminal deoxynucleoside triphosphates (dNTPs) via a complex set of interactions. Computer experiments reveal that our model exhibits adequate and reasonable sensitivity in terms of dNTP pool levels and rates of DNA synthesis when inputs to the system are varied. These simulation experiments reveal that in order to achieve maximum DNA synthesis (in terms of purine metabolism), a proper balance is required in guanine and adenine input into this metabolic system. Excessive inputs will become inhibitory to DNA synthesis. In addition, studies are carried out on rates of DNA synthesis when various parameters are changed quantitatively. The current system is formulated by 110 differential equations.

Frequency-difference MIT imaging of cerebral haemorrhage with a hemispherical coil array: numerical modelling.

PubMed

Zolgharni, M; Griffiths, H; Ledger, P D

2010-08-01

The feasibility of detecting a cerebral haemorrhage with a hemispherical MIT coil array consisting of 56 exciter/sensor coils of 10 mm radius and operating at 1 and 10 MHz was investigated. A finite difference method combined with an anatomically realistic head model comprising 12 tissue types was used to simulate the strokes. Frequency-difference images were reconstructed from the modelled data with different levels of the added phase noise and two types of a priori boundary errors: a displacement of the head and a size scaling error. The results revealed that a noise level of 3 m degrees (standard deviation) was adequate for obtaining good visualization of a peripheral stroke (volume approximately 49 ml). The simulations further showed that the displacement error had to be within 3-4 mm and the scaling error within 3-4% so as not to cause unacceptably large artefacts on the images.

Simulation of high temperature thermal energy storage system based on coupled metal hydrides for solar driven steam power plants

DOE PAGES

d'Entremont, Anna; Corgnale, Claudio; Hardy, Bruce; ...

2018-01-11

Concentrating solar power plants can achieve low cost and efficient renewable electricity production if equipped with adequate thermal energy storage systems. Metal hydride based thermal energy storage systems are appealing candidates due to their demonstrated potential for very high volumetric energy densities, high exergetic efficiencies, and low costs. The feasibility and performance of a thermal energy storage system based on NaMgH 2F hydride paired with TiCr 1.6Mn 0.2 is examined, discussing its integration with a solar-driven ultra-supercritical steam power plant. The simulated storage system is based on a laboratory-scale experimental apparatus. It is analyzed using a detailed transport model accountingmore » for the thermochemical hydrogen absorption and desorption reactions, including kinetics expressions adequate for the current metal hydride system. The results show that the proposed metal hydride pair can suitably be integrated with a high temperature steam power plant. The thermal energy storage system achieves output energy densities of 226 kWh/m 3, 9 times the DOE SunShot target, with moderate temperature and pressure swings. Also, simulations indicate that there is significant scope for performance improvement via heat-transfer enhancement strategies.« less

Simulation of high temperature thermal energy storage system based on coupled metal hydrides for solar driven steam power plants

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

d'Entremont, Anna; Corgnale, Claudio; Hardy, Bruce

Concentrating solar power plants can achieve low cost and efficient renewable electricity production if equipped with adequate thermal energy storage systems. Metal hydride based thermal energy storage systems are appealing candidates due to their demonstrated potential for very high volumetric energy densities, high exergetic efficiencies, and low costs. The feasibility and performance of a thermal energy storage system based on NaMgH 2F hydride paired with TiCr 1.6Mn 0.2 is examined, discussing its integration with a solar-driven ultra-supercritical steam power plant. The simulated storage system is based on a laboratory-scale experimental apparatus. It is analyzed using a detailed transport model accountingmore » for the thermochemical hydrogen absorption and desorption reactions, including kinetics expressions adequate for the current metal hydride system. The results show that the proposed metal hydride pair can suitably be integrated with a high temperature steam power plant. The thermal energy storage system achieves output energy densities of 226 kWh/m 3, 9 times the DOE SunShot target, with moderate temperature and pressure swings. Also, simulations indicate that there is significant scope for performance improvement via heat-transfer enhancement strategies.« less

Comparison of PARASOL Observations with Polarized Reflectances Simulated Using Different Ice Habit Mixtures

NASA Technical Reports Server (NTRS)

Cole, Benjamin H.; Yang, Ping; Baum, Bryan A.; Riedi, Jerome; Labonnote, Laurent C.; Thieuleux, Francois; Platnick, Steven

2012-01-01

Insufficient knowledge of the habit distribution and the degree of surface roughness of ice crystals within ice clouds is a source of uncertainty in the forward light scattering and radiative transfer simulations required in downstream applications involving these clouds. The widely used MODerate Resolution Imaging Spectroradiometer (MODIS) Collection 5 ice microphysical model assumes a mixture of various ice crystal shapes with smooth-facets except aggregates of columns for which a moderately rough condition is assumed. When compared with PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) polarized reflection data, simulations of polarized reflectance using smooth particles show a poor fit to the measurements, whereas very rough-faceted particles provide an improved fit to the polarized reflectance. In this study a new microphysical model based on a mixture of 9 different ice crystal habits with severely roughened facets is developed. Simulated polarized reflectance using the new ice habit distribution is calculated using a vector adding-doubling radiative transfer model, and the simulations closely agree with the polarized reflectance observed by PARASOL. The new general habit mixture is also tested using a spherical albedo differences analysis, and surface roughening is found to improve the consistency of multi-angular observations. It is suggested that an ice model incorporating an ensemble of different habits with severely roughened surfaces would potentially be an adequate choice for global ice cloud retrievals.

Spacecraft Jitter Attenuation Using Embedded Piezoelectric Actuators

NASA Technical Reports Server (NTRS)

Belvin, W. Keith

1995-01-01

Remote sensing from spacecraft requires precise pointing of measurement devices in order to achieve adequate spatial resolution. Unfortunately, various spacecraft disturbances induce vibrational jitter in the remote sensing instruments. The NASA Langley Research Center has performed analysis, simulations, and ground tests to identify the more promising technologies for minimizing spacecraft pointing jitter. These studies have shown that the use of smart materials to reduce spacecraft jitter is an excellent match between a maturing technology and an operational need. This paper describes the use of embedding piezoelectric actuators for vibration control and payload isolation. In addition, recent advances in modeling, simulation, and testing of spacecraft pointing jitter are discussed.

Dynamic Deployment Simulations of Inflatable Space Structures

NASA Technical Reports Server (NTRS)

Wang, John T.

2005-01-01

The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.

Thermal Structure Analysis of SIRCA Tile for X-34 Wing Leading Edge TPS

NASA Technical Reports Server (NTRS)

Milos, Frank S.; Squire, Thomas H.; Rasky, Daniel J. (Technical Monitor)

1997-01-01

This paper will describe in detail thermal/structural analyses of SIRCA tiles which were performed at NASA Ames under the The Tile Analysis Task of the X-34 Program. The analyses used the COSMOS/M finite element software to simulate the material response in arc-jet tests, mechanical deflection tests, and the performance of candidate designs for the TPS system. Purposes of the analysis were to verify thermal and structural models for the SIRCA tiles, to establish failure criteria for stressed tiles, to simulate the TPS response under flight aerothermal and mechanical load, and to confirm that adequate safety margins exist for the actual TPS design.

Modeling of detachment experiments at DIII-D

DOE PAGES

Canik, John M.; Briesemeister, Alexis R.; Lasnier, C. J.; ...

2014-11-26

Edge fluid–plasma/kinetic–neutral modeling of well-diagnosed DIII-D experiments is performed in order to document in detail how well certain aspects of experimental measurements are reproduced within the model as the transition to detachment is approached. Results indicate, that at high densities near detachment onset, the poloidal temperature profile produced in the simulations agrees well with that measured in experiment. However, matching the heat flux in the model requires a significant increase in the radiated power compared to what is predicted using standard chemical sputtering rates. Lastly, these results suggest that the model is adequate to predict the divertor temperature, provided thatmore » the discrepancy in radiated power level can be resolved.« less

Modelling the morphology of migrating bacterial colonies

NASA Astrophysics Data System (ADS)

Nishiyama, A.; Tokihiro, T.; Badoual, M.; Grammaticos, B.

2010-08-01

We present a model which aims at describing the morphology of colonies of Proteus mirabilis and Bacillus subtilis. Our model is based on a cellular automaton which is obtained by the adequate discretisation of a diffusion-like equation, describing the migration of the bacteria, to which we have added rules simulating the consolidation process. Our basic assumption, following the findings of the group of Chuo University, is that the migration and consolidation processes are controlled by the local density of the bacteria. We show that it is possible within our model to reproduce the morphological diagrams of both bacteria species. Moreover, we model some detailed experiments done by the Chuo University group, obtaining a fine agreement.

Numerical model for the uptake of groundwater contaminants by phreatophytes

USGS Publications Warehouse

Widdowson, M.A.; El-Sayed, A.; Landmeyer, J.E.

2008-01-01

Conventional solute transport models do not adequately account for the effects of phreatophytic plant systems on contaminant concentrations in shallow groundwater systems. A numerical model was developed and tested to simulate threedimensional reactive solute transport in a heterogeneous porous medium. Advective-dispersive transport is coupled to biodegradation, sorption, and plantbased attenuation processes including plant uptake and sorption by plant roots. The latter effects are a function of the physical-chemical properties of the individual solutes and plant species. Models for plant uptake were tested and evaluated using the experimental data collected at a field site comprised of hybrid poplar trees. A non-linear equilibrium isotherm model best represented site conditions.

Exploring the large-scale structure of Taylor–Couette turbulence through Large-Eddy Simulations

NASA Astrophysics Data System (ADS)

Ostilla-Mónico, Rodolfo; Zhu, Xiaojue; Verzicco, Roberto

2018-04-01

Large eddy simulations (LES) of Taylor-Couette (TC) flow, the flow between two co-axial and independently rotating cylinders are performed in an attempt to explore the large-scale axially-pinned structures seen in experiments and simulations. Both static and dynamic LES models are used. The Reynolds number is kept fixed at Re = 3.4 · 104, and the radius ratio η = ri /ro is set to η = 0.909, limiting the effects of curvature and resulting in frictional Reynolds numbers of around Re τ ≈ 500. Four rotation ratios from Rot = ‑0.0909 to Rot = 0.3 are simulated. First, the LES of TC is benchmarked for different rotation ratios. Both the Smagorinsky model with a constant of cs = 0.1 and the dynamic model are found to produce reasonable results for no mean rotation and cyclonic rotation, but deviations increase for increasing rotation. This is attributed to the increasing anisotropic character of the fluctuations. Second, “over-damped” LES, i.e. LES with a large Smagorinsky constant is performed and is shown to reproduce some features of the large-scale structures, even when the near-wall region is not adequately modeled. This shows the potential for using over-damped LES for fast explorations of the parameter space where large-scale structures are found.

Comprehensive Modeling and Visualization of Cardiac Anatomy and Physiology from CT Imaging and Computer Simulations

PubMed Central

Sun, Peng; Zhou, Haoyin; Ha, Seongmin; Hartaigh, Bríain ó; Truong, Quynh A.; Min, James K.

2016-01-01

In clinical cardiology, both anatomy and physiology are needed to diagnose cardiac pathologies. CT imaging and computer simulations provide valuable and complementary data for this purpose. However, it remains challenging to gain useful information from the large amount of high-dimensional diverse data. The current tools are not adequately integrated to visualize anatomic and physiologic data from a complete yet focused perspective. We introduce a new computer-aided diagnosis framework, which allows for comprehensive modeling and visualization of cardiac anatomy and physiology from CT imaging data and computer simulations, with a primary focus on ischemic heart disease. The following visual information is presented: (1) Anatomy from CT imaging: geometric modeling and visualization of cardiac anatomy, including four heart chambers, left and right ventricular outflow tracts, and coronary arteries; (2) Function from CT imaging: motion modeling, strain calculation, and visualization of four heart chambers; (3) Physiology from CT imaging: quantification and visualization of myocardial perfusion and contextual integration with coronary artery anatomy; (4) Physiology from computer simulation: computation and visualization of hemodynamics (e.g., coronary blood velocity, pressure, shear stress, and fluid forces on the vessel wall). Substantially, feedback from cardiologists have confirmed the practical utility of integrating these features for the purpose of computer-aided diagnosis of ischemic heart disease. PMID:26863663

Human-robot interaction modeling and simulation of supervisory control and situational awareness during field experimentation with military manned and unmanned ground vehicles

NASA Astrophysics Data System (ADS)

Johnson, Tony; Metcalfe, Jason; Brewster, Benjamin; Manteuffel, Christopher; Jaswa, Matthew; Tierney, Terrance

2010-04-01

The proliferation of intelligent systems in today's military demands increased focus on the optimization of human-robot interactions. Traditional studies in this domain involve large-scale field tests that require humans to operate semiautomated systems under varying conditions within military-relevant scenarios. However, provided that adequate constraints are employed, modeling and simulation can be a cost-effective alternative and supplement. The current presentation discusses a simulation effort that was executed in parallel with a field test with Soldiers operating military vehicles in an environment that represented key elements of the true operational context. In this study, "constructive" human operators were designed to represent average Soldiers executing supervisory control over an intelligent ground system. The constructive Soldiers were simulated performing the same tasks as those performed by real Soldiers during a directly analogous field test. Exercising the models in a high-fidelity virtual environment provided predictive results that represented actual performance in certain aspects, such as situational awareness, but diverged in others. These findings largely reflected the quality of modeling assumptions used to design behaviors and the quality of information available on which to articulate principles of operation. Ultimately, predictive analyses partially supported expectations, with deficiencies explicable via Soldier surveys, experimenter observations, and previously-identified knowledge gaps.

Improving Simulations of Extreme Flows by Coupling a Physically-based Hydrologic Model with a Machine Learning Model

NASA Astrophysics Data System (ADS)

Mohammed, K.; Islam, A. S.; Khan, M. J. U.; Das, M. K.

2017-12-01

With the large number of hydrologic models presently available along with the global weather and geographic datasets, streamflows of almost any river in the world can be easily modeled. And if a reasonable amount of observed data from that river is available, then simulations of high accuracy can sometimes be performed after calibrating the model parameters against those observed data through inverse modeling. Although such calibrated models can succeed in simulating the general trend or mean of the observed flows very well, more often than not they fail to adequately simulate the extreme flows. This causes difficulty in tasks such as generating reliable projections of future changes in extreme flows due to climate change, which is obviously an important task due to floods and droughts being closely connected to people's lives and livelihoods. We propose an approach where the outputs of a physically-based hydrologic model are used as an input to a machine learning model to try and better simulate the extreme flows. To demonstrate this offline-coupling approach, the Soil and Water Assessment Tool (SWAT) was selected as the physically-based hydrologic model, the Artificial Neural Network (ANN) as the machine learning model and the Ganges-Brahmaputra-Meghna (GBM) river system as the study area. The GBM river system, located in South Asia, is the third largest in the world in terms of freshwater generated and forms the largest delta in the world. The flows of the GBM rivers were simulated separately in order to test the performance of this proposed approach in accurately simulating the extreme flows generated by different basins that vary in size, climate, hydrology and anthropogenic intervention on stream networks. Results show that by post-processing the simulated flows of the SWAT models with ANN models, simulations of extreme flows can be significantly improved. The mean absolute errors in simulating annual maximum/minimum daily flows were minimized from 4967 cusecs to 1294 cusecs for Ganges, from 5695 cusecs to 2115 cusecs for Brahmaputra and from 689 cusecs to 321 cusecs for Meghna. Using this approach, simulations of hydrologic variables other than streamflow can also be improved given that a decent amount of observed data for that variable is available.

Processes influencing model-data mismatch in drought-stressed, fire-disturbed eddy flux sites

NASA Astrophysics Data System (ADS)

Mitchell, Stephen; Beven, Keith; Freer, Jim; Law, Beverly

2011-06-01

Semiarid forests are very sensitive to climatic change and among the most difficult ecosystems to accurately model. We tested the performance of the Biome-BGC model against eddy flux data taken from young (years 2004-2008), mature (years 2002-2008), and old-growth (year 2000) ponderosa pine stands at Metolius, Oregon, and subsequently examined several potential causes for model-data mismatch. We used the Generalized Likelihood Uncertainty Estimation methodology, which involved 500,000 model runs for each stand (1,500,000 total). Each simulation was run with randomly generated parameter values from a uniform distribution based on published parameter ranges, resulting in modeled estimates of net ecosystem CO2 exchange (NEE) that were compared to measured eddy flux data. Simulations for the young stand exhibited the highest level of performance, though they overestimated ecosystem C accumulation (-NEE) 99% of the time. Among the simulations for the mature and old-growth stands, 100% and 99% of the simulations underestimated ecosystem C accumulation. One obvious area of model-data mismatch is soil moisture, which was overestimated by the model in the young and old-growth stands yet underestimated in the mature stand. However, modeled estimates of soil water content and associated water deficits did not appear to be the primary cause of model-data mismatch; our analysis indicated that gross primary production can be accurately modeled even if soil moisture content is not. Instead, difficulties in adequately modeling ecosystem respiration, mainly autotrophic respiration, appeared to be the fundamental cause of model-data mismatch.

Processes influencing model-data mismatch in drought-stressed, fire-disturbed eddy flux sites

NASA Astrophysics Data System (ADS)

Mitchell, S. R.; Beven, K.; Freer, J. E.; Law, B. E.

2010-12-01

Semi-arid forests are very sensitive to climatic change and among the most difficult ecosystems to accurately model. We tested the performance of the Biome-BGC model against eddy flux data taken from young (years 2004-2008), mature (years 2002-2008), and old-growth (year 2000) Ponderosa pine stands at Metolius, Oregon, and subsequently examined several potential causes for model-data mismatch. We used the generalized likelihood uncertainty estimation (GLUE) methodology, which involved 500,000 model runs for each stand (1,500,000 total). Each simulation was run with randomly generated parameter values from a uniform distribution based on published parameter ranges, resulting in modeled estimates of net ecosystem CO2 exchange (NEE) that were compared to measured eddy flux data. Simulations for the young stand exhibited the highest level of performance, though they over-estimated ecosystem C accumulation (-NEE) 99% of the time. Among the simulations for the mature and old-growth stands, 100% and 99% of the simulations under-estimated ecosystem C accumulation. One obvious area of model-data mismatch is soil moisture, which was overestimated by the model in the young and old-growth stands yet underestimated in the mature stand. However, modeled estimates of soil water content and associated water deficits did not appear to be the primary cause of model-data mismatch; our analysis indicated that gross primary production can be accurately modeled even if soil moisture content is not. Instead, difficulties in adequately modeling ecosystem respiration, both autotrophic and heterotrophic, appeared to be fundamental causes of model-data mismatch.

A Film Depositional Model of Permeability for Mineral Reactions in Unsaturated Media.

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

Freedman, Vicky L.; Saripalli, Prasad; Bacon, Diana H.

2004-11-15

A new modeling approach based on the biofilm models of Taylor et al. (1990, Water Resources Research, 26, 2153-2159) has been developed for modeling changes in porosity and permeability in saturated porous media and implemented in an inorganic reactive transport code. Application of the film depositional models to mineral precipitation and dissolution reactions requires that calculations of mineral films be dynamically changing as a function of time dependent reaction processes. Since calculations of film thicknesses do not consider mineral density, results show that the film porosity model does not adequately describe volumetric changes in the porous medium. These effects canmore » be included in permeability calculations by coupling the film permeability models (Mualem and Childs and Collis-George) to a volumetric model that incorporates both mineral density and reactive surface area. Model simulations demonstrate that an important difference between the biofilm and mineral film models is in the translation of changes in mineral radii to changes in pore space. Including the effect of tortuosity on pore radii changes improves the performance of the Mualem permeability model for both precipitation and dissolution. Results from simulation of simultaneous dissolution and secondary mineral precipitation provides reasonable estimates of porosity and permeability. Moreover, a comparison of experimental and simulated data show that the model yields qualitatively reasonable results for permeability changes due to solid-aqueous phase reactions.« less

Influence of mass transfer resistance on overall nitrate removal rate in upflow sludge bed reactors.

PubMed

Ting, Wen-Huei; Huang, Ju-Sheng

2006-09-01

A kinetic model with intrinsic reaction kinetics and a simplified model with apparent reaction kinetics for denitrification in upflow sludge bed (USB) reactors were proposed. USB-reactor performance data with and without sludge wasting were also obtained for model verification. An independent batch study showed that the apparent kinetic constants k' did not differ from the intrinsic k but the apparent Ks' was significantly larger than the intrinsic Ks suggesting that the intra-granule mass transfer resistance can be modeled by changes in Ks. Calculations of the overall effectiveness factor, Thiele modulus, and Biot number combined with parametric sensitivity analysis showed that the influence of internal mass transfer resistance on the overall nitrate removal rate in USB reactors is more significant than the external mass transfer resistance. The simulated residual nitrate concentrations using the simplified model were in good agreement with the experimental data; the simulated results using the simplified model were also close to those using the kinetic model. Accordingly, the simplified model adequately described the overall nitrate removal rate and can be used for process design.

Experimental Investigations And Numerical Modelling of 210CR12 Steel in Semi-Solid State

NASA Astrophysics Data System (ADS)

Macioł, Piotr; Zalecki, Władysław; Kuziak, Roman; Jakubowicz, Aleksandra; Weglarczyk, Stanisław

2011-05-01

Experimental investigation, including hot compression and simple closed die filling was performed. Temperature range of tests was between 1225 °C and 1320 °C. Temperature selection was adequate with liquid fraction between 20 and 60%, which is typical for thixoforming processes. In the die filling test, steel dies with ceramic layer was used (highly refractory air-setting mortar JM 3300 manufactured by Thermal Ceramics). Experiments were carried out on the Gleeble 3800 physical simulator with MCU unit. In the paper, methodology of experimental investigation is described. Dependency of forming forces on temperature and forming velocities is analysed. Obtained results are discussed. The second part of the paper concerns numerical modelling of semi-solid forming. Numerical models for both sets of test were developed. Structural and Computational Fluid Dynamics models are compared. Initial works in microstructural modelling of 210CR12 steel behaviour are described. Lattice Boltzman Method model for thixotropic flows is introduced. Microscale and macroscale models were integrated into multiscale simulation of semi-solid forming. Some fundamental issues related to multiscale modelling of thixoforming are discussed.

Modeling of the flow stress for AISI H13 Tool Steel during Hard Machining Processes

NASA Astrophysics Data System (ADS)

Umbrello, Domenico; Rizzuti, Stefania; Outeiro, José C.; Shivpuri, Rajiv

2007-04-01

In general, the flow stress models used in computer simulation of machining processes are a function of effective strain, effective strain rate and temperature developed during the cutting process. However, these models do not adequately describe the material behavior in hard machining, where a range of material hardness between 45 and 60 HRC are used. Thus, depending on the specific material hardness different material models must be used in modeling the cutting process. This paper describes the development of a hardness-based flow stress and fracture models for the AISI H13 tool steel, which can be applied for range of material hardness mentioned above. These models were implemented in a non-isothermal viscoplastic numerical model to simulate the machining process for AISI H13 with various hardness values and applying different cutting regime parameters. Predicted results are validated by comparing them with experimental results found in the literature. They are found to predict reasonably well the cutting forces as well as the change in chip morphology from continuous to segmented chip as the material hardness change.

Simulating Space Capsule Water Landing with Explicit Finite Element Method

NASA Technical Reports Server (NTRS)

Wang, John T.; Lyle, Karen H.

2007-01-01

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Simulating smoke transport from wildland fires with a regional-scale air quality model: sensitivity to spatiotemporal allocation of fire emissions.

PubMed

Garcia-Menendez, Fernando; Hu, Yongtao; Odman, Mehmet T

2014-09-15

Air quality forecasts generated with chemical transport models can provide valuable information about the potential impacts of fires on pollutant levels. However, significant uncertainties are associated with fire-related emission estimates as well as their distribution on gridded modeling domains. In this study, we explore the sensitivity of fine particulate matter concentrations predicted by a regional-scale air quality model to the spatial and temporal allocation of fire emissions. The assessment was completed by simulating a fire-related smoke episode in which air quality throughout the Atlanta metropolitan area was affected on February 28, 2007. Sensitivity analyses were carried out to evaluate the significance of emission distribution among the model's vertical layers, along the horizontal plane, and into hourly inputs. Predicted PM2.5 concentrations were highly sensitive to emission injection altitude relative to planetary boundary layer height. Simulations were also responsive to the horizontal allocation of fire emissions and their distribution into single or multiple grid cells. Additionally, modeled concentrations were greatly sensitive to the temporal distribution of fire-related emissions. The analyses demonstrate that, in addition to adequate estimates of emitted mass, successfully modeling the impacts of fires on air quality depends on an accurate spatiotemporal allocation of emissions. Copyright © 2014 Elsevier B.V. All rights reserved.

Simulating the Risk of Liver Fluke Infection using a Mechanistic Hydro-epidemiological Model

NASA Astrophysics Data System (ADS)

Beltrame, Ludovica; Dunne, Toby; Rose, Hannah; Walker, Josephine; Morgan, Eric; Vickerman, Peter; Wagener, Thorsten

2016-04-01

Liver Fluke (Fasciola hepatica) is a common parasite found in livestock and responsible for considerable economic losses throughout the world. Risk of infection is strongly influenced by climatic and hydrological conditions, which characterise the host environment for parasite development and transmission. Despite on-going control efforts, increases in fluke outbreaks have been reported in recent years in the UK, and have been often attributed to climate change. Currently used fluke risk models are based on empirical relationships derived between historical climate and incidence data. However, hydro-climate conditions are becoming increasingly non-stationary due to climate change and direct anthropogenic impacts such as land use change, making empirical models unsuitable for simulating future risk. In this study we introduce a mechanistic hydro-epidemiological model for Liver Fluke, which explicitly simulates habitat suitability for disease development in space and time, representing the parasite life cycle in connection with key environmental conditions. The model is used to assess patterns of Liver Fluke risk for two catchments in the UK under current and potential future climate conditions. Comparisons are made with a widely used empirical model employing different datasets, including data from regional veterinary laboratories. Results suggest that mechanistic models can achieve adequate predictive ability and support adaptive fluke control strategies under climate change scenarios.

On the use of three hydrological models as hypotheses to investigate the behaviour of a small Mediterranean catchment

NASA Astrophysics Data System (ADS)

Ruiz Pérez, Guiomar; Latron, Jérôme; Llorens, Pilar; Gallart, Francesc; Francés, Félix

2017-04-01

Selecting an adequate hydrological model is the first step to carry out a rainfall-runoff modelling exercise. A hydrological model is a hypothesis of catchment functioning, encompassing a description of dominant hydrological processes and predicting how these processes interact to produce the catchment's response to external forcing. Current research lines emphasize the importance of multiple working hypotheses for hydrological modelling instead of only using a single model. In line with this philosophy, here different hypotheses were considered and analysed to simulate the nonlinear response of a small Mediterranean catchment and to progress in the analysis of its hydrological behaviour. In particular, three hydrological models were considered representing different potential hypotheses: two lumped models called LU3 and LU4, and one distributed model called TETIS. To determine how well each specific model performed and to assess whether a model was more adequate than another, we raised three complementary tests: one based on the analysis of residual errors series, another based on a sensitivity analysis and the last one based on using multiple evaluation criteria associated to the concept of Pareto frontier. This modelling approach, based on multiple working hypotheses, helped to improve our perceptual model of the catchment behaviour and, furthermore, could be used as a guidance to improve the performance of other environmental models.

Numerical Study of a Convective Turbulence Encounter

NASA Technical Reports Server (NTRS)

Proctor, Fred H.; Hamilton, David W.; Bowles, Roland L.

2002-01-01

A numerical simulation of a convective turbulence event is investigated and compared with observational data. The specific case was encountered during one of NASA's flight tests and was characterized by severe turbulence. The event was associated with overshooting convective turrets that contained low to moderate radar reflectivity. Model comparisons with observations are quite favorable. Turbulence hazard metrics are proposed and applied to the numerical data set. Issues such as adequate grid size are examined.

Doppler lidar wind measurement on Eos

NASA Technical Reports Server (NTRS)

Fitzjarrald, D.; Bilbro, J.; Beranek, R.; Mabry, J.

1985-01-01

A polar-orbiting platform segment of the Earth Observing System (EOS) could carry a CO2-laser based Doppler lidar for recording global wind profiles. Development goals would include the manufacture of a 10 J laser with a 2 yr operational life, space-rating the optics and associated software, and the definition of models for global aerosol distributions. Techniques will be needed for optimal scanning and generating computer simulations which will provide adequately accurate weather predictions.

A model for phosphorus transformation and runoff loss for surface-applied manures.

PubMed

Vadas, P A; Gburek, W J; Sharpley, A N; Kleinman, P J A; Moore, P A; Cabrera, M L; Harmel, R D

2007-01-01

Agricultural P transport in runoff is an environmental concern. An important source of P runoff is surface-applied, unincorporated manures, but computer models used to assess P transport do not adequately simulate P release and transport from surface manures. We developed a model to address this limitation. The model operates on a daily basis and simulates manure application to the soil surface, letting 60% of manure P infiltrate into soil if manure slurry with less than 15% solids is applied. The model divides manure P into four pools, water-extractable inorganic and organic P, and stable inorganic and organic P. The model simulates manure dry matter decomposition, and manure stable P transformation to water-extractable P. Manure dry matter and P are assimilated into soil to simulate bioturbation. Water-extractable P is leached from manure when it rains, and a portion of leached P can be transferred to surface runoff. Eighty percent of manure P leached into soil by rain remains in the top 2 cm, while 20% leaches deeper. This 2-cm soil layer contributes P to runoff via desorption. We used data from field studies in Texas, Pennsylvania, Georgia, and Arkansas to build and validate the model. Validation results show the model accurately predicted cumulative P loads in runoff, reflecting successful simulation of the dynamics of manure dry matter, manure and soil P pools, and storm-event runoff P concentrations. Predicted runoff P concentrations were significantly related to (r2=0.57) but slightly less than measured concentrations. Our model thus represents an important modification for field or watershed scale models that assess P loss from manured soils.

Regional climate simulations over South America: sensitivity to model physics and to the treatment of lateral boundary conditions using the MM5 model

NASA Astrophysics Data System (ADS)

Solman, Silvina A.; Pessacg, Natalia L.

2012-01-01

In this study the capability of the MM5 model in simulating the main mode of intraseasonal variability during the warm season over South America is evaluated through a series of sensitivity experiments. Several 3-month simulations nested into ERA40 reanalysis were carried out using different cumulus schemes and planetary boundary layer schemes in an attempt to define the optimal combination of physical parameterizations for simulating alternating wet and dry conditions over La Plata Basin (LPB) and the South Atlantic Convergence Zone regions, respectively. The results were compared with different observational datasets and model evaluation was performed taking into account the spatial distribution of monthly precipitation and daily statistics of precipitation over the target regions. Though every experiment was able to capture the contrasting behavior of the precipitation during the simulated period, precipitation was largely underestimated particularly over the LPB region, mainly due to a misrepresentation in the moisture flux convergence. Experiments using grid nudging of the winds above the planetary boundary layer showed a better performance compared with those in which no constrains were imposed to the regional circulation within the model domain. Overall, no single experiment was found to perform the best over the entire domain and during the two contrasting months. The experiment that outperforms depends on the area of interest, being the simulation using the Grell (Kain-Fritsch) cumulus scheme in combination with the MRF planetary boundary layer scheme more adequate for subtropical (tropical) latitudes. The ensemble of the sensitivity experiments showed a better performance compared with any individual experiment.

Effect of lift-to-drag ratio in pilot rating of the HL-20 landing task

NASA Technical Reports Server (NTRS)

Jackson, E. B.; Rivers, Robert A.; Bailey, Melvin L.

1993-01-01

A man-in-the-loop simulation study of the handling qualities of the HL-20 lifting-body vehicle was made in a fixed-base simulation cockpit at NASA Langley Research Center. The purpose of the study was to identify and substantiate opportunities for improving the original design of the vehicle from a handling qualities and landing performance perspective. Using preliminary wind-tunnel data, a subsonic aerodynamic model of the HL-20 was developed. This model was adequate to simulate the last 75-90 s of the approach and landing. A simple flight-control system was designed and implemented. Using this aerodynamic model as a baseline, visual approaches and landings were made at several vehicle lift-to-drag ratios. Pilots rated the handling characteristics of each configuration using a conventional numerical pilot-rating scale. Results from the study showed a high degree of correlation between the lift-to-drag ratio and pilot rating. Level 1 pilot ratings were obtained when the L/D ratio was approximately 3.8 or higher.

Effect of lift-to-drag ratio in pilot rating of the HL-20 landing task

NASA Astrophysics Data System (ADS)

Jackson, E. B.; Rivers, Robert A.; Bailey, Melvin L.

1993-10-01

A man-in-the-loop simulation study of the handling qualities of the HL-20 lifting-body vehicle was made in a fixed-base simulation cockpit at NASA Langley Research Center. The purpose of the study was to identify and substantiate opportunities for improving the original design of the vehicle from a handling qualities and landing performance perspective. Using preliminary wind-tunnel data, a subsonic aerodynamic model of the HL-20 was developed. This model was adequate to simulate the last 75-90 s of the approach and landing. A simple flight-control system was designed and implemented. Using this aerodynamic model as a baseline, visual approaches and landings were made at several vehicle lift-to-drag ratios. Pilots rated the handling characteristics of each configuration using a conventional numerical pilot-rating scale. Results from the study showed a high degree of correlation between the lift-to-drag ratio and pilot rating. Level 1 pilot ratings were obtained when the L/D ratio was approximately 3.8 or higher.

Independent Review of Simulation of Net Infiltration for Present-Day and Potential Future Climates

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

Review Panel: Soroosh Sorooshian, Ph.D., Panel Chairperson, University of California, Irvine; Jan M. H. Hendrickx, Ph.D., New Mexico Institute of Mining and Technology; Binayak P. Mohanty, Ph.D., Texas A&M University

The DOE Office of Civilian Radioactive Waste Management (OCRWM) tasked Oak Ridge Institute for Science and Education (ORISE) with providing an independent expert review of the documented model and prediction results for net infiltration of water into the unsaturated zone at Yucca Mountain. The specific purpose of the model, as documented in the report MDL-NBS-HS-000023, Rev. 01, is “to provide a spatial representation, including epistemic and aleatory uncertainty, of the predicted mean annual net infiltration at the Yucca Mountain site ...” (p. 1-1) The expert review panel assembled by ORISE concluded that the model report does not provide a technicallymore » credible spatial representation of net infiltration at Yucca Mountain. Specifically, the ORISE Review Panel found that: • A critical lack of site-specific meteorological, surface, and subsurface information prevents verification of (i) the net infiltration estimates, (ii) the uncertainty estimates of parameters caused by their spatial variability, and (iii) the assumptions used by the modelers (ranges and distributions) for the characterization of parameters. The paucity of site-specific data used by the modeling team for model implementation and validation is a major deficiency in this effort. • The model does not incorporate at least one potentially important hydrologic process. Subsurface lateral flow is not accounted for by the model, and the assumption that the effect of subsurface lateral flow is negligible is not adequately justified. This issue is especially critical for the wetter climate periods. This omission may be one reason the model results appear to underestimate net infiltration beneath wash environments and therefore imprecisely represent the spatial variability of net infiltration. • While the model uses assumptions consistently, such as uniform soil depths and a constant vegetation rooting depth, such assumptions may not be appropriate for this net infiltration simulation because they oversimplify a complex landscape and associated hydrologic processes, especially since the model assumptions have not been adequately corroborated by field and laboratory observations at Yucca Mountain.« less

Technical tips and advancements in pediatric minimally invasive surgical training on porcine based simulations.

PubMed

Narayanan, Sarath Kumar; Cohen, Ralph Clinton; Shun, Albert

2014-06-01

Minimal access techniques have transformed the way pediatric surgery is practiced. Due to various constraints, surgical residency programs have not been able to tutor adequate training skills in the routine setting. The advent of new technology and methods in minimally invasive surgery (MIS), has similarly contributed to the need for systematic skills' training in a safe, simulated environment. To enable the training of the proper technique among pediatric surgery trainees, we have advanced a porcine non-survival model for endoscopic surgery. The technical advancements over the past 3 years and a subjective validation of the porcine model from 114 participating trainees using a standard questionnaire and a 5-point Likert scale have been described here. Mean attitude scores and analysis of variance (ANOVA) were used for statistical analysis of the data. Almost all trainees agreed or strongly agreed that the animal-based model was appropriate (98.35%) and also acknowledged that such workshops provided adequate practical experience before attempting on human subjects (96.6%). Mean attitude score for respondents was 19.08 (SD 3.4, range 4-20). Attitude scores showed no statistical association with years of experience or the level of seniority, indicating a positive attitude among all groups of respondents. Structured porcine-based MIS training should be an integral part of skill acquisition for pediatric surgery trainees and the experience gained can be transferred into clinical practice. We advocate that laparoscopic training should begin in a controlled workshop setting before procedures are attempted on human patients.

3D SPH numerical simulation of the wave generated by the Vajont rockslide

NASA Astrophysics Data System (ADS)

Vacondio, R.; Mignosa, P.; Pagani, S.

2013-09-01

A 3D numerical modeling of the wave generated by the Vajont slide, one of the most destructive ever occurred, is presented in this paper. A meshless Lagrangian Smoothed Particle Hydrodynamics (SPH) technique was adopted to simulate the highly fragmented violent flow generated by the falling slide in the artificial reservoir. The speed-up achievable via General Purpose Graphic Processing Units (GP-GPU) allowed to adopt the adequate resolution to describe the phenomenon. The comparison with the data available in literature showed that the results of the numerical simulation reproduce satisfactorily the maximum run-up, also the water surface elevation in the residual lake after the event. Moreover, the 3D velocity field of the flow during the event and the discharge hydrograph which overtopped the dam, were obtained.

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

Ewy, Ann; Heim, Kenneth J.; McGonigal, Sean T.

A comparative groundwater hydrogeologic modeling analysis is presented herein to simulate potential contaminant migration pathways in a sole source aquifer in Nassau County, Long Island, New York. The source of contamination is related to historical operations at the Sylvania Corning Plant ('Site'), a 9.49- acre facility located at 70, 100 and 140 Cantiague Rock Road, Town of Oyster Bay in the westernmost portion of Hicksville, Long Island. The Site had historically been utilized as a nuclear materials manufacturing facility (e.g., cores, slug, and fuel elements) for reactors used in both research and electric power generation in early 1950's until latemore » 1960's. The Site is contaminated with various volatile organic and inorganic compounds, as well as radionuclides. The major contaminants of concern at the Site are tetrachloroethene (PCE), trichloroethene (TCE), nickel, uranium, and thorium. These compounds are present in soil and groundwater underlying the Site and have migrated off-site. The Site is currently being investigated as part of the Formerly Utilized Sites Remedial Action Program (FUSRAP). The main objective of the current study is to simulate the complex hydrogeologic features in the region, such as numerous current and historic production well fields; large, localized recharge basins; and, multiple aquifers, and to assess potential contaminant migration pathways originating from the Site. For this purpose, the focus of attention was given to the underlying Magothy formation, which has been impacted by the contaminants of concern. This aquifer provides more than 90% of potable water supply in the region. Nassau and Suffolk Counties jointly developed a three-dimensional regional groundwater flow model to help understand the factors affecting groundwater flow regime in the region, to determine adequate water supply for public consumption, to investigate salt water intrusion in localized areas, to evaluate the impacts of regional pumping activity, and to better understand the contaminant transport and fate mechanisms through the underlying aquifers. This regional model, developed for the N.Y. State Department of Environmental Conservation (NYSDEC) by Camp Dresser and McKee (CDM), uses the finite element model DYNFLOW developed by CDM, Cambridge, Massachusetts. The coarseness of the regional model, however, could not adequately capture the hydrogeologic heterogeneity of the aquifer. Specifically, the regional model did not adequately capture the interbedded nature of the Magothy aquifer and, as such, simulated particles tended to track down-gradient from the Site in relatively straight lines while the movement of groundwater in such a heterogeneous aquifer is expected to proceed along a more tortuous path. This paper presents a qualitative comparison of site-specific groundwater flow modeling results with results obtained from the regional model. In order to assess the potential contaminant migration pathways, a particle tracking method was employed. Available site-specific and regional hydraulic conductivity data measured in-situ with respect to depth and location were incorporated into the T-PROG module in GMS model to define statistical variation to better represent the actual stratigraphy and layer heterogeneity. The groundwater flow characteristics in the Magothy aquifer were simulated with the stochastic hydraulic conductivity variation as opposed to constant values as employed in the regional model. Contaminant sources and their exact locations have been fully delineated at the Site during the Remedial Investigation (RI) phase of the project. Contaminant migration pathways originating from these source locations at the Site are qualitatively traced within the sole source aquifer utilizing particles introduced at source locations. Contaminant transport mechanism modeled in the current study is based on pure advection (i.e., plug flow) and mechanical dispersion while molecular diffusion effects are neglected due to relatively high groundwater velocities encountered in the aquifer. In addition, fate of contaminants is ignored hereby to simulate the worst-case scenario, which considers the contaminants of concern as tracer-like compounds for modeling purposes. The results of the modeling analysis are qualitatively compared with the County's regional model, and patterns of contaminant migration in the region are presented. (authors)« less

Land Use Management in the Panama Canal Watershed to Maximize Hydrologic Ecosystem Services Benefits: Explicit Simulation of Preferential Flow Paths in an HPC Environment

NASA Astrophysics Data System (ADS)

Regina, J. A.; Ogden, F. L.; Steinke, R. C.; Frazier, N.; Cheng, Y.; Zhu, J.

2017-12-01

Preferential flow paths (PFP) resulting from biotic and abiotic factors contribute significantly to the generation of runoff in moist lowland tropical watersheds. Flow through PFPs represents the dominant mechanism by which land use choices affect hydrological behavior. The relative influence of PFP varies depending upon land-use management practices. Assessing the possible effects of land-use and landcover change on flows, and other ecosystem services, in the humid tropics partially depends on adequate simulation of PFP across different land-uses. Currently, 5% of global trade passes through the Panama Canal, which is supplied with fresh water from the Panama Canal Watershed. A third set of locks, recently constructed, are expected to double the capacity of the Canal. We incorporated explicit simulation of PFPs in to the ADHydro HPC distributed hydrological model to simulate the effects of land-use and landcover change due to land management incentives on water resources availability in the Panama Canal Watershed. These simulations help to test hypotheses related to the effectiveness of various proposed payments for ecosystem services schemes. This presentation will focus on hydrological model formulation and performance in an HPC environment.

Short-term responses of leaf growth rate to water deficit scale up to whole-plant and crop levels: an integrated modelling approach in maize.

PubMed

Chenu, Karine; Chapman, Scott C; Hammer, Graeme L; McLean, Greg; Salah, Halim Ben Haj; Tardieu, François

2008-03-01

Physiological and genetic studies of leaf growth often focus on short-term responses, leaving a gap to whole-plant models that predict biomass accumulation, transpiration and yield at crop scale. To bridge this gap, we developed a model that combines an existing model of leaf 6 expansion in response to short-term environmental variations with a model coordinating the development of all leaves of a plant. The latter was based on: (1) rates of leaf initiation, appearance and end of elongation measured in field experiments; and (2) the hypothesis of an independence of the growth between leaves. The resulting whole-plant leaf model was integrated into the generic crop model APSIM which provided dynamic feedback of environmental conditions to the leaf model and allowed simulation of crop growth at canopy level. The model was tested in 12 field situations with contrasting temperature, evaporative demand and soil water status. In observed and simulated data, high evaporative demand reduced leaf area at the whole-plant level, and short water deficits affected only leaves developing during the stress, either visible or still hidden in the whorl. The model adequately simulated whole-plant profiles of leaf area with a single set of parameters that applied to the same hybrid in all experiments. It was also suitable to predict biomass accumulation and yield of a similar hybrid grown in different conditions. This model extends to field conditions existing knowledge of the environmental controls of leaf elongation, and can be used to simulate how their genetic controls flow through to yield.

THE SIMULATION OF FINE SCALE NOCTURNAL BOUNDARY LAYER MOTIONS WITH A MESO-SCALE ATMOSPHERIC MODEL

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

Werth, D.; Kurzeja, R.; Parker, M.

A field project over the Atmospheric Radiation Measurement-Clouds and Radiation Testbed (ARM-CART) site during a period of several nights in September, 2007 was conducted to explore the evolution of the low-level jet (LLJ). Data was collected from a tower and a sodar and analyzed for turbulent behavior. To study the full range of nocturnal boundary layer (NBL) behavior, the Regional Atmospheric Modeling System (RAMS) was used to simulate the ARM-CART NBL field experiment and validated against the data collected from the site. This model was run at high resolution, and is ideal for calculating the interactions among the various motionsmore » within the boundary layer and their influence on the surface. The model reproduces adequately the synoptic situation and the formation and dissolution cycles of the low-level jet, although it suffers from insufficient cloud production and excessive nocturnal cooling. The authors suggest that observed heat flux data may further improve the realism of the simulations both in the cloud formation and in the jet characteristics. In a higher resolution simulation, the NBL experiences motion on a range of timescales as revealed by a wavelet analysis, and these are affected by the presence of the LLJ. The model can therefore be used to provide information on activity throughout the depth of the NBL.« less

Evaluation of the Surface Representation of the Greenland Ice Sheet in a General Circulation Model

NASA Technical Reports Server (NTRS)

Cullather, Richard I.; Nowicki, Sophie M. J.; Zhao, Bin; Suarez, Max J.

2014-01-01

Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS 5) atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980-2008 GrIS SMB average is 24.7+/-4.5 cm yr(- 1) water-equivalent (w.e.) at.5 degree model grid spacing, and 18.2+/-3.3 cm yr(- 1) w.e. for 2 degree grid spacing. The spatial variability and seasonal cycle of the simulation compare favorably to recent studies using regional climate models, while results from 2 degree integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.

Atomistic simulation of the coupled adsorption and unfolding of protein GB1 on the polystyrenes nanoparticle surface

NASA Astrophysics Data System (ADS)

Xiao, HuiFang; Huang, Bin; Yao, Ge; Kang, WenBin; Gong, Sheng; Pan, Hai; Cao, Yi; Wang, Jun; Zhang, Jian; Wang, Wei

2018-03-01

Understanding the processes of protein adsorption/desorption on nanoparticles' surfaces is important for the development of new nanotechnology involving biomaterials; however, an atomistic resolution picture for these processes and for the simultaneous protein conformational change is missing. Here, we report the adsorption of protein GB1 on a polystyrene nanoparticle surface using atomistic molecular dynamic simulations. Enabled by metadynamics, we explored the relevant phase space and identified three protein states, each involving both the adsorbed and desorbed modes. We also studied the change of the secondary and tertiary structures of GB1 during adsorption and the dominant interactions between the protein and surface in different adsorption stages. The results we obtained from simulation were found to be more adequate and complete than the previous one. We believe the model presented in this paper, in comparison with the previous ones, is a better theoretical model to understand and explain the experimental results.

A simple statistical model for geomagnetic reversals

NASA Technical Reports Server (NTRS)

Constable, Catherine

1990-01-01

The diversity of paleomagnetic records of geomagnetic reversals now available indicate that the field configuration during transitions cannot be adequately described by simple zonal or standing field models. A new model described here is based on statistical properties inferred from the present field and is capable of simulating field transitions like those observed. Some insight is obtained into what one can hope to learn from paleomagnetic records. In particular, it is crucial that the effects of smoothing in the remanence acquisition process be separated from true geomagnetic field behavior. This might enable us to determine the time constants associated with the dominant field configuration during a reversal.

Momentum loss in proton-nucleus and nucleus-nucleus collisions

NASA Technical Reports Server (NTRS)

Khan, Ferdous; Townsend, Lawrence W.

1993-01-01

An optical model description, based on multiple scattering theory, of longitudinal momentum loss in proton-nucleus and nucleus-nucleus collisions is presented. The crucial role of the imaginary component of the nucleon-nucleon transition matrix in accounting for longitudinal momentum transfer is demonstrated. Results obtained with this model are compared with Intranuclear Cascade (INC) calculations, as well as with predictions from Vlasov-Uehling-Uhlenbeck (VUU) and quantum molecular dynamics (QMD) simulations. Comparisons are also made with experimental data where available. These indicate that the present model is adequate to account for longitudinal momentum transfer in both proton-nucleus and nucleus-nucleus collisions over a wide range of energies.

Ventilation inhomogeneity in obstructive lung diseases measured by electrical impedance tomography: a simulation study.

PubMed

Schullcke, B; Krueger-Ziolek, S; Gong, B; Jörres, R A; Mueller-Lisse, U; Moeller, K

2017-10-10

Electrical impedance tomography (EIT) has mostly been used in the Intensive Care Unit (ICU) to monitor ventilation distribution but is also promising for the diagnosis in spontaneously breathing patients with obstructive lung diseases. Beside tomographic images, several numerical measures have been proposed to quantitatively assess the lung state. In this study two common measures, the 'Global Inhomogeneity Index' and the 'Coefficient of Variation' were compared regarding their capability to reflect the severity of lung obstruction. A three-dimensional simulation model was used to simulate obstructed lungs, whereby images were reconstructed on a two-dimensional domain. Simulations revealed that minor obstructions are not adequately recognized in the reconstructed images and that obstruction above and below the electrode plane may result in misleading values of inhomogeneity measures. EIT measurements on several electrode planes are necessary to apply these measures in patients with obstructive lung diseases in a promising manner.

Large eddy simulation of turbulent premixed combustion using tabulated detailed chemistry and presumed probability density function

NASA Astrophysics Data System (ADS)

Zhang, Hongda; Han, Chao; Ye, Taohong; Ren, Zhuyin

2016-03-01

A method of chemistry tabulation combined with presumed probability density function (PDF) is applied to simulate piloted premixed jet burner flames with high Karlovitz number using large eddy simulation. Thermo-chemistry states are tabulated by the combination of auto-ignition and extended auto-ignition model. To evaluate the predictive capability of the proposed tabulation method to represent the thermo-chemistry states under the condition of different fresh gases temperature, a-priori study is conducted by performing idealised transient one-dimensional premixed flame simulations. Presumed PDF is used to involve the interaction of turbulence and flame with beta PDF to model the reaction progress variable distribution. Two presumed PDF models, Dirichlet distribution and independent beta distribution, respectively, are applied for representing the interaction between two mixture fractions that are associated with three inlet streams. Comparisons of statistical results show that two presumed PDF models for the two mixture fractions are both capable of predicting temperature and major species profiles, however, they are shown to have a significant effect on the predictions for intermediate species. An analysis of the thermo-chemical state-space representation of the sub-grid scale (SGS) combustion model is performed by comparing correlations between the carbon monoxide mass fraction and temperature. The SGS combustion model based on the proposed chemistry tabulation can reasonably capture the peak value and change trend of intermediate species. Aspects regarding model extensions to adequately predict the peak location of intermediate species are discussed.

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models. [probability density function

NASA Technical Reports Server (NTRS)

Hsu, Andrew T.

1992-01-01

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Theoretical models for coronary vascular biomechanics: Progress & challenges

PubMed Central

Waters, Sarah L.; Alastruey, Jordi; Beard, Daniel A.; Bovendeerd, Peter H.M.; Davies, Peter F.; Jayaraman, Girija; Jensen, Oliver E.; Lee, Jack; Parker, Kim H.; Popel, Aleksander S.; Secomb, Timothy W.; Siebes, Maria; Sherwin, Spencer J.; Shipley, Rebecca J.; Smith, Nicolas P.; van de Vosse, Frans N.

2013-01-01

A key aim of the cardiac Physiome Project is to develop theoretical models to simulate the functional behaviour of the heart under physiological and pathophysiological conditions. Heart function is critically dependent on the delivery of an adequate blood supply to the myocardium via the coronary vasculature. Key to this critical function of the coronary vasculature is system dynamics that emerge via the interactions of the numerous constituent components at a range of spatial and temporal scales. Here, we focus on several components for which theoretical approaches can be applied, including vascular structure and mechanics, blood flow and mass transport, flow regulation, angiogenesis and vascular remodelling, and vascular cellular mechanics. For each component, we summarise the current state of the art in model development, and discuss areas requiring further research. We highlight the major challenges associated with integrating the component models to develop a computational tool that can ultimately be used to simulate the responses of the coronary vascular system to changing demands and to diseases and therapies. PMID:21040741

PUFoam : A novel open-source CFD solver for the simulation of polyurethane foams

NASA Astrophysics Data System (ADS)

Karimi, M.; Droghetti, H.; Marchisio, D. L.

2017-08-01

In this work a transient three-dimensional mathematical model is formulated and validated for the simulation of polyurethane (PU) foams. The model is based on computational fluid dynamics (CFD) and is coupled with a population balance equation (PBE) to describe the evolution of the gas bubbles/cells within the PU foam. The front face of the expanding foam is monitored on the basis of the volume-of-fluid (VOF) method using a compressible solver available in OpenFOAM version 3.0.1. The solver is additionally supplemented to include the PBE, solved with the quadrature method of moments (QMOM), the polymerization kinetics, an adequate rheological model and a simple model for the foam thermal conductivity. The new solver is labelled as PUFoam and is, for the first time in this work, validated for 12 different mixing-cup experiments. Comparison of the time evolution of the predicted and experimentally measured density and temperature of the PU foam shows the potentials and limitations of the approach.

Towards understanding of magnetization reversal in Nd-Fe-B nanocomposites: analysis by high-throughput micromagnetic simulations

NASA Astrophysics Data System (ADS)

Erokhin, Sergey; Berkov, Dmitry; Ito, Masaaki; Kato, Akira; Yano, Masao; Michels, Andreas

2018-03-01

We demonstrate how micromagnetic simulations can be employed in order to characterize and analyze the magnetic microstructure of nanocomposites. For the example of nanocrystalline Nd-Fe-B, which is a potential material for future permanent-magnet applications, we have compared three different models for the micromagnetic analysis of this material class: (i) a description of the nanocomposite microstructure in terms of Stoner-Wohlfarth particles with and without the magnetodipolar interaction; (ii) a model based on the core-shell representation of the nanograins; (iii) the latter model including a contribution of superparamagnetic clusters. The relevant parameter spaces have been systematically scanned with the aim to establish which micromagnetic approach can most adequately describe experimental data for this material. According to our results, only the last, most sophisticated model is able to provide an excellent agreement with the measured hysteresis loop. The presented methodology is generally applicable to multiphase magnetic nanocomposites and it highligths the complex interrelationship between the microstructure, magnetic interactions, and the macroscopic magnetic properties.

Numerical modeling of subsurface radioactive solute transport from waste seepage ponds at the Idaho National Engineering Laboratory

USGS Publications Warehouse

Robertson, John B.

1976-01-01

Aqueous chemical and low-level radioactive effluents have been disposed to seepage ponds since 1952 at the Idaho National Engineering Laboratory. The solutions percolate toward the Snake River Plain aquifer (135 m below) through interlayered basalts and unconsolidated sediments and an extensive zone of ground water perched on a sedimentary layer about 40 m beneath the ponds. A three-segment numerical model was developed to simulate the system, including effects of convection, hydrodynamic dispersion, radioactive decay, and adsorption. Simulated hydraulics and solute migration patterns for all segments agree adequately with the available field data. The model can be used to project subsurface distributions of waste solutes under a variety of assumed conditions for the future. Although chloride and tritium reached the aquifer several years ago, the model analysis suggests that the more easily sorbed solutes, such as cesium-137 and strontium-90, would not reach the aquifer in detectable concentrations within 150 years for the conditions assumed. (Woodard-USGS)

Link-prediction to tackle the boundary specification problem in social network surveys

PubMed Central

De Wilde, Philippe; Buarque de Lima-Neto, Fernando

2017-01-01

Diffusion processes in social networks often cause the emergence of global phenomena from individual behavior within a society. The study of those global phenomena and the simulation of those diffusion processes frequently require a good model of the global network. However, survey data and data from online sources are often restricted to single social groups or features, such as age groups, single schools, companies, or interest groups. Hence, a modeling approach is required that extrapolates the locally restricted data to a global network model. We tackle this Missing Data Problem using Link-Prediction techniques from social network research, network generation techniques from the area of Social Simulation, as well as a combination of both. We found that techniques employing less information may be more adequate to solve this problem, especially when data granularity is an issue. We validated the network models created with our techniques on a number of real-world networks, investigating degree distributions as well as the likelihood of links given the geographical distance between two nodes. PMID:28426826

Comparison of penumbra regions produced by ancient Gamma knife model C and Gamma ART 6000 using Monte Carlo MCNP6 simulation.

PubMed

Banaee, Nooshin; Asgari, Sepideh; Nedaie, Hassan Ali

2018-07-01

The accuracy of penumbral measurements in radiotherapy is pivotal because dose planning computers require accurate data to adequately modeling the beams, which in turn are used to calculate patient dose distributions. Gamma knife is a non-invasive intracranial technique based on principles of the Leksell stereotactic system for open deep brain surgeries, invented and developed by Professor Lars Leksell. The aim of this study is to compare the penumbra widths of Leksell Gamma Knife model C and Gamma ART 6000. Initially, the structure of both systems were simulated by using Monte Carlo MCNP6 code and after validating the accuracy of simulation, beam profiles of different collimators were plotted. MCNP6 beam profile calculations showed that the penumbra values of Leksell Gamma knife model C and Gamma ART 6000 for 18, 14, 8 and 4 mm collimators are 9.7, 7.9, 4.3, 2.6 and 8.2, 6.9, 3.6, 2.4, respectively. The results of this study showed that since Gamma ART 6000 has larger solid angle in comparison with Gamma Knife model C, it produces better beam profile penumbras than Gamma Knife model C in the direct plane. Copyright © 2017 Elsevier Ltd. All rights reserved.

Development of a Non-Invasive Biomonitoring Approach to Determine Exposure to the Organophosphorus Insecticide Chlorpyrifos in Rat Saliva

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

Timchalk, Chuck; Campbell, James A.; Liu, Guodong

2007-03-01

Abstract Non-invasive biomonitoring approaches are being developed using reliable portable analytical systems to quantify dosimetry utilizing readily obtainable body fluids, such as saliva. In the current study, rats were given single oral gavage doses (1, 10 or 50 mg/kg) of the insecticide chlorpyrifos (CPF), saliva and blood were collected from groups of animals (4/time-point) at 3, 6, and 12 hr post-dosing, and the samples were analyzed for the CPF metabolite trichlorpyridinol (TCP). Trichlorpyridinol was detected in both blood and saliva at all doses and the TCP concentration in blood exceeded saliva, although the kinetics in blood and saliva were comparable.more » A physiologically based pharmacokinetic and pharmacodynamic (PBPK/PD) model for CPF incorporated a compartment model to describe the time-course of TCP in blood and saliva. The model adequately simulated the experimental results over the dose ranges evaluated. A rapid and sensitive sequential injection (SI) electrochemical immunoassay was developed to monitor TCP, and the reported detection limit for TCP in water was 6 ng/L. Computer model simulation in the range of the Allowable Daily Intake (ADI) or Reference Dose (RfD) for CPF (0.01-0.003 mg/kg/day) suggest that the electrochemical immunoassay had adequate sensitivity to detect and quantify TCP in saliva at these low exposure levels. To validate this approach further studies are needed to more fully understand the pharmacokinetics of CPF and TCP excretion in saliva. The utilization of saliva as a biomonitoring matrix, coupled to real-time quantitation and PBPK/PD modeling represents a novel approach with broad application for evaluating both occupational and environmental exposures to insecticides.« less

14 CFR 142.15 - Facilities.

Code of Federal Regulations, 2014 CFR

2014-01-01

... adequate periods of time and at a location approved by the Administrator, adequate flight training equipment and courseware, including at least one flight simulator or advanced flight training device. [Doc... significant distractions caused by flight operations and maintenance operations at the airport. (b) An...

14 CFR 142.15 - Facilities.

Code of Federal Regulations, 2011 CFR

2011-01-01

... adequate periods of time and at a location approved by the Administrator, adequate flight training equipment and courseware, including at least one flight simulator or advanced flight training device. [Doc... significant distractions caused by flight operations and maintenance operations at the airport. (b) An...

14 CFR 142.15 - Facilities.

Code of Federal Regulations, 2013 CFR

2013-01-01

... adequate periods of time and at a location approved by the Administrator, adequate flight training equipment and courseware, including at least one flight simulator or advanced flight training device. [Doc... significant distractions caused by flight operations and maintenance operations at the airport. (b) An...

14 CFR 142.15 - Facilities.

Code of Federal Regulations, 2012 CFR

2012-01-01

... adequate periods of time and at a location approved by the Administrator, adequate flight training equipment and courseware, including at least one flight simulator or advanced flight training device. [Doc... significant distractions caused by flight operations and maintenance operations at the airport. (b) An...

Differentiability of simulated MEG hippocampal, medial temporal and neocortical temporal epileptic spike activity.

PubMed

Stephen, Julia M; Ranken, Doug M; Aine, Cheryl J; Weisend, Michael P; Shih, Jerry J

2005-12-01

Previous studies have shown that magnetoencephalography (MEG) can measure hippocampal activity, despite the cylindrical shape and deep location in the brain. The current study extended this work by examining the ability to differentiate the hippocampal subfields, parahippocampal cortex, and neocortical temporal sources using simulated interictal epileptic activity. A model of the hippocampus was generated on the MRIs of five subjects. CA1, CA3, and dentate gyrus of the hippocampus were activated as well as entorhinal cortex, presubiculum, and neocortical temporal cortex. In addition, pairs of sources were activated sequentially to emulate various hypotheses of mesial temporal lobe seizure generation. The simulated MEG activity was added to real background brain activity from the five subjects and modeled using a multidipole spatiotemporal modeling technique. The waveforms and source locations/orientations for hippocampal and parahippocampal sources were differentiable from neocortical temporal sources. In addition, hippocampal and parahippocampal sources were differentiated to varying degrees depending on source. The sequential activation of hippocampal and parahippocampal sources was adequately modeled by a single source; however, these sources were not resolvable when they overlapped in time. These results suggest that MEG has the sensitivity to distinguish parahippocampal and hippocampal spike generators in mesial temporal lobe epilepsy.

Modeling of rock friction 2. Simulation of preseismic slip

USGS Publications Warehouse

Dieterich, J.H.

1979-01-01

The constitutive relations developed in the companion paper are used to model detailed observations of preseismic slip and the onset of unstable slip in biaxial laboratory experiments. The simulations employ a deterministic plane strain finite element model to represent the interactions both within the sliding blocks and between the blocks and the loading apparatus. Both experiments and simulations show that preseismic slip is controlled by initial inhomogeneity of shear stress along the sliding surface relative to the frictional strength. As a consequence of the inhomogeneity, stable slip begins at a point on the surface and the area of slip slowly expands as the external loading increases. A previously proposed correlation between accelerating rates of stable slip and growth of the area of slip is supported by the simulations. In the simulations and in the experiments, unstable slip occurs shortly after a propagating slip event traverses the sliding surface and breaks out at the ends of the sample. In the model the breakout of stable slip causes a sudden acceleration of slip rates. Because of velocity dependency of the constitutive relationship for friction, the rapid acceleration of slip causes a decrease in frictional strength. Instability occurs when the frictional strength decreases with displacement at a rate that exceeds the intrinsic unloading characteristics of the sample and test machine. A simple slider-spring model that does not consider preseismic slip appears to approximate the transition adequately from stable sliding to unstable slip as a function of normal stress, machine stiffness, and surface roughness for small samples. However, for large samples and for natural faults the simulations suggest that the simple model may be inaccurate because it does not take into account potentially large preseismic displacements that will alter the friction parameters prior to instability. Copyright ?? 1979 by the American Geophysical Union.

Linear Combinations of Multiple Outcome Measures to Improve the Power of Efficacy Analysis ---Application to Clinical Trials on Early Stage Alzheimer Disease

PubMed Central

Xiong, Chengjie; Luo, Jingqin; Morris, John C; Bateman, Randall

2018-01-01

Modern clinical trials on Alzheimer disease (AD) focus on the early symptomatic stage or even the preclinical stage. Subtle disease progression at the early stages, however, poses a major challenge in designing such clinical trials. We propose a multivariate mixed model on repeated measures to model the disease progression over time on multiple efficacy outcomes, and derive the optimum weights to combine multiple outcome measures by minimizing the sample sizes to adequately power the clinical trials. A cross-validation simulation study is conducted to assess the accuracy for the estimated weights as well as the improvement in reducing the sample sizes for such trials. The proposed methodology is applied to the multiple cognitive tests from the ongoing observational study of the Dominantly Inherited Alzheimer Network (DIAN) to power future clinical trials in the DIAN with a cognitive endpoint. Our results show that the optimum weights to combine multiple outcome measures can be accurately estimated, and that compared to the individual outcomes, the combined efficacy outcome with these weights significantly reduces the sample size required to adequately power clinical trials. When applied to the clinical trial in the DIAN, the estimated linear combination of six cognitive tests can adequately power the clinical trial. PMID:29546251

Construction of anthropomorphic hybrid, dual-lattice voxel models for optimizing image quality and dose in radiography

NASA Astrophysics Data System (ADS)

Petoussi-Henss, Nina; Becker, Janine; Greiter, Matthias; Schlattl, Helmut; Zankl, Maria; Hoeschen, Christoph

2014-03-01

In radiography there is generally a conflict between the best image quality and the lowest possible patient dose. A proven method of dosimetry is the simulation of radiation transport in virtual human models (i.e. phantoms). However, while the resolution of these voxel models is adequate for most dosimetric purposes, they cannot provide the required organ fine structures necessary for the assessment of the imaging quality. The aim of this work is to develop hybrid/dual-lattice voxel models (called also phantoms) as well as simulation methods by which patient dose and image quality for typical radiographic procedures can be determined. The results will provide a basis to investigate by means of simulations the relationships between patient dose and image quality for various imaging parameters and develop methods for their optimization. A hybrid model, based on NURBS (Non Linear Uniform Rational B-Spline) and PM (Polygon Mesh) surfaces, was constructed from an existing voxel model of a female patient. The organs of the hybrid model can be then scaled and deformed in a non-uniform way i.e. organ by organ; they can be, thus, adapted to patient characteristics without losing their anatomical realism. Furthermore, the left lobe of the lung was substituted by a high resolution lung voxel model, resulting in a dual-lattice geometry model. "Dual lattice" means in this context the combination of voxel models with different resolution. Monte Carlo simulations of radiographic imaging were performed with the code EGS4nrc, modified such as to perform dual lattice transport. Results are presented for a thorax examination.

Simulating flight boundary conditions for orbiter payload modal survey

NASA Technical Reports Server (NTRS)

Chung, Y. T.; Sernaker, M. L.; Peebles, J. H.

1993-01-01

An approach to simulate the characteristics of the payload/orbiter interfaces for the payload modal survey was developed. The flexure designed for this approach is required to provide adequate stiffness separation in the free and constrained interface degrees of freedom to closely resemble the flight boundary condition. Payloads will behave linearly and demonstrate similar modal effective mass distribution and load path as the flight if the flexure fixture is used for the payload modal survey. The potential non-linearities caused by the trunnion slippage during the conventional fixed base modal survey may be eliminated. Consequently, the effort to correlate the test and analysis models can be significantly reduced. An example is given to illustrate the selection and the sensitivity of the flexure stiffness. The advantages of using flexure fixtures for the modal survey and for the analytical model verification are also demonstrated.

The role of the antecedent soil moisture condition on the distributed hydrologic modelling of the Toce alpine basin floods.

NASA Astrophysics Data System (ADS)

Ravazzani, G.; Montaldo, N.; Mancini, M.; Rosso, R.

2003-04-01

Event-based hydrologic models need the antecedent soil moisture condition, as critical boundary initial condition for flood simulation. Land-surface models (LSMs) have been developed to simulate mass and energy transfers, and to update the soil moisture condition through time from the solution of water and energy balance equations. They are recently used in distributed hydrologic modeling for flood prediction systems. Recent developments have made LSMs more complex by inclusion of more processes and controlling variables, increasing parameter number and uncertainty of their estimates. This also led to increasing of computational burden and parameterization of the distributed hydrologic models. In this study we investigate: 1) the role of soil moisture initial conditions in the modeling of Alpine basin floods; 2) the adequate complexity level of LSMs for the distributed hydrologic modeling of Alpine basin floods. The Toce basin is the case study; it is located in the North Piedmont (Italian Alps), and it has a total drainage area of 1534 km2 at Candoglia section. Three distributed hydrologic models of different level of complexity are developed and compared: two (TDLSM and SDLSM) are continuous models, one (FEST02) is an event model based on the simplified SCS-CN method for rainfall abstractions. In the TDLSM model a two-layer LSM computes both saturation and infiltration excess runoff, and simulates the evolution of the water table spatial distribution using the topographic index; in the SDLSM model a simplified one-layer distributed LSM only computes hortonian runoff, and doesn’t simulate the water table dynamic. All the three hydrologic models simulate the surface runoff propagation through the Muskingum-Cunge method. TDLSM and SDLSM models have been applied for the two-year (1996 and 1997) simulation period, during which two major floods occurred in the November 1996 and in the June 1997. The models have been calibrated and tested comparing simulated and observed hydrographs at Candoglia. Sensitivity analysis of the models to significant LSM parameters were also performed. The performances of the three models in the simulation of the two major floods are compared. Interestingly, the results indicate that the SDLSM model is able to sufficiently well predict the major floods of this Alpine basin; indeed, this model is a good compromise between the over-parameterized and too complex TDLSM model and the over-simplified FEST02 model.

Material Characterization and Computer Model Simulation of Low Density Polyurethane Foam Used in a Rodent Traumatic Brain Injury Model

PubMed Central

Zhang, Liying; Gurao, Manish; Yang, King H.; King, Albert I.

2011-01-01

Computer models of the head can be used to simulate the events associated with traumatic brain injury (TBI) and quantify biomechanical response within the brain. Marmarou’s impact acceleration rodent model is a widely used experimental model of TBI mirroring axonal pathology in humans. The mechanical properties of the low density polyurethane (PU) foam, an essential piece of energy management used in Marmarou’s impact device, has not been fully characterized. The foam used in Marmarou’s device was tested at seven strain rates ranging from quasi-static to dynamic (0.014 ~ 42.86 s−1) to quantify the stress-strain relationships in compression. Recovery rate of the foam after cyclic compression was also determined through the periods of recovery up to three weeks. The experimentally determined stress-strain curves were incorporated into a material model in an explicit Finite Element (FE) solver to validate the strain rate dependency of the FE foam model. Compression test results have shown that the foam used in the rodent impact acceleration model is strain rate dependent. The foam has been found to be reusable for multiple impacts. However the stress resistance of used foam is reduced to 70% of the new foam. The FU_CHANG_FOAM material model in an FE solver has been found to be adequate to simulate this rate sensitive foam. PMID:21459114

Material characterization and computer model simulation of low density polyurethane foam used in a rodent traumatic brain injury model.

PubMed

Zhang, Liying; Gurao, Manish; Yang, King H; King, Albert I

2011-05-15

Computer models of the head can be used to simulate the events associated with traumatic brain injury (TBI) and quantify biomechanical response within the brain. Marmarou's impact acceleration rodent model is a widely used experimental model of TBI mirroring axonal pathology in humans. The mechanical properties of the low density polyurethane (PU) foam, an essential piece of energy management used in Marmarou's impact device, has not been fully characterized. The foam used in Marmarou's device was tested at seven strain rates ranging from quasi-static to dynamic (0.014-42.86 s⁻¹) to quantify the stress-strain relationships in compression. Recovery rate of the foam after cyclic compression was also determined through the periods of recovery up to three weeks. The experimentally determined stress-strain curves were incorporated into a material model in an explicit Finite Element (FE) solver to validate the strain rate dependency of the FE foam model. Compression test results have shown that the foam used in the rodent impact acceleration model is strain rate dependent. The foam has been found to be reusable for multiple impacts. However the stress resistance of used foam is reduced to 70% of the new foam. The FU_CHANG_FOAM material model in an FE solver has been found to be adequate to simulate this rate sensitive foam. Copyright © 2011 Elsevier B.V. All rights reserved.

Structure of dysprosium(111) dl-tartrate dimer in aqueous solution

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

Chevela, V.V.; Vulfson, S.G.; Salnikov, Y.I.

1994-10-01

The paramagnetic birefringence method was supplemented by numerical simulation to determine the molar paramagnetic-birefringence constant of the dysprosium dl-tartrate dimer Dy{sub 2}(d-L)(l-L){sup 2-} (I), where d-L{sup 4-} and l-L{sup 4-} are the deprotonated d- and l-tartaric acid molecules, respectively. The structure of the ligand and hydration surroundings of I was modeled by molecular mechanic calculations (the Dashevskii-Pylamovatyi model). It is shown that adequate results can be obtained only if one takes into account the coordination of I to the Na{sup +} ion.

Report of NPSAT1 Battery Thermal Contact Resistance Testing, Modeling and Simulation

DTIC Science & Technology

2012-10-01

lithium ion battery is the spacecraft component with the smallest temperature range of 0?C to 45?C during operation. Thermal analysis results, however, can only provide adequate results if there is sufficient fidelity in thermal modeling. Arguably, the values used in defining thermal coupling for components are the most difficult to estimate because of the many variables that define them. This document describes the work performed by the authors starting in the 2012 winter quarter as part of the SS3900 directed study course. The objectives of the study were to

Simulation of the summer circulation over South America by two regional climate models. Part I: Mean climatology

NASA Astrophysics Data System (ADS)

Fernandez, J. P. R.; Franchito, S. H.; Rao, V. B.

2006-09-01

This study investigates the capabilities of two regional models (the ICTP RegCM3 and the climate version of the CPTEC Eta model - EtaClim) in simulating the mean climatological features of the summer quasi-stationary circulations over South America. Comparing the results with the NCEP/DOE reanalysis II data it is seen that the RegCM3 simulates a weaker and southward shifted Bolivian high (BH). But, the Nordeste low (NL) is located close to its climatological position. In the EtaClim the position of the BH is reproduced well, but the NL is shifted towards the interior of the continent. To the east of Andes, the RegCM3 simulates a weaker low level jet and a weaker basic flow from the tropical Atlantic to Amazonia while they are stronger in the EtaClim. In general, the RegCM3 and EtaClim show, respectively a negative and positive bias in the surface temperature in almost all regions of South America. For both models, the correlation coefficients between the simulated precipitation and the GPCP data are high over most of South America. Although the RegCM3 and EtaClim overestimate the precipitation in the Andes region they show a negative bias in general over the entire South America. The simulations of upper and lower level circulations and precipitation fields in EtaClim were better than that of the RegCM3. In central Amazonia both models were unable to simulate the precipitation correctly. The results showed that although the RegCM3 and EtaClim are capable of simulating the main climatological features of the summer climate over South America, there are areas which need improvement. This indicates that the models must be more adequately tuned in order to give reliable predictions in the different regions of South America.

Calculation for simulation of archery goal value using a web camera and ultrasonic sensor

NASA Astrophysics Data System (ADS)

Rusjdi, Darma; Abdurrasyid, Wulandari, Dewi Arianti

2017-08-01

Development of the device simulator digital indoor archery-based embedded systems as a solution to the limitations of the field or open space is adequate, especially in big cities. Development of the device requires simulations to calculate the value of achieving the target based on the approach defined by the parabolic motion variable initial velocity and direction of motion of the arrow reaches the target. The simulator device should be complemented with an initial velocity measuring device using ultrasonic sensors and measuring direction of the target using a digital camera. The methodology uses research and development of application software from modeling and simulation approach. The research objective to create simulation applications calculating the value of the achievement of the target arrows. Benefits as a preliminary stage for the development of the simulator device of archery. Implementation of calculating the value of the target arrows into the application program generates a simulation game of archery that can be used as a reference development of the digital archery simulator in a room with embedded systems using ultrasonic sensors and web cameras. Applications developed with the simulation calculation comparing the outer radius of the circle produced a camera from a distance of three meters.

Population pharmacokinetics of recombinant coagulation factor VIII-SingleChain in patients with severe hemophilia A.

PubMed

Zhang, Y; Roberts, J; Tortorici, M; Veldman, A; St Ledger, K; Feussner, A; Sidhu, J

2017-06-01

Essentials rVIII-SingleChain is a unique recombinant factor VIII (FVIII) molecule. A population pharmacokinetic model was based on FVIII activity of severe hemophilia A patients. The model was used to simulate factor VIII activity-time profiles for various dosing scenarios. The model supports prolonged dosing of rVIII-SingleChain with intervals of up to twice per week. Background Single-chain recombinant coagulation factor VIII (rVIII-SingleChain) is a unique recombinant coagulation factor VIII molecule. Objectives To: (i) characterize the population pharmacokinetics (PK) of rVIII-SingleChain in patients with severe hemophilia A; (ii) identify correlates of variability in rVIII-SingleChain PK; and (iii) simulate various dosing scenarios of rVIII-SingleChain. Patients/Methods A population PK model was developed, based on FVIII activity levels of 130 patients with severe hemophilia A (n = 91 for ≥ 12-65 years; n = 39 for < 12 years) who had participated in a single-dose PK investigation with rVIII-SingleChain 50 IU kg -1 . PK sampling was performed for up to 96 h. Results A two-compartment population PK model with first-order elimination adequately described FVIII activity. Body weight and predose level of von Willebrand factor were significant covariates on clearance, and body weight was a significant covariate on the central distribution volume. Simulations using the model with various dosing scenarios estimated that > 85% and > 93% of patients were predicted to maintain FVIII activity level above 1 IU dL -1 , at all times with three-times-weekly dosing (given on days 0, 2, and 4.5) at the lowest (20 IU kg -1 ) and highest (50 IU kg -1 ) doses, respectively. For twice weekly dosing (days 0 and 3.5) of 50 IU kg -1 rVIII-SingleChain, 62-80% of patients across all ages were predicted to maintain a FVIII activity level above 1 IU dL -1 at day 7. Conclusions The population PK model adequately characterized rVIII-SingleChain PK, and the model can be utilized to simulate FVIII activity-time profiles for various dosing scenarios. © 2017 The Authors. Journal of Thrombosis and Haemostasis published by Wiley Periodicals, Inc. on behalf of International Society on Thrombosis and Haemostasis.

A computer model of solar panel-plasma interactions

NASA Technical Reports Server (NTRS)

Cooke, D. L.; Freeman, J. W.

1980-01-01

High power solar arrays for satellite power systems are presently being planned with dimensions of kilometers, and with tens of kilovolts distributed over their surface. Such systems face many plasma interaction problems, such as power leakage to the plasma, particle focusing, and anomalous arcing. These effects cannot be adequately modeled without detailed knowledge of the plasma sheath structure and space charge effects. Laboratory studies of 1 by 10 meter solar array in a simulated low Earth orbit plasma are discussed. The plasma screening process is discussed, program theory is outlined, and a series of calibration models is presented. These models are designed to demonstrate that PANEL is capable of accurate self consistant space charge calculations. Such models include PANEL predictions for the Child-Langmuir diode problem.

Cell population modelling of yeast glycolytic oscillations.

PubMed Central

Henson, Michael A; Müller, Dirk; Reuss, Matthias

2002-01-01

We investigated a cell-population modelling technique in which the population is constructed from an ensemble of individual cell models. The average value or the number distribution of any intracellular property captured by the individual cell model can be calculated by simulation of a sufficient number of individual cells. The proposed method is applied to a simple model of yeast glycolytic oscillations where synchronization of the cell population is mediated by the action of an excreted metabolite. We show that smooth one-dimensional distributions can be obtained with ensembles comprising 1000 individual cells. Random variations in the state and/or structure of individual cells are shown to produce complex dynamic behaviours which cannot be adequately captured by small ensembles. PMID:12206713

Development of a model for on-line control of crystal growth by the AHP method

NASA Astrophysics Data System (ADS)

Gonik, M. A.; Lomokhova, A. V.; Gonik, M. M.; Kuliev, A. T.; Smirnov, A. D.

2007-05-01

The possibility to apply a simplified 2D model for heat transfer calculations in crystal growth by the axial heat close to phase interface (AHP) method is discussed in this paper. A comparison with global heat transfer calculations with the CGSim software was performed to confirm the accuracy of this model. The simplified model was shown to provide adequate results for the shape of the melt-crystal interface and temperature field in an opaque (Ge) and a transparent crystal (CsI:Tl). The model proposed is used for identification of the growth setup as a control object, for synthesis of a digital controller (PID controller at the present stage) and, finally, in on-line simulations of crystal growth control.

Implementing Dynamic Root Optimization in Noah-MP for Simulating Phreatophytic Root Water Uptake

NASA Astrophysics Data System (ADS)

Wang, Ping; Niu, Guo-Yue; Fang, Yuan-Hao; Wu, Run-Jian; Yu, Jing-Jie; Yuan, Guo-Fu; Pozdniakov, Sergey P.; Scott, Russell L.

2018-03-01

Widely distributed in arid and semiarid regions, phreatophytic roots extend into the saturated zone and extract water directly from groundwater. In this paper, we implemented a vegetation optimality model of root dynamics (VOM-ROOT) in the Noah land surface model with multiparameterization options (Noah-MP LSM) to model the extraction of groundwater through phreatophytic roots at a riparian site with a hyperarid climate (with precipitation of 35 mm/yr) in northwestern China. VOM-ROOT numerically describes the natural optimization of the root profile in response to changes in subsurface water conditions. The coupled Noah-MP/VOM-ROOT model substantially improves the simulation of surface energy and water fluxes, particularly during the growing season, compared to the prescribed static root profile in the default Noah-MP. In the coupled model, more roots are required to grow into the saturated zone to meet transpiration demand when the groundwater level declines over the growing season. The modeling results indicate that at the study site, the modeled annual transpiration is 472 mm, accounting for 92.3% of the total evapotranspiration. Direct root water uptake from the capillary fringe and groundwater, which is supplied by lateral groundwater flow, accounts for approximately 84% of the total transpiration. This study demonstrates the importance of implementing a dynamic root scheme in a land surface model for adequately simulating phreatophytic root water uptake and the associated latent heat flux.

Pronounced differences between observed and CMIP5-simulated multidecadal climate variability in the twentieth century

NASA Astrophysics Data System (ADS)

Kravtsov, Sergey

2017-06-01

Identification and dynamical attribution of multidecadal climate undulations to either variations in external forcings or to internal sources is one of the most important topics of modern climate science, especially in conjunction with the issue of human-induced global warming. Here we utilize ensembles of twentieth century climate simulations to isolate the forced signal and residual internal variability in a network of observed and modeled climate indices. The observed internal variability so estimated exhibits a pronounced multidecadal mode with a distinctive spatiotemporal signature, which is altogether absent in model simulations. This single mode explains a major fraction of model-data differences over the entire climate index network considered; it may reflect either biases in the models' forced response or models' lack of requisite internal dynamics, or a combination of both.