Long-term in vitro degradation of PDLLA/bioglass bone scaffolds in acellular simulated body fluid.
Blaker, J J; Nazhat, S N; Maquet, V; Boccaccini, A R
2011-02-01
The long-term (600days) in vitro degradation of highly porous poly(D,L-lactide) (PDLLA)/Bioglass-filled composite foams developed for bone tissue engineering scaffolds has been investigated in simulated body fluid (SBF). Foams of ∼93% porosity were produced by thermally induced phase separation (TIPS). The degradation profile for foams of neat PDLLA and the influence of Bioglass addition were comprehensively assessed in terms of changes in dimensional stability, pore morphology, weight loss, molecular weight and mechanical properties (dry and wet states). It is shown that the degradation process proceeded in several stages: (a) a quasi-stable stage, where water absorption and plasticization occurred together with weight loss due to Bioglass particle loss and dissolution, resulting in decreased wet mechanical properties; (b) a stage showing a slight increase in the wet mechanical properties and a moderate decrease in dimensions, with the properties remaining moderately constant until the onset of significant weight loss, whilst molecular weight continued to decrease; (c) an end stage of massive weight loss, disruption of the pore structure and the formation of blisters and embrittlement of the scaffold (evident on handling). The findings from this long-term in vitro degradation investigation underpin studies that have been and continue to be performed on highly porous poly(α-hydroxyesters) scaffolds filled with bioactive glasses for bone tissue engineering applications. PMID:20849987
Electronic Falling Body Simulator
ERIC Educational Resources Information Center
Goodman, John M.
1975-01-01
Describes an analog electronic simulator for the motion of an object projected horizontally in a vertical acceleration field. The device features adjustable values for the horizontal speeds and "gravity." Trajectories are displayed on an oscilloscope. (Author/CP)
Cosmological N-body Simulation
NASA Astrophysics Data System (ADS)
Lake, George
1994-05-01
.90ex> }}} The ``N'' in N-body calculations has doubled every year for the last two decades. To continue this trend, the UW N-body group is working on algorithms for the fast evaluation of gravitational forces on parallel computers and establishing rigorous standards for the computations. In these algorithms, the computational cost per time step is ~ 10(3) pairwise forces per particle. A new adaptive time integrator enables us to perform high quality integrations that are fully temporally and spatially adaptive. SPH--smoothed particle hydrodynamics will be added to simulate the effects of dissipating gas and magnetic fields. The importance of these calculations is two-fold. First, they determine the nonlinear consequences of theories for the structure of the Universe. Second, they are essential for the interpretation of observations. Every galaxy has six coordinates of velocity and position. Observations determine two sky coordinates and a line of sight velocity that bundles universal expansion (distance) together with a random velocity created by the mass distribution. Simulations are needed to determine the underlying structure and masses. The importance of simulations has moved from ex post facto explanation to an integral part of planning large observational programs. I will show why high quality simulations with ``large N'' are essential to accomplish our scientific goals. This year, our simulations have N >~ 10(7) . This is sufficient to tackle some niche problems, but well short of our 5 year goal--simulating The Sloan Digital Sky Survey using a few Billion particles (a Teraflop-year simulation). Extrapolating past trends, we would have to ``wait'' 7 years for this hundred-fold improvement. Like past gains, significant changes in the computational methods are required for these advances. I will describe new algorithms, algorithmic hacks and a dedicated computer to perform Billion particle simulations. Finally, I will describe research that can be enabled by
Sato, Y; Sato, H
1991-01-01
There is no suitable animal model for pertussis encephalopathy in humans. In this study, we have compared the toxicity of acellular pertussis vaccine with whole cell pertussis vaccine in mice or guinea pigs. Two lots of acellular and two lots of whole cell vaccine produced in different countries were assayed in the test. 1. There was no statistical difference in mouse protective potency between these acellular or whole cell pertussis vaccines. 2. There were no differences in chemical ingredients between acellular and whole cell pertussis vaccines except for protein nitrogen content. The protein nitrogen content of whole cell vaccine was at least three times higher than that of the acellular product. 3. Anti-PT antibody productivity of the acellular vaccine was higher than that of the whole cell vaccine. 4. Anti-agglutinogen antibody productivity of the whole cell vaccine was higher than that of the acellular vaccine. 5. There was no pyrogenic activity with the acellular vaccine, but high pyrogenicity was seen with whole cell vaccine. 6. There was high body-weight decreasing toxicity in mice and guinea pigs by the whole cell vaccine. 7. The mice died when they received whole cell pertussis vaccine iv, but no deaths occurred in the mice which received acellular pertussis vaccine. PMID:1778317
Laboratory Simulation Of Upper-Body Work
NASA Technical Reports Server (NTRS)
Lantz, R.; Vykukal, H.; Webbon, B.
1992-01-01
Paper describes exercise method evoking muscular, cardiovascular, respiratory, and thermoregulatory responses typical of activity involving upper body, to simulate effects of working in zero gravity in space suit. Used in research on thermoregulatory subsystem of advanced portable life-support system for space suit. Discusses requirements for simulation of zero-gravity work in space suit and describes evolution of method through three versions. In preferred version, one performs arm-crank ergometry with gimbaled body-support mechanism that react forces at feet.
A simulation method for the fruitage body
NASA Astrophysics Data System (ADS)
Lu, Ling; Song, Weng-lin; Wang, Lei
2009-07-01
An effective visual modeling for creating the fruitage body has been present. According to the geometry shape character of fruitage, we build up its face model base on ellipsoid deformation. The face model is relation with radius. We consider different radius become a face in the fruitage, and uses same method to simulate the shape of fruitage inside. The body model is formed by combine face model and radius direction. Our method can simulate virtual inter and outer structure for fruitage body. The method decreases a lot of data and increases display speed. Another, the texture model of fruitage is defined by sum of different base function. This kind of method is simple and speed. We show the feasibility of our method by creating a winter-jujube and an apricot. They include exocorp, mesocorp and endocarp. It is useful that develop virtual plant.
Robust three-body water simulation model
NASA Astrophysics Data System (ADS)
Tainter, C. J.; Pieniazek, P. A.; Lin, Y.-S.; Skinner, J. L.
2011-05-01
The most common potentials used in classical simulations of liquid water assume a pairwise additive form. Although these models have been very successful in reproducing many properties of liquid water at ambient conditions, none is able to describe accurately water throughout its complicated phase diagram. The primary reason for this is the neglect of many-body interactions. To this end, a simulation model with explicit three-body interactions was introduced recently [R. Kumar and J. L. Skinner, J. Phys. Chem. B 112, 8311 (2008), 10.1021/jp8009468]. This model was parameterized to fit the experimental O-O radial distribution function and diffusion constant. Herein we reparameterize the model, fitting to a wider range of experimental properties (diffusion constant, rotational correlation time, density for the liquid, liquid/vapor surface tension, melting point, and the ice Ih density). The robustness of the model is then verified by comparing simulation to experiment for a number of other quantities (enthalpy of vaporization, dielectric constant, Debye relaxation time, temperature of maximum density, and the temperature-dependent second and third virial coefficients), with good agreement.
Realistic Simulation for Body Area and Body-To-Body Networks
Alam, Muhammad Mahtab; Ben Hamida, Elyes; Ben Arbia, Dhafer; Maman, Mickael; Mani, Francesco; Denis, Benoit; D’Errico, Raffaele
2016-01-01
In this paper, we present an accurate and realistic simulation for body area networks (BAN) and body-to-body networks (BBN) using deterministic and semi-deterministic approaches. First, in the semi-deterministic approach, a real-time measurement campaign is performed, which is further characterized through statistical analysis. It is able to generate link-correlated and time-varying realistic traces (i.e., with consistent mobility patterns) for on-body and body-to-body shadowing and fading, including body orientations and rotations, by means of stochastic channel models. The full deterministic approach is particularly targeted to enhance IEEE 802.15.6 proposed channel models by introducing space and time variations (i.e., dynamic distances) through biomechanical modeling. In addition, it helps to accurately model the radio link by identifying the link types and corresponding path loss factors for line of sight (LOS) and non-line of sight (NLOS). This approach is particularly important for links that vary over time due to mobility. It is also important to add that the communication and protocol stack, including the physical (PHY), medium access control (MAC) and networking models, is developed for BAN and BBN, and the IEEE 802.15.6 compliance standard is provided as a benchmark for future research works of the community. Finally, the two approaches are compared in terms of the successful packet delivery ratio, packet delay and energy efficiency. The results show that the semi-deterministic approach is the best option; however, for the diversity of the mobility patterns and scenarios applicable, biomechanical modeling and the deterministic approach are better choices. PMID:27104537
Realistic Simulation for Body Area and Body-To-Body Networks.
Alam, Muhammad Mahtab; Ben Hamida, Elyes; Ben Arbia, Dhafer; Maman, Mickael; Mani, Francesco; Denis, Benoit; D'Errico, Raffaele
2016-01-01
In this paper, we present an accurate and realistic simulation for body area networks (BAN) and body-to-body networks (BBN) using deterministic and semi-deterministic approaches. First, in the semi-deterministic approach, a real-time measurement campaign is performed, which is further characterized through statistical analysis. It is able to generate link-correlated and time-varying realistic traces (i.e., with consistent mobility patterns) for on-body and body-to-body shadowing and fading, including body orientations and rotations, by means of stochastic channel models. The full deterministic approach is particularly targeted to enhance IEEE 802.15.6 proposed channel models by introducing space and time variations (i.e., dynamic distances) through biomechanical modeling. In addition, it helps to accurately model the radio link by identifying the link types and corresponding path loss factors for line of sight (LOS) and non-line of sight (NLOS). This approach is particularly important for links that vary over time due to mobility. It is also important to add that the communication and protocol stack, including the physical (PHY), medium access control (MAC) and networking models, is developed for BAN and BBN, and the IEEE 802.15.6 compliance standard is provided as a benchmark for future research works of the community. Finally, the two approaches are compared in terms of the successful packet delivery ratio, packet delay and energy efficiency. The results show that the semi-deterministic approach is the best option; however, for the diversity of the mobility patterns and scenarios applicable, biomechanical modeling and the deterministic approach are better choices. PMID:27104537
N-body simulations of γ gravity
NASA Astrophysics Data System (ADS)
Vargas dos Santos, Marcelo; Winther, Hans A.; Mota, David F.; Waga, Ioav
2016-03-01
We have investigated structure formation in the γ gravity f(R) model with N-body simulations. The γ gravity model is a proposal which, unlike other viable f(R) models, not only changes the gravitational dynamics, but can in principle also have signatures at the background level that are different from those obtained in ΛCDM (Cosmological constant, Cold Dark Matter). The aim of this paper is to study the nonlinear regime of the model in the case where, at late times, the background differs from ΛCDM. We quantify the signatures produced on the power spectrum, the halo mass function, and the density and velocity profiles. To appreciate the features of the model, we have compared it to ΛCDM and the Hu-Sawicki f(R) models. For the considered set of parameters we find that the screening mechanism is ineffective, which gives rise to deviations in the halo mass function that disagree with observations. This does not rule out the model per se, but requires choices of parameters such that | fR0 | is much smaller, which would imply that its cosmic expansion history cannot be distinguished from ΛCDM at the background level.
Whole body measurement systems. [for weightlessness simulation
NASA Technical Reports Server (NTRS)
Ogle, J. S. (Inventor)
1973-01-01
A system for measuring the volume and volume variations of a human body under zero gravity conditions is disclosed. An enclosed chamber having a defined volume and arranged for receiving a human body is provided with means for infrasonically varying the volume of the chamber. The changes in volume produce resultant changes in pressure, and under substantially isentropic conditions, an isentropic relationship permits a determination of gas volume which, in turn, when related to total chamber volume permits a determination of the body volume. By comparison techniques, volume changes of a human independent of gravity conditions can be determined.
In vitro assessment of biodurability: acellular systems.
de Meringo, A; Morscheidt, C; Thélohan, S; Tiesler, H
1994-01-01
The assessment of biodurability of man-made vitreous fibers is essential to the limitation of health hazards associated with human exposure to environments in which respirable fibers are present. In vitro acellular systems provide effective test methods of measuring fiber solubility provided care is taken to select the most suitable solvent and test conditions for the specific fiber type and dimension. PMID:7882955
Human acellular dermal wound matrix: evidence and experience.
Kirsner, Robert S; Bohn, Greg; Driver, Vickie R; Mills, Joseph L; Nanney, Lillian B; Williams, Marie L; Wu, Stephanie C
2015-12-01
A chronic wound fails to complete an orderly and timely reparative process and places patients at increased risk for wound complications that negatively impact quality of life and require greater health care expenditure. The role of extracellular matrix (ECM) is critical in normal and chronic wound repair. Not only is ECM the largest component of the dermal skin layer, but also ECM proteins provide structure and cell signalling that are necessary for successful tissue repair. Chronic wounds are characterised by their inflammatory and proteolytic environment, which degrades the ECM. Human acellular dermal matrices, which provide an ECM scaffold, therefore, are being used to treat chronic wounds. The ideal human acellular dermal wound matrix (HADWM) would support regenerative healing, providing a structure that could be repopulated by the body's cells. Experienced wound care investigators and clinicians discussed the function of ECM, the evidence related to a specific HADWM (Graftjacket(®) regenerative tissue matrix, Wright Medical Technology, Inc., licensed by KCI USA, Inc., San Antonio, TX), and their clinical experience with this scaffold. This article distills these discussions into an evidence-based and practical overview for treating chronic lower extremity wounds with this HADWM. PMID:24283346
Wagner, Darcy E.; Fenn, Spencer L.; Bonenfant, Nicholas R.; Marks, Elliot R.; Borg, Zachary; Saunders, Patrick; Oldinski, Rachael A.; Weiss, Daniel J.
2015-01-01
Whole organ decellularization of complex organs, such as lungs, presents a unique opportunity for use of acellular scaffolds for ex vivo tissue engineering or for studying cell-extracellular matrix interactions ex vivo. A growing body of literature investigating decellularizing and recellularizing rodent lungs has provided important proof of concept models and rodent lungs are readily available for high throughput studies. In contrast, comparable progress in large animal and human lungs has been impeded owing to more limited availability and difficulties in handling larger tissue. While the use of smaller segments of acellular large animal or human lungs would maximize usage from a single lung, excision of small acellular segments compromises the integrity of the pleural layer, leaving the terminal ends of blood vessels and airways exposed. We have developed a novel pleural coating using non-toxic ionically crosslinked alginate or photocrosslinked methacrylated alginate which can be applied to excised acellular lung segments, permits inflation of small segments, and significantly enhances retention of cells inoculated through cannulated airways or blood vessels. Further, photocrosslinking methacrylated alginate, using eosin Y and triethanolamine (TEOA) at 530nm wavelength, results in a mechanically stable pleural coating that permits effective cyclic 3-dimensional stretch, i.e. mechanical ventilation, of individual segments. PMID:25750684
A General Simulation Method for Multiple Bodies in Proximate Flight
NASA Technical Reports Server (NTRS)
Meakin, Robert L.
2003-01-01
Methods of unsteady aerodynamic simulation for an arbitrary number of independent bodies flying in close proximity are considered. A novel method to efficiently detect collision contact points is described. A method to compute body trajectories in response to aerodynamic loads, applied loads, and inter-body collisions is also given. The physical correctness of the methods are verified by comparison to a set of analytic solutions. The methods, combined with a Navier-Stokes solver, are used to demonstrate the possibility of predicting the unsteady aerodynamics and flight trajectories of moving bodies that involve rigid-body collisions.
Cartilage oligomeric matrix protein enhances the vascularization of acellular nerves
Cui, Wei-ling; Qiu, Long-hai; Lian, Jia-yan; Li, Jia-chun; Hu, Jun; Liu, Xiao-lin
2016-01-01
Vascularization of acellular nerves has been shown to contribute to nerve bridging. In this study, we used a 10-mm sciatic nerve defect model in rats to determine whether cartilage oligomeric matrix protein enhances the vascularization of injured acellular nerves. The rat nerve defects were treated with acellular nerve grafting (control group) alone or acellular nerve grafting combined with intraperitoneal injection of cartilage oligomeric matrix protein (experimental group). As shown through two-dimensional imaging, the vessels began to invade into the acellular nerve graft from both anastomotic ends at day 7 post-operation, and gradually covered the entire graft at day 21. The vascular density, vascular area, and the velocity of revascularization in the experimental group were all higher than those in the control group. These results indicate that cartilage oligomeric matrix protein enhances the vascularization of acellular nerves. PMID:27127495
N-body simulations of disks. [of stellar systems
NASA Technical Reports Server (NTRS)
Hohl, F.
1975-01-01
The methods used in large-scale n-body simulations are discussed. However, the present review concentrates on results already obtained in n-body simulations using systems containing up to 200,000 simulation stars. Results are presented which show that the stability criterion developed for flattened systems applies only to truly axisymmetric instabilities. Purely stellar disks acquire rather large velocity dispersions, generally two or more times the velocity dispersion required by Toomre (1964) for axisymmetric stability. In computer simulations, the bar-forming instability can be prevented only by comparatively large velocity dispersions. However, simulations including the effects of the galactic halo and core as a fixed background field show that bar formation can be prevented for fixed halo components as large or larger than the self-consistent disk component. Experiments performed to determine the collisional relaxation time for large-scale gravitational n-body calculations show that these models are indeed 'collisionless'.
Data from acellular human heart matrix.
Sánchez, Pedro L; Fernández-Santos, M Eugenia; Espinosa, M Angeles; González-Nicolas, M Angeles; Acebes, Judith R; Costanza, Salvatore; Moscoso, Isabel; Rodríguez, Hugo; García, Julio; Romero, Jesús; Kren, Stefan M; Bermejo, Javier; Yotti, Raquel; Del Villar, Candelas Pérez; Sanz-Ruiz, Ricardo; Elizaga, Jaime; Taylor, Doris A; Fernández-Avilés, Francisco
2016-09-01
Perfusion decellularization of cadaveric hearts removes cells and generates a cell-free extracellular matrix scaffold containing acellular vascular conduits, which are theoretically sufficient to perfuse and support tissue-engineered heart constructs. This article contains additional data of our experience decellularizing and testing structural integrity and composition of a large series of human hearts, "Acellular human heart matrix: a critical step toward whole heat grafts" (Sanchez et al., 2015) [1]. Here we provide the information about the heart decellularization technique, the valve competence evaluation of the decellularized scaffolds, the integrity evaluation of epicardial and myocardial coronary circulation, the pressure volume measurements, the primers used to assess cardiac muscle gene expression and, the characteristics of donors, donor hearts, scaffolds and perfusion decellularization process. PMID:27331090
Post-Newtonian N-body simulations
NASA Astrophysics Data System (ADS)
Aarseth, Sverre J.
2007-06-01
We report on the first fully consistent conventional cluster simulation which includes terms up to the third-order post-Newtonian approximation. Numerical problems for treating extremely energetic binaries orbiting a single massive object are circumvented by employing the special `wheel-spoke' regularization method of Zare which has not been used in large-N simulations before. Idealized models containing N = 1 × 105 particles of mass 1Msolar with a central black hole (BH) of 300Msolar have been studied on GRAPE-type computers. An initial half-mass radius of rh ~= 0.1 pc is sufficiently small to yield examples of relativistic coalescence. This is achieved by significant binary shrinkage within a density cusp environment, followed by the generation of extremely high eccentricities which are induced by Kozai cycles and/or resonant relaxation. More realistic models with white dwarfs and 10 times larger half-mass radii also show evidence of general relativity effects before disruption. An experimentation with the post-Newtonian terms suggests that reducing the time-scales for activating the different orders progressively may be justified for obtaining qualitatively correct solutions without aiming for precise predictions of the final gravitational radiation wave form. The results obtained suggest that the standard loss-cone arguments underestimate the swallowing rate in globular clusters containing a central BH.
Computer simulation of multigrid body dynamics and control
NASA Technical Reports Server (NTRS)
Swaminadham, M.; Moon, Young I.; Venkayya, V. B.
1990-01-01
The objective is to set up and analyze benchmark problems on multibody dynamics and to verify the predictions of two multibody computer simulation codes. TREETOPS and DISCOS have been used to run three example problems - one degree-of-freedom spring mass dashpot system, an inverted pendulum system, and a triple pendulum. To study the dynamics and control interaction, an inverted planar pendulum with an external body force and a torsional control spring was modeled as a hinge connected two-rigid body system. TREETOPS and DISCOS affected the time history simulation of this problem. System state space variables and their time derivatives from two simulation codes were compared.
Extended Eden model reproduces growth of an acellular slime mold
NASA Astrophysics Data System (ADS)
Wagner, Geri; Halvorsrud, Ragnhild; Meakin, Paul
1999-11-01
A stochastic growth model was used to simulate the growth of the acellular slime mold Physarum polycephalum on substrates where the nutrients were confined in separate drops. Growth of Physarum on such substrates was previously studied experimentally and found to produce a range of different growth patterns [Phys. Rev. E 57, 941 (1998)]. The model represented the aging of cluster sites and differed from the original Eden model in that the occupation probability of perimeter sites depended on the time of occupation of adjacent cluster sites. This feature led to a bias in the selection of growth directions. A moderate degree of persistence was found to be crucial to reproduce the biological growth patterns under various conditions. Persistence in growth combined quick propagation in heterogeneous environments with a high probability of locating sources of nutrients.
Whole-body mathematical model for simulating intracranial pressure dynamics
NASA Technical Reports Server (NTRS)
Lakin, William D. (Inventor); Penar, Paul L. (Inventor); Stevens, Scott A. (Inventor); Tranmer, Bruce I. (Inventor)
2007-01-01
A whole-body mathematical model (10) for simulating intracranial pressure dynamics. In one embodiment, model (10) includes 17 interacting compartments, of which nine lie entirely outside of intracranial vault (14). Compartments (F) and (T) are defined to distinguish ventricular from extraventricular CSF. The vasculature of the intracranial system within cranial vault (14) is also subdivided into five compartments (A, C, P, V, and S, respectively) representing the intracranial arteries, capillaries, choroid plexus, veins, and venous sinus. The body's extracranial systemic vasculature is divided into six compartments (I, J, O, Z, D, and X, respectively) representing the arteries, capillaries, and veins of the central body and the lower body. Compartments (G) and (B) include tissue and the associated interstitial fluid in the intracranial and lower regions. Compartment (Y) is a composite involving the tissues, organs, and pulmonary circulation of the central body and compartment (M) represents the external environment.
Robotic Simulation of Flexible-Body Spacecraft Dynamics
NASA Technical Reports Server (NTRS)
Brannan, Justin C.; Carignan, Craig R.
2016-01-01
A robotic testbed has been developed to conduct hardware-in-the-loop simulations of a robotic servicer interacting with a client satellite on-orbit. By creating an analytical model of a satellite with flexible appendages, it is possible to simulate the system response to external force and torque inputs and compare the predicted system motion to a robot mass simulator outfitted with physical appendages. This validation effort includes multiple test cases that encompass the types of interaction forces a satellite might experience during a nominal on-orbit servicing mission and aims to show the simulation's ability to capture the physical system response. After incorporating the flexible-body dynamics into the robotic mass simulator at NASA Goddard Space Flight Center (GSFC), a hardware-in-the-loop simulation can be used to characterize the potential impact of structural flexibility on an end-to-end satellite servicing mission.
Discreteness noise versus force errors in N-body simulations
NASA Technical Reports Server (NTRS)
Hernquist, Lars; Hut, Piet; Makino, Jun
1993-01-01
A low accuracy in the force calculation per time step of a few percent for each particle pair is sufficient for collisionless N-body simulations. Higher accuracy is made meaningless by the dominant discreteness noise in the form of two-body relaxation, which can be reduced only by increasing the number of particles. Since an N-body simulation is a Monte Carlo procedure in which each particle-particle force is essentially random, i.e., carries an error of about 1000 percent, the only requirement is a systematic averaging-out of these intrinsic errors. We illustrate these assertions with two specific examples in which individual pairwise forces are deliberately allowed to carry significant errors: tree-codes on supercomputers and algorithms on special-purpose machines with low-precision hardware.
N-Body Simulations of Disk Galaxy Interactions
NASA Astrophysics Data System (ADS)
Myers, Jeannette; Messick, Garrett; The, Lih-Sin
2004-11-01
We present results of N-body simulations for two interacting disk galaxies. The simulated galaxies are initially modeled with parameters fitting observations of the Milky Way Galaxy. We use a four-component model of a disk, dark matter halo, spheroid, and central concentration. In our present simulations, we neglect the gas component and associated star formation. Among our analysis, we address the question of member transference from one galaxy to another. We also examine the rate of member ejections from each disk system and final configurations after the disk interactions.
Acellular organ scaffolds for tumor tissue engineering
NASA Astrophysics Data System (ADS)
Guller, Anna; Trusova, Inna; Petersen, Elena; Shekhter, Anatoly; Kurkov, Alexander; Qian, Yi; Zvyagin, Andrei
2015-12-01
Rationale: Tissue engineering (TE) is an emerging alternative approach to create models of human malignant tumors for experimental oncology, personalized medicine and drug discovery studies. Being the bottom-up strategy, TE provides an opportunity to control and explore the role of every component of the model system, including cellular populations, supportive scaffolds and signalling molecules. Objectives: As an initial step to create a new ex vivo TE model of cancer, we optimized protocols to obtain organ-specific acellular matrices and evaluated their potential as TE scaffolds for culture of normal and tumor cells. Methods and results: Effective decellularization of animals' kidneys, ureter, lungs, heart, and liver has been achieved by detergent-based processing. The obtained scaffolds demonstrated biocompatibility and growthsupporting potential in combination with normal (Vero, MDCK) and tumor cell lines (C26, B16). Acellular scaffolds and TE constructs have been characterized and compared with morphological methods. Conclusions: The proposed methodology allows creation of sustainable 3D tumor TE constructs to explore the role of organ-specific cell-matrix interaction in tumorigenesis.
Long-lived spiral waves in N-body simulations
NASA Technical Reports Server (NTRS)
Comins, Neil F.; Schroeder, Michael C.
1989-01-01
Results of N-body simulations of disc galaxies using a two-dimensional Cartesian N-body code are presented. Both trailing arm spirals (TAS) and leading arm spirals (LAS) were used with varieties of pitch angles and pattern speeds. LAS perturbations transferred their energy to TAS via swing amplification; TAS perturbations led to TAS arms. In both cases the spiral arms persisted for more than 5 rotation periods, but the maximum amplitude 2-armed spirals were generated by LAS perturbations. The persistence of the trailing arm spiral waves is thought to be caused by the kinematic spiral arm mechanism described by Kalnajs (1973).
Model reduction in the simulation of interconnected flexible bodies
NASA Technical Reports Server (NTRS)
Eke, Fidelis O.; Man, Guy K.
1988-01-01
Given the control system specifications for a system of interconnected rigid and flexible bodies, methods now exist for determining the system modes that do not interact 'strongly' with the controller. Once these important system modes are known, there still remains the problem of determining the modes of individual bodies that should be retained, since, in the final analysis, it is the modal information at the component level that must be fed into any multibody simulation code. Systematic identification of these component modes is achieved through a two-phase matrix diagonalization process starting with judiciously chosen submatrices of the system modal matrix.
Direct simulation of hypersonic flows over blunt slender bodies
NASA Technical Reports Server (NTRS)
Moss, J. N.; Cuda, V., Jr.
1986-01-01
Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.
n-body simulations using message passing parallel computers.
NASA Astrophysics Data System (ADS)
Grama, A. Y.; Kumar, V.; Sameh, A.
The authors present new parallel formulations of the Barnes-Hut method for n-body simulations on message passing computers. These parallel formulations partition the domain efficiently incurring minimal communication overhead. This is in contrast to existing schemes that are based on sorting a large number of keys or on the use of global data structures. The new formulations are augmented by alternate communication strategies which serve to minimize communication overhead. The impact of these communication strategies is experimentally studied. The authors report on experimental results obtained from an astrophysical simulation on an nCUBE2 parallel computer.
Simulation of Wearable Antennas for Body Centric Wireless Communication
NASA Astrophysics Data System (ADS)
Cousin, Richard; Rütschlin, Marc; Wittig, Tilmann; Bhattacharya, Arnab
2015-11-01
The performance of a body area network (BAN) is strongly dependent on several parameters which make wireless communication quite challenging. For instance, the performance of the antenna itself could be affected by its geometric deformation when the structure is directly integrated into clothes. Operation of the antenna close to the human body necessitates adjusting its design for the intended applications whereas the maximum SAR value estimated in such conditions has to respect the standards. In this context, simulation tools that can take into account specific biological models offer a range of possibilities for investigating and optimizing the performance of BAN devices. Two different applications are presented here: the case of an RFID tag operating at 870 MHz, and a UWB antenna working in a frequency range between 3 and 6 GHz. The simulation tools developed by CST are used in this context to optimize the implementation of BAN devices shown in this paper.
Acellular Biomaterials: An Evolving Alternative to Cell-Based Therapies
Burdick, Jason A.; Mauck, Robert L.; Gorman, Joseph H.; Gorman, Robert C.
2014-01-01
Acellular biomaterials can stimulate the local environment to repair tissues without the regulatory and scientific challenges of cell-based therapies. A greater understanding of the mechanisms of such endogenous tissue repair is furthering the design and application of these biomaterials. We discuss recent progress in acellular materials for tissue repair, using cartilage and cardiac tissues as examples of applications with substantial intrinsic hurdles, but where human translation is now occurring. PMID:23486777
Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct N-body
NASA Astrophysics Data System (ADS)
Rodriguez, Carl L.; Morscher, Meagan; Wang, Long; Chatterjee, Sourav; Rasio, Frederic A.; Spurzem, Rainer
2016-08-01
We present the first detailed comparison between million-body globular cluster simulations computed with a Hénon-type Monte Carlo code, CMC, and a direct N-body code, NBODY6++GPU. Both simulations start from an identical cluster model with 106 particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find good agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (>1000) are retained for 12 Gyr. Thus, the very accurate direct N-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-N dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct N-body, is capable of producing an accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes.
Spherical indentation of free-standing acellular extracellular matrix membranes.
Cloonan, Aidan J; O'Donnell, Michael R; Lee, William T; Walsh, Michael T; De Barra, Eamonn; McGloughlin, Tim M
2012-01-01
Numerous scaffold materials have been developed for tissue engineering and regenerative medicine applications to replace or repair damaged tissues and organs. Naturally occurring scaffold materials derived from acellular xenogeneic and autologous extracellular matrix (ECM) are currently in clinical use. These biological scaffold materials possess inherent variations in mechanical properties. Spherical indentation or ball burst testing has commonly been used to evaluate ECM and harvested tissue due to its ease of use and simulation of physiological biaxial loading, but has been limited by complex material deformation profiles. An analytical methodology has been developed and applied to experimental load-deflection data of a model hyperelastic material and lyophilized ECM scaffolds. An optimum rehydration protocol was developed based on water absorption, hydration relaxation and dynamic mechanical analysis. The analytical methodology was compared with finite element simulations of the tests and excellent correlation was seen between the computed biaxial stress resultants and geometry deformations. A minimum rehydration period of 5 min at 37°C was sufficient for the evaluated multilaminated ECM materials. The proposed approach may be implemented for convenient comparative analysis of ECM materials and source tissues, process optimization or during lot release testing. PMID:21864728
Simulation of Strongly Correlated Quantum Many-Body Systems
NASA Astrophysics Data System (ADS)
Bilgin, Ersen
In this thesis, we address the problem of solving for the properties of interacting quantum many-body systems in thermal equilibrium. The complexity of this problem increases exponentially with system size, limiting exact numerical simulations to very small systems. To tackle more complex systems, one needs to use heuristic algorithms that approximate solutions to these systems. Belief propagation is one such algorithm that we discuss in chapters 2 and 3. Using belief propagation, we demonstrate that it is possible to solve for static properties of highly correlated quantum many-body systems for certain geometries at all temperatures. In chapter 4, we generalize the multiscale renormalization ansatz to the anyonic setting to solve for the ground state properties of anyonic quantum many-body systems. The algorithms we present in chapters 2, 3, and 4 are very successful in certain settings, but they are not applicable to the most general quantum mechanical systems. For this, we propose using quantum computers as we discuss in chapter 5. The dimension reduction algorithm we consider in chapter 5 enables us to prepare thermal states of any quantum many-body system on a quantum computer faster than any previously known algorithm. Using these thermal states as the initialization of a quantum computer, one can study both static and dynamic properties of quantum systems without any memory overhead.
Simulating typical entanglement with many-body Hamiltonian dynamics
Nakata, Yoshifumi; Murao, Mio
2011-11-15
We study the time evolution of the amount of entanglement generated by one-dimensional spin-1/2 Ising-type Hamiltonians composed of many-body interactions. We investigate sets of states randomly selected during the time evolution generated by several types of time-independent Hamiltonians by analyzing the distributions of the amount of entanglement of the sets. We compare such entanglement distributions with that of typical entanglement, entanglement of a set of states randomly selected from a Hilbert space with respect to the unitarily invariant measure. We show that the entanglement distribution obtained by a time-independent Hamiltonian can simulate the average and standard deviation of the typical entanglement, if the Hamiltonian contains suitable many-body interactions. We also show that the time required to achieve such a distribution is polynomial in the system size for certain types of Hamiltonians.
Many-body localization: Entanglement and efficient numerical simulations
NASA Astrophysics Data System (ADS)
Pollmann, Frank
Many-body localization (MBL) occurs in isolated quantum systems when Anderson localization persists in the presence of finite interactions. To understand this phenomenon, the development of new efficient numerical methods to find highly excited many-body eigenstates is essential. In this talk, we will discuss two complimentary approaches to simulate MBL systems: First, we introduce a variant of the density-matrix renormalization group (DMRG) method that obtains individual highly excited eigenstates of MBL systems to machine precision accuracy at moderate to large disorder. This method explicitly takes advantage of the local spatial structure and the low entanglement which is characteristic for MBL eigenstates. Second, we propose an approach to directly find an approximate compact representation of the diagonalizing unitary by using a variational unitary matrix-product operator.
Simulation capability for dynamics of two-body flexible satellites
NASA Technical Reports Server (NTRS)
Austin, F.; Zetkov, G.
1973-01-01
An analysis and computer program were prepared to realistically simulate the dynamic behavior of a class of satellites consisting of two end bodies separated by a connecting structure. The shape and mass distribution of the flexible end bodies are arbitrary; the connecting structure is flexible but massless and is capable of deployment and retraction. Fluid flowing in a piping system and rigid moving masses, representing a cargo elevator or crew members, have been modeled. Connecting structure characteristics, control systems, and externally applied loads are modeled in easily replaced subroutines. Subroutines currently available include a telescopic beam-type connecting structure as well as attitude, deployment, spin and wobble control. In addition, a unique mass balance control system was developed to sense and balance mass shifts due to the motion of a cargo elevator. The mass of the cargo may vary through a large range. Numerical results are discussed for various types of runs.
Dynamic mixing in magma bodies - Theory, simulations, and implications
NASA Technical Reports Server (NTRS)
Oldenburg, Curtis M.; Spera, Frank J.; Yuen, David A.; Sewell, Granville
1989-01-01
The magma-mixing process is different from the mantle mixing process in that the mixing components of magma are dynamically active, with the melt density depending strongly on composition. This paper describes simulations of time-dependent variable-viscosity double-diffusive convection which were carried out to investigate quantitatively the mixing dynamics of magma in melt-dominated magma bodies. Results show that the dynamics of double-diffusive convection can impart complex patterns of composition, through time and space. The mixing time depends nonlinearly on many factors, including heat flux driving convection, the rate of diffusion of chemical species, the relative importance of thermal and chemical buoyancy, the viscosities of the mixing components, and the shape of the magma body.
CFD Approaches for Simulation of Wing-Body Stage Separation
NASA Technical Reports Server (NTRS)
Buning, Pieter G.; Gomez, Reynaldo J.; Scallion, William I.
2004-01-01
A collection of computational fluid dynamics tools and techniques are being developed and tested for application to stage separation and abort simulation for next-generation launch vehicles. In this work, an overset grid Navier-Stokes flow solver has been enhanced and demonstrated on a matrix of proximity cases and on a dynamic separation simulation of a belly-to-belly wing-body configuration. Steady cases show excellent agreement between Navier-Stokes results, Cartesian grid Euler solutions, and wind tunnel data at Mach 3. Good agreement has been obtained between Navier-Stokes, Euler, and wind tunnel results at Mach 6. An analysis of a dynamic separation at Mach 3 demonstrates that unsteady aerodynamic effects are not important for this scenario. Results provide an illustration of the relative applicability of Euler and Navier-Stokes methods to these types of problems.
Flow Simulation of N2B Hybrid Wing Body Configuration
NASA Technical Reports Server (NTRS)
Kim, Hyoungjin; Liou, Meng-Sing
2012-01-01
The N2B hybrid wing body aircraft was conceptually designed to meet environmental and performance goals for the N+2 generation transport set by the subsonic fixed wing project. In this study, flow fields around the N2B configuration is simulated using a Reynolds-averaged Navier-Stokes flow solver using unstructured meshes. Boundary conditions at engine fan face and nozzle exhaust planes are provided by response surfaces of the NPSS thermodynamic engine cycle model. The present flow simulations reveal challenging design issues arising from boundary layer ingestion offset inlet and nacelle-airframe interference. The N2B configuration can be a good test bed for application of multidisciplinary design optimization technology.
Parallelizing N-Body Simulations on a Heterogeneous Cluster
NASA Astrophysics Data System (ADS)
Stenborg, T. N.
2009-10-01
This thesis evaluates quantitatively the effectiveness of a new technique for parallelising direct gravitational N-body simulations on a heterogeneous computing cluster. In addition to being an investigation into how a specific computational physics task can be optimally load balanced across the heterogeneity factors of a distributed computing cluster, it is also, more generally, a case study in effective heterogeneous parallelisation of an all-pairs programming task. If high-performance computing clusters are not designed to be heterogeneous initially, they tend to become so over time as new nodes are added, or existing nodes are replaced or upgraded. As a result, effective techniques for application parallelisation on heterogeneous clusters are needed if maximum cluster utilisation is to be achieved and is an active area of research. A custom C/MPI parallel particle-particle N-body simulator was developed, validated and deployed for this evaluation. Simulation communication proceeds over cluster nodes arranged in a logical ring and employs nonblocking message passing to encourage overlap of communication with computation. Redundant calculations arising from force symmetry given by Newton's third law are removed by combining chordal data transfer of accumulated forces with ring passing data transfer. Heterogeneity in node computation speed is addressed by decomposing system data across nodes in proportion to node computation speed, in conjunction with use of evenly sized communication buffers. This scheme is shown experimentally to have some potential in improving simulation performance in comparison with an even decomposition of data across nodes. Techniques for further heterogeneous cluster load balancing are discussed and remain an opportunity for further work.
N-body simulations for coupled scalar-field cosmology
Li Baojiu; Barrow, John D.
2011-01-15
We describe in detail the general methodology and numerical implementation of consistent N-body simulations for coupled-scalar-field models, including background cosmology and the generation of initial conditions (with the different couplings to different matter species taken into account). We perform fully consistent simulations for a class of coupled-scalar-field models with an inverse power-law potential and negative coupling constant, for which the chameleon mechanism does not work. We find that in such cosmological models the scalar-field potential plays a negligible role except in the background expansion, and the fifth force that is produced is proportional to gravity in magnitude, justifying the use of a rescaled gravitational constant G in some earlier N-body simulation works for similar models. We then study the effects of the scalar coupling on the nonlinear matter power spectra and compare with linear perturbation calculations to see the agreement and places where the nonlinear treatment deviates from the linear approximation. We also propose an algorithm to identify gravitationally virialized matter halos, trying to take account of the fact that the virialization itself is also modified by the scalar-field coupling. We use the algorithm to measure the mass function and study the properties of dark-matter halos. We find that the net effect of the scalar coupling helps produce more heavy halos in our simulation boxes and suppresses the inner (but not the outer) density profile of halos compared with the {Lambda}CDM prediction, while the suppression weakens as the coupling between the scalar field and dark-matter particles increases in strength.
Acellular Dermal Matrix in Rotator Cuff Surgery.
Cooper, Joseph; Mirzayan, Raffy
2016-01-01
The success of rotator cuff repair (RCR) surgery can be measured clinically (validated outcome scores, range of motion) as well as structurally (re-tear rates using imaging studies). Regardless of repair type or technique, most studies have shown that patients do well clinically. However, multiple studies have also shown that structurally, the failure rate can be very high. A variety of factors, including poor tendon quality, age over 63 years, smoking, advanced fatty infiltration into the muscle, and the inability of the tendon to heal to bone, have been implicated as the cause of the high re-tear rate in RCRs. The suture-tendon interface is felt to be the weakest link in the RCR construct, and suture pullout through the tendon is believed to be the most common method of failure. This review of the published literature seeks to determine if there is support for augmentation of RCR with acellular dermal matrices to strengthen the suture-tendon interface and reduce the re-tear rate. PMID:27552454
Local stochastic non-Gaussianity and N-body simulations
Smith, Kendrick M.; LoVerde, Marilena E-mail: marilena@ias.edu
2011-11-01
Large-scale clustering of highly biased tracers of large-scale structure has emerged as one of the best observational probes of primordial non-Gaussianity of the local type (i.e. f{sub NL}{sup local}). This type of non-Gaussianity can be generated in multifield models of inflation such as the curvaton model. Recently, Tseliakhovich, Hirata, and Slosar showed that the clustering statistics depend qualitatively on the ratio of inflaton to curvaton power ξ after reheating, a free parameter of the model. If ξ is significantly different from zero, so that the inflaton makes a non-negligible contribution to the primordial adiabatic curvature, then the peak-background split ansatz predicts that the halo bias will be stochastic on large scales. In this paper, we test this prediction in N-body simulations. We find that large-scale stochasticity is generated, in qualitative agreement with the prediction, but that the level of stochasticity is overpredicted by ≈ 30%. Other predictions, such as ξ independence of the halo bias, are confirmed by the simulations. Surprisingly, even in the Gaussian case we do not find that halo model predictions for stochasticity agree consistently with simulations, suggesting that semi-analytic modeling of stochasticity is generally more difficult than modeling halo bias.
Numerical simulation of separated flows past bluff bodies
NASA Astrophysics Data System (ADS)
Alexandrou, Andreas Neophytou
1986-12-01
The steady two-dimensional flow past bluff bodies is simulated numerically using a hybrid Eulerian-Lagrangian model. The boundary layer effects, such as the location of the separation points and the rate of the generation of vorticity, are determined by a boundary layer solver. This solver uses Prandtl's boundary layer equations transformed by the Falkner-Skan transformation and then solved using a cubic spline approximation and a mean weighted residual technique. The vorticity generated at separation is discretized into elemental point vortices and convected downstream into the wake in a Lagrangian manner. The wake is modeled in a finite Eulerian computational domain using a modified Cloud-in-Cell (CIC) method. The velocity field at each time step is obtained as a solution to the rotationality condition using the finite element method in a cartesian mesh with nine-node elements and biquadratic shape functions. The biquadratic shape functions introduce a higher order interpolation scheme for the distribution of the vorticity at the nodal points than the bilinear (area) interpolation used in the original CIC method. The higher order interpolation as used in the CIC formulation performs better than the bilinear interpolation of the original method. This is demonstrated by the simulation of an isolated Rankine vortex. The ability of the CIC method to simulate the dynamics of vortex structures is also tested for the cases of flow past a flat plate and a circular cylinder.
Challenges to acellular biological scaffold mediated skeletal muscle tissue regeneration.
Corona, Benjamin T; Greising, Sarah M
2016-10-01
Volumetric muscle loss (VML) injuries present a complex and heterogeneous clinical problem that results in a chronic loss of muscle tissue and strength. The primary limitation to muscle tissue regeneration after VML injury is the frank loss of all native muscle constituents in the defect, especially satellite cells and the basal lamina. Recent advancements in regenerative medicine have set forth encouraging and emerging translational and therapeutic options for these devastating injuries including the surgical implantation of acellular biological scaffolds. While these biomaterials can modulate the wound environment, the existing data do not support their capacity to promote appreciable muscle fiber regeneration that can contribute to skeletal muscle tissue functional improvements. An apparent restriction of endogenous satellite cell (i.e., pax7(+)) migration to acellular biological scaffolds likely underlies this deficiency. This work critically evaluates the role of an acellular biological scaffold in orchestrating skeletal muscle tissue regeneration, specifically when used as a regenerative medicine approach for VML injury. PMID:27472161
The biasing scheme in N-body simulations
NASA Technical Reports Server (NTRS)
Park, Changbom
1991-01-01
A popular prescription for the galaxy formation, i.e., identifying peaks in the linear density field as sites of galaxy formation, is numerically justified under the Cold Dark Matter cosmogony. In particular, the peak-background scheme, originally developed to simplify analytic calculations of peak properties in a Gaussian random density field, is proven to be an excellent prescription for assigning peak tracers in low-resolution N-body simulations through many statistical comparisons. It is found that statistical properties of peak tracers allocated by the peak-background scheme closely resemble those of true peaks when the background smoothing length is about three times the galaxy scale and the correspondence improves as they evolve gravitationally. It is also shown that these peak particles are indeed strongly associated with nonlinear objects that are presently collapsed.
N-BODY SIMULATION OF THE STEPHAN'S QUINTET
Renaud, Florent; Appleton, Philip N.; Xu, C. Kevin
2010-11-20
The evolution of compact groups of galaxies may represent one of the few places in the nearby universe in which massive galaxies are being forged through a complex set of processes involving tidal interaction, ram-pressure stripping, and perhaps finally 'dry mergers' of galaxies stripped of their cool gas. Using collisionless N-body simulations, we propose a possible scenario for the formation of one of the best-studied compact groups: Stephan's Quintet. We define a serial approach which allows us to consider the history of the group as a sequence of galaxy-galaxy interactions seen as relatively separate events in time, but chained together in such a way as to provide a plausible scenario that ends in the current configuration of the galaxies. By covering a large set of parameters, we claim that it is very unlikely that both major tidal tails of the group have been created by the interaction between the main galaxy and a single intruder. We propose instead a scenario based on two satellites orbiting the main disk, plus the recent involvement of an additional interloper, coming from the background at high speed. This purely N-body study of the quintet will provide a parameter-space exploration of the basic dynamics of the group that can be used as a basis for a more sophisticated N-body/hydrodynamic study of the group that is necessary to explain the giant shock structure and other purely gaseous phenomena observed in both the cold, warm and hot gas in the group.
Second-order variational equations for N-body simulations
NASA Astrophysics Data System (ADS)
Rein, Hanno; Tamayo, Daniel
2016-07-01
First-order variational equations are widely used in N-body simulations to study how nearby trajectories diverge from one another. These allow for efficient and reliable determinations of chaos indicators such as the Maximal Lyapunov characteristic Exponent (MLE) and the Mean Exponential Growth factor of Nearby Orbits (MEGNO). In this paper we lay out the theoretical framework to extend the idea of variational equations to higher order. We explicitly derive the differential equations that govern the evolution of second-order variations in the N-body problem. Going to second order opens the door to new applications, including optimization algorithms that require the first and second derivatives of the solution, like the classical Newton's method. Typically, these methods have faster convergence rates than derivative-free methods. Derivatives are also required for Riemann manifold Langevin and Hamiltonian Monte Carlo methods which provide significantly shorter correlation times than standard methods. Such improved optimization methods can be applied to anything from radial-velocity/transit-timing-variation fitting to spacecraft trajectory optimization to asteroid deflection. We provide an implementation of first- and second-order variational equations for the publicly available REBOUND integrator package. Our implementation allows the simultaneous integration of any number of first- and second-order variational equations with the high-accuracy IAS15 integrator. We also provide routines to generate consistent and accurate initial conditions without the need for finite differencing.
Numerical techniques for large cosmological N-body simulations
NASA Technical Reports Server (NTRS)
Efstathiou, G.; Davis, M.; White, S. D. M.; Frenk, C. S.
1985-01-01
Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.
New material model for simulating large impacts on rocky bodies
NASA Astrophysics Data System (ADS)
Tonge, A.; Barnouin, O.; Ramesh, K.
2014-07-01
Large impact craters on an asteroid can provide insights into its internal structure. These craters can expose material from the interior of the body at the impact site [e.g., 1]; additionally, the impact sends stress waves throughout the body, which interrogate the asteroid's interior. Through a complex interplay of processes, such impacts can result in a variety of motions, the consequence of which may appear as lineaments that are exposed over all or portions of the asteroid's surface [e.g., 2,3]. While analytic, scaling, and heuristic arguments can provide some insight into general phenomena on asteroids, interpreting the results of a specific impact event, or series of events, on a specific asteroid geometry generally necessitates the use of computational approaches that can solve for the stress and displacement history resulting from an impact event. These computational approaches require a constitutive model for the material, which relates the deformation history of a small material volume to the average force on the boundary of that material volume. In this work, we present a new material model that is suitable for simulating the failure of rocky materials during impact events. This material model is similar to the model discussed in [4]. The new material model incorporates dynamic sub-scale crack interactions through a micro-mechanics-based damage model, thermodynamic effects through the use of a Mie-Gruneisen equation of state, and granular flow of the fully damaged material. The granular flow model includes dilatation resulting from the mutual interaction of small fragments of material (grains) as they are forced to slide and roll over each other and includes a P-α type porosity model to account for compaction of the granular material in a subsequent impact event. The micro-mechanics-based damage model provides a direct connection between the flaw (crack) distribution in the material and the rate-dependent strength. By connecting the rate
New material model for simulating large impacts on rocky bodies
NASA Astrophysics Data System (ADS)
Tonge, A.; Barnouin, O.; Ramesh, K.
2014-07-01
Large impact craters on an asteroid can provide insights into its internal structure. These craters can expose material from the interior of the body at the impact site [e.g., 1]; additionally, the impact sends stress waves throughout the body, which interrogate the asteroid's interior. Through a complex interplay of processes, such impacts can result in a variety of motions, the consequence of which may appear as lineaments that are exposed over all or portions of the asteroid's surface [e.g., 2,3]. While analytic, scaling, and heuristic arguments can provide some insight into general phenomena on asteroids, interpreting the results of a specific impact event, or series of events, on a specific asteroid geometry generally necessitates the use of computational approaches that can solve for the stress and displacement history resulting from an impact event. These computational approaches require a constitutive model for the material, which relates the deformation history of a small material volume to the average force on the boundary of that material volume. In this work, we present a new material model that is suitable for simulating the failure of rocky materials during impact events. This material model is similar to the model discussed in [4]. The new material model incorporates dynamic sub-scale crack interactions through a micro-mechanics-based damage model, thermodynamic effects through the use of a Mie-Gruneisen equation of state, and granular flow of the fully damaged material. The granular flow model includes dilatation resulting from the mutual interaction of small fragments of material (grains) as they are forced to slide and roll over each other and includes a P-α type porosity model to account for compaction of the granular material in a subsequent impact event. The micro-mechanics-based damage model provides a direct connection between the flaw (crack) distribution in the material and the rate-dependent strength. By connecting the rate
Direct numerical simulation of the sea flows around blunt bodies
NASA Astrophysics Data System (ADS)
Matyushin, Pavel V.; Gushchin, Valentin A.
2015-11-01
The aim of the present paper is the demonstration of the opportunities of the mathematical modeling of the separated flows of the sea water around blunt bodies on the basis of the Navier-Stokes equations (NSE) in the Boussinesq approximation. The 3D density stratified incompressible viscous fluid flows around a sphere have been investigated by means of the direct numerical simulation (DNS) on supercomputers and the visualization of the 3D vortex structures in the wake. For solving of NSE the Splitting on physical factors Method for Incompressible Fluid flows (SMIF) with hybrid explicit finite difference scheme (second-order accuracy in space, minimum scheme viscosity and dispersion, capable for work in wide range of the Reynolds (Re) and the internal Froude (Fr) numbers and monotonous) has been developed and successfully applied. The different transitions in sphere wakes with increasing of Re (10 < Re < 500) and decreasing of Fr (0.005 < Fr < 100) have been investigated in details. Thus the classifications of the viscous fluid flow regimes around a sphere have been refined.
Relaxation in N-body simulations of spherical systems
NASA Astrophysics Data System (ADS)
Sellwood, J. A.
2015-11-01
I present empirical measurements of the rate of relaxation in N-body simulations of stable spherical systems and distinguish two separate types of relaxation: energy diffusion that is largely independent of particle mass, and energy exchange between particles of differing masses. While diffusion is generally regarded as a Fokker-Planck process, it can equivalently be viewed as the consequence of collective oscillations that are driven by shot noise. Empirical diffusion rates scale as N-1 in inhomogeneous models, in agreement with Fokker-Planck predictions, but collective effects cause relaxation to scale more nearly as N-1/2 in the special case of a uniform sphere. I use four different methods to compute the gravitational field, and a 100-fold range in the numbers of particles in each case. I find the rate at which energy is exchanged between particles of differing masses does not depend at all on the force determination method, but I do find the energy diffusion rate is marginally lower when a field method is used. The relaxation rate in 3D is virtually independent of the method used because it is dominated by distant encounters; any method to estimate the gravitational field that correctly captures the contributions from distant particles must also capture their statistical fluctuations and the collective modes they drive.
FORMING CIRCUMBINARY PLANETS: N-BODY SIMULATIONS OF KEPLER-34
Lines, S.; Leinhardt, Z. M.; Paardekooper, S.; Baruteau, C.; Thebault, P.
2014-02-10
Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors to high energy, potentially destructive collisions that inhibit the growth of aspiring planets. In this work, we conduct high-resolution, inter-particle gravity enabled N-body simulations to investigate the feasibility of planetesimal growth in the Kepler-34 system. We improve upon previous work by including planetesimal disk self-gravity and an extensive collision model to accurately handle inter-planetesimal interactions. We find that super-catastrophic erosion events are the dominant mechanism up to and including the orbital radius of Kepler-34(AB)b, making in situ growth unlikely. It is more plausible that Kepler-34(AB)b migrated from a region beyond 1.5 AU. Based on the conclusions that we have made for Kepler-34, it seems likely that all of the currently known circumbinary planets have also migrated significantly from their formation location with the possible exception of Kepler-47(AB)c.
SKID: Finding Gravitationally Bound Groups in N-body Simulations
NASA Astrophysics Data System (ADS)
N-Body Shop
2011-02-01
SKID finds gravitationally bound groups in N-body simulations. The SKID program will group different types of particles depending on the type of input binary file. This could be either dark matter particles, gas particles, star particles or gas and star particles depending on what is in the input tipsy binary file. Once groups with at least a certain minimum number of members have been determined, SKID will remove particles which are not bound to the group. SKID must use the original positions of all the particles to determine whether or not particles are bound. This procedure which we call unbinding, is again dependent on the type of grouping we are dealing with. There are two cases, one for dark matter only or star particles only (case 1 unbinding), the other for inputs including gas (also stars in a dark matter environment this is case 2 unbinding). Skid version 1.3 is a much improved version of the old denmax-1.1 version. The new name was given to avoid confusion with the DENMAX program of Gelb & Bertschinger, and although it is based on the same idea it represents a substantial evolution in the method.
Corrosion and tribocorrosion of hafnium in simulated body fluids.
Rituerto Sin, J; Neville, A; Emami, N
2014-08-01
Hafnium is a passive metal with good biocompatibility and osteogenesis, however, little is known about its resistance to wear and corrosion in biological environments. The corrosion and tribocorrosion behavior of hafnium and commercially pure (CP) titanium in simulated body fluids were investigated using electrochemical techniques. Cyclic polarization scans and open circuit potential measurements were performed in 0.9% NaCl solution and 25% bovine calf serum solution to assess the effect of organic species on the corrosion behavior of the metal. A pin-on-plate configuration tribometer and a three electrode electrochemical cell were integrated to investigate the tribocorrosion performance of the studied materials. The results showed that hafnium has good corrosion resistance. The corrosion density currents measured in its passive state were lower than those measured in the case of CP titanium; however, it showed a higher tendency to suffer from localized corrosion, which was more acute when imperfections were present on the surface. The electrochemical breakdown of the oxide layer was retarded in the presence of proteins. Tribocorrosion tests showed that hafnium has the ability to quickly repassivate after the oxide layer was damaged; however, it showed higher volumetric loss than CP titanium in equivalent wear-corrosion conditions. © 2013 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 102B: 1157-1164, 2014. PMID:24376175
N-body simulations of viscous instability of planetary rings
NASA Astrophysics Data System (ADS)
Salo, Heikki; Schmidt, Jürgen
2010-04-01
We study viscous instability of planetary rings in terms of N-body simulations. We show that for rings composed of fairly elastic particles (e.g. as in Hatzes et al. [Hatzes, A., Bridges, F.G., Lin, D.N.C., 1988. Collisional properties of ice spheres at low impact velocities. Mon. Not. R. Astron. Soc. 231, 1091-1115]) the instability may lead to the spontaneous formation of dense ringlets in a background of lower density. In most parts of Saturn's rings the particle collisions are probably much more dissipative, as suggested by the presence of self-gravity wakes, and classic viscous instability should be suppressed. However, our results demonstrate that the mechanism of viscous instability itself is valid. The dynamical effects of size-dependent elasticity in a system with a size distribution have never been studied before. We show that this may in principle lead to a size-selective viscous instability, small particles concentrating on ringlets against the more uniform background of large particles.
Classical simulation of quantum many-body systems
NASA Astrophysics Data System (ADS)
Huang, Yichen
Classical simulation of quantum many-body systems is in general a challenging problem for the simple reason that the dimension of the Hilbert space grows exponentially with the system size. In particular, merely encoding a generic quantum many-body state requires an exponential number of bits. However, condensed matter physicists are mostly interested in local Hamiltonians and especially their ground states, which are highly non-generic. Thus, we might hope that at least some physical systems allow efficient classical simulation. Starting with one-dimensional (1D) quantum systems (i.e., the simplest nontrivial case), the first basic question is: Which classes of states have efficient classical representations? It turns out that this question is quantitatively related to the amount of entanglement in the state, for states with "little entanglement'' are well approximated by matrix product states (a data structure that can be manipulated efficiently on a classical computer). At a technical level, the mathematical notion for "little entanglement'' is area law, which has been proved for unique ground states in 1D gapped systems. We establish an area law for constant-fold degenerate ground states in 1D gapped systems and thus explain the effectiveness of matrix-product-state methods in (e.g.) symmetry breaking phases. This result might not be intuitively trivial as degenerate ground states in gapped systems can be long-range correlated. Suppose an efficient classical representation exists. How can one find it efficiently? The density matrix renormalization group is the leading numerical method for computing ground states in 1D quantum systems. However, it is a heuristic algorithm and the possibility that it may fail in some cases cannot be completely ruled out. Recently, a provably efficient variant of the density matrix renormalization group has been developed for frustration-free 1D gapped systems. We generalize this algorithm to all (i.e., possibly frustrated) 1D
Takano, Yoshiro; Sakai, Hideo; Watanabe, Eiko; Ideguchi-Ohma, Noriko; Jayawardena, Chantha K; Arai, Kazumi; Asawa, Yukiyo; Nakano, Yukiko; Shuda, Yoko; Sakamoto, Yujiro; Terashima, Tatsuo
2003-01-01
The mechanism whereby a region-specific deposition of the two types of cementum (cellular cementum and acellular extrinsic fibre cementum) is regulated on the growing root surface was tested using bisphosphonate-affected teeth of young rats and guinea pigs. The animals were injected subcutaneously with 8 or 10 mg P x kg body weight(-1) x day(-1) of 1-hydroxyethylidene-1,1-bisphosphonate (HEBP) for 1 or 2 weeks. In rat molars, HEBP prevented mineralization of newly formed root dentin matrix and totally inhibited de novo deposition of acellular extrinsic fibre cementum. Instead, thick cellular cementum was induced on the non-mineralized root dentin surface, irrespective of the position of the root. In both animals, cellular cementum was also induced on the non-mineralized surface of root analogue dentin in HEBP-affected incisors, where only acellular extrinsic fibre cementum is deposited under normal conditions. In normal rat molars, dentin sialoprotein (DSP) was concentrated along the dentin-cellular cementum border, but not that of dentin and acellular extrinsic fibre cementum. In HEBP-affected rat incisors, DSP was shown to penetrate through the non-mineralized dentin into the surrounding tissues, but not through the mineralized portions. These data suggest that, at the site of cellular cementum formation, putative inducing factors for cellular cementum might diffuse into the periodontal space through the newly deposited mantle dentin matrix before it is mineralized. At earlier stages of root formation, mantle dentin might mineralize more promptly not to allow such diffusion. The timing of mineralization of mantle dentin matrix might be the key determinant of the types of the cementum deposited on the growing root surface. PMID:14756246
Smooth Potential Chaos and N-Body Simulations
NASA Astrophysics Data System (ADS)
Kandrup, Henry E.; Sideris, Ioannis V.
2003-03-01
Integrations in fixed N-body realizations of smooth density distributions corresponding to a chaotic galactic potential can be used to derive reliable estimates of the largest (finite-time) Lyapunov exponent χS associated with an orbit in the smooth potential generated from the same initial condition, even though the N-body orbit is typically characterized by an N-body exponent χN>>χS. This can be accomplished by either comparing initially nearby orbits in a single N-body system or tracking orbits with the same initial condition evolved in two different N-body realizations of the same smooth density.
Huang, Haitao; Xiao, Hongxi; Liu, Huawei; Niu, Yu; Yan, Rongzeng; Hu, Min
2015-10-01
Acellular nerves are composed of a basal lamina tube, which retains sufficient bioactivity to promote axon regeneration, thereby repairing peripheral nerve gaps. However, the clinical application of acellular allografts has been restricted due to its limited availability. To investigate whether xenografts, a substitute to allograft acellular nerves in abundant supply, could efficiently promote nerve regeneration, rabbit and rat acellular nerve grafts were used to reconstruct 1 cm defects in Wistar rat facial nerves. Autologous peroneal nerve grafts served as a positive control group. A total of 12 weeks following the surgical procedure, the axon number, myelinated axon number, myelin sheath thickness, and nerve conduction velocity of the rabbit and rat‑derived acellular nerve grafts were similar, whereas the fiber diameter of the rabbit‑derived acellular xenografts decreased, as compared with those of rat‑derived acellular allografts. Autografts exerted superior effects on nerve regeneration; however, no significant difference was observed between the axon number in the autograft group, as compared with the two acellular groups. These results suggested that autografts perform better than acellular nerve grafts, and chemically extracted acellular allografts and xenografts have similar effects on the regeneration of short facial nerve defects. PMID:26239906
Numerical Simulation of Floating Bodies in Extreme Free Surface Waves
NASA Astrophysics Data System (ADS)
Hu, Zheng Zheng; Causon, Derek; Mingham, Clive; Qiang, Ling
2010-05-01
and efficient. Firstly, extreme design wave conditions are generated in an empty NWT and compared with physical experiments as a precursor to calculations to investigate the survivability of the Bobber device operating in a challenging wave climate. Secondly, we consider a bench-mark test case involving in a first order regular wave maker acting on a fixed cylinder and Pelamis. Finally, a floating Bobber has been simulated under extreme wave conditions. These results will be reported at the meeting. Causon D.M., Ingram D.M., Mingham C.G., Yang G. Pearson R.V. (2000). Calculation of shallow water flows using a Cartesian cut cell approach. Advances in Water resources, 23: 545-562. Causon D.M., Ingram D.M., Mingham C.G. (2000). A Cartesian cut cell method for shallow water flows with moving boundaries. Advances in Water resources, 24: 899-911. Dalzell J.F. 1999 A note on finite depth second-order wave-wave interactions. Appl. Ocean Res. 21, 105-111. Ning D.Z., Zang J., Liu S.X. Eatock Taylor R. Teng B. & Taylor P.H. 2009 Free surface and wave kinematics for nonlinear focused wave groups. J. Ocean Engineering. Accepted. Hu Z.Z., Causon D.M., Mingham C.M. and Qian L.(2009). Numerical wave tank study of a wave energy converter in heave. Proceedlings 19th ISOPE conference, Osaka, Japan Qian L., Causon D.M. & Mingham C.G., Ingram D.M. 2006 A free-surface capturing method for two fluid flows with moving bodies. Proc. Roy. Soc. London, Vol. A 462 21-42.
Bluff Body Flow Simulation Using a Vortex Element Method
Anthony Leonard; Phillippe Chatelain; Michael Rebel
2004-09-30
Heavy ground vehicles, especially those involved in long-haul freight transportation, consume a significant part of our nation's energy supply. it is therefore of utmost importance to improve their efficiency, both to reduce emissions and to decrease reliance on imported oil. At highway speeds, more than half of the power consumed by a typical semi truck goes into overcoming aerodynamic drag, a fraction which increases with speed and crosswind. Thanks to better tools and increased awareness, recent years have seen substantial aerodynamic improvements by the truck industry, such as tractor/trailer height matching, radiator area reduction, and swept fairings. However, there remains substantial room for improvement as understanding of turbulent fluid dynamics grows. The group's research effort focused on vortex particle methods, a novel approach for computational fluid dynamics (CFD). Where common CFD methods solve or model the Navier-Stokes equations on a grid which stretches from the truck surface outward, vortex particle methods solve the vorticity equation on a Lagrangian basis of smooth particles and do not require a grid. They worked to advance the state of the art in vortex particle methods, improving their ability to handle the complicated, high Reynolds number flow around heavy vehicles. Specific challenges that they have addressed include finding strategies to accurate capture vorticity generation and resultant forces at the truck wall, handling the aerodynamics of spinning bodies such as tires, application of the method to the GTS model, computation time reduction through improved integration methods, a closest point transform for particle method in complex geometrics, and work on large eddy simulation (LES) turbulence modeling.
Modal reduction strategies for interconnected flexible bodies simulation
NASA Technical Reports Server (NTRS)
Eke, F. O.; Man, G. K.
1989-01-01
Multi-body dynamics programs require characterization of each body. The Galileo spacecraft system modes to be retained were determined using available criteria, modal influence coefficients, and bode. The descent to component level was achieved via a two-phase diagonalization process starting with submatrices of truncated augmented system modal matrix.
Rigid Body Motion in Stereo 3D Simulation
ERIC Educational Resources Information Center
Zabunov, Svetoslav
2010-01-01
This paper addresses the difficulties experienced by first-grade students studying rigid body motion at Sofia University. Most quantities describing the rigid body are in relations that the students find hard to visualize and understand. They also lose the notion of cause-result relations between vector quantities, such as the relation between…
Body Constraints on Motor Simulation in Autism Spectrum Disorders
ERIC Educational Resources Information Center
Conson, Massimiliano; Hamilton, Antonia; De Bellis, Francesco; Errico, Domenico; Improta, Ilaria; Mazzarella, Elisabetta; Trojano, Luigi; Frolli, Alessandro
2016-01-01
Developmental data suggested that mental simulation skills become progressively dissociated from overt motor activity across development. Thus, efficient simulation is rather independent from current sensorimotor information. Here, we tested the impact of bodily (sensorimotor) information on simulation skills of adolescents with Autism Spectrum…
Infection in the Nasal Tip Caused by Acellular Dermal Matrix.
Lee, Kun Hee
2015-12-01
A 19-year-old female patient visited our clinic for rhinoplasty. She complained about her low take-off point, which was apparent in profile view, and wanted slight tip projection. She refused additional cartilage harvesting from ears or ribs but consented to the use of homologous tissue, including acellular dermal matrix, for her dorsum and tip. Septoturbinoplasty was performed, and only a very small amount of septal cartilage could be harvested. It was used as both the columellar strut and the alar rim graft. Nasal dorsum and tip were augmented with acellular dermal matrix. Three months postoperatively, she experienced a few episodes of edema and redness on her nasal tip, followed by pus exudation from the nasal skin. Six months postoperatively, she underwent revision rhinoplasty for removal of inflamed grafts, and onlay tip graft with homologous rib cartilage was performed. Nasal dorsum or tip grafts are an integral part of Asian rhinoplasty. Autogenous tissue is the gold standard for grafting materials. However, the limited availability of autogenous tissue and the preference of patients and surgeons for artificial surgical implants make Asian rhinoplasty challenging. Unavailability of autogenous cartilage and patient refusal of artificial implants led to the use of acellular dermal matrix (ADM) in the nasal dorsum and tip for this case. This is the first report of postoperative complication because of infection rather than absorption after ADM use. PMID:26894006
NASA Technical Reports Server (NTRS)
Kremic, Tibor; Vento, Dan; Lalli, Nick; Palinski, Timothy
2014-01-01
Science, technology, and planetary mission communities have a growing interest in components and systems that are capable of working in extreme (high) temperature and pressure conditions. Terrestrial applications range from scientific research, aerospace, defense, automotive systems, energy storage and power distribution, deep mining and others. As the target environments get increasingly extreme, capabilities to develop and test the sensors and systems designed to operate in such environments will be required. An application of particular importance to the planetary science community is the ability for a robotic lander to survive on the Venus surface where pressures are nearly 100 times that of Earth and temperatures approach 500C. The scientific importance and relevance of Venus missions are stated in the current Planetary Decadal Survey. Further, several missions to Venus were proposed in the most recent Discovery call. Despite this interest, the ability to accurately simulate Venus conditions at a scale that can test and validate instruments and spacecraft systems and accurately simulate the Venus atmosphere has been lacking. This paper discusses and compares the capabilities that are known to exist within and outside the United States to simulate the extreme environmental conditions found in terrestrial or planetary surfaces including the Venus atmosphere and surface. The paper then focuses on discussing the recent additional capability found in the NASA Glenn Extreme Environment Rig (GEER). The GEER, located at the NASA Glenn Research Center in Cleveland, Ohio, is designed to simulate not only the temperature and pressure extremes described, but can also accurately reproduce the atmospheric compositions of bodies in the solar system including those with acidic and hazardous elements. GEER capabilities and characteristics are described along with operational considerations relevant to potential users. The paper presents initial operating results and concludes
Miri, Amir K; Muja, Naser; Kamranpour, Neysan O; Lepry, William C; Boccaccini, Aldo R; Clarke, Susan A; Nazhat, Showan N
2016-04-01
Gel aspiration-ejection (GAE) has recently been introduced as an effective technique for the rapid production of injectable dense collagen (IDC) gel scaffolds with tunable collagen fibrillar densities (CFDs) and microstructures. Herein, a GAE system was applied for the advanced production and delivery of IDC and IDC-Bioglass(®) (IDC-BG) hybrid gel scaffolds for potential bone tissue engineering applications. The efficacy of GAE in generating mineralizable IDC-BG gels (from an initial 75-25 collagen-BG ratio) produced through needle gauge numbers 8G (3.4 mm diameter and 6 wt% CFD) and 14G (1.6 mm diameter and 14 wt% CFD) was investigated. Second harmonic generation (SHG) imaging of as-made gels revealed an increase in collagen fibril alignment with needle gauge number. In vitro mineralization of IDC-BG gels was confirmed where carbonated hydroxyapatite was detected as early as day 1 in simulated body fluid, which progressively increased up to day 14. In vivo mineralization of, and host response to, acellular IDC and IDC-BG gel scaffolds were further investigated following subcutaneous injection in adult rats. Mineralization, neovascularization and cell infiltration into the scaffolds was enhanced by the addition of BG and at day 21 post injection, there was evidence of remodelling of granulation tissue into woven bone-like tissue in IDC-BG. SHG imaging of explanted scaffolds indicated collagen fibril remodelling through cell infiltration and mineralization over time. In sum, the results suggest that IDC-BG hybrid gels have osteoinductive properties and potentially offer a novel therapeutic approach for procedures requiring the injectable delivery of a malleable and dynamic bone graft that mineralizes under physiological conditions. PMID:26871889
Message in the "Body": Effects of Simulation in Sentence Production
ERIC Educational Resources Information Center
Sato, Manami
2010-01-01
This study investigates the role of mental simulation in message formulation and grammatical encoding in two typologically distinct languages, English and Japanese. It examines relationships among physical motion, mental simulation, and sentence production, following the claims of Perceptual Symbol Systems (Barsalou, 1999) that people understand…
Efficiency and capabilities of multi-body simulations
NASA Technical Reports Server (NTRS)
Vandervoort, R. J.
1989-01-01
Simulation efficiency and capability go hand in hand. The more capability you have the lower the efficiency will be. The efficiency and capabilities are discussed. The lesson learned about generic simulation is: Don't rule out any capabilities at the beginning, but keep each one on a switch so it can be bypassed when warranted by a specific application.
NASA Astrophysics Data System (ADS)
Jin, Yuzhen; Li, Jun; Zhu, Linhang; Du, Jiayou; Jin, Yingzi; Lin, Peifeng
2014-12-01
For the large deformation of the flexible body may cause the fluid grid distortion, which will make the numerical calculation tedious, even to end, the numerical simulation of the flexible body coupling with the fluid is always a tough problem. In this paper, the flexible body is under two kinds of constrained conditions and the ratio of length-diameter is 1:30. The Reynolds number of the airflow is 513, belonging to the area of low Reynolds number. The control equations of the coupling of flexible body with airflow are built and the adaptive grid control method is adopted to conduct the three-dimensional numerical simulation of the movement of the flexible body. The numerical results show that it is possible to simulate the characteristics of the flexible body's movement in the low Reynolds number airflow when the appropriate control equations are modeled and suitable equation-solving method is adopted. Unconstrained flexible body would turn over forward along the airflow's diffusion direction, while constrained flexible body in the flow field will make periodic rotation motion along the axis of the flexible body, and the bending deformation is more obvious than that of unconstrained flexible body. The preliminary three-dimensional numerical simulation can provide references for further research on the characteristics of the yarn movement in high Reynolds number airflow.
Visualization of N-body Simulations in Virtual Worlds
NASA Astrophysics Data System (ADS)
Knop, Robert A.; Ames, J.; Djorgovski, G.; Farr, W.; Hut, P.; Johnson, A.; McMillan, S.; Nakasone, A.; Vesperini, E.
2010-01-01
We report on work to use virtual worlds for visualizing the results of N-body calculations, on three levels. First, we have written a demonstration 3-body solver entirely in the scripting language of the popularly used virtual world Second Life. Second, we have written a physics module for the open source virtual world OpenSim that performs N-body calculations as the physics engine for the server, allowing natural 3-d visualization of the solution as the solution is being performed. Finally, we give an initial report on the potential use of virtual worlds to visualize calculations which have previously been performed, or which are being performed in other processes and reported to the virtual world server. This work has been performed as part of the Meta-Institute of Computational Astrophysics (MICA). http://www.mica-vw.org
The direct numerical simulations of the turbulent wakes of axisymmetric bodies
NASA Technical Reports Server (NTRS)
Riley, J. J.; Metcalfe, R. W.
1978-01-01
Results of direct numerical simulations of turbulence are compared with both laboratory data and self-similarity theory for the case of the turbulent wakes of towed, axisymmetric bodies. In general, the agreement of the simulation results with both the laboratory data and the self-similarity theory is good, although the comparisons are hampered by inadequate procedures for initializing the numerical simulations.
Computer simulation of plasma and N-body problems
NASA Technical Reports Server (NTRS)
Harries, W. L.; Miller, J. B.
1975-01-01
The following FORTRAN language computer codes are presented: (1) efficient two- and three-dimensional central force potential solvers; (2) a three-dimensional simulator of an isolated galaxy which incorporates the potential solver; (3) a two-dimensional particle-in-cell simulator of the Jeans instability in an infinite self-gravitating compressible gas; and (4) a two-dimensional particle-in-cell simulator of a rotating self-gravitating compressible gaseous system of which rectangular coordinate and superior polar coordinate versions were written.
G-Guidance Interface Design for Small Body Mission Simulation
NASA Technical Reports Server (NTRS)
Acikmese, Behcet; Carson, John; Phan, Linh
2008-01-01
The G-Guidance software implements a guidance and control (G and C) algorithm for small-body, autonomous proximity operations, developed under the Small Body GN and C task at JPL. The software is written in Matlab and interfaces with G-OPT, a JPL-developed optimization package written in C that provides G-Guidance with guaranteed convergence to a solution in a finite computation time with a prescribed accuracy. The resulting program is computationally efficient and is a prototype of an onboard, real-time algorithm for autonomous guidance and control. Two thruster firing schemes are available in G-Guidance, allowing tailoring of the software for specific mission maneuvers. For example, descent, landing, or rendezvous benefit from a thruster firing at the maneuver termination to mitigate velocity errors. Conversely, ascent or separation maneuvers benefit from an immediate firing to avoid potential drift toward a second body. The guidance portion of this software explicitly enforces user-defined control constraints and thruster silence times while minimizing total fuel usage. This program is currently specialized to small-body proximity operations, but the underlying method can be generalized to other applications.
Porcine acellular lung matrix for wound healing and abdominal wall reconstruction: A pilot study
Fernandez-Moure, Joseph S; Van Eps, Jeffrey L; Rhudy, Jessica R; Cabrera, Fernando J; Acharya, Ghanashyam S; Tasciotti, Ennio; Sakamoto, Jason; Nichols, Joan E
2016-01-01
Surgical wound healing applications require bioprosthetics that promote cellular infiltration and vessel formation, metrics associated with increased mechanical strength and resistance to infection. Porcine acellular lung matrix is a novel tissue scaffold known to promote cell adherence while minimizing inflammatory reactions. In this study, we evaluate the capacity of porcine acellular lung matrix to sustain cellularization and neovascularization in a rat model of subcutaneous implantation and chronic hernia repair. We hypothesize that, compared to human acellular dermal matrix, porcine acellular lung matrix would promote greater cell infiltration and vessel formation. Following pneumonectomy, porcine lungs were processed and characterized histologically and by scanning electron microscopy to demonstrate efficacy of the decellularization. Using a rat model of subcutaneou implantation, porcine acellular lung matrices (n = 8) and human acellular dermal matrices (n = 8) were incubated in vivo for 6 weeks. To evaluate performance under mechanically stressed conditions, porcine acellular lung matrices (n = 7) and human acellular dermal matrices (n = 7) were implanted in a rat model of chronic ventral incisional hernia repair for 6 weeks. After 6 weeks, tissues were evaluated using hematoxylin and eosin and Masson’s trichrome staining to quantify cell infiltration and vessel formation. Porcine acellular lung matrices were shown to be successfully decellularized. Following subcutaneous implantation, macroscopic vessel formation was evident. Porcine acellular lung matrices demonstrated sufficient incorporation and showed no evidence of mechanical failure after ventral hernia repair. Porcine acellular lung matrices demonstrated significantly greater cellular density and vessel formation when compared to human acellular dermal matrix. Vessel sizes were similar across all groups. Cell infiltration and vessel formation are well-characterized metrics of incorporation
Porcine acellular lung matrix for wound healing and abdominal wall reconstruction: A pilot study.
Fernandez-Moure, Joseph S; Van Eps, Jeffrey L; Rhudy, Jessica R; Cabrera, Fernando J; Acharya, Ghanashyam S; Tasciotti, Ennio; Sakamoto, Jason; Nichols, Joan E
2016-01-01
Surgical wound healing applications require bioprosthetics that promote cellular infiltration and vessel formation, metrics associated with increased mechanical strength and resistance to infection. Porcine acellular lung matrix is a novel tissue scaffold known to promote cell adherence while minimizing inflammatory reactions. In this study, we evaluate the capacity of porcine acellular lung matrix to sustain cellularization and neovascularization in a rat model of subcutaneous implantation and chronic hernia repair. We hypothesize that, compared to human acellular dermal matrix, porcine acellular lung matrix would promote greater cell infiltration and vessel formation. Following pneumonectomy, porcine lungs were processed and characterized histologically and by scanning electron microscopy to demonstrate efficacy of the decellularization. Using a rat model of subcutaneou implantation, porcine acellular lung matrices (n = 8) and human acellular dermal matrices (n = 8) were incubated in vivo for 6 weeks. To evaluate performance under mechanically stressed conditions, porcine acellular lung matrices (n = 7) and human acellular dermal matrices (n = 7) were implanted in a rat model of chronic ventral incisional hernia repair for 6 weeks. After 6 weeks, tissues were evaluated using hematoxylin and eosin and Masson's trichrome staining to quantify cell infiltration and vessel formation. Porcine acellular lung matrices were shown to be successfully decellularized. Following subcutaneous implantation, macroscopic vessel formation was evident. Porcine acellular lung matrices demonstrated sufficient incorporation and showed no evidence of mechanical failure after ventral hernia repair. Porcine acellular lung matrices demonstrated significantly greater cellular density and vessel formation when compared to human acellular dermal matrix. Vessel sizes were similar across all groups. Cell infiltration and vessel formation are well-characterized metrics of incorporation
Forward dynamics simulation of human body under tilting perturbations
NASA Astrophysics Data System (ADS)
Naderi, D.; Pasha Zanoosi, A. A.; Sadeghi-Mehr, M.
2012-02-01
Human body uses different strategies to maintain its stability and these strategies vary from fixed-foot strategies to strategies which foot is moved in order to increase the support base. Tilting movement of foot is one type of the perturbations usually is exposed to human body. In the presence of such perturbations human body must employ appropriate reactions to prevent threats like falling. But it is not clear that how human body maintains its stability by central nervous system (CNS). At present study it is tried that by presenting a musculoskeletal model of human lower extremity with four links, three degrees of freedom (DOF) and eight skeletal muscles, the level of muscle activations causes the maintenance of stability, be investigated. Using forward dynamics solution, leads to a more general problem, rather than inverse dynamics. Hence, forward dynamics solution by forward optimization has been used for solving this highly nonlinear problem. To this end, first the system's equations of motion has been derived using lagrangian dynamics. Eight Hill-type muscles as actuators of the system were modeled. Because determination of muscle forces considering their number is an undetermined problem, optimization of an appropriate goal function should be practiced. For optimization problem, the characteristics of genetic algorithms as a method based on direct search, and the direct collocation method, has been profited. Also by considering requirements of problem, some constraints such as conservation of model stability are entered into optimization procedure. Finally to investigate validation of model, the results from optimization and experimental data are compared and good agreements are obtained.
NASA Technical Reports Server (NTRS)
Croston, R. C.; Fitzjerrell, D. G.
1974-01-01
A mathematical model and digital computer simulation of the human cardiovascular system and its controls have been developed to simulate pulsatile dynamic responses to the cardiovascular experiments of the Skylab missions and to selected physiological stresses of manned space flight. Specific model simulations of the bicycle ergometry, lower body negative pressure, and tilt experiments have been developed and verified for 1-g response by comparison with available experimental data. The zero-g simulations of two Skylab experiments are discussed.
The Million-Body Problem: Particle Simulations in Astrophysics
Rasio, Fred [Northwestern University
2009-09-01
Computer simulations using particles play a key role in astrophysics. They are widely used to study problems across the entire range of astrophysical scales, from the dynamics of stars, gaseous nebulae, and galaxies, to the formation of the largest-scale structures in the universe. The 'particles' can be anything from elementary particles to macroscopic fluid elements, entire stars, or even entire galaxies. Using particle simulations as a common thread, this talk will present an overview of computational astrophysics research currently done in our theory group at Northwestern. Topics will include stellar collisions and the gravothermal catastrophe in dense star clusters.
Kokkalis, Zinon T; Zanaros, George; Sotereanos, Dean G
2009-03-01
Ligament reconstruction tendon interposition arthroplasty is currently the preferred technique for carpometacarpal joint arthritis of the thumb by most surgeons. Despite its efficacy, morbidity has been associated with the harvest of the flexor carpi radialis tendon. Using an allograft as material for arthroplasty, donor site morbidity is avoided. In this report, we present our surgical technique to perform ligament reconstruction tendon interposition arthroplasty using an acellular dermal matrix allograft (GraftJacket) in patients with Eaton stages II, III, and IV symptomatic first carpometacarpal arthritis.One hundred thumbs with trapeziometacarpal osteoarthritis underwent surgical treatment using GraftJacket allograft instead of the flexor carpi radialis tendon autograft. Each patient was followed for a minimum of 12 months. The surgical procedure included trapezial excision and identification of the flexor carpi radialis. The allograft was cut to create a 15-cm strip. The ligament reconstruction was performed by passing the strip around the flexor carpi radialis tendon and suturing it to the base of the thumb metacarpal base through an intramedullary drill hole. The remaining portion of the allograft was fashioned as an interposition mass (anchovy) and interposed between the scaphoid and the base of the first metacarpal.All but 1 patient experienced significant improvement in his or her pain scale rating and grip and pinch strengths. Outcomes from this study compare very favorably with those of other series. No patients experienced a foreign body reaction or infection in this series. We believe that the use of an acellular dermal allograft for both ligament reconstruction and tendon interposition provides a safe and an effective alternative technique for the treatment of advanced first carpometacarpal arthritis. PMID:19276927
Zhang, Yanru; Zhang, Hui; Katiella, Kaka; Huang, Wenhua
2014-01-01
A chemically extracted acellular allogeneic nerve graft can reduce postoperative immune rejection, similar to an autologous nerve graft, and can guide neural regeneration. However, it remains poorly understood whether a chemically extracted acellular allogeneic nerve graft combined with neurotrophic factors provides a good local environment for neural regeneration. This study investigated the repair of injured rat sciatic nerve using a chemically extracted acellular allogeneic nerve graft combined with ciliary neurotrophic factor. An autologous nerve anastomosis group and a chemical acellular allogeneic nerve bridging group were prepared as controls. At 8 weeks after repair, sciatic functional index, evoked potential amplitude of the soleus muscle, triceps wet weight recovery rate, total number of myelinated nerve fibers and myelin sheath thickness were measured. For these indices, values in the three groups showed the autologous nerve anastomosis group > chemically extracted acellular nerve graft + ciliary neurotrophic factor group > chemical acellular allogeneic nerve bridging group. These results suggest that chemically extracted acellular nerve grafts combined with ciliary neurotrophic factor can repair sciatic nerve defects, and that this repair is inferior to autologous nerve anastomosis, but superior to chemically extracted acellular allogeneic nerve bridging alone. PMID:25221592
Zhang, Yanru; Zhang, Hui; Katiella, Kaka; Huang, Wenhua
2014-07-15
A chemically extracted acellular allogeneic nerve graft can reduce postoperative immune rejection, similar to an autologous nerve graft, and can guide neural regeneration. However, it remains poorly understood whether a chemically extracted acellular allogeneic nerve graft combined with neurotrophic factors provides a good local environment for neural regeneration. This study investigated the repair of injured rat sciatic nerve using a chemically extracted acellular allogeneic nerve graft combined with ciliary neurotrophic factor. An autologous nerve anastomosis group and a chemical acellular allogeneic nerve bridging group were prepared as controls. At 8 weeks after repair, sciatic functional index, evoked potential amplitude of the soleus muscle, triceps wet weight recovery rate, total number of myelinated nerve fibers and myelin sheath thickness were measured. For these indices, values in the three groups showed the autologous nerve anastomosis group > chemically extracted acellular nerve graft + ciliary neurotrophic factor group > chemical acellular allogeneic nerve bridging group. These results suggest that chemically extracted acellular nerve grafts combined with ciliary neurotrophic factor can repair sciatic nerve defects, and that this repair is inferior to autologous nerve anastomosis, but superior to chemically extracted acellular allogeneic nerve bridging alone. PMID:25221592
Bluff-body flow simulations using hybrid RANS/LES.
Roy, Christopher John; DeChant, Lawrence Justin.; Payne, Jeffrey L.; Blottner, Frederick G.
2003-06-01
The Detached Eddy Simulation (DES) and steady-state Reynolds-Averaged Navier-Stokes (RANS) turbulence modeling approaches are examined for the incompressible flow over a square cross-section cylinder at a Reynolds number of 21,400. A compressible flow code is used which employes a second-order Roe upwind spatial discretization. Efforts are made to assess the numerical accuracy of the DES predictions with regards to statistical convergence, iterative convergence, and temporal and spatial discretization error. Three-dimensional DES simulations compared well with two-dimensional DES simulations, suggesting that the dominant vortex shedding mechanism is effectively two-dimensional. The two-dimensional simulations are validated via comparison to experimental data for mean and RMS velocities as well as Reynolds stress in the cylinder wake. The steady-state RANS models significantly overpredict the size of the recirculation zone, thus underpredicting the drag coefficient relative to the experimental value. The DES model is found to give good agreement with the experimental velocity data in the wake, drag coefficient, and recirculation zone length.
Analytical stability and simulation response study for a coupled two-body system
NASA Technical Reports Server (NTRS)
Tao, K. M.; Roberts, J. R.
1975-01-01
An analytical stability study and a digital simulation response study of two connected rigid bodies are documented. Relative rotation of the bodies at the connection is allowed, thereby providing a model suitable for studying system stability and response during a soft-dock regime. Provisions are made of a docking port axes alignment torque and a despin torque capability for encountering spinning payloads. Although the stability analysis is based on linearized equations, the digital simulation is based on nonlinear models.
Phase II trial of whole-cell pertussis vaccine vs an acellular vaccine containing agglutinogens.
Miller, E; Ashworth, L A; Robinson, A; Waight, P A; Irons, L I
1991-01-12
An acellular pertussis vaccine containing agglutinogens 2 and 3, pertussis toxin, and filamentous haemagglutinin was developed by the Centre for Applied Microbiology and Research in the UK. 188 infants were entered into a randomised blind trial and received either the acellular or a whole-cell vaccine, combined with diphtheria and tetanus toxoids, in a 3, 5, and 8-10 month schedule. Local reactions were similar in the two groups but significantly fewer infants had systemic symptoms after the acellular vaccine. Mean log-antibody titres to the agglutinogen and toxin components were higher with the acellular than with the whole-cell vaccine. Persistence of antibodies one year after the third dose was also better in the acellular group. PMID:1670725
Use of an Acellular Regenerative Tissue Matrix Over Chronic Wounds
Stacey, D. Heath
2013-01-01
Objectives: Bioengineered skin grafts, including acellular dermal matrices, may be effective in treating lower extremity and trunk wounds that are not responsive to traditional wound management. Acellular dermal wound matrix is derived from human acellular dermal wound matrix (HADWM) tissue and provides a scaffold that supports cellular repopulation and revascularization. The major structural components of the dermis are retained during processing, and a single application has been shown to help achieve wound closure. Methods: This patient case series examined the use of HADWM on lower extremity and trunk wounds in 11 patients (6 male and 5 female) with a mean age of 55 years (range: 31–83 years). Wounds were debrided 1 to 2 times, followed by placement of HADWM (range: 4–330 cm2) on wounds that varied from the dorsal surface of the foot, lower abdomen, and lower extremity to the Achilles flap. A nonadherent layer in conjunction with bacitracin was placed over HADWM. Negative pressure wound therapy (NPWT) was placed over the HADWM and initiated continuously at −125 mm Hg for 1 to 2 weeks. After the application of NPWT, HADWM was covered with various gauze dressings using mineral oil. Results: All patients completed their treatment successfully, and follow-up ranged from 1 week to 6 months. One patient experienced an infection, which resulted in partial graft loss that required replacement with HADWM and NPWT. No additional complications occurred in the other patients. Conclusions: This patient case series demonstrated successful use of HADWM and NPWT, which further supports published studies documenting HADWM success in chronic wounds. PMID:24324850
Whooping cough, twenty years from acellular vaccines introduction.
Greco, D; Esposito, S; Tozzi, A; Pandolfi, E; Icardi, G; Giammanco, A
2015-01-01
Clinical pertussis resulting from infection with B. pertussis is a significant medical and public health problem, despite the huge success of vaccination that has greatly reduced its incidence. The whole cell vaccine had an undeniable success over the last 50 years, but its acceptance was strongly inhibited by fear, only partially justified, of severe side effects, but also, in the Western world, by the difficulty to enter in combination with other vaccines: today multi-vaccine formulations are essential to maintain a high vaccination coverage. The advent of acellular vaccines was greeted with enthusiasm by the public health world: in the Nineties, several controlled vaccine trials were carried out: they demonstrated a high safety and good efficacy of new vaccines. In fact, in the Western world, the acellular vaccines completely replaced the whole cells ones. In the last years, ample evidence on the variety of protection of these vaccines linked to the presence of different antigens of Bordetella pertussis was collected. It also became clear that the protection provided, on average around 80%, leaves every year a significant cohort of vaccinated susceptible even in countries with a vaccination coverage of 95%, such as Italy. Finally, it was shown that, as for the pertussis disease, protection decreases over time, to leave a proportion of adolescents and adults unprotected. Waiting for improved pertussis vaccines, the disease control today requires a different strategy that includes a booster at 5 years for infants, but also boosters for teenagers and young adults, re-vaccination of health care personnel, and possibly of pregnant women and of those who are in contact with infants (cocooning). Finally, the quest for better vaccines inevitably tends towards pertussis acellular vaccines with at least three components, which have demonstrated superior effectiveness and have been largely in use in Italy for fifteen years. PMID:26051141
Finite element simulation of food transport through the esophageal body
Yang, Wei; Fung, Tat Ching; Chian, Kerm Sim; Chong, Chuh Khiun
2007-01-01
The peristaltic transport of swallowed material in the esophagus is a neuro-muscular function involving the nerve control, bolus-structure interaction, and structure-mechanics relationship of the tissue. In this study, a finite element model (FEM) was developed to simulate food transport through the esophagus. The FEM consists of three components, i.e., tissue, food bolus and peristaltic wave, as well as the interactions between them. The transport process was simulated as three stages, i.e., the filling of fluid, contraction of circular muscle and traveling of peristaltic wave. It was found that the maximal passive intraluminal pressure due to bolus expansion was in the range of 0.8-10 kPa and it increased with bolus volume and fluid viscosity. It was found that the highest normal and shear stresses were at the inner surface of muscle layer. In addition, the peak pressure required for the fluid flow was predicted to be 1-15 kPa at the bolus tail. The diseases of systemic sclerosis or osteogenesis imperfecta, with the remodeled microstructures and mechanical properties, might induce the malfunction of esophageal transport. In conclusion, the current simulation was demonstrated to be able to capture the main characteristics in the intraluminal pressure and bolus geometry as measured experimentally. Therefore, the finite element model established in this study could be used to further explore the mechanism of esophageal transport in various clinical applications. PMID:17457965
Modeling and Simulation of Anchoring Processess for Small Body Exploration
NASA Technical Reports Server (NTRS)
Quadrelli, Marco B.; Mazahar, Hammad; Negrut, Dan
2012-01-01
This paper describes recent work done in modeling and simulation of anchoring processes in granular media, with applications to anchoring on a Near Earth Object (NEO), where the forces due to interactions with the regolith are much stronger than the local surface gravity field. This effort is part of a larger systems engineering capability developed at JPL to answer key questions, validate requirements, conduct key system and mission trades,and evaluate performance and risk related to NEO operations for any proposed human or robotic missions to a NEO.
TOXIWASP: A DYNAMIC MODEL FOR SIMULATING THE TRANSPORT AND FATE OF TOXIC CHEMICALS IN WATER BODIES
TOXIWASP is a dynamic model for simulating the transport and fate of toxic chemicals in water bodies. Two state variables are simulated: organic chemical and total sediment. The generalized chemical model can be used for problems requiring dynamic transport and loading capabiliti...
Simulations of Bluff Body Flow Interaction for Noise Source Modeling
NASA Technical Reports Server (NTRS)
Khorrami, Medi R.; Lockard David P.; Choudhari, Meelan M.; Jenkins, Luther N.; Neuhart, Dan H.; McGinley, Catherine B.
2006-01-01
The current study is a continuation of our effort to characterize the details of flow interaction between two cylinders in a tandem configuration. This configuration is viewed to possess many of the pertinent flow features of the highly interactive unsteady flow field associated with the main landing gear of large civil transports. The present effort extends our previous two-dimensional, unsteady, Reynolds Averaged Navier-Stokes computations to three dimensions using a quasilaminar, zonal approach, in conjunction with a two-equation turbulence model. Two distinct separation length-to-diameter ratios of L/D = 3.7 and 1.435, representing intermediate and short separation distances between the two cylinders, are simulated. The Mach 0.166 simulations are performed at a Reynolds number of Re = 1.66 105 to match the companion experiments at NASA Langley Research Center. Extensive comparisons with the measured steady and unsteady surface pressure and off-surface particle image velocimetry data show encouraging agreement. Both prominent and some of the more subtle trends in the mean and fluctuating flow fields are correctly predicted. Both computations and the measured data reveal a more robust and energetic shedding process at L/D = 3.7 in comparison with the weaker shedding in the shorter separation case of L/D = 1.435. The vortex shedding frequency based on the computed surface pressure spectra is in reasonable agreement with the measured Strouhal frequency.
Direct numerical simulation of rigid bodies in multiphase flow within an Eulerian framework
NASA Astrophysics Data System (ADS)
Rauschenberger, P.; Weigand, B.
2015-06-01
A new method is presented to simulate rigid body motion in the Volume-of-Fluid based multiphase code Free Surface 3D. The specific feature of the new method is that it works within an Eulerian framework without the need for a Lagrangian representation of rigid bodies. Several test cases are shown to prove the validity of the numerical scheme. The technique is able to conserve the shape of arbitrarily shaped rigid bodies and predict terminal velocities of rigid spheres. The instability of a falling ellipsoid is captured. Multiple rigid bodies including collisions may be considered using only one Volume-of-Fluid variable which allows to simulate the drafting, kissing and tumbling phenomena of two rigid spheres. The method can easily be extended to rigid bodies undergoing phase change processes.
Numerical Simulation of the Sedimentation of a Tripole-like Body in an Incompressible Viscous Fluid
Juarez, L H.; Glowinski, R; Pettitt, Bernard M.
2002-08-01
In this note, we discuss the application of a methodology combining distributed Lagrange multiplier based fictitious domain techniques, finite-element approximations and operator splitting, to the numerical simulation of the motion of a tripole-like rigid body falling in a Newtonian incompressible viscous fluid. The motion of the body is driven by the hydrodynamical forces and gravity. The numerical simulation shows that the distribution of mass of this rigid body and added moment of inertia compared to a simple cylinder (circular or elliptic) plays a significant role on the particle-fluid interaction. Apparently, for the parameters examined, the action of the moving rigid body on the fluid is stronger than the hydrodynamic forces acting on the rigid body.
Fliszkiewicz, Monika; Giejdasz, Karol; Wasielewski, Oskar; Krishnan, Natraj
2012-12-01
The influence of simulated climate change on body weight and depletion of fat body reserves was studied during diapause in the European solitary bee Osmia rufa L. (Hymenoptera: Megachilidae). Insects (females) were reared and collected from outdoor nests from September to March. One cohort of females was weighed and dissected immediately for analyses, whereas another cohort was subjected to simulated warmer temperature (15°C for 7 d) before analyses. A gradual decline in body mass and fat body content was recorded with declining temperatures from September to January in female bees from natural conditions. Temperature increased gradually from January to March with a further decline in body mass and fat body content. The fat body development index dropped from five in September-October (≈ 89% individuals) to four for the period from November to February (≈ 84% individuals) and further to three in March (95% individuals) before emergence. Simulated warmer winter temperature also resulted in a similar decline in body weight and fat body content; however, body weight and fat body content declined faster. The fat body development index dropped to three in December in the majority of individuals and continued at this level until March just before emergence. Taken together, our data indicate an earlier depletion of fat body reserves under simulated climate change conditions that may impact ovarian development and reproductive fitness in O. rufa. PMID:23321111
A simple model simulating a fan as a source of axial and circumferential body forces
Energy Science and Technology Software Center (ESTSC)
2002-07-01
This software can be used in a computational fluids dynamics (CFD) code to represent a fan as a source of axial and circumferential body forces. The combined software can be used effectively in car design analyses that involve many underhood thermal management simulations. FANMOD uses as input the rotational speed of the fan, geometric fan data, and the lift and drag coefficients of the blades, and predicts the body forces generated by the fan inmore » the axial and circumferential directions. These forces can be used as momentum forces in a CFD code to simulate the effect of the fan in an underhood thermal management simulation.« less
A Special-Purpose Computer forN-Body Simulations: GRAPE-2A
NASA Astrophysics Data System (ADS)
Ito, Tomoyoshi; Makino, Junichiro; Fukushige, Toshiyuki; Ebisuzaki, Toshikazu; Okumura, Sachiko K.; Sugimoto, Daiichiro
1993-06-01
We have developed GRAPE-2A, which is a back-end processor used to accelerate simulations of gravitational N-body systems, such as stellar clusters, a proto planetary system, and the structure formation of the universe. GRAPE-2A calculates the forces exerted on one particle from all other particles. The host computer, which is connected to GRAPE-2A through a VME bus, performs other calculations, such as time integration. In the simulation of gravitational N-body systems, almost all of the computing time is consumed in calculating the forces between particles. GRAPE-2A performs this force calculation with a speed that is much faster than that of a general-purpose computer. GRAPE-2A can be used for cosmological N-body simulations with periodic boundary conditions using the Ewald method, and for molecular-dynamics simulations of proteins and crystals. The computational speed of GRAPE-2A is 180 Mflops.
Acellular dermal matrices in breast reconstructions - a literature review.
Skovsted Yde, Simon; Brunbjerg, Mette Eline; Damsgaard, Tine Engberg
2016-08-01
During the last two decades, acellular dermal matrices (ADM) have been more widely used in reconstructive procedures i.e. breast reconstructions. Several, both synthetic and biologic products derived from human, porcine and bovine tissue, have been introduced. Until this point postoperative complications for the acellular dermal matrices, as a group, have been the main focus. The purpose of this literature review is to summarize the current knowledge on the each biologic product used in breast reconstructions, including product specific complication frequencies. A systematic search of the literature was performed in the PubMed and EMBASE databases, identifying 55 relevant articles, mainly evidence level III. AlloDerm seems to be associated with severe complicating matters in the reconstructive process compared to other products. This could be due to the higher number of investigating studies relative to the others. The surgical area faces certain challenges comparing results, due to surgical variance, the data collection and follow-up. More well-defined guidelines and more high-evidence randomized studies could increase the overall level of evidence in this area. PMID:26881927
Unsteady transonic flow simulation on a full-span-wing-body configuration
NASA Technical Reports Server (NTRS)
Guruswamy, Guru P.; Goorjian, Peter M.
1987-01-01
The presence of a body influences both the aerodynamic and aeroelastic performance of wings. Such effects are more pronounced in the transonic regime. To accurately account for the effect of the body, particularly when the wings are experiencing asymmetric modal motions, it is necessary to model the full configuration in the nonlinear transonic regime. In this study, full-span-wing-body configurations are simulated for the first time by a theoretical method that uses the unsteady potential equations based on the small-disturbance theory. The body geometry is modeled exactly as the physical shape, instead of as a rectangular box, which has been done in the past. Steady pressure computations for wing-body configurations compare well with the available experimental data. Unsteady pressure computations when the wings are oscillating in asymmetric modes show significant influence of the body.
Simulation and study of stratified flows around finite bodies
NASA Astrophysics Data System (ADS)
Gushchin, V. A.; Matyushin, P. V.
2016-06-01
The flows past a sphere and a square cylinder of diameter d moving horizontally at the velocity U in a linearly density-stratified viscous incompressible fluid are studied. The flows are described by the Navier-Stokes equations in the Boussinesq approximation. Variations in the spatial vortex structure of the flows are analyzed in detail in a wide range of dimensionless parameters (such as the Reynolds number Re = Ud/ ν and the internal Froude number Fr = U/( Nd), where ν is the kinematic viscosity and N is the buoyancy frequency) by applying mathematical simulation (on supercomputers of Joint Supercomputer Center of the Russian Academy of Sciences) and three-dimensional flow visualization. At 0.005 < Fr < 100, the classification of flow regimes for the sphere (for 1 < Re < 500) and for the cylinder (for 1 < Re < 200) is improved. At Fr = 0 (i.e., at U = 0), the problem of diffusion-induced flow past a sphere leading to the formation of horizontal density layers near the sphere's upper and lower poles is considered. At Fr = 0.1 and Re = 50, the formation of a steady flow past a square cylinder with wavy hanging density layers in the wake is studied in detail.
NVision: A 3D Visualization Environment for N-Body Simulations
NASA Astrophysics Data System (ADS)
Markiel, J. A.
2000-05-01
NVision: A 3D Visualization Environment for N-Body Simulations We are developing a set of packages for 3D visualization and analysis of our numerical N-body simulations. These tools are intended to be generalizable to a wide range of related problems including cosmological, planetary dynamics, and molecular dynamics simulations. The applications and source code will be fully available to the community. To prototype this project we have adopted the Java platform with the newly released Java3D extension to take advantage of its portability, object-oriented environment, and availability of extensive documentation and class libraries. We will describe the goals and design principles of the project and demo the currently implemented features, including visualization of cosmological simulations and the simulated collision of two rubble-pile asteroids. This research is supported by NSF grants AST99-73209 and AST99-79891.
Numerical simulation of tidal evolution of a viscoelastic body modelled with a mass-spring network
NASA Astrophysics Data System (ADS)
Frouard, Julien; Quillen, Alice C.; Efroimsky, Michael; Giannella, David
2016-05-01
We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number k2 divided by the tidal quality factor Q) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of k2/Q, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalized to inhomogeneous and/or non-spherical bodies.
Monte Carlo simulation of efficient data acquisition for an entire-body PET scanner
NASA Astrophysics Data System (ADS)
Isnaini, Ismet; Obi, Takashi; Yoshida, Eiji; Yamaya, Taiga
2014-07-01
Conventional PET scanners can image the whole body using many bed positions. On the other hand, an entire-body PET scanner with an extended axial FOV, which can trace whole-body uptake images at the same time and improve sensitivity dynamically, has been desired. The entire-body PET scanner would have to process a large amount of data effectively. As a result, the entire-body PET scanner has high dead time at a multiplex detector grouping process. Also, the entire-body PET scanner has many oblique line-of-responses. In this work, we study an efficient data acquisition for the entire-body PET scanner using the Monte Carlo simulation. The simulated entire-body PET scanner based on depth-of-interaction detectors has a 2016-mm axial field-of-view (FOV) and an 80-cm ring diameter. Since the entire-body PET scanner has higher single data loss than a conventional PET scanner at grouping circuits, the NECR of the entire-body PET scanner decreases. But, single data loss is mitigated by separating the axially arranged detector into multiple parts. Our choice of 3 groups of axially-arranged detectors has shown to increase the peak NECR by 41%. An appropriate choice of maximum ring difference (MRD) will also maintain the same high performance of sensitivity and high peak NECR while at the same time reduces the data size. The extremely-oblique line of response for large axial FOV does not contribute much to the performance of the scanner. The total sensitivity with full MRD increased only 15% than that with about half MRD. The peak NECR was saturated at about half MRD. The entire-body PET scanner promises to provide a large axial FOV and to have sufficient performance values without using the full data.
The simulation of electrostatic coupling intra-body communication based on finite-element models
NASA Astrophysics Data System (ADS)
Song, Yong; Yang, Guang; Hao, Qun; Wang, Ming
2008-12-01
Intra-body Body Communication (IBC) is a communication technology in which human body is used as a signal transmission medium. Due to its unique characters, IBC technology is proposed as a novel and promising technology for personal area network (PAN), computer network access, implant biomedical monitoring, human energy transmission, etc. In this paper, investigation has been done in the computer simulation of the electrostatic coupling IBC by using the developed finite-element models, in which (1) the incidence and reflection of electronic signal in the upper arm model were analyzed by using the theory of electromagnetic wave, (2) the finite-element models of electrostatic coupling IBC were developed by using the electromagnetic analysis package of ANSYS software, (3) the signal attenuation of electrostatic coupling IBC were simulated under the conditions of different signal frequency, electrodes direction, electrodes size and transmission distance. Finally, some important conclusions are deduced on the basis of simulation results.
Multibody Simulation Software Testbed for Small-Body Exploration and Sampling
NASA Technical Reports Server (NTRS)
Acikmese, Behcet; Blackmore, James C.; Mandic, Milan
2011-01-01
G-TAG is a software tool for the multibody simulation of a spacecraft with a robotic arm and a sampling mechanism, which performs a touch-and-go (TAG) maneuver for sampling from the surface of a small celestial body. G-TAG utilizes G-DYN, a multi-body simulation engine described in the previous article, and interfaces to controllers, estimators, and environmental forces that affect the spacecraft. G-TAG can easily be adapted for the analysis of the mission stress cases to support the design of a TAG system, as well as for comprehensive Monte Carlo simulations to analyze and evaluate a particular TAG system design. Any future small-body mission will benefit from using G-TAG, which has already been extensively used in Comet Odyssey and Galahad Asteroid New Frontiers proposals.
2HOT: An Improved Parallel Hashed Oct-Tree N-Body Algorithm for Cosmological Simulation
Warren, Michael S.
2014-01-01
We report on improvements made over the past two decades to our adaptive treecode N-body method (HOT). A mathematical and computational approach to the cosmological N-body problem is described, with performance and scalability measured up to 256k (2 18 ) processors. We present error analysis and scientific application results from a series of more than ten 69 billion (4096 3 ) particle cosmological simulations, accounting for 4×10 20 floating point operations. These results include the first simulations using the new constraints on the standard model of cosmology from the Planck satellite. Our simulations set a new standard for accuracymore » and scientific throughput, while meeting or exceeding the computational efficiency of the latest generation of hybrid TreePM N-body methods.« less
Simulations of dipolar fluids using effective many-body isotropic interactions.
Sindt, Julien O; Camp, Philip J
2015-07-14
The partition function of a system with pairwise-additive anisotropic dipole-dipole interactions is equal to that of a hypothetical system with many-body isotropic interactions [G. Stell, Phys. Rev. Lett. 32, 286 (1974)]. The effective many-body interactions contain n-body contributions of all orders. Each contribution is known as an expansion in terms of the particle-particle distances r, and the coefficients are temperature dependent. The leading-order two-body term is the familiar -r(-6) attraction, and the leading-order three-body term is equivalent to the Axilrod-Teller interaction. In this work, a fluid of particles with the leading-order two-body and three-body interactions is compared to an equivalent dipolar soft-sphere fluid. Molecular simulations are used to determine the conditions under which the effective many-body interactions reproduce the fluid-phase structures of the dipolar system. The effective many-body interaction works well at moderately high temperatures but fails at low temperatures where particle chaining is expected to occur. It is shown that an adjustment of the coefficients of the two-body and three-body terms leads to a good description of the structure of the dipolar fluid even in the chaining regime, due primarily to the ground-state linear configuration of the three-body Axilrod-Teller interaction. The vapor-liquid phase diagrams of systems with different Axilrod-Teller contributions are determined. As the strength of the three-body interaction is increased, the critical temperature and density both decrease and disappear completely above a threshold strength, where chaining eventually suppresses the condensation transition. PMID:26178112
Simulation of SAR in the Human Body to Determine Effects of RF Heating
NASA Astrophysics Data System (ADS)
Michiyama, Tetsuyuki; Nikawa, Yoshio
The body area network (BAN) has attracted attention because of its potential for high-grade wireless communication technology and its safety and high durability. Also, human area transmission of a BAN propagating at an ultra-wide band (UWB) has been demonstrated recently. When considering the efficiency of electromagnetic (EM) propagation inside the human body for BAN and hyperthermia treatment using RF, it is important to determine the mechanism of EM dissipation in the human body. A body heating system for hyperthermia must deposit EM energy deep inside the body. Also, it is important that the EM field generated by the implant system is sufficiently strong. In this study, the specific absorption rate (SAR) distribution is simulated using an EM simulator to consider the biological transmission mechanism and its effects. To utilize the EM field distribution using an implant system for hyperthermia treatment, the SAR distribution inside the human body is simulated. As a result, the SAR distribution is concentrated on the surface of human tissue, the muscle-bolus interface, the pancreas, the stomach, the spleen and the regions around bones. It can also be concentrated in bone marrow and cartilage. From these results, the appropriate location for the implant system is revealed on the basis of the current distribution and differences in the wave impedance of interfacing tissues. The possibility of accurate data transmission and suitable treatment planning is confirmed.
A generic multi-flex-body dynamics, controls simulation tool for space station
NASA Technical Reports Server (NTRS)
London, Ken W.; Lee, John F.; Singh, Ramen P.; Schubele, Buddy
1991-01-01
An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.
NASA Astrophysics Data System (ADS)
Song, Yong; Chu, Yingfang; Zhang, Kai; Kang, Bangzhi; Hao, Qun
2010-11-01
The simulation based on the finite-element (FE) method plays an important role in the investigation of the intra-body communication (IBC). In this paper, the method for modeling the whole human body based on the finite-element method is proposed, while a finite-element model of the whole human body used for the simulations of the waveguide intra-body communication has been developed. Finally, the simulations of the waveguide IBC with different signal transmission paths have been achieved by using the developed finite-element model. Moreover, both the potential distributions and the signal attenuations of the simulation results are discussed in detail, which indicate that the proposed method and model offer the significant advantages in the theoretical analysis and the system design of the waveguide intra-body communication.
Robotic simulation of flexible-body spacecraft dynamics in a satellite servicing testbed
NASA Astrophysics Data System (ADS)
Brannan, Justin Cory
Satellite failures that once led to end-of-life may eventually be addressed using robotic servicing platforms. The ability to model and simulate the physical interaction between two free-floating spacecraft is a key aspect of robotic servicing, and understanding how large appendages such as solar panels, antenna arrays and booms affect the combined system dynamics may be critical to mission operations. This research presents a model of the coupled rigid- and flexible-body satellite dynamics that can be implemented on a robotic satellite simulator. The coupled dynamics are validated against a commercially available dynamics software package, and robot hardware-in-the-loop tests are conducted to demonstrate how the dynamics model is able to predict the response of a robot mass simulator outfitted with physical appendages. Through both validation efforts, a flexible-body simulation is developed to observe the resulting dynamics of a given satellite system on-orbit.
Application of acellular dermis and autograft on burns and scars.
Ramos Duron, L E; Martínez Pardo, M E; Olivera Zavaleta, V; Silva Diaz, T; Reyes Frías, M L; Luna Zaragoza, D
1999-01-01
The cases of two patients with burns treated with dermis allograft and of one patient for lip reconstructive aesthetic filling treated with less than one mm3 of radiosterilised acellular dermis are presented. This paper emphasizes the treatment with radiosterilised dermal grafts with a permanent character so far. Hospitals, therefore, can satisfy the demand for this kind of tissue in the case of disaster and patients with serious injuries. In the cases cited, histocompatibility analysis was not required, thus having the advantage of long-time storage of the radiosterilised dermis used on these patients. Neither inflammatory reaction nor acute phase re-absorption were observed. Moreover, shrink (contract) healing was diminished. After two years, the results are still satisfactory. PMID:10853787
Fast Generation of Ensembles of Cosmological N-Body Simulations via Mode-Resampling
Schneider, M D; Cole, S; Frenk, C S; Szapudi, I
2011-02-14
We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the large-scale (Gaussian-distributed) Fourier modes in a periodic N-body simulation box in a manner that properly accounts for the nonlinear mode-coupling between large and small scales. We find that our method for adding new large-scale mode realizations recovers the nonlinear power spectrum to sub-percent accuracy on scales larger than about half the Nyquist frequency of the simulation box. Using 20 N-body simulations, we obtain a power spectrum covariance matrix estimate that matches the estimator from Takahashi et al. (from 5000 simulations) with < 20% errors in all matrix elements. Comparing the rates of convergence, we determine that our algorithm requires {approx}8 times fewer simulations to achieve a given error tolerance in estimates of the power spectrum covariance matrix. The degree of success of our algorithm indicates that we understand the main physical processes that give rise to the correlations in the matter power spectrum. Namely, the large-scale Fourier modes modulate both the degree of structure growth through the variation in the effective local matter density and also the spatial frequency of small-scale perturbations through large-scale displacements. We expect our algorithm to be useful for noise modeling when constraining cosmological parameters from weak lensing (cosmic shear) and galaxy surveys, rescaling summary statistics of N-body simulations for new cosmological parameter values, and any applications where the influence of Fourier modes larger than the simulation size must be accounted for.
Kato, T; Goshima, T; Nakajima, N; Kaku, H; Arimoto, Y; Hayashi, F
1989-12-01
To evaluate the vaccine efficacy of an acellular pertussis vaccine which has been in clinical use in Japan since 1981, a retrospective study was performed by a questionnaire survey of secondary pertussis attacks through family contact in 146 children with pertussis diagnosed in the period from January 1981 through May 1988. In this study, acellular vaccine made by Takeda Pharmaceutical Company, which contains a high level of FHA (filamentous hemagglutinin), a low level of PT (pertussis toxin) and a small amount of agglutinogen, was evaluated. Secondary pertussis attacks through family contact were found in 17 of 29 siblings (58.6%) not immunized with pertussis vaccine. On the other hand, 27 siblings immunized with Takeda's acellular vaccine were exposed to pertussis through family contact and a secondary attack was seen in only one of them (3.7%). The present study revealed an efficacy rate of 93.7% for Takeda's acellular pertussis vaccine. PMID:2516396
Kato, T; Matsuyoshi, S; Goshima, T; Nakajima, N; Yamamoto, H; Arimoto, Y; Kaku, H; Hayashi, F
1989-09-01
To evaluate the vaccine efficacy of acellular pertussis vaccine which has been in clinical use in Japan since 1981, a retrospective study was made by a questionnaire from secondary pertussis attack through family contact in 149 children with pertussis diagnosed in the period from January 1981 through May 1988. In this study, Takeda's acellular vaccine which contains a high level of FHA, low level of PT and a small amount of agglutinogen, was evaluated. Secondary pertussis attacks through family contact were found in 17 of 29 siblings (58.6%) not immunized with pertussis vaccine. On the other hand of the siblings immunized with Takeda's acellular vaccine 27 were exposed to pertussis through family contact and a secondary attack was seen in only one of them (3.4%). The present study revealed an efficacy rate of 94.2% for the Takeda's acellular pertussis vaccine. PMID:2509597
Kato, T; Kaku, H; Arimoto, Y
1988-01-01
To evaluate the vaccine efficacy of an acellular pertussis vaccine which has been in clinical use in Japan since 1981, a retrospective study was performed by a questionnaire survey of secondary pertussis attacks through family contact in 146 children with pertussis diagnosed in the period from January 1981 through May 1988. In this study, Takeda's acellular vaccine which contains a high level of FHA, low level of PT and a small amount of agglutinogen, was evaluated. Secondary pertussis attacks through family contact were found in 17 of 27 siblings (62.9%) not immunized with pertussis vaccine. On the other hand, 26 siblings immunized with Takeda's acellular vaccine were exposed to pertussis through family contact and a secondary attack was seen in only one of them (3.8%). The present study revealed an efficacy rate of 93.9% for Takeda's acellular pertussis vaccine. PMID:3078808
Rigid body dynamics approach to Stokesian dynamics simulations of nonspherical particles
NASA Astrophysics Data System (ADS)
Kutteh, Ramzi
2010-05-01
We describe an algorithm for performing Stokesian dynamics (SD) simulations of suspensions of arbitrary shape rigid particles with hydrodynamic interactions, modeled as rigid groups of spheres, the hydrodynamic mobility matrix of which is accurately computable by several established schemes for spheres. The algorithm is based on Stokesian rigid body equations of translational and rotational motion, which we have derived by an approach formally analogous to that of Newtonian rigid body dynamics. Particle orientation is represented in terms of Euler parameters (quaternion of rotation). This rigid body SD algorithm (RBSDA) complements recently described constraint SD algorithms [R. Kutteh, J. Chem. Phys. 119, 9280 (2003); R. Kutteh, Phys. Rev. E 69, 011406 (2004)], over which it offers the same computational advantages in imposing total rigidity that the basic rigid body molecular dynamics (MD) algorithm offers over constraint MD algorithms. We show that SD simulation results generated with the RBSDA, in bounded and unbounded geometries, agree very well with those from experiment and other SD and non-SD methods, and are numerically identical to those from a constraint SD algorithm, HSHAKE. Finally, for completeness we also describe a third (additional to the constraint SD and rigid body SD approaches) more traditional approach for SD simulations of arbitrary shape rigid particles modeled as rigid groups of spheres.
Rigid body dynamics approach to Stokesian dynamics simulations of nonspherical particles.
Kutteh, Ramzi
2010-05-01
We describe an algorithm for performing Stokesian dynamics (SD) simulations of suspensions of arbitrary shape rigid particles with hydrodynamic interactions, modeled as rigid groups of spheres, the hydrodynamic mobility matrix of which is accurately computable by several established schemes for spheres. The algorithm is based on Stokesian rigid body equations of translational and rotational motion, which we have derived by an approach formally analogous to that of Newtonian rigid body dynamics. Particle orientation is represented in terms of Euler parameters (quaternion of rotation). This rigid body SD algorithm (RBSDA) complements recently described constraint SD algorithms [R. Kutteh, J. Chem. Phys. 119, 9280 (2003); R. Kutteh, Phys. Rev. E 69, 011406 (2004)], over which it offers the same computational advantages in imposing total rigidity that the basic rigid body molecular dynamics (MD) algorithm offers over constraint MD algorithms. We show that SD simulation results generated with the RBSDA, in bounded and unbounded geometries, agree very well with those from experiment and other SD and non-SD methods, and are numerically identical to those from a constraint SD algorithm, HSHAKE. Finally, for completeness we also describe a third (additional to the constraint SD and rigid body SD approaches) more traditional approach for SD simulations of arbitrary shape rigid particles modeled as rigid groups of spheres. PMID:20459156
The Rufous Hummingbird in hovering flight -- full-body 3D immersed boundary simulation
NASA Astrophysics Data System (ADS)
Ferreira de Sousa, Paulo; Luo, Haoxiang; Bocanegra Evans, Humberto
2009-11-01
Hummingbirds are an interesting case study for the development of micro-air vehicles since they combine the high flight stability of insects with the low metabolic power per unit of body mass of bats, during hovering flight. In this study, simulations of a full-body hummingbird in hovering flight were performed at a Reynolds number around 3600. The simulations employ a versatile sharp-interface immersed boundary method recently enhanced at our lab that can treat thin membranes and solid bodies alike. Implemented on a Cartesian mesh, the numerical method allows us to capture the vortex dynamics of the wake accurately and efficiently. The whole-body simulation will allow us to clearly identify the three general patterns of flow velocity around the body of the hummingbird referred in Altshuler et al. (Exp Fluids 46 (5), 2009). One focus of the current study is to understand the interaction between the wakes of the two wings at the end of the upstroke, and how the tail actively defects the flow to contribute to pitch stability. Another focus of the study will be to identify the pair of unconnected loops underneath each wing.
Advances in Chimera Grid Tools for Multi-Body Dynamics Simulations and Script Creation
NASA Technical Reports Server (NTRS)
Chan, William M.
2004-01-01
This viewgraph presentation contains information about (1) Framework for multi-body dynamics - Geometry Manipulation Protocol (GMP), (2) Simulation procedure using Chimera Grid Tools (CGT) and OVERFLOW-2 (3) Further recent developments in Chimera Grid Tools OVERGRID, Grid modules, Script library and (4) Future work.
Dynamic Simulation and Static Matching for Action Prediction: Evidence from Body Part Priming
ERIC Educational Resources Information Center
Springer, Anne; Brandstadter, Simone; Prinz, Wolfgang
2013-01-01
Accurately predicting other people's actions may involve two processes: internal real-time simulation (dynamic updating) and matching recently perceived action images (static matching). Using a priming of body parts, this study aimed to differentiate the two processes. Specifically, participants played a motion-controlled video game with…
RAY-RAMSES: a code for ray tracing on the fly in N-body simulations
NASA Astrophysics Data System (ADS)
Barreira, Alexandre; Llinares, Claudio; Bose, Sownak; Li, Baojiu
2016-05-01
We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.
Phantom-GRAPE: SIMD accelerated numerical library for N-body simulations
NASA Astrophysics Data System (ADS)
Tanikawa, Ataru; Yoshikawa, Kohji; Nitadori, Keigo; Okamoto, Takashi
2012-09-01
Phantom-GRAPE is a numerical software library to accelerate collisionless N-body simulation with SIMD instruction set on x86 architecture. The Newton's forces and also central forces with an arbitrary shape f(r), which have a finite cutoff radius r_cut (i.e. f(r)=0 at r>r_cut), can be quickly computed.
NASA Astrophysics Data System (ADS)
Graham, Hannah Robyn
In order to be able to qualify and quantify radiation exposure in terms of dose, a Fastscan whole body counter must be calibrated correctly. Current calibration methods do not take the full range of body types into consideration when creating efficiency curve calibrations. The goal of this work is the creation of a Monte Carlo (MCNP) model, that allows the simulation of efficiency curves for a diverse population of subjects. Models were created for both the Darlington and the Pickering Fastscan WBCs, and the simulations were benchmarked against experimental results with good agreement. The Pickering Fastscan was found to have agreement to within +/-9%, and the Darlington Fastscan had agreement to within +/-11%. Further simulations were conducted to investigate the effects of increased body fat on the detected activity, as well as locating the position of external contamination using front/back ratios of activity. Simulations were also conducted to create efficiency calibrations that had good agreement with the manufacturer's efficiency curves. The work completed in this thesis can be used to create efficiency calibration curves for unique body compositions in the future.
Dynamically coupled fluid body interactions in vorticity-based numerical simulations
NASA Astrophysics Data System (ADS)
Eldredge, Jeff D.
2008-11-01
A novel method is presented for robustly simulating coupled dynamics in fluid-body interactions with vorticity-based flow solvers. In this work, the fluid dynamics are simulated with a viscous vortex particle method. In the first substep of each time increment, the fluid convective and diffusive processes are treated, while a predictor is used to independently advance the body configuration. An iterative corrector is then used to simultaneously remove the spurious slip - via vorticity flux - and compute the end-of-step body configuration. Fluid inertial forces are isolated and combined with body inertial terms to ensure robust treatment of dynamics for bodies of arbitrary mass. The method is demonstrated for dynamics of articulated rigid bodies, including a falling cylinder, flow-induced vibration of a circular cylinder and free swimming of a three-link 'fish'. The error and momentum conservation properties of the method are explored. In the case of the vibrating cylinder, comparison with previous work demonstrates good agreement.
Simulation and Analyses of Multi-Body Separation in Launch Vehicle Staging Environment
NASA Technical Reports Server (NTRS)
Pamadi, Bandu N.; Hotchko, Nathaniel J.; Samareh, Jamshid; Covell, Peter F.; Tartabini, Paul V.
2006-01-01
The development of methodologies, techniques, and tools for analysis and simulation of multi-body separation is critically needed for successful design and operation of next generation launch vehicles. As a part of this activity, ConSep simulation tool is being developed. ConSep is a generic MATLAB-based front-and-back-end to the commercially available ADAMS. solver, an industry standard package for solving multi-body dynamic problems. This paper discusses the 3-body separation capability in ConSep and its application to the separation of the Shuttle Solid Rocket Boosters (SRBs) from the External Tank (ET) and the Orbiter. The results are compared with STS-1 flight data.
Stellar evolution in N-body simulations of disk galaxies. I
NASA Technical Reports Server (NTRS)
Comins, N. F.
1983-01-01
The Kalnajs (1972, 1976) Omega models of global mass and velocity distributions are employed in the present two-dimensional N-body simulation, which allows for a spectrum of particle masses, stellar explosions, explosion remnant interactions with an interstellar medium, and the creation of new stars from the gas. Two sequences of runs using the Omega values of 0.8 and 0.9 examine the separate and combined effects of particle mass distribution, the gravitational influence of an interstellar gas distribution on the N-body particles, and stellar evolution, allowing for stellar explosions and star formation from the gas. It is found that both Omega values' nonequilibrium results dramatically change when evolution is allowed to occur. These results call for more realistic coupled N-body and evolution simulations in order to improve the understanding of disk galaxy evolution.
Three-body interactions in complex fluids: Virial coefficients from simulation finite-size effects
Ashton, Douglas J.; Wilding, Nigel B.
2014-06-28
A simulation technique is described for quantifying the contribution of three-body interactions to the thermodynamical properties of coarse-grained representations of complex fluids. The method is based on a new approach for determining virial coefficients from the measured volume-dependent asymptote of a certain structural function. By comparing the third virial coefficient B{sub 3} for a complex fluid with that of an approximate coarse-grained model described by a pair potential, three body effects can be quantified. The strategy is applicable to both Molecular Dynamics and Monte Carlo simulation. Its utility is illustrated via measurements of three-body effects in models of star polymers and in highly size-asymmetrical colloid-polymer mixtures.
Cosmological N-body Simulation of Galaxy and Large-Scale Structure Formation: The Gravity Frontier
NASA Astrophysics Data System (ADS)
Klypin, Anatoly
2015-04-01
One of the first N-body simulations done almost 50 years ago had only 200 self-gravitating particles. Even this first baby step made substantial impact on understanding how astronomical objects should form. Now powerful supercomputers and new algorithms allow astronomers produce N-body simulations that employ up to a trillion dark matter particles and produce vital theoretical predictions regarding formation, evolution, structure and statistics of objects ranging from dwarf galaxies to clusters and superclusters of galaxies. With only gravity involved in these theoretical models, one would naively expect that by now we should know everything we need about N-body dynamics of cosmological fluctuations. Not the case. It appears that the Universe was not cooperative and gave us divergencies in the initial conditions generated during the Inflation epoch and subsequent expansion of the Universe - the infinite phase-space density and divergent density fluctuations. Ever increasing observational demands on statistics and accuracy of theoretical predictions is another driving force for more realistic and larger N-body simulations. Large current and new planned observational projects such as BOSS, eBOSS, Euclid, LSST will bring information on spatial distribution, motion, and properties of millions of galaxies at different redshifts. Direct simulations of evolution of gas and formation of stars for millions of forming galaxies will not be available for years leaving astronomers with the only option - to develop methods to combine large N-body simulations with models of galaxy formation to produce accurate theoretical predictions. I will discuss the current status of the field and directions of its development.
Global High-resolution N-body Simulation of Planet Formation. I. Planetesimal-driven Migration
NASA Astrophysics Data System (ADS)
Kominami, J. D.; Daisaka, H.; Makino, J.; Fujimoto, M.
2016-03-01
We investigated whether outward planetesimal-driven migration (PDM) takes place or not in simulations when the self-gravity of planetesimals is included. We performed N-body simulations of planetesimal disks with a large width (0.7-4 au) that ranges over the ice line. The simulations consisted of two stages. The first-stage simulations were carried out to see the runaway growth phase using the planetesimals of initially the same mass. The runaway growth took place both at the inner edge of the disk and at the region just outside the ice line. This result was utilized for the initial setup of the second-stage simulations, in which the runaway bodies just outside the ice line were replaced by the protoplanets with about the isolation mass. In the second-stage simulations, the outward migration of the protoplanet was followed by the stopping of the migration due to the increase of the random velocity of the planetesimals. Owing to this increase of random velocities, one of the PDM criteria derived in Minton & Levison was broken. In the current simulations, the effect of the gas disk is not considered. It is likely that the gas disk plays an important role in PDM, and we plan to study its effect in future papers.
Development of a Multi-body Dynamics Simulation Tool for Tracked Vehicles
NASA Astrophysics Data System (ADS)
Huh, Kunsoo; Choi, Jinhwan; Yoo, Honghee
Multi-body Dynamics simulation of tracked vehicles is very useful not only for the analysis of dynamic behaviors, but also for the performance evaluation of the chassis controllers. The track tension is closely related to the maneuverability of tracked vehicles and the durability of tracks and suspension systems. In order to minimize the excessive load on the tracks and to prevent the peal-off of tracks from the sprocket, the track tension needs to be maintained at the optimum level throughout the maneuver. In this paper, a co-simulation tool is developed such that the performance of the track tension control system can be investigated for various maneuvering tasks. The MBD (multi-body dynamics) vehicle model for tracked vehicles is very complicated (189 bodies, 36 revolute joints, 152 bushing elements and 954 degrees of freedom) and cannot be easily implemented in commercial software. Besides, the track tension controller based on fuzzy logic can be easily constructed in the commercial control software. Therefore, co-simulation methodology is proposed so that the designed tension controller is interfaced into the MBD simulation software. The performance of the tension control system is verified through the proposed co-simulation tool.
On the simulation of indistinguishable fermions in the many-body Wigner formalism
Sellier, J.M. Dimov, I.
2015-01-01
The simulation of quantum systems consisting of interacting, indistinguishable fermions is an incredible mathematical problem which poses formidable numerical challenges. Many sophisticated methods addressing this problem are available which are based on the many-body Schrödinger formalism. Recently a Monte Carlo technique for the resolution of the many-body Wigner equation has been introduced and successfully applied to the simulation of distinguishable, spinless particles. This numerical approach presents several advantages over other methods. Indeed, it is based on an intuitive formalism in which quantum systems are described in terms of a quasi-distribution function, and highly scalable due to its Monte Carlo nature. In this work, we extend the many-body Wigner Monte Carlo method to the simulation of indistinguishable fermions. To this end, we first show how fermions are incorporated into the Wigner formalism. Then we demonstrate that the Pauli exclusion principle is intrinsic to the formalism. As a matter of fact, a numerical simulation of two strongly interacting fermions (electrons) is performed which clearly shows the appearance of a Fermi (or exchange–correlation) hole in the phase-space, a clear signature of the presence of the Pauli principle. To conclude, we simulate 4, 8 and 16 non-interacting fermions, isolated in a closed box, and show that, as the number of fermions increases, we gradually recover the Fermi–Dirac statistics, a clear proof of the reliability of our proposed method for the treatment of indistinguishable particles.
GalevNB: a conversion from N-body simulations to observations
NASA Astrophysics Data System (ADS)
Pang, Xiao-Ying; Olczak, Christoph; Guo, Di-Feng; Spurzem, Rainer; Kotulla, Ralf
2016-03-01
We present GalevNB (Galev for N-body simulations), a utility that converts fundamental stellar properties of N-body simulations into observational properties using the GALEV (GAlaxy EVolutionary synthesis models) package, and allowing direct comparisons between observations and N-body simulations. It works by converting fundamental stellar properties, such as stellar mass, temperature, luminosity and metallicity into observational magnitudes for a variety of filters used by mainstream instruments/telescopes, such as HST, ESO, SDSS, 2MASS, etc., and into spectra that span the range from far-UV (90 Å) to near-IR (160 μm). As an application, we use GalevNB to investigate the secular evolution of the spectral energy distribution (SED) and color magnitude diagram (CMD) of a simulated star cluster over a few hundred million years. With the results given by GalevNB we discover a UV-excess in the SED of the cluster over the whole simulation time. We also identify four candidates that contribute to the FUV peak: core helium burning stars, second asymptotic giant branch (AGB) stars, white dwarfs and naked helium stars.
On the simulation of indistinguishable fermions in the many-body Wigner formalism
NASA Astrophysics Data System (ADS)
Sellier, J. M.; Dimov, I.
2015-01-01
The simulation of quantum systems consisting of interacting, indistinguishable fermions is an incredible mathematical problem which poses formidable numerical challenges. Many sophisticated methods addressing this problem are available which are based on the many-body Schrödinger formalism. Recently a Monte Carlo technique for the resolution of the many-body Wigner equation has been introduced and successfully applied to the simulation of distinguishable, spinless particles. This numerical approach presents several advantages over other methods. Indeed, it is based on an intuitive formalism in which quantum systems are described in terms of a quasi-distribution function, and highly scalable due to its Monte Carlo nature. In this work, we extend the many-body Wigner Monte Carlo method to the simulation of indistinguishable fermions. To this end, we first show how fermions are incorporated into the Wigner formalism. Then we demonstrate that the Pauli exclusion principle is intrinsic to the formalism. As a matter of fact, a numerical simulation of two strongly interacting fermions (electrons) is performed which clearly shows the appearance of a Fermi (or exchange-correlation) hole in the phase-space, a clear signature of the presence of the Pauli principle. To conclude, we simulate 4, 8 and 16 non-interacting fermions, isolated in a closed box, and show that, as the number of fermions increases, we gradually recover the Fermi-Dirac statistics, a clear proof of the reliability of our proposed method for the treatment of indistinguishable particles.
Fu, Yao E-mail: jhsong@cec.sc.edu; Song, Jeong-Hoon E-mail: jhsong@cec.sc.edu
2014-08-07
Hardy stress definition has been restricted to pair potentials and embedded-atom method potentials due to the basic assumptions in the derivation of a symmetric microscopic stress tensor. Force decomposition required in the Hardy stress expression becomes obscure for multi-body potentials. In this work, we demonstrate the invariance of the Hardy stress expression for a polymer system modeled with multi-body interatomic potentials including up to four atoms interaction, by applying central force decomposition of the atomic force. The balance of momentum has been demonstrated to be valid theoretically and tested under various numerical simulation conditions. The validity of momentum conservation justifies the extension of Hardy stress expression to multi-body potential systems. Computed Hardy stress has been observed to converge to the virial stress of the system with increasing spatial averaging volume. This work provides a feasible and reliable linkage between the atomistic and continuum scales for multi-body potential systems.
Flow Simulation of N3-X Hybrid Wing-Body Configuration
NASA Technical Reports Server (NTRS)
Kim, Hyoungjin; Liou, Meng-Sing
2013-01-01
System studies show that a N3-X hybrid wing-body aircraft with a turboelectric distributed propulsion system using a mail-slot inlet/nozzle nacelle can meet the environmental and performance goals for N+3 generation transports (three generations beyond the current air transport technology level) set by NASA s Subsonic Fixed Wing Project. In this study, a Navier-Stokes flow simulation of N3-X on hybrid unstructured meshes was conducted, including the mail-slot propulsor. The geometry of the mail-slot propulsor was generated by a CAD (Computer-Aided Design)-free shape parameterization. A body force approach was used for a more realistic and efficient simulation of the turning and loss effects of the fan blades and the inlet-fan interactions. Flow simulation results of the N3-X demonstrates the validity of the present approach.
Particle-Based Simulation of Shock-Induced Deformation of Elastic Bodies
NASA Astrophysics Data System (ADS)
Sakamura, Y.; Sugimoto, T.; Nakayama, K.
Shock-induced deformations of solid bodies are of practical interest to those who are concerned with explosive processing of materials, demolition of buildings, precautions against accidental explosions, etc. In order to simulate the shock-induced deformations of solid bodies, a large number of numerical codes based on continuum mechanics, which are called hydrocodes, have been developed so far [1, 2]. When the amount of deformation is relatively small, Lagrangian hydrocodes have been used to simulate the dynamic response of shock-loaded materials. When the deformation is large, Eulerian hydrocodes have been utilized instead. This is because the computational grids distorted along with the deformation of materials in the Lagrangian approach make the simulations either inaccurate or unstable, while the Eulerian approach where grids are fixed in space can handle such large deformations of materials. On the contrary, material interfaces that are precisely defined in the Lagrangian approach are not traced exactly in the Eulerian one.
NASA Astrophysics Data System (ADS)
Inácio, O.; Antunes, J.; Wright, M. C. M.
2008-02-01
Most theoretical studies of bowed-string instruments deal with isolated strings, pinned on fixed supports. In others, the instrument body dynamics have been accounted by using extremely simplified models of the string-body interaction through the instrument bridge. Such models have, nevertheless, been instrumental to the understanding of a very common and musically undesirable phenomenon known as the wolf note—a strong beating interplay between string and body vibrations. Cellos, bad and good, are particularly prone to this problem. In previous work, a computational method that allows efficient time-domain modelling of bowed strings based on a modal approach has been introduced. This has been extended to incorporate the complex dynamics of real-life instrument bodies, and their coupling to the string motions, using experimental dynamical body data. The string is modelled using its unconstrained modes, assuming pinned-pinned boundary conditions at the tailpiece and the nut. At the intermediary bridge location, the string-body coupling is enforced using the body impulse-response or modal data, as measured at the instrument bridge. In the present paper, this computational approach is applied to a specific cello, which provided experimental wolf-behaviour data under several bowing conditions, as well as laboratory measurements of the bridge impulse responses on which the numerical simulations were based. Interesting aspects of the string-body dynamical responses are highlighted by numerical simulations and the corresponding sounds and animations produced. Finally, a qualitative (and, when possible, quantitative) comparison of the experimental and numerical results is presented.
Single-shot simulations of dynamic quantum many-body systems
NASA Astrophysics Data System (ADS)
Sakmann, Kaspar; Kasevich, Mark
2016-05-01
Single experimental shots of ultracold quantum gases sample the many-particle probability distribution. In a few cases such single shots could be successfully simulated from a given many-body wavefunction, but for realistic time-dependent many-body dynamics this has been difficult to achieve. Here, we show how single shots can be simulated from numerical solutions of the time-dependent many-body Schrödinger equation. Using this approach, we provide first-principle explanations for fluctuations in the collision of attractive Bose-Einstein condensates (BECs), for the appearance of randomly fluctuating vortices and for the centre-of-mass fluctuations of attractive BECs in a harmonic trap. We also show how such simulations provide full counting distributions and correlation functions of any order. Such calculations have not been previously possible and our method is broadly applicable to many-body systems whose phenomenology is driven by information beyond what is typically available in low-order correlation functions.
RANS Simulation of the Heave Response of a Two-Body Floating Point Wave Absorber: Preprint
Yu, Y.; Li, Y.
2011-03-01
A preliminary study on a two-body floating wave absorbers is presented in this paper. A Reynolds-Averaged Navier-Stokes computational method is applied for analyzing the hydrodynamic heave response of the absorber in operational wave conditions. The two-body floating wave absorber contains a float section and a submerged reaction section. For validation purposes, our model is first assumed to be locked. The two sections are forced to move together with each other. The locked single body model is used in a heave decay test, where the RANS result is validated with the experimental measurement. For the two-body floating point absorber simulation, the two sections are connected through a mass-spring-damper system, which is applied to simulate the power take-off mechanism under design wave conditions. Overall, the details of the flow around the absorber and its nonlinear interaction with waves are investigated, and the power absorption efficiency of the two-body floating wave absorber in waves with a constant value spring-damper system is examined.
Numerical simulation of a moving rigid body in a rarefied gas
NASA Astrophysics Data System (ADS)
Shrestha, Samir; Tiwari, Sudarshan; Klar, Axel; Hardt, Steffen
2015-07-01
In this paper we present a numerical scheme to simulate a moving rigid body with arbitrary shape suspended in a rarefied gas. The rarefied gas is simulated by solving the Boltzmann equation using a DSMC particle method. The motion of the rigid body is governed by the Newton-Euler equations, where the force and the torque on the rigid body are computed from the momentum transfer of the gas molecules colliding with the body. On the other hand, the motion of the rigid body influences the gas flow in its surroundings. We validate the numerical scheme by considering a moving piston problem in 1D and the Einstein relation for Brownian motion of the suspended particle in 3D. In the piston problem it is shown that the equilibrium position of the moving piston converges to the analytical solution for a wide range of Knudsen numbers. In the case of Brownian motion the translational as well as the rotational degrees of freedom are taken into account. In this case it is shown that the numerically computed translational and rotational diffusion coefficients converge to the theoretical values. Finally, the motion of an object of complex shape under the influence of a thermophoretic force is investigated.
NASA Astrophysics Data System (ADS)
Hoffmann, H.; Seiß, M.; Salo, H.; Spahn, F.
2014-04-01
Small moonlets in Saturn's rings induce propeller called structures into the surrounding ring material. Images of Saturn's rings, taken by the Cassini spacecraft near Saturn's equinox in 2009, show shadows cast by these propellers [1], offering the opportunity to study their vertical structure. We compare results from an extended hydrodynamical propeller model with results from local N-body box simulations of propeller structures. In the hydrodynamical model, maximal propeller heights are determined from the gravitational scattering of the ring particles by the moonlet. Afterwards the disturbed balance of viscous heating and collisional cooling is considered as main mechanism of the propeller height relaxation [2]. For the N-body box simulations we use the code by Salo [3], which was also applied in the propeller simulations of [4] and [5]. We find that the exponential height relaxation predicted by the hydrodynamical modelling is confirmed by N-body simulations of non-self gravitating ring particles. By projecting the propeller height evolution of the hydrodynamical model into observations of the shadows cast by the Earhart propeller, we determine the exponential cooling constant of the height relaxation. With this cooling constant we estimate collision frequencies of about 6 collisions per particle per orbit in the propeller gap region or about 11 collisions per particle per orbit in the propeller wake region of the Earhart propeller. The N-body simulations lead to maximal propeller heights between 60 to 70 percent of the Hill radius of the corresponding moonlet. Moonlet sizes estimated by this relation are in fair agreement with size estimates from radial propeller scalings [5, 6] for propeller structures with observed shadows.
Shi, Xiaobo; Li, Wei; Song, Jeungeun; Hossain, M Shamim; Mizanur Rahman, Sk Md; Alelaiwi, Abdulhameed
2016-10-01
With the development of IoT (Internet of Thing), big data analysis and cloud computing, traditional medical information system integrates with these new technologies. The establishment of cloud-based smart healthcare application gets more and more attention. In this paper, semi-physical simulation technology is applied to cloud-based smart healthcare system. The Body sensor network (BSN) of system transmit has two ways of data collection and transmission. The one is using practical BSN to collect data and transmitting it to the data center. The other is transmitting real medical data to practical data center by simulating BSN. In order to transmit real medical data to practical data center by simulating BSN under semi-physical simulation environment, this paper designs an OPNET packet structure, defines a gateway node model between simulating BSN and practical data center and builds a custom protocol stack. Moreover, this paper conducts a large amount of simulation on the real data transmission through simulation network connecting with practical network. The simulation result can provides a reference for parameter settings of fully practical network and reduces the cost of devices and personnel involved. PMID:27562482
Vlasov simulation of the interaction between the solar wind and a dielectric body
NASA Astrophysics Data System (ADS)
Umeda, Takayuki; Kimura, Tetsuya; Togano, Kentaro; Fukazawa, Keiichiro; Matsumoto, Yosuke; Miyoshi, Takahiro; Terada, Naoki; Nakamura, Takuma K. M.; Ogino, Tatsuki
2011-01-01
The global structure of wake field behind an unmagnetized object in the solar wind is studied by means of a 2.5-dimensional full-electromagnetic Vlasov simulation. The interaction of a plasma flow with an unmagnetized object is quite different from that with a magnetized object such as the Earth. Due to the absence of the global magnetic field, the unmagnetized object absorbs plasma particles that reach the surface, generating a plasma cavity called "wake" on the antisolar side of the object. For numerical simulations of electromagnetic structures around the wake, it is important to include the charging effect in global-scale simulations. The present study is one of the first attempts to study the formation of wake fields via a full-kinetic Vlasov simulation. It has been confirmed that the spatial structures of wake fields depend on the direction of interplanetary magnetic fields as well as the distance from the body.
Vlasov simulation of the interaction between the solar wind and a dielectric body
Umeda, Takayuki; Kimura, Tetsuya; Togano, Kentaro; Matsumoto, Yosuke; Ogino, Tatsuki; Fukazawa, Keiichiro; Miyoshi, Takahiro; Terada, Naoki; Nakamura, Takuma K. M.
2011-01-15
The global structure of wake field behind an unmagnetized object in the solar wind is studied by means of a 2.5-dimensional full-electromagnetic Vlasov simulation. The interaction of a plasma flow with an unmagnetized object is quite different from that with a magnetized object such as the Earth. Due to the absence of the global magnetic field, the unmagnetized object absorbs plasma particles that reach the surface, generating a plasma cavity called ''wake'' on the antisolar side of the object. For numerical simulations of electromagnetic structures around the wake, it is important to include the charging effect in global-scale simulations. The present study is one of the first attempts to study the formation of wake fields via a full-kinetic Vlasov simulation. It has been confirmed that the spatial structures of wake fields depend on the direction of interplanetary magnetic fields as well as the distance from the body.
Ingram, Joanne; Fisher, John; Ingham, Eileen
2011-01-01
Previously, we have described the development of an acellular porcine meniscal scaffold. The aims of this study were to determine the immunocompatibility of the scaffold and capacity for cellular attachment and infiltration to gain insight into its potential for meniscal repair and replacement. Porcine menisci were decellularized by exposing the tissue to freeze–thaw cycles, incubation in hypotonic tris buffer, 0.1% (w/v) sodium dodecyl sulfate in hypotonic buffer plus protease inhibitors, nucleases, hypertonic buffer followed by disinfection using 0.1% (v/v) peracetic, and final washing in phosphate-buffered saline. In vivo immunocompatibility was assessed after implantation of the acellular meniscal scaffold subcutaneously into galactosyltransferase knockout mice for 3 months in comparison to fresh and acellular tissue treated with α-galactosidase (negative control). The cellular infiltrates in the explants were assessed by histology and characterized using monoclonal antibodies against: CD3, CD4, CD34, F4/80, and C3c. Static culture was used to assess the potential of acellular porcine meniscal scaffold to support the attachment and infiltration of primary human dermal fibroblasts and primary porcine meniscal cells in vitro. The explants were surrounded by capsules that were more pronounced for the fresh meniscal tissue compared to the acellular tissues. Cellular infiltrates compromised mononuclear phagocytes, CD34-positive cells, and nonlabeled fibroblastic cells. T-lymphocytes were sparse in all explanted tissue types and there was no evidence of C3c deposition. The analysis revealed an absence of a specific immune response to all of the implanted tissues. Acellular porcine meniscus was shown to be capable of supporting the attachment and infiltration of primary human fibroblasts and primary porcine meniscal cells. In conclusion, acellular porcine meniscal tissue exhibits excellent immunocompatibility and potential for cellular regeneration in the longer
Turbine Powered Simulator Calibration and Testing for Hybrid Wing Body Powered Airframe Integration
NASA Technical Reports Server (NTRS)
Shea, Patrick R.; Flamm, Jeffrey D.; Long, Kurtis R.; James, Kevin D.; Tompkins, Daniel M.; Beyar, Michael D.
2016-01-01
Propulsion airframe integration testing on a 5.75% scale hybrid wing body model us- ing turbine powered simulators was completed at the National Full-Scale Aerodynamics Complex 40- by 80-foot test section. Four rear control surface con gurations including a no control surface de ection con guration were tested with the turbine powered simulator units to investigate how the jet exhaust in uenced the control surface performance as re- lated to the resultant forces and moments on the model. Compared to ow-through nacelle testing on the same hybrid wing body model, the control surface e ectiveness was found to increase with the turbine powered simulator units operating. This was true for pitching moment, lift, and drag although pitching moment was the parameter of greatest interest for this project. With the turbine powered simulator units operating, the model pitching moment was seen to increase when compared to the ow-through nacelle con guration indicating that the center elevon and vertical tail control authority increased with the jet exhaust from the turbine powered simulator units.
Large eddy simulation of flows after a bluff body: Coherent structures and mixing properties
NASA Astrophysics Data System (ADS)
Zhang, Pei; Han, Chao; Chen, Yiliang
2013-10-01
This paper performs large eddy simulations (LES) to investigate coherent structures in the flows after the Sydney bluff-body burner, a circular bluff body with an orifice at its center. The simulations are validated by comparison to existing experimental data. The Q function method is used to visualize the instantaneous vortex structures. Three kinds of structures are found, a cylindrical shell structure in the outer shear layer, a ring structure and some hairpin-like structures in the inner shear layer. An eduction scheme is employed to investigate the coherent structures in this flow. Some large streaks constituted by counter-rotating vortices are found in the outer shear layer and some well-organized strong structures are found in the inner shear layer. Finally, the influences of coherent structures on scalar mixing are studied and it is shown that scalar in the recirculation region is transported outward by coherent structures.
Clouds-in-clouds, clouds-in-cells physics for many-body plasma simulation
Birdsall, C.K.; Fuss, D.
1997-08-01
A clouds-interacting-with-clouds, clouds-in-cells method (CIC) is presented for many-body nonlinear plasma problems. Density and force are obtained by assuming that the particles have finite size, are tenuous, and may pass through one another; the particles are thus called clouds. They obey a Coulomb force ({approx} 1/r or 1/r{sup 3}) when separated and a linear force ({approx}r) when overlapping, allowing simple harmonic oscillations at small separation. CIC is contrasted with the zero-size particle and nearest-grid-point approach, ZSP-NGP. CIC appears to have substantially less unwanted noise than ZSP-NGP and should be more useful in simulating dense plasmas. Initial runs have been encouraging. The methods may find use in other many-body simulations, such as with stars, or with particles in phase space. 9 refs., 13 figs.
Flight-Simulated Launch-Pad-Abort-to-Landing Maneuvers for a Lifting Body
NASA Technical Reports Server (NTRS)
Jackson, E. Bruce; Rivers, Robert A.
1998-01-01
The results of an in-flight investigation of the feasibility of conducting a successful landing following a launch-pad abort of a vertically-launched lifting body are presented. The study attempted to duplicate the abort-to-land-ing trajectory from the point of apogee through final flare and included the steep glide and a required high-speed, low-altitude turn to the runway heading. The steep glide was flown by reference to ground-provided guidance. The low-altitude turn was flown visually with a reduced field- of-view duplicating that of the simulated lifting body. Results from the in-flight experiment are shown to agree with ground-based simulation results; however, these tests should not be regarded as a definitive due to performance and control law dissimilarities between the two vehicles.
Simulation of the many-body dynamical quantum Hall effect in an optical lattice
NASA Astrophysics Data System (ADS)
Zhang, Dan-Wei; Yang, Xu-Chen
2016-05-01
We propose an experimental scheme to simulate the many-body dynamical quantum Hall effect with ultra-cold bosonic atoms in a one-dimensional optical lattice. We first show that the required model Hamiltonian of a spin-1/2 Heisenberg chain with an effective magnetic field and tunable parameters can be realized in this system. For dynamical response to ramping the external fields, the quantized plateaus emerge in the Berry curvature of the interacting atomic spin chain as a function of the effective spin-exchange interaction. The quantization of this response in the parameter space with the interaction-induced topological transition characterizes the many-body dynamical quantum Hall effect. Furthermore, we demonstrate that this phenomenon can be observed in practical cold atom experiments with numerical simulations.
Wavelet-based adaptive numerical simulation of unsteady 3D flow around a bluff body
NASA Astrophysics Data System (ADS)
de Stefano, Giuliano; Vasilyev, Oleg
2012-11-01
The unsteady three-dimensional flow past a two-dimensional bluff body is numerically simulated using a wavelet-based method. The body is modeled by exploiting the Brinkman volume-penalization method, which results in modifying the governing equations with the addition of an appropriate forcing term inside the spatial region occupied by the obstacle. The volume-penalized incompressible Navier-Stokes equations are numerically solved by means of the adaptive wavelet collocation method, where the non-uniform spatial grid is dynamically adapted to the flow evolution. The combined approach is successfully applied to the simulation of vortex shedding flow behind a stationary prism with square cross-section. The computation is conducted at transitional Reynolds numbers, where fundamental unstable three-dimensional vortical structures exist, by well-predicting the unsteady forces arising from fluid-structure interaction.
Comparative study of the biodegradability of porous silicon films in simulated body fluid.
Peckham, J; Andrews, G T
2015-01-01
The biodegradability of oxidized microporous, mesoporous and macroporous silicon films in a simulated body fluid with ion concentrations similar to those found in human blood plasma were studied using gravimetry. Film dissolution rates were determined by periodically weighing the samples after removal from the fluid. The dissolution rates for microporous silicon were found to be higher than those for mesoporous silicon of comparable porosity. The dissolution rate of macroporous silicon was much lower than that for either microporous or mesoporous silicon. This is attributed to the fact that its specific surface area is much lower than that of microporous and mesoporous silicon. Using an equation adapted from [Surf. Sci. Lett. 306 (1994), L550-L554], the dissolution rate of porous silicon in simulated body fluid can be estimated if the film thickness and specific surface area are known. PMID:25585985
Cartesian Off-Body Grid Adaption for Viscous Time- Accurate Flow Simulation
NASA Technical Reports Server (NTRS)
Buning, Pieter G.; Pulliam, Thomas H.
2011-01-01
An improved solution adaption capability has been implemented in the OVERFLOW overset grid CFD code. Building on the Cartesian off-body approach inherent in OVERFLOW and the original adaptive refinement method developed by Meakin, the new scheme provides for automated creation of multiple levels of finer Cartesian grids. Refinement can be based on the undivided second-difference of the flow solution variables, or on a specific flow quantity such as vorticity. Coupled with load-balancing and an inmemory solution interpolation procedure, the adaption process provides very good performance for time-accurate simulations on parallel compute platforms. A method of using refined, thin body-fitted grids combined with adaption in the off-body grids is presented, which maximizes the part of the domain subject to adaption. Two- and three-dimensional examples are used to illustrate the effectiveness and performance of the adaption scheme.
A Numerical Simulation of a Fishlike Body's Self-propelled C-start
NASA Astrophysics Data System (ADS)
Liu, G.; Yu, Y. L.; Tong, B. G.
2011-09-01
This paper presents a numerical method to deal with a two-dimensional deformable fishlike body's large deformation self-propelled swimming. Overset grids are employed to discretize the flow domain around the large deforming body which is simulated by a foil. The kinematics, energetics and the flow structures of a typical C-start are predicted by a coupling solution of the two-dimensional incompressible fluid dynamics and the deforming body dynamics. As a typical practice, the foil performs a C-start-like motion in stationary water based on the prescribed deforming mode. It is found that the locomotion of the foil is similar to a real fish's C-start and the hydrodynamic efficiency of this C-start model is about 29% which is close to the value calculated by the previous theoretical estimation. Particularly, a structure of three significant concentrated vortices is discovered in the wake.
Attitude dynamics simulation subroutines for systems of hinge-connected rigid bodies
NASA Technical Reports Server (NTRS)
Fleischer, G. E.; Likins, P. W.
1974-01-01
Several computer subroutines are designed to provide the solution to minimum-dimension sets of discrete-coordinate equations of motion for systems consisting of an arbitrary number of hinge-connected rigid bodies assembled in a tree topology. In particular, these routines may be applied to: (1) the case of completely unrestricted hinge rotations, (2) the totally linearized case (all system rotations are small), and (3) the mixed, or partially linearized, case. The use of the programs in each case is demonstrated using a five-body spacecraft and attitude control system configuration. The ability of the subroutines to accommodate prescribed motions of system bodies is also demonstrated. Complete listings and user instructions are included for these routines (written in FORTRAN V) which are intended as multi- and general-purpose tools in the simulation of spacecraft and other complex electromechanical systems.
Positive Tensor Network Approach for Simulating Open Quantum Many-Body Systems
NASA Astrophysics Data System (ADS)
Werner, A. H.; Jaschke, D.; Silvi, P.; Kliesch, M.; Calarco, T.; Eisert, J.; Montangero, S.
2016-06-01
Open quantum many-body systems play an important role in quantum optics and condensed matter physics, and capture phenomena like transport, the interplay between Hamiltonian and incoherent dynamics, and topological order generated by dissipation. We introduce a versatile and practical method to numerically simulate one-dimensional open quantum many-body dynamics using tensor networks. It is based on representing mixed quantum states in a locally purified form, which guarantees that positivity is preserved at all times. Moreover, the approximation error is controlled with respect to the trace norm. Hence, this scheme overcomes various obstacles of the known numerical open-system evolution schemes. To exemplify the functioning of the approach, we study both stationary states and transient dissipative behavior, for various open quantum systems ranging from few to many bodies.
Positive Tensor Network Approach for Simulating Open Quantum Many-Body Systems.
Werner, A H; Jaschke, D; Silvi, P; Kliesch, M; Calarco, T; Eisert, J; Montangero, S
2016-06-10
Open quantum many-body systems play an important role in quantum optics and condensed matter physics, and capture phenomena like transport, the interplay between Hamiltonian and incoherent dynamics, and topological order generated by dissipation. We introduce a versatile and practical method to numerically simulate one-dimensional open quantum many-body dynamics using tensor networks. It is based on representing mixed quantum states in a locally purified form, which guarantees that positivity is preserved at all times. Moreover, the approximation error is controlled with respect to the trace norm. Hence, this scheme overcomes various obstacles of the known numerical open-system evolution schemes. To exemplify the functioning of the approach, we study both stationary states and transient dissipative behavior, for various open quantum systems ranging from few to many bodies. PMID:27341253
Assessment of simulated high-dose partial-body irradiation by PCC-R assay
Romero, Ivonne; García, Omar; Lamadrid, Ana I.; Gregoire, Eric; González, Jorge E.; Morales, Wilfredo; Martin, Cécile; Barquinero, Joan-Francesc; Voisin, Philippe
2013-01-01
The estimation of the dose and the irradiated fraction of the body is important information in the primary medical response in case of a radiological accident. The PCC-R assay has been developed for high-dose estimations, but little attention has been given to its applicability for partial-body irradiations. In the present work we estimated the doses and the percentage of the irradiated fraction in simulated partial-body radiation exposures at high doses using the PCC-R assay. Peripheral whole blood of three healthy donors was exposed to doses from 0–20 Gy, with 60Co gamma radiation. To simulate partial body irradiations, irradiated and non-irradiated blood was mixed to obtain proportions of irradiated blood from 10–90%. Lymphocyte cultures were treated with Colcemid and Calyculin-A before harvest. Conventional and triage scores were performed for each dose, proportion of irradiated blood and donor. The Papworth's u test was used to evaluate the PCC-R distribution per cell. A dose-response relationship was fitted according to the maximum likelihood method using the frequencies of PCC-R obtained from 100% irradiated blood. The dose to the partially irradiated blood was estimated using the Contaminated Poisson method. A new D0 value of 10.9 Gy was calculated and used to estimate the initial fraction of irradiated cells. The results presented here indicate that by PCC-R it is possible to distinguish between simulated partial- and whole-body irradiations by the u-test, and to accurately estimate the dose from 10–20 Gy, and the initial fraction of irradiated cells in the interval from 10–90%. PMID:23596200
NOVel Adaptive softening for collisionless N-body simulations: eliminating spurious haloes
NASA Astrophysics Data System (ADS)
Hobbs, Alexander; Read, Justin I.; Agertz, Oscar; Iannuzzi, Francesca; Power, Chris
2016-05-01
We describe a NOVel form of Adaptive softening (NOVA) for collisionless N-body simulations, implemented in the RAMSES adaptive mesh refinement code. In RAMSES - that we refer to as a `standard N-body method' - cells are only split if they contain more than eight particles (a mass refinement criterion). Here, we introduce an additional criterion that the particle distribution within each cell be sufficiently isotropic, as measured by the ratio of the maximum to minimum eigenvalues of its moment of inertia tensor: η = λmax/λmin. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening ɛ is always of order twice the largest interparticle spacing in a cell. This more conservative force softening criterion is designed to minimize spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NOVA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that - unlike the standard N-body method - NOVA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NOVA converges significantly more rapidly than standard N-body, producing little or no spurious haloes on small scales. We will use NOVA in forthcoming papers to study the issue of halo formation below the free-streaming mass Mfs; filament stability; and to obtain new constraints on the temperature of dark matter.
Unsteady aerodynamic simulation of multiple bodies in relative motion: A prototype method
NASA Technical Reports Server (NTRS)
Meakin, Robert L.
1989-01-01
A prototype method for time-accurate simulation of multiple aerodynamic bodies in relative motion is presented. The method is general and features unsteady chimera domain decomposition techniques and an implicit approximately factored finite-difference procedure to solve the time-dependent thin-layer Navier-Stokes equations. The method is applied to a set of two- and three- dimensional test problems to establish spatial and temporal accuracy, quantify computational efficiency, and begin to test overall code robustness.
Molina, John J; Yamamoto, Ryoichi
2013-12-21
An improved formulation of the "Smoothed Profile" method is introduced to perform direct numerical simulations of arbitrary rigid body dispersions in a Newtonian host solvent. Previous implementations of the method were restricted to spherical particles, severely limiting the types of systems that could be studied. The validity of the method is carefully examined by computing the friction/mobility tensors for a wide variety of geometries and comparing them to reference values obtained from accurate solutions to the Stokes-Equation. PMID:24359350
Three-dimensional finite element simulations of vertebral body thermal treatment (Invited Paper)
NASA Astrophysics Data System (ADS)
Ryan, Thomas P.; Patel, Samit J.; Morris, Ronit; Hoopes, P. J.; Bergeron, Jeffrey A.; Mahajan, Roop
2005-04-01
Lower back pain affects a large group of people worldwide and when in its early stages, has no viable interventional treatment. In order to avoid the eventuality of an invasive surgical procedure, which is further down the Care Pathway, an interventional treatment that is minimally invasive and arrests the patient's pain would be of tremendous clinical benefit. There is a hypothesis that if the basivertebral nerve in the vertebral body is defunctionalized, lower back pain may be lessened. To further investigate creating a means to provide localized thermal therapy, bench and animal studies were planned, but to help select the applicator configuration and placement, numerical modeling studies were undertaken. A 3D finite element model was utilized to predict the electric field pattern and power deposition pattern of radiofrequency (RF) based electrodes. Three types of tissues were modeled: 1) porcine (ex-vivo), ovine (in-vivo preclinical), and 3) human (ex-vivo, in-vivo). Two types of RF devices were simulated: 1) a pair of converging, hollow electrodes, and 2) an in-line pair of spaced-apart electrodes. Temperature distributions over time were plotted using the electric field results and the bioheat equation. Since the thermal and electrical properties of the vertebral bodies of porcine, ovine, and human tissue were not available, measurements were undertaken to capture these data to input into the model. The measurements of electrical and thermal properties of cancellous and cortical vertebral body were made over a range of temperatures. The simulation temperature results agreed with live animal and human cadaver studies. In addition, the lesion shapes predicted in the simulations matched CT and MRI studies done during the chronic ovine study, as well as histology results. In conclusion, the simulations aided in shaping and sizing the RF electrodes, as well as positioning them in the vertebral body structures to assure that the basivertebral nerve was ablated, but
WINE-1: Special-Purpose Computer forN-Body Simulations with a Periodic Boundary Condition
NASA Astrophysics Data System (ADS)
Fukushige, Toshiyuki; Makino, Junichiro; Ito, Tomoyoshi; Okumura, Sachiko K.; Ebisuzaki, Toshikazu; Sugimoto, Daiichiro
1993-06-01
We have developed WINE-1 (Wave space INtegrator for Ewald method), a special-purpose computer for N-body simulations with a periodic boundary condition. In N-body simulations with a periodic boundary condition such as cosmological N-body simulations, we use the Ewald method to calculate the gravitational interaction. With the Ewald method, we can calculate the interaction more accurately than a calculation with other methods, such as the PM method, the P(3) M method, or the tree algorithm. In the Ewald method, the total force exerted on a particle is divided into contributions from real space and wave-number space so that the infinite sum can converge exponentially in both spaces. WINE is a special-purpose computer used to calculate the interaction in wave-number space. WINE is connected to a host computer via the VME bus. We have developed the first machine, WINE-1. It is made of one board having a size of 38 cm by 40 cm, on which 31 LSI chips and 46 IC chips are wire-wrapped. The peak speed of WINE-1 is equivalent to 480 Mflops. The summation in real space is calculated using a GRAPE system, another special-purpose computer for the direct calculation of the interparticle force. For example, we can perform a cosmological N-body simulation for N=80,000 (500 steps) within a week if we use GRAPE-2A for the summation in real space and WINE-1 for that in wave-number space.
Comparative biological activities of acellular pertussis vaccines produced by Kitasato.
Watanabe, M; Izumiya, K; Sato, T; Yoshino, K; Nakagawa, N; Ohoishi, M; Hoshino, M
1991-04-01
The quality of 14 lots of acellular pertussis-diphtheria-tetanus (AC-PDT) vaccines manufactured by the Kitasato Institute during the period 1987-1990 were investigated. The geometric means of HSU, LPU, and BWDU were 0.078, 0.257, and 7.33 per ml respectively. The potency was higher than 14 IU per ml. These results indicated the consistency of the Kitasato AC-PDT vaccines. The antibody response to the AC-PDT vaccines was measured in primary and secondary vaccinated mice by ELISA. IgG antibody response to FHA and PT was obtained in all immunized mice (P less than 0.001) after the primary injection. In contrast, IgG antibody response to fimbriae 2 showed a significant titer rise (P less than 0.001) after the booster injection. The results indicated that the Kitasato AC-P vaccines consisted of protein, PT and FHA as the major antigens, and a little agglutinogen as the minor antigen. PMID:1798236
Egorova, N B; Efremova, V N; Kurbatova, E A; Gruber, I M
2008-01-01
Results of experimental, clinical and immunological effects of acellular dry staphylococcal vaccine "Staphylovac" developed in Mechnikov Research Institute of Vaccines and Sera are presented. Original mildly virulent strains of Staphylococcus aureus having high immunogenicity, and intra- and interspecies protective activity against different representatives of opportunistic microflora were used for construction of the preparation. Low-toxicity and weak anapylactogenicity of the vaccine were established. In experiments on mice, guinea pigs and rabbits significant protective, antigenic and immunomodulate activity of the preparation was revealed with low sensitization of animals. Clinical trials performed in different centers showed that inclusion of vaccinotherapy in complex treatment of chronic staphylococcal infections (chronic pyodermia, lung abscess etc.) resulted in prolonged pathologic locus, decrease of number and severity of exacerbations, prolongation of remission, and complete recovery in significant number of patients. Activation of innate and adaptive immunity was revealed in the same patients. It was shown on the large group of athletes that administration of the vaccine by aerosol route prevents disruption of immunologic adaptation occurring due to excess physical activity and stress situations during competitions. PMID:19186558
A Complication Analysis of 2 Acellular Dermal Matrices in Prosthetic-based Breast Reconstruction
Page, Eugenia K.; Hart, Alexandra; Rudderman, Randall; Carlson, Grant W.; Losken, Albert
2016-01-01
Background: Acellular dermal matrices (ADM) are now routine in postmastectomy prosthetic-based breast reconstruction. The goal of the current study was to compare the complications of 2 ADM products—AlloDerm and Cortiva. Methods: A retrospective analysis of prosthetic-based breast reconstruction in Atlanta, Ga., over 5 years. Inclusion criteria were the use of the ADM types (AlloDerm or Cortiva) and use of a tissue expander or implant. Statistical analysis compared group demographics, risk factors, and early complications. Results: Of the 298 breast reconstructions, 174 (58.4%) used AlloDerm and 124 (41.6%) used Cortiva. There was no difference in overall complication frequency (16 AlloDerm and 18 Cortiva; P = 0.195). Within specific categories, there was a difference in mastectomy skin flap necrosis, but, based on further regression analysis, this was attributable to differences in body mass index (P = 0.036). Furthermore, there were no differences in the rates of infection (6 AlloDerm and 5 Cortiva; P = 1.0), seroma/hematoma (9 AlloDerm and 7 Cortiva; P = 1.0), or drain duration (13.2 day AlloDerm and 14.2 day Cortiva, P = 0.2). By using a general estimating equation for binomial logistical regression, it was found that only current tobacco use (P = 0.033) was a significant predictor for a complication. Trending predictors were body mass index (P = 0.074) and age (P = 0.093). The type of matrix was not a significant predictor for any of the recorded complication (P = 0.160). Conclusions: Although AlloDerm is well established, we have shown that Cortiva has an equivalent complication frequency. Future work will focus on long-term outcome measures and histological evaluation of vascularization and integration. PMID:27536479
Modeling, simulation and optimization approaches for design of lightweight car body structures
NASA Astrophysics Data System (ADS)
Kiani, Morteza
Simulation-based design optimization and finite element method are used in this research to investigate weight reduction of car body structures made of metallic and composite materials under different design criteria. Besides crashworthiness in full frontal, offset frontal, and side impact scenarios, vibration frequencies, static stiffness, and joint rigidity are also considered. Energy absorption at the component level is used to study the effectiveness of carbon fiber reinforced polymer (CFRP) composite material with consideration of different failure criteria. A global-local design strategy is introduced and applied to multi-objective optimization of car body structures with CFRP components. Multiple example problems involving the analysis of full-vehicle crash and body-in-white models are used to examine the effect of material substitution and the choice of design criteria on weight reduction. The results of this study show that car body structures that are optimized for crashworthiness alone may not meet the vibration criterion. Moreover, optimized car body structures with CFRP components can be lighter with superior crashworthiness than the baseline and optimized metallic structures.
General relativistic N-body simulations in the weak field limit
NASA Astrophysics Data System (ADS)
Adamek, Julian; Daverio, David; Durrer, Ruth; Kunz, Martin
2013-11-01
We develop a formalism for general relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives to second order, and second derivatives to all orders, thus taking into account the most important nonlinear effects of Einstein gravity. It is also expected that any significant “backreaction” should appear at this order. We show that the simulation scheme is feasible in practice by implementing it for a plane-symmetric situation and running two test cases, one with only cold dark matter, and one which also includes a cosmological constant. For these plane-symmetric situations, the deviations from the usual Newtonian N-body simulations remain small and, apart from a nontrivial correction to the background, can be accurately estimated within the Newtonian framework. The correction to the background scale factor, which is a genuine backreaction effect, can be robustly obtained with our algorithm. Our numerical approach is also naturally suited for the inclusion of extra relativistic fields and thus for dark energy or modified gravity simulations.
Acceleration of hybrid MPI parallel NBODY6++ for large N-body globular cluster simulations
NASA Astrophysics Data System (ADS)
Wang, Long; Spurzem, Rainer; Aarseth, Sverre; Nitadori, Keigo; Berczik, Peter; Kouwenhoven, M. B. N.; Naab, Thorsten
2016-02-01
Previous research on globular clusters (GCs) dynamics is mostly based on semi-analytic, Fokker-Planck, Monte-Carlo methods and on direct N-body (NB) simulations. These works have great advantages but also limits since GCs are massive and compact and close encounters and binaries play very important roles in their dynamics. The former three methods make approximations and assumptions, while expensive computing time and number of stars limit the latter method. The current largest direct NB simulation has ~ 500k stars (Heggie 2014). Here, we accelerate the direct NB code NBODY6++ (which extends NBODY6 to supercomputers by using MPI) with new parallel computing technologies (GPU, OpenMP + SSE/AVX). Our aim is to handle large N (up to 106) direct NB simulations to obtain better understanding of the dynamical evolution of GCs.
Real-time simulation model of the HL-20 lifting body
NASA Technical Reports Server (NTRS)
Jackson, E. Bruce; Cruz, Christopher I.; Ragsdale, W. A.
1992-01-01
A proposed manned spacecraft design, designated the HL-20, has been under investigation at Langley Research Center. Included in that investigation are flight control design and flying qualities studies utilizing a man-in-the-loop real-time simulator. This report documents the current real-time simulation model of the HL-20 lifting body vehicle, known as version 2.0, presently in use at NASA Langley Research Center. Included are data on vehicle aerodynamics, inertias, geometries, guidance and control laws, and cockpit displays and controllers. In addition, trim case and dynamic check case data is provided. The intent of this document is to provide the reader with sufficient information to develop and validate an equivalent simulation of the HL-20 for use in real-time or analytical studies.
Hamai, Ryo; Shirosaki, Yuki; Miyazaki, Toshiki
2016-10-01
Phosphate groups on materials surfaces are known to contribute to apatite formation upon exposure of the materials in simulated body fluid and improved affinity of the materials for osteoblast-like cells. Typically, polymers containing phosphate groups are organic matrices consisting of apatite-polymer composites prepared by biomimetic process using simulated body fluid. Ca(2+) incorporation into the polymer accelerates apatite formation in simulated body fluid owing because of increase in the supersaturation degree, with respect to apatite in simulated body fluid, owing to Ca(2+) release from the polymer. However, the effects of phosphate content on the Ca(2+) release and apatite-forming abilities of copolymers in simulated body fluid are rather elusive. In this study, a phosphate-containing copolymer prepared from vinylphosphonic acid, 2-hydroxyethyl methacrylate, and triethylene glycol dimethacrylate was examined. The release of Ca(2+) in Tris-NaCl buffer and simulated body fluid increased as the additive amount of vinylphosphonic acid increased. However, apatite formation was suppressed as the phosphate groups content increased despite the enhanced release of Ca(2+) from the polymer. This phenomenon was reflected by changes in the surface zeta potential. Thus, it was concluded that the apatite-forming ability of vinylphosphonic acid-2-hydroxyethyl methacrylate-triethylene glycol dimethacrylate copolymer treated with CaCl2 solution was governed by surface state rather than Ca(2+) release in simulated body fluid. PMID:27585911
Large-scale hydrodynamical and N-body simulations of viscous overstability in Saturn's rings
NASA Astrophysics Data System (ADS)
Latter, H.; Rein, H.; Ogilvie, G.
2012-09-01
We aim to understand axisymmetric structure formation in Saturn's A and B-rings on scales of 100 m to several kms through nonlinear hydrodynamical and N-body simulations of the viscous overstability. The viscous overstability is a robust generator of structure on short scales, as witnessed by previous hydrodynamical and N-body simulations (Schmidt and Tscharnuter 1999, Salo et al. 2001), and is hence the most likely candidate responsible for recently observed periodic microstructure (Colwell et al. 2007, Thomson et al. 2007). It is also possible that during its nonlinear saturation the instability gives rise to axisymmetric patterns on slightly longer scales that may correspond to observed irregular structure on 1-10 km (Porco et al. 2005, Latter and Ogilvie 2009, 2010). Our hydrodynamical and N-body simulations are undertaken in local Cartesian domains that can extend over 10 km in radius and can be evolved forward in time for more than 1000 orbits. These hence provide the scope to fully describe the nonlinear saturation of the overstability and to mainfest the full range of its dynamics. Self-gravity is omitted at this stage, but will be included in future work. Nonlinear wavetrains dominate all the simulations, and we associate them with the observed periodic microstructure. The preferred lengthscale of these waves (~200 m) is set by secondary modulational instabilities. These wavetrains undergo small chaotic fluctuations in their phases and amplitudes, and may be punctuated by more formidable `wave-defects', that are distributed on longer scales (~ 1-5 km). It is possible that the defects are connected to the irregular larger-scale variations observed. We also speculate on the azimuthal extent of the waves and the influence of self-gravity wakes on their dynamics.
Drift trajectories of a floating human body simulated in a hydraulic model of Puget Sound.
Ebbesmeyer, C C; Haglund, W D
1994-01-01
After a young man jumped off a 221-foot (67 meters) high bridge, the drift of the body that beached 20 miles (32 km) away at Alki Point in Seattle, Washington was simulated with a hydraulic model. Simulations for the appropriate time period were performed using a small floating bead to represent the body in the hydraulic model at the University of Washington. Bead movements were videotaped and transferred to Computer Aided Drafting (AutoCAD) charts on a personal computer. Because of strong tidal currents in the narrow passage under the bridge (The Narrows near Tacoma, WA), small changes in the time of the jump (+/- 30 minutes) made large differences in the distance the body traveled (30 miles; 48 km). Hydraulic and other types of oceanographic models may be located by contacting technical experts known as physical oceanographers at local universities, and can be utilized to demonstrate trajectories of floating objects and the time required to arrive at selected locations. Potential applications for forensic death investigators include: to be able to set geographic and time limits for searches; determine potential origin of remains found floating or beached; and confirm and correlate information regarding entry into the water and sightings of remains. PMID:8113703
Numerical simulations of impacts involving porous bodies. II. Comparison with laboratory experiments
NASA Astrophysics Data System (ADS)
Jutzi, Martin; Michel, Patrick; Hiraoka, Kensuke; Nakamura, Akiko M.; Benz, Willy
2009-06-01
In this paper, we compare the outcome of high-velocity impact experiments on porous targets, composed of pumice, with the results of simulations by a 3D SPH hydrocode in which a porosity model has been implemented. The different populations of small bodies of our Solar System are believed to be composed, at least partially, of objects with a high degree of porosity. To describe the fragmentation of such porous objects, a different model is needed than that used for non-porous bodies. In the case of porous bodies, the impact process is not only driven by the presence of cracks which propagate when a stress threshold is reached, it is also influenced by the crushing of pores and compaction. Such processes can greatly affect the whole body's response to an impact. Therefore, another physical model is necessary to improve our understanding of the collisional process involving porous bodies. Such a model has been developed recently and introduced successfully in a 3D SPH hydrocode [Jutzi, M., Benz, W., Michel, P., 2008. Icarus 198, 242-255]. Basic tests have been performed which already showed that it is implemented in a consistent way and that theoretical solutions are well reproduced. However, its full validation requires that it is also capable of reproducing the results of real laboratory impact experiments. Here we present simulations of laboratory experiments on pumice targets for which several of the main material properties have been measured. We show that using the measured material properties and keeping the remaining free parameters fixed, our numerical model is able to reproduce the outcome of these experiments carried out under different impact conditions. This first complete validation of our model, which will be tested for other porous materials in the future, allows us to start addressing problems at larger scale related to small bodies of our Solar System, such as collisions in the Kuiper Belt or the formation of a family by the disruption of a porous
Three-Body Abrasion Testing Using Lunar Dust Simulants to Evaluate Surface System Materials
NASA Technical Reports Server (NTRS)
Kobrick, Ryan L.; Budinski, Kenneth G.; Street, Kenneth W., Jr.; Klaus, David M.
2010-01-01
Numerous unexpected operational issues relating to the abrasive nature of lunar dust, such as scratched visors and spacesuit pressure seal leaks, were encountered during the Apollo missions. To avoid reoccurrence of these unexpected detrimental equipment problems on future missions to the Moon, a series of two- and three-body abrasion tests were developed and conducted in order to begin rigorously characterizing the effect of lunar dust abrasiveness on candidate surface system materials. Two-body scratch tests were initially performed to examine fundamental interactions of a single particle on a flat surface. These simple and robust tests were used to establish standardized measurement techniques for quantifying controlled volumetric wear. Subsequent efforts described in the paper involved three-body abrasion testing designed to be more representative of actual lunar interactions. For these tests, a new tribotester was developed to expose samples to a variety of industrial abrasives and lunar simulants. The work discussed in this paper describes the three-body hardware setup consisting of a rotating rubber wheel that applies a load on a specimen as a loose abrasive is fed into the system. The test methodology is based on ASTM International (ASTM) B611, except it does not mix water with the abrasive. All tests were run under identical conditions. Abraded material specimens included poly(methyl methacrylate) (PMMA), hardened 1045 steel, 6061-T6 aluminum (Al) and 1018 steel. Abrasives included lunar mare simulant JSC- 1A-F (nominal size distribution), sieved JSC-1A-F (<25 m particle diameter), lunar highland simulant NU-LHT-2M, alumina (average diameter of 50 m used per ASTM G76), and silica (50/70 mesh used per ASTM G65). The measured mass loss from each specimen was converted using standard densities to determine total wear volume in cm3. Abrasion was dominated by the alumina and the simulants were only similar to the silica (i.e., sand) on the softer materials of
Simulation, Modeling and Feedback Control of the flow around a Square-Back Bluff Body
NASA Astrophysics Data System (ADS)
Dalla Longa, Laurent; Morgans, Aimee; Imperial College London-Flow Control Team
2015-11-01
Because of capacity, aesthetic and comfort requirements, most road vehicles are not streamlined but blunt bluff bodies. The flow exhibits a large wake recirculation area leading to high pressure drag, which at highway speeds, represents the main source of energy loss. In this work, Large Eddy Simulations of the flow past a square-back bluff body with interacting shear layers are performed with the aim of reducing aerodynamic drag. A linear feedback control strategy is applied to increase the back face pressure and therefore obtain drag reduction. Synthetic jets located along the perimeter of the back face are used for actuation while body mounted sensors record the base pressure. System identification, via harmonic actuator forcing, is used to characterize the flow response to actuation, which is assumed to be dynamically linear. Based on the identified frequency response, a feedback controller is designed in the frequency domain which aims to either attenuate or amplify base pressure fluctuations by shaping of the sensitivity transfer function. This is first done for a D-shaped body. Current work extends this strategy to a simplified lorry geometry on which experiments were carried out recently.
Shi, Xiangnan; Cao, Libo; Reed, Matthew P; Rupp, Jonathan D; Hu, Jingwen
2015-01-01
The objective of this study is to investigate the effects of obesity on occupant responses in frontal crashes using whole-body human finite element (FE) models representing occupants with different obesity levels. In this study, the geometry of THUMS 4 midsize male model was varied using mesh morphing techniques with target geometries defined by statistical models of external body contour and exterior ribcage geometry. Models with different body mass indices (BMIs) were calibrated against cadaver test data under high-speed abdomen loading and frontal crash conditions. A parametric analysis was performed to investigate the effects of BMI on occupant injuries in frontal crashes based on the Taguchi method while controlling for several vehicle design parameters. Simulations of obese occupants predicted significantly higher risks of injuries to the thorax and lower extremities in frontal crashes compared with non-obese occupants, which is consistent with previous field data analyses. These higher injury risks are mainly due to the increased body mass and relatively poor belt fit caused by soft tissues for obese occupants. This study demonstrated the feasibility of using a parametric human FE model to investigate the obesity effects on occupant responses in frontal crashes. PMID:24666169
NASA Astrophysics Data System (ADS)
Thirouin, Audrey; Bagati, A. C.; Ortiz, J.; Duffard, R.; Benavidez, P.; Richardson, D.
2010-10-01
Haumea is one of the most studied and probably one of the most interesting Trans-Neptunian Objects (TNOs) and a fast spinning dwarf planet (Rabinowitz et al., 2006; Thirouin et al. 2010) that has at least two satellites and whose orbital elements are related to a group/family of bodies. A catastrophic collision at high relative velocity (around 3 km/s) between two bodies in the 1000-1500 km size range and a mass ratio 0.2 has been suggested for the formation of the 'family' and the presence of satellites (Brown et al. 2007). Leinhardt et al. (2010) proposed another mechanism in which two 1300 km size bodies undergo a grazing collision with reacculumation of part of the mass and dispersion of the rest, partly into satellites. The likelihood of both scenarios is small when tested against collisional and dynamical evolution model predictions of collisional probabilities. Furthermore, these mechanisms have trouble in explaining the velocity dispersion of the family members and the fast spin of Haumea. Schlichting and Sari (2009) proposed that a former collision at low speed formed a proto-satellite that should lately undergo a final impact finally leading to the current observed system. In this work, we propose alternative mechanisms for the formation of the Haumea system ant test them by means of N-body numerical simulations (PKDGRAV code. Richardson, 1994).
2010-01-01
In vitro bioaccessibility (IVBA) studies were carried out on samples of mercury (Hg) mine-waste calcine (roasted Hg ore) by leaching with simulated human body fluids. The objective was to estimate potential human exposure to Hg due to inhalation of airborne calcine particulates and hand-to-mouth ingestion of Hg-bearing calcines. Mine waste calcines collected from Hg mines at Almadén, Spain, and Terlingua, Texas, contain Hg sulfide, elemental Hg, and soluble Hg compounds, which constitute primary ore or compounds formed during Hg retorting. Elevated leachate Hg concentrations were found during calcine leaching using a simulated gastric fluid (as much as 6200 μg of Hg leached/g sample). Elevated Hg concentrations were also found in calcine leachates using a simulated lung fluid (as much as 9200 μg of Hg leached/g), serum-based fluid (as much as 1600 μg of Hg leached/g), and water of pH 5 (as much as 880 μg of Hg leached/g). The leaching capacity of Hg is controlled by calcine mineralogy; thus, calcines containing soluble Hg compounds contain higher leachate Hg concentrations. Results indicate that ingestion or inhalation of Hg mine-waste calcine may lead to increased Hg concentrations in the human body, especially through the ingestion pathway. PMID:20491469
Hierarchical tree algorithm for collisional N-body simulations on GRAPE
NASA Astrophysics Data System (ADS)
Fukushige, Toshiyuki; Kawai, Atsushi
2016-06-01
We present an implementation of the hierarchical tree algorithm on the individual timestep algorithm (the Hermite scheme) for collisional N-body simulations, running on the GRAPE-9 system, a special-purpose hardware accelerator for gravitational many-body simulations. Such a combination of the tree algorithm and the individual timestep algorithm was not easy on the previous GRAPE system mainly because its memory addressing scheme was limited only to sequential access to a full set of particle data. The present GRAPE-9 system has an indirect memory addressing unit and a particle memory large enough to store all the particle data and also the tree node data. The indirect memory addressing unit stores interaction lists for the tree algorithm, which is constructed on the host computer, and, according to the interaction lists, force pipelines calculate only the interactions necessary. In our implementation, the interaction calculations are significantly reduced compared to direct N2 summation in the original Hermite scheme. For example, we can achieve about a factor 30 of speedup (equivalent to about 17 teraflops) against the Hermite scheme for a simulation of an N = 106 system, using hardware of a peak speed of 0.6 teraflops for the Hermite scheme.
Hierarchical tree algorithm for collisional N-body simulations on GRAPE
NASA Astrophysics Data System (ADS)
Fukushige, Toshiyuki; Kawai, Atsushi
2016-03-01
We present an implementation of the hierarchical tree algorithm on the individual timestep algorithm (the Hermite scheme) for collisional N-body simulations, running on the GRAPE-9 system, a special-purpose hardware accelerator for gravitational many-body simulations. Such a combination of the tree algorithm and the individual timestep algorithm was not easy on the previous GRAPE system mainly because its memory addressing scheme was limited only to sequential access to a full set of particle data. The present GRAPE-9 system has an indirect memory addressing unit and a particle memory large enough to store all the particle data and also the tree node data. The indirect memory addressing unit stores interaction lists for the tree algorithm, which is constructed on the host computer, and, according to the interaction lists, force pipelines calculate only the interactions necessary. In our implementation, the interaction calculations are significantly reduced compared to direct N2 summation in the original Hermite scheme. For example, we can achieve about a factor 30 of speedup (equivalent to about 17 teraflops) against the Hermite scheme for a simulation of an N = 106 system, using hardware of a peak speed of 0.6 teraflops for the Hermite scheme.
NASA Astrophysics Data System (ADS)
Sato, Masanori; Matsubara, Takahiko
2013-06-01
It is crucial to understand and model a behavior of galaxy biasing for future ambitious galaxy redshift surveys. Using 40 large cosmological N-body simulations for a standard ΛCDM cosmology, we study the cross-correlation coefficient between matter and the halo density field, which is an indicator of the stochasticity of bias, over a wide redshift range 0≤z≤3. The cross-correlation coefficient is important to extract information on the matter density field, e.g., by combining galaxy clustering and galaxy-galaxy lensing measurements. We compare the simulation results with integrated perturbation theory (iPT) proposed by one of the present authors and standard perturbation theory combined with a phenomenological model of local bias. The cross-correlation coefficient derived from the iPT agrees with N-body simulation results down to r˜15(10)h-1Mpc within 0.5 (1.0)% for all redshifts and halo masses we consider. The standard perturbation theory with local bias does not explain complicated behaviors on quasilinear scales at low redshifts, while roughly reproduces the general behavior of the cross-correlation coefficient on fully nonlinear scales. The iPT is powerful to predict the cross-correlation coefficient down to quasilinear regimes with a high precision.
Quantification of discreteness effects in cosmological N-body simulations: Initial conditions
Joyce, M.; Marcos, B.
2007-03-15
The relation between the results of cosmological N-body simulations, and the continuum theoretical models they simulate, is currently not understood in a way which allows a quantification of N dependent effects. In this first of a series of papers on this issue, we consider the quantification of such effects in the initial conditions of such simulations. A general formalism developed in [A. Gabrielli, Phys. Rev. E 70, 066131 (2004).] allows us to write down an exact expression for the power spectrum of the point distributions generated by the standard algorithm for generating such initial conditions. Expanded perturbatively in the amplitude of the input (i.e. theoretical, continuum) power spectrum, we obtain at linear order the input power spectrum, plus two terms which arise from discreteness and contribute at large wave numbers. For cosmological type power spectra, one obtains as expected, the input spectrum for wave numbers k smaller than that characteristic of the discreteness. The comparison of real space correlation properties is more subtle because the discreteness corrections are not as strongly localized in real space. For cosmological type spectra the theoretical mass variance in spheres and two-point correlation function are well approximated above a finite distance. For typical initial amplitudes this distance is a few times the interparticle distance, but it diverges as this amplitude (or, equivalently, the initial redshift of the cosmological simulation) goes to zero, at fixed particle density. We discuss briefly the physical significance of these discreteness terms in the initial conditions, in particular, with respect to the definition of the continuum limit of N-body simulations.
Rahman Prize Talk: Pushing the frontier in the simulation of correlated quantum many body systems
NASA Astrophysics Data System (ADS)
Troyer, Matthias
Amazing progress in the simulation of correlated quantum many body systems has been achieved in the past two decades by combining significant advances in new algorithms with efficient implementations on ever faster supercomputers. This has enabled the accurate simulation of an increasing number of problems and helped settle many open questions. I will review a selection of results that my collaborators and I have worked on, from quantum phase transitions in quantum magnets, over supersolidity of bosons in lattice models and Helium-4 to recent simulations of correlated fermions and quantum gases. I will then provide an outlook to the future and discuss how in the short term analog quantum simulators can help tackle problems for which no efficient simulation algorithms exist and how in the longer term quantum computers can be used to solve many of the still open questions in the field. I will finally connect to the topic of the remainder of this symposium by touching on how the design of new topological materials will help in the construction of these quantum computers.
Halo Statistics Analysis Within Medium Volume Cosmological N-Body Simulation
NASA Astrophysics Data System (ADS)
Martinovic, N.
2015-06-01
In this paper we present halo statistics analysis of a Λ CDM N-body cosmological simulation (from first halo formation until z = 0). We study mean major merger rate as a function of time, where for time we consider both per redshift and per Gyr dependence. For latter we find that it scales as the well known power law (1+z)^n for which we obtain n = 2.4. The halo mass function and halo growth function are derived and compared both with analytical and empirical fits. We analyse halo growth through out entire simulation, making it possible to continuously monitor evolution of halo number density within given mass ranges. The halo formation redshift is studied exploring possibility for a new simple preliminary analysis during the simulation run. Visualization of the simulation is portrayed as well. At redshifts z = 0-7 halos from simulation have good statistics for further analysis especially in mass range of 10^{11} - 10^{14} M_{⊙}/h.
Characterization of a lower-body exoskeleton for simulation of space-suited locomotion
NASA Astrophysics Data System (ADS)
Carr, Christopher E.; Newman, Dava J.
2008-02-01
In a previous analysis of suited and unsuited locomotion energetics, we found evidence that space suits act as springs during running. Video images from the lunar surface suggest that knee torques create, in large part, this spring effect. We hypothesized that a lower-body exoskeleton, properly constructed, could be used to simulate the knee torques of a range of space suits. Here we report characterization of a lower-body exoskeleton. Equivalent spring stiffness of each exoskeleton leg varies as a function of exoskeleton knee angle and load, and the exoskeleton joint-torque relationship closely matches the current NASA space suit, or Extravehicular Mobility Unit, knee torques in form and magnitude. We have built an exoskeleton with two physical non-linear springs, which achieve space-suit like joint-torques. Therefore space-suit legs act as springs, with this effect most pronounced when locomotion requires large changes in knee flexion such as during running.
MHD simulations of large conducting bodies moving through a planetary magnetosphere
NASA Astrophysics Data System (ADS)
Kopp, A.; Schröer, A.
1998-01-01
The interaction between a large body with finite conductivity moving relatively to a planetary magnetic field and the magnetospheric plasma is studied by means of numerical simulations in the framework of resistive MHD. Applications discussed are spacecrafts (tethered satellites) in the Earth's ionosphere and the interaction of Jupiter with its satellite Io. The interaction excites MHD waves in which consequence an electrical current flow through the conductor is generated which extends into the surrounding plasma and propagates along the Alfvén characteristics. Thus, kinetic energy is converted into electrical energy which may be utilized for satellite projects or explain Jupiter's decametric radiation. The linear and nonlinear evolution of the plasma environment and the current system is investigated with special care taken of the principal phenomena occuring in the vicinity of the conducting body.
Barred Galaxy Photometry: Comparing results from the Cananea sample with N-body simulations
NASA Astrophysics Data System (ADS)
Athanassoula, E.; Gadotti, D. A.; Carrasco, L.; Bosma, A.; de Souza, R. E.; Recillas, E.
2009-11-01
We compare the results of the photometrical analysis of barred galaxies with those of a similar analysis from N-body simulations. The photometry is for a sample of nine barred galaxies observed in the J and K[s] bands with the CANICA near infrared (NIR) camera at the 2.1 m telescope of the Observatorio Astrofísico Guillermo Haro (OAGH) in Cananea, Sonora, Mexico. The comparison includes radial ellipticity profiles and surface brightness (density for the N-body galaxies) profiles along the bar major and minor axes. We find very good agreement, arguing that the exchange of angular momentum within the galaxy plays a determinant role in the evolution of barred galaxies.
[Pertussis vaccines: acellular versus whole cell. Perhaps a return to the past?].
Cofré, José
2015-10-01
The resurgence of pertussis in the world and in our country has questioned the effectiveness of cellular and acellular vaccines. The reason why pertussis has not been controlled or eliminated after 70 years of implementation of the vaccination is probably multifactorial. This article, on the basis of questions and answers, describes the benefits and limitations of both cellular and acellular vaccines and suggests new strategies of vaccination in childhood. It is a fact that the currently applied vaccination does not eliminate the circulation of Bordetella pertussis in the community. Perhaps the introduction of vaccines with live B. pertussis, inhalation, will be able to eliminate the disease around the world. PMID:26633113
Plastic Surgery and Acellular Dermal Matrix: Highlighting Trends from 1999 to 2013
Daar, David A; Gandy, Jessica R; Clark, Emily G; Mowlds, Donald S; Paydar, Keyianoosh Z; Wirth, Garrett A
2016-01-01
The last decade has ushered in a rapidly expanding global discussion regarding acellular dermal matrix (ADM) applications, economic analyses, technical considerations, benefits, and risks, with recent emphasis on ADM use in breast surgery. This study aims to evaluate global trends in ADM research using bibliometric analysis. The top nine Plastic Surgery journals were determined by impact factor (IF). Each issue of the nine journals between 1999 and 2013 was accessed to compile a database of articles discussing ADM. Publications were further classified by IF, authors’ geographic location, study design, and level of evidence (LOE, I-V). Productivity index and productivity share were calculated for each region. In total, 256 ADM articles were accessed. The annual global publication volume increased significantly by 4.2 (0.87) articles per year (p<0.001), with a mean productivity index of 36.3 (59.0). The mean impact factor of the nine journals increased significantly from 0.61 (0.11) to 2.47 (0.99) from 1993 to 2013 (p<0.001). Despite this increase in the global ADM literature, the majority of research was of weaker LOE (level I: 2.29% and level II: 9.17%). USA contributed the most research (87%), followed by Asia (4.76%) and Western Europe (4.71%). USA contributed the greatest volume of research. Regarding clinical application of ADM, the majority of publications focused on ADM use in breast surgery, specifically breast reconstruction (154 articles, 60.2%). The majority of research was of lower LOE; thus, efforts should be made to strengthen the body of literature, particularly with regard to cost analysis. PMID:27579264
Plastic Surgery and Acellular Dermal Matrix: Highlighting Trends from 1999 to 2013.
Daar, David A; Gandy, Jessica R; Clark, Emily G; Mowlds, Donald S; Paydar, Keyianoosh Z; Wirth, Garrett A
2016-05-01
The last decade has ushered in a rapidly expanding global discussion regarding acellular dermal matrix (ADM) applications, economic analyses, technical considerations, benefits, and risks, with recent emphasis on ADM use in breast surgery. This study aims to evaluate global trends in ADM research using bibliometric analysis. The top nine Plastic Surgery journals were determined by impact factor (IF). Each issue of the nine journals between 1999 and 2013 was accessed to compile a database of articles discussing ADM. Publications were further classified by IF, authors' geographic location, study design, and level of evidence (LOE, I-V). Productivity index and productivity share were calculated for each region. In total, 256 ADM articles were accessed. The annual global publication volume increased significantly by 4.2 (0.87) articles per year (p<0.001), with a mean productivity index of 36.3 (59.0). The mean impact factor of the nine journals increased significantly from 0.61 (0.11) to 2.47 (0.99) from 1993 to 2013 (p<0.001). Despite this increase in the global ADM literature, the majority of research was of weaker LOE (level I: 2.29% and level II: 9.17%). USA contributed the most research (87%), followed by Asia (4.76%) and Western Europe (4.71%). USA contributed the greatest volume of research. Regarding clinical application of ADM, the majority of publications focused on ADM use in breast surgery, specifically breast reconstruction (154 articles, 60.2%). The majority of research was of lower LOE; thus, efforts should be made to strengthen the body of literature, particularly with regard to cost analysis. PMID:27579264
Carruthers, Christopher A.; Dearth, Christopher L.; Reing, Janet E.; Kramer, Caroline R.; Gagne, Darcy H.; Crapo, Peter M.; Garcia, Onelio; Badhwar, Amit; Scott, Jeffrey R.
2015-01-01
Background: Acellular dermal matrices (ADMs) have been commonly used in expander-based breast reconstruction to provide inferolateral prosthesis coverage. Although the clinical performance of these biologic scaffold materials varies depending on a number of factors, an in-depth systematic characterization of the host response is yet to be performed. The present study evaluates the biochemical composition and structure of two ADMs, AlloDerm® Regenerative Tissue Matrix and AlloMax™ Surgical Graft, and provides a comprehensive spatiotemporal characterization in a porcine model of tissue expander breast reconstruction. Methods: Each ADM was characterized with regard to thickness, permeability, donor nucleic acid content, (residual double-stranded DNA [dsDNA]), and growth factors (basic fibroblast growth factor [bFGF], vascular endothelial growth factor [VEGF], and transforming growth factor-beta 1 [TGF-β1]). Cytocompatibility was evaluated by in vitro cell culture on the ADMs. The host response was evaluated at 4 and 12 weeks at various locations within the ADMs using established metrics of the inflammatory and tissue remodeling response: cell infiltration, multinucleate giant cell formation, extent of ADM remodeling, and neovascularization. Results: AlloMax incorporated more readily with surrounding host tissue as measured by earlier and greater cell infiltration, fewer foreign body giant cells, and faster remodeling of ADM. These findings correlated with the in vitro composition and cytocompatibility analysis, which showed AlloMax to more readily support in vitro cell growth. Conclusions: AlloMax and AlloDerm demonstrated distinct remodeling characteristics in a porcine model of tissue expander breast reconstruction. PMID:24941900
NASA Astrophysics Data System (ADS)
Ning, Chengyun; Cheng, Haimei; Zhu, Wenjun; Yin, Zhaoyi; Chen, Hao; Zheng, Huade; Lei, Shumei; Yin, Shiheng; Tan, Guoxin
2008-11-01
Porous polycaprolactone (PCL) scaffolds were fabricated by combination of porogen-leaching and freeze-drying processes. Ice particulates were used as porogen materials. The porous PCL scaffolds were modified by potassium hydroxide solution with concentration of 1 mol/L at room temperature for 8 h, subsequently biomineralized in simulated body fluid for 2 h and 8 h, respectively. The microstructure and characteristics of the PCL scaffolds were investigated by scanning electron microscope (SEM) and EDS. The results showed (1) PCL scaffolds had high degree of connectivity and different pore sizes. (2) Plate-like apatite was observed on the surface of the scaffolds after being immersed into SBF for 8 h.
Molecular dynamics simulation of interparticle spacing and many-body effect in gold supracrystals.
Liu, X P; Ni, Y; He, L H
2016-04-01
Interparticle spacing in supracrystals is a crucial parameter for photoelectric applications as it dominates the transport rates between neighboring nanoparticles (NPs). Based on large-scale molecular dynamics simulations, we calculate interparticle spacing in alkylthiol-stabilized gold supracrystals as a function of the NP size, ligand length and external pressure. The repulsive many-body interactions in the supracrystals are also quantified by comparing the interparticle spacing with that between two individual NPs at equilibrium. Our results are consistent with available experiments, and are expected to help precise control of interparticle spacing in supracrystal devices. PMID:26909856
Efficient time-symmetric simulation of torqued rigid bodies using Jacobi elliptic functions
NASA Astrophysics Data System (ADS)
Celledoni, E.; Säfström, N.
2006-05-01
If the three moments of inertia are distinct, the solution to the Euler equations for the free rigid body is given in terms of Jacobi elliptic functions. Using the arithmetic-geometric mean algorithm (Abramowitz and Stegun 1992 Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (New York: Dover)), these functions can be calculated efficiently and accurately. Compared to standard numerical ODE and Lie-Poisson solvers, the overall approach yields a faster and more accurate numerical solution to the Euler equations. This approach is designed for mass asymmetric rigid bodies. In the case of symmetric bodies, the exact solution is available in terms of trigonometric functions, see Dullweber et al (1997 J. Chem. Phys. 107 5840-51), Reich (1996 Fields Inst. Commun. 10 181-91) and Benettin et al (2001 SIAM J. Sci. Comp. 23 1189-203) for details. In this paper, we consider the case of asymmetric rigid bodies subject to external forces. We consider a strategy similar to the symplectic splitting method proposed in Reich (1996 Fields Inst. Commun. 10 181-91) and Dullweber et al (1997 J. Chem. Phys. 107 5840-51). The method proposed here is time-symmetric. We decompose the vector field of our problem into a free rigid body (FRB) problem and another completely integrable vector field. The FRB problem consists of the Euler equations and a differential equation for the 3 × 3 orientation matrix. The Euler equations are integrated exactly while the matrix equation is approximated using a truncated Magnus series. In our experiments, we observe that the overall numerical solution benefits greatly from the very accurate solution of the Euler equations. We apply the method to the heavy top and the simulation of artificial satellite attitude dynamics.
Higher moments of primordial non-Gaussianity and N-body simulations
Adhikari, Saroj; Shandera, Sarah; Dalal, Neal E-mail: shandera@gravity.psu.edu
2014-06-01
We perform cosmological N-body simulations with non-Gaussian initial conditions generated from two independent fields. The dominant contribution to the perturbations comes from a purely Gaussian field, but we allow the second field to have local non-Gaussianity that need not be weak. This scenario allows us to adjust the relative importance of non-Gaussian contributions beyond the skewness, producing a scaling of the higher moments different from (and stronger than) the scaling in the usual single field local ansatz. We compare semi-analytic prescriptions for the non-Gaussian mass function, large scale halo bias, and stochastic bias against the simulation results. We discuss applications of this work to large scale structure measurements that can test a wider range of models for the primordial fluctuations than is usually explored.
Relationship between trawl selectivity and fish body size in a simulated population
NASA Astrophysics Data System (ADS)
Sun, Peng; Liang, Zhenlin; Huang, Liuyi; Tang, Yanli; He, Xin
2013-03-01
Trawl is a main fishing gear in Chinese fishery, capturing large fish and letting small ones at large. However, long-term use of trawl would result in changes of phenotypic traits of the fish stocks, such as smaller size-at-age and earlier age-at-maturation. In this study, we simulated a fish population with size characteristics of trawl fishing and the population produces one generation of offspring and lives for one year, used trawl to exploit the simulated fish population, and captured individuals by body size. We evaluated the impact of the changes on selectivity parameters, such as selective range and the length at 50% retention. Under fishing pressure, we specified the selectivity parameters, and determined that smaller selection rates and greater length at 50% retention were associated with an increased tendency towards miniaturization.
NASA Technical Reports Server (NTRS)
Otto, John C.; Paraschivoiu, Marius; Yesilyurt, Serhat; Patera, Anthony T.
1995-01-01
Engineering design and optimization efforts using computational systems rapidly become resource intensive. The goal of the surrogate-based approach is to perform a complete optimization with limited resources. In this paper we present a Bayesian-validated approach that informs the designer as to how well the surrogate performs; in particular, our surrogate framework provides precise (albeit probabilistic) bounds on the errors incurred in the surrogate-for-simulation substitution. The theory and algorithms of our computer{simulation surrogate framework are first described. The utility of the framework is then demonstrated through two illustrative examples: maximization of the flowrate of fully developed ow in trapezoidal ducts; and design of an axisymmetric body that achieves a target Stokes drag.
Recent advances in numerical simulation and control of asymmetric flows around slender bodies
NASA Technical Reports Server (NTRS)
Kandil, Osama A.; Wong, Tin-Chee; Sharaf, Hazem H.; Liu, C. H.
1992-01-01
The problems of asymmetric flow around slender bodies and its control are formulated using the unsteady, compressible, thin-layer or full Navier-Stokes equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. The problem is numerically simulated for both locally-conical and three-dimensional flows. The numerical applications include studies of the effects of relative incidence, Mach number and Reynolds number on the flow asymmetry. For the control of flow asymmetry, the numerical simulation cover passive and active control methods. For the passive control, the effectiveness of vertical fins placed in the leeward plane of geometric symmetry and side strakes with different orientations is studied. For the active control, the effectiveness of normal and tangential flow injection and surface heating and a combination of these methods is studied.
NASA Astrophysics Data System (ADS)
Waye, A. B.; Krygiel, R. G.; Susin, T. B.; Baptista, R.; Rehnberg, L.; Heidner, G. S.; de Campos, F.; Falcão, F. P.; Russomano, T.
2013-09-01
Performance of efficient single-person cardiopulmonary resuscitation (CPR) is vital to maintain cardiac and cerebral perfusion during the 2-4 min it takes for deployment of advanced life support during a space mission. The aim of the present study was to investigate potential differences in upper body muscle activity during CPR performance at terrestrial gravity (+1Gz) and in simulated microgravity (μG). Muscle activity of the triceps brachii, erector spinae, rectus abdominis and pectoralis major was measured via superficial electromyography in 20 healthy male volunteers. Four sets of 30 external chest compressions (ECCs) were performed on a mannequin. Microgravity was simulated using a body suspension device and harness; the Evetts-Russomano (ER) method was adopted for CPR performance in simulated microgravity. Heart rate and perceived exertion via Borg scores were also measured. While a significantly lower depth of ECCs was observed in simulated microgravity, compared with +1Gz, it was still within the target range of 40-50 mm. There was a 7.7% decrease of the mean (±SEM) ECC depth from 48 ± 0.3 mm at +1Gz, to 44.3 ± 0.5 mm during microgravity simulation (p < 0.001). No significant difference in number or rate of compressions was found between the two conditions. Heart rate displayed a significantly larger increase during CPR in simulated microgravity than at +1Gz, the former presenting a mean (±SEM) of 23.6 ± 2.91 bpm and the latter, 76.6 ± 3.8 bpm (p < 0.001). Borg scores were 70% higher post-microgravity compressions (17 ± 1) than post +1Gz compressions (10 ± 1) (p < 0.001). Intermuscular comparisons showed the triceps brachii to have significantly lower muscle activity than each of the other three tested muscles, in both +1Gz and microgravity. As shown by greater Borg scores and heart rate increases, CPR performance in simulated microgravity is more fatiguing than at +1Gz. Nevertheless, no significant difference in muscle activity between conditions
Maps of CMB lensing deflection from N-body simulations in Coupled Dark Energy Cosmologies
Carbone, Carmelita; Baldi, Marco; Baccigalupi, Carlo E-mail: marco.baldi5@unibo.it E-mail: bacci@sissa.it
2013-09-01
We produce lensing potential and deflection-angle maps in order to simulate the weak gravitational lensing of the Cosmic Microwave Background (CMB) via ray-tracing through the COupled Dark Energy Cosmological Simulations (CoDECS), the largest suite of N-body simulations to date for interacting Dark Energy cosmologies. The constructed maps faithfully reflect the N-body cosmic structures on a range of scales going from the arcminute to the degree scale, limited only by the resolution and extension of the simulations. We investigate the variation of the lensing pattern due to the underlying Dark Energy (DE) dynamics, characterised by different background and perturbation behaviours as a consequence of the interaction between the DE field and Cold Dark Matter (CDM). In particular, we study in detail the results from three cosmological models differing in the background and perturbations evolution at the epoch in which the lensing cross section is most effective, corresponding to a redshift of ∼ 1, with the purpose to isolate their imprints in the lensing observables, regardless of the compatibility of these models with present constraints. The scenarios investigated here include a reference ΛCDM cosmology, a standard coupled DE (cDE) scenario, and a ''bouncing'' cDE scenario. For the standard cDE scenario, we find that typical differences in the lensing potential result from two effects: the enhanced growth of linear CDM density fluctuations with respect to the ΛCDM case, and the modified nonlinear dynamics of collapsed structures induced by the DE-CDM interaction. As a consequence, CMB lensing highlights the DE impact in the cosmological expansion, even in the degenerate case where the amplitude of the linear matter density perturbations, parametrised through σ{sub 8}, is the same in both the standard cDE and ΛCDM cosmologies. For the ''bouncing'' scenario, we find that the two opposite behaviours of the lens density contrast and of the matter abundance lead to
N-body simulations with generic non-Gaussian initial conditions II: halo bias
NASA Astrophysics Data System (ADS)
Wagner, Christian; Verde, Licia
2012-03-01
We present N-body simulations for generic non-Gaussian initial conditions with the aim of exploring and modelling the scale-dependent halo bias. This effect is evident on very large scales requiring large simulation boxes. In addition, the previously available prescription to implement generic non-Gaussian initial conditions has been improved to keep under control higher-order terms which were spoiling the power spectrum on large scales. We pay particular attention to the differences between physical, inflation-motivated primordial bispectra and their factorizable templates, and to the operational definition of the non-Gaussian halo bias (which has both a scale-dependent and an approximately scale-independent contributions). We find that analytic predictions for both the non-Gaussian halo mass function and halo bias work well once a fudge factor (which was introduced before but still lacks convincing physical explanation) is calibrated on simulations. The halo bias remains therefore an extremely promising tool to probe primordial non-Gaussianity and thus to give insights into the physical mechanism that generated the primordial perturbations. The simulation outputs and tables of the analytic predictions will be made publicly available via the non-Gaussian comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html.
Algorithm for simulation of quantum many-body dynamics using dynamical coarse-graining
NASA Astrophysics Data System (ADS)
Khasin, M.; Kosloff, R.
2010-04-01
An algorithm for simulation of quantum many-body dynamics having su(2) spectrum-generating algebra is developed. The algorithm is based on the idea of dynamical coarse-graining. The original unitary dynamics of the target observables—the elements of the spectrum-generating algebra—is simulated by a surrogate open-system dynamics, which can be interpreted as weak measurement of the target observables, performed on the evolving system. The open-system state can be represented by a mixture of pure states, localized in the phase space. The localization reduces the scaling of the computational resources with the Hilbert-space dimension n by factor n3/2(lnn)-1 compared to conventional sparse-matrix methods. The guidelines for the choice of parameters for the simulation are presented and the scaling of the computational resources with the Hilbert-space dimension of the system is estimated. The algorithm is applied to the simulation of the dynamics of systems of 2×104 and 2×106 cold atoms in a double-well trap, described by the two-site Bose-Hubbard model.
Stehr, K; Heininger, U; Uhlenbusch, R; Angersbach, P; Hackell, J; Eckhardt, T
1995-03-01
Immunization against pertussis has been re-recommended for healthy children in Germany in 1991. In addition the former restriction of immunizing only in the first 2 years of life was abolished. In children born before 1991 immunization rates against pertussis were 15% or less. With the new recommendations physicians are now faced with an increasing demand of parents for catch-up vaccinations in these children. Since they were immunized against diphtheria and tetanus previously monovalent pertussis vaccines are needed for this indication. Therefore a monovalent, multicomponent acellular pertussis vaccine was studied in 249 German children 15 months to 6 years of age. Three doses were administered at 6-10 week intervals. Reactogenicity and antibody responses against the vaccine antigens pertussis toxin (PT), filamentous haemagglutinin (FHA), 69-kd antigen (pertactin) and fimbriae-2 (agglutinogen) were investigated. Local and systemic reactions were minimal in frequency and severity. Antibody responses against all vaccine antigens were pronounced with 93%-100% of vaccinees demonstrating at least four fold titre rises above pre-immunization after the third dose. These findings indicate that this monovalent, multicomponent acellular pertussis vaccine with excellent immunogenicity and low reactogenicity is an appropriate candidate for closing immunization gaps in older children in countries with previously low vaccination rates against pertussis. Based on the results of this study the monovalent acellular pertussis vaccine was licensed in Germany in January 1994. PMID:7758519
Potential sites for the perception of gravity in the acellular slime mold Physarum polycephalum
NASA Astrophysics Data System (ADS)
Block, I.; Briegleb, W.
Recently a gravisensitivity of the acellular slime mold Physarum polycephalum, which possesses no specialized gravireceptor, could be established by conducting experiments under simulated and under real near weightlessness. In these experiments macroplasmodia showed a modulation of their contraction rhythm followed by regulation phenomena. Until now the perception mechanism for the gravistimulus is unknown, but several findings indicate the involvement of mitochondria: A) During the impediment of respiration the Og-reaction is inhibited and the regulation is reduced. B) The response to a light stimulus and the following regulation phenomena strongly resemble the behavior during exposure to Og, the only difference is that the two reactions are directed into opposite directions. In the blue-light reaction a flavin of the mitochondrial matrix seems to be involved in the light perception. C) The contraction rhythm as well as its modulations are coupled to rhythmic changes in the levels of ATP and calcium ions, involving the mitochondria as sites of energy production and of Ca++-storage. - So the mitochondria could be the site of the regulation and they possibly are the receptor sites for the light and gravity stimuli. - Also the observation of a morphologic polarity of the slime mold's plasmodial strands has to be considered: Cross-sections reveal that the ectoplasmic wall surrounding the streaming endoplasm is much thinner on the physically lower side than on the upper side of the strand - this applies to strands lying on or hanging on a horizontal surface. So, in addition to the mitochondria, also the morphologic polarity may be involved in the perception mechanism of the observed gravisensitivity and of the recently established geotaxis. - The potential role of the nuclei and of the contractile elements in the perception of gravity is also discussed.
When clusters collide - A numerical Hydro/N-body simulation of merging galaxy clusters
NASA Technical Reports Server (NTRS)
Roettiger, Kurt; Burns, Jack; Loken, Chris
1993-01-01
A 3D numerical simulation of two merging clusters of galaxies, using a hybrid Hydro/N-body code, is presented. The hydrodynamics of the code is solved by an Eulerian finite difference method. Initial results disclose that the X-ray emission of the dominant cluster becomes elongated and broadened; heating occurs at the core of the dominant cluster as a result of multiple shocks, and high velocity gas motions within the intracluster medium. It is predicted that clusters which have undergone recent mergers and do not have cooling flows will have high peculiar gas velocities and that the shocks and turbulence generated during the merger may power cluster-wide radio halos. Prolonged high-velocity gas motions through the dominant cluster core possibly play a major role in the formation and shaping of wide-angle tailed radio sources associated with central dominant galaxies. The N-body component of the simulation reveals the subcluster to be dispersed as it passes through the dominant cluster.
Robust human body model injury prediction in simulated side impact crashes.
Golman, Adam J; Danelson, Kerry A; Stitzel, Joel D
2016-01-01
This study developed a parametric methodology to robustly predict occupant injuries sustained in real-world crashes using a finite element (FE) human body model (HBM). One hundred and twenty near-side impact motor vehicle crashes were simulated over a range of parameters using a Toyota RAV4 (bullet vehicle), Ford Taurus (struck vehicle) FE models and a validated human body model (HBM) Total HUman Model for Safety (THUMS). Three bullet vehicle crash parameters (speed, location and angle) and two occupant parameters (seat position and age) were varied using a Latin hypercube design of Experiments. Four injury metrics (head injury criterion, half deflection, thoracic trauma index and pelvic force) were used to calculate injury risk. Rib fracture prediction and lung strain metrics were also analysed. As hypothesized, bullet speed had the greatest effect on each injury measure. Injury risk was reduced when bullet location was further from the B-pillar or when the bullet angle was more oblique. Age had strong correlation to rib fractures frequency and lung strain severity. The injuries from a real-world crash were predicted using two different methods by (1) subsampling the injury predictors from the 12 best crush profile matching simulations and (2) using regression models. Both injury prediction methods successfully predicted the case occupant's low risk for pelvic injury, high risk for thoracic injury, rib fractures and high lung strains with tight confidence intervals. This parametric methodology was successfully used to explore crash parameter interactions and to robustly predict real-world injuries. PMID:26158552
Kaur, Balwinder; Srivastava, Rajendra; Satpati, Biswarup; Kondepudi, Kanthi Kiran; Bishnoi, Mahendra
2015-11-01
Silver ion-exchanged nanocrystalline zeolite (Ag-Nano-ZSM-5) and silver ion-exchanged conventional zeolite (Ag-ZSM-5) were synthesized. Zeolites were incubated in simulated body fluid at 310K for different time periods to grow hydroxyapatite in their matrixes. Significant large amount of hydroxyapatite was grown in Ag-Nano-ZSM-5 matrix after incubation in simulated body fluid when compared to Ag-ZSM-5. The resultant material was characterized using X-ray diffraction, N2-adsorption, scanning/transmission electron microscopy, energy dispersive X-ray, and inductively coupled plasma analysis. Mechanical properties such as compressive modulus, compressive strength, and strain at failure of the parent materials were evaluated. Biocompatibility assays suggested that Ag-Nano-ZSM-5 and hydroxyapatite grown in Ag-Nano-ZSM-5 were compatible and did not impose any toxicity to RAW 264.7 cells macrophase and Caco2 cells suggesting considerable potential for biomedical applications such as bone implants. PMID:26255163
Peculiar velocities in redshift space: formalism, N-body simulations and perturbation theory
Okumura, Teppei; Seljak, Uroš; Vlah, Zvonimir; Desjacques, Vincent E-mail: useljak@berkeley.edu E-mail: Vincent.Desjacques@unige.ch
2014-05-01
Direct measurements of peculiar velocities of galaxies and clusters of galaxies can in principle provide explicit information on the three dimensional mass distribution, but this information is modulated by the fact that velocity field is sampled at galaxy positions, and is thus probing galaxy momentum. We derive expressions for the cross power spectrum between the density and momentum field and the auto spectrum of the momentum field in redshift space, by extending the distribution function method to these statistics. The resulting momentum cross and auto power spectra in redshift space are expressed as infinite sums over velocity moment correlators in real space, as is the case for the density power spectrum in redshift space. We compute each correlator using Eulerian perturbation theory (PT) and halo biasing model and compare the resulting redshift-space velocity statistics to those measured from N-body simulations for both dark matter and halos. We find that in redshift space linear theory predictions for the density-momentum cross power spectrum as well as for the momentum auto spectrum fail to predict the N-body results at very large scales. On the other hand, our nonlinear PT prediction for these velocity statistics, together with real-space power spectrum for dark matter from simulations, improves the accuracy for both dark matter and halos. We also present the same analysis in configuration space, computing the redshift-space pairwise mean infall velocities and velocity correlation function and compare to nonlinear PT.
SVD-GFD scheme to simulate complex moving body problems in 3D space
NASA Astrophysics Data System (ADS)
Wang, X. Y.; Yu, P.; Yeo, K. S.; Khoo, B. C.
2010-03-01
The present paper presents a hybrid meshfree-and-Cartesian grid method for simulating moving body incompressible viscous flow problems in 3D space. The method combines the merits of cost-efficient and accurate conventional finite difference approximations on Cartesian grids with the geometric freedom of generalized finite difference (GFD) approximations on meshfree grids. Error minimization in GFD is carried out by singular value decomposition (SVD). The Arbitrary Lagrangian-Eulerian (ALE) form of the Navier-Stokes equations on convecting nodes is integrated by a fractional-step projection method. The present hybrid grid method employs a relatively simple mode of nodal administration. Nevertheless, it has the geometrical flexibility of unstructured mesh-based finite-volume and finite element methods. Boundary conditions are precisely implemented on boundary nodes without interpolation. The present scheme is validated by a moving patch consistency test as well as against published results for 3D moving body problems. Finally, the method is applied on low-Reynolds number flapping wing applications, where large boundary motions are involved. The present study demonstrates the potential of the present hybrid meshfree-and-Cartesian grid scheme for solving complex moving body problems in 3D.
Borazjani, Iman; Ge, Liang; Sotiropoulos, Fotis
2008-08-10
The sharp-interface CURVIB approach of Ge and Sotiropoulos [L. Ge, F. Sotiropoulos, A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries, Journal of Computational Physics 225 (2007) 1782-1809] is extended to simulate fluid structure interaction (FSI) problems involving complex 3D rigid bodies undergoing large structural displacements. The FSI solver adopts the partitioned FSI solution approach and both loose and strong coupling strategies are implemented. The interfaces between immersed bodies and the fluid are discretized with a Lagrangian grid and tracked with an explicit front-tracking approach. An efficient ray-tracing algorithm is developed to quickly identify the relationship between the background grid and the moving bodies. Numerical experiments are carried out for two FSI problems: vortex induced vibration of elastically mounted cylinders and flow through a bileaflet mechanical heart valve at physiologic conditions. For both cases the computed results are in excellent agreement with benchmark simulations and experimental measurements. The numerical experiments suggest that both the properties of the structure (mass, geometry) and the local flow conditions can play an important role in determining the stability of the FSI algorithm. Under certain conditions unconditionally unstable iteration schemes result even when strong coupling FSI is employed. For such cases, however, combining the strong-coupling iteration with under-relaxation in conjunction with the Aitken's acceleration technique is shown to effectively resolve the stability problems. A theoretical analysis is presented to explain the findings of the numerical experiments. It is shown that the ratio of the added mass to the mass of the structure as well as the sign of the local time rate of change of the force or moment imparted on the structure by the fluid determine the stability and convergence of the FSI
Borazjani, Iman; Ge, Liang; Sotiropoulos, Fotis
2010-01-01
The sharp-interface CURVIB approach of Ge and Sotiropoulos [L. Ge, F. Sotiropoulos, A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries, Journal of Computational Physics 225 (2007) 1782–1809] is extended to simulate fluid structure interaction (FSI) problems involving complex 3D rigid bodies undergoing large structural displacements. The FSI solver adopts the partitioned FSI solution approach and both loose and strong coupling strategies are implemented. The interfaces between immersed bodies and the fluid are discretized with a Lagrangian grid and tracked with an explicit front-tracking approach. An efficient ray-tracing algorithm is developed to quickly identify the relationship between the background grid and the moving bodies. Numerical experiments are carried out for two FSI problems: vortex induced vibration of elastically mounted cylinders and flow through a bileaflet mechanical heart valve at physiologic conditions. For both cases the computed results are in excellent agreement with benchmark simulations and experimental measurements. The numerical experiments suggest that both the properties of the structure (mass, geometry) and the local flow conditions can play an important role in determining the stability of the FSI algorithm. Under certain conditions unconditionally unstable iteration schemes result even when strong coupling FSI is employed. For such cases, however, combining the strong-coupling iteration with under-relaxation in conjunction with the Aitken’s acceleration technique is shown to effectively resolve the stability problems. A theoretical analysis is presented to explain the findings of the numerical experiments. It is shown that the ratio of the added mass to the mass of the structure as well as the sign of the local time rate of change of the force or moment imparted on the structure by the fluid determine the stability and convergence of the
Dynamic modeling and simulation of multi-body systems using the Udwadia-Kalaba theory
NASA Astrophysics Data System (ADS)
Zhao, Han; Zhen, Shengchao; Chen, Ye-Hwa
2013-09-01
Laboratory experiments were conducted for falling U-chain, but explicit analytic form of the general equations of motion was not presented. Several modeling methods were developed for fish robots, however they just focused on the whole fish’s locomotion which does little favor to understand the detailed swimming behavior of fish. Udwadia-Kalaba theory is used to model these two multi-body systems and obtain explicit analytic equations of motion. For falling U-chain, the mass matrix is non-singular. Second-order constraints are used to get the constraint force and equations of motion and the numerical simulation is conducted. Simulation results show that the chain tip falls faster than the freely falling body. For fish robot, two-joint Carangiform fish robot is focused on. Quasi-steady wing theory is used to approximately calculate fluid lift force acting on the caudal fin. Based on the obtained explicit analytic equations of motion (the mass matrix is singular), propulsive characteristics of each part of the fish robot are obtained. Through these two cases of U chain and fish robot, how to use Udwadia-Kalaba equation to obtain the dynamical model is shown and the modeling methodology for multi-body systems is presented. It is also shown that Udwadia-Kalaba theory is applicable to systems whether or not their mass matrices are singular. In the whole process of applying Udwadia-Kalaba equation, Lagrangian multipliers and quasi-coordinates are not used. Udwadia-Kalaba theory is creatively applied to dynamical modeling of falling U-chain and fish robot problems and explicit analytic equations of motion are obtained.
Corneal Stroma Regeneration with Acellular Corneal Stroma Sheets and Keratocytes in a Rabbit Model
Ma, Xiao Yun; Zhang, Yun; Zhu, Dan; Lu, Yang; Zhou, Guangdong; Liu, Wei; Cao, Yilin; Zhang, Wen Jie
2015-01-01
Acellular corneal stroma matrix has been used for corneal stroma engineering. However, because of its compact tissue structure, regrowth of keratocytes into the scaffold is difficult. Previously, we developed a sandwich model for cartilage engineering using acellular cartilage sheets. In the present study, we tested this model for corneal stroma regeneration using acellular porcine corneal stroma (APCS) sheets and keratocytes. Porcine corneas were decellularized by NaCl treatment, and the APCS was cut into 20-μm-thick sheets. A rabbit corneal stroma defect model was created by lamellar keratoplasty and repaired by transplantation of five pieces of APCS sheets with keratocytes. Six months after transplantation, transparent corneas were present in the experimental group, which were confirmed by anterior segment optical coherence tomography examination and transmittance examination. The biomechanical properties in the experimental group were similar to those of normal cornea. Histological analyses showed an even distribution of keratocytes and well-oriented matrix in the stroma layer in the experimental group. Together, these results demonstrated that the sandwich model using acellular corneal stroma sheets and keratocytes could be potentially useful for corneal stroma regeneration. PMID:26167895
Repair of a Gingival Fenestration Using an Acellular Dermal Matrix Allograft.
Breault, Lawrence G; Brentson, Raquel C; Fowler, Edward B; Bisch, Frederick C
2016-01-01
A case report illustrating the successful treatment of a gingival fenestration with an acellular dermal matrix (ADM) allograft. After 2½ months of healing, the ADM was completely integrated into the soft tissues of the mandibular anterior gingiva with complete resolution of the gingival fenestration, resulting in excellent gingival esthetics. PMID:26874103
NASA Technical Reports Server (NTRS)
Chan, William M.
1992-01-01
The following papers are presented: (1) numerical methods for the simulation of complex multi-body flows with applications for the Integrated Space Shuttle vehicle; (2) a generalized scheme for 3-D hyperbolic grid generation; (3) collar grids for intersecting geometric components within the Chimera overlapped grid scheme; and (4) application of the Chimera overlapped grid scheme to simulation of Space Shuttle ascent flows.
Yamin, Stephanie; Stinchcombe, Arne; Gagnon, Sylvain
2015-01-01
Driving is a multifactorial behaviour drawing on multiple cognitive, sensory, and physical systems. Dementia is a progressive and degenerative neurological condition that impacts the cognitive processes necessary for safe driving. While a number of studies have examined driving among individuals with Alzheimer's disease, less is known about the impact of Dementia with Lewy Bodies (DLB) on driving safety. The present study compared simulated driving performance of 15 older drivers with mild DLB with that of 21 neurologically healthy control drivers. DLB drivers showed poorer performance on all indicators of simulated driving including an increased number of collisions in the simulator and poorer composite indicators of overall driving performance. A measure of global cognitive function (i.e., the Mini Mental State Exam) was found to be related to the overall driving performance. In addition, measures of attention (i.e., Useful Field of View, UFOV) and space processing (Visual Object and Space Perception, VOSP, Test) correlated significantly with a rater's assessment of driving performance. PMID:26713169
Experimental quantum simulations of many-body physics with trapped ions.
Schneider, Ch; Porras, Diego; Schaetz, Tobias
2012-02-01
Direct experimental access to some of the most intriguing quantum phenomena is not granted due to the lack of precise control of the relevant parameters in their naturally intricate environment. Their simulation on conventional computers is impossible, since quantum behaviour arising with superposition states or entanglement is not efficiently translatable into the classical language. However, one could gain deeper insight into complex quantum dynamics by experimentally simulating the quantum behaviour of interest in another quantum system, where the relevant parameters and interactions can be controlled and robust effects detected sufficiently well. Systems of trapped ions provide unique control of both the internal (electronic) and external (motional) degrees of freedom. The mutual Coulomb interaction between the ions allows for large interaction strengths at comparatively large mutual ion distances enabling individual control and readout. Systems of trapped ions therefore exhibit a prominent system in several physical disciplines, for example, quantum information processing or metrology. Here, we will give an overview of different trapping techniques of ions as well as implementations for coherent manipulation of their quantum states and discuss the related theoretical basics. We then report on the experimental and theoretical progress in simulating quantum many-body physics with trapped ions and present current approaches for scaling up to more ions and more-dimensional systems. PMID:22790343
ATOMISTIC SIMULATIONS OF DIFFUSIONAL CREEP IN A NANOCRYSTALLINE BODY-CENTERED CUBIC MATERIAL
Paul C. Millett; Tapan Desai; Vesselin Yamakov; Dieter Wolf
2008-08-01
Molecular dynamics (MD) simulations are used to study diffusion-accommodated creep deformation in nanocrystalline molybdenum, a body-centered cubic metal. In our simulations, the microstructures are subjected to constant-stress loading at levels below the dislocation nucleation threshold and at high temperatures (i.e., T > 0.75Tmelt), thereby ensuring that the overall deformation is indeed attributable to atomic self-diffusion. The initial microstructures were designed to consist of hexagonally shaped columnar grains bounded by high-energy asymmetric tilt grain boundaries (GBs). Remarkably the creep rates, which exhibit a double-exponential dependence on temperature and a double power-law dependence on grain size, indicate that both GB diffusion in the form of Coble creep and lattice diffusion in the form of Nabarro–Herring creep contribute to the overall deformation. For the first time in an MD simulation, we observe the formation and emission of vacancies from high-angle GBs into the grain interiors, thus enabling bulk diffusion.
`Grandeur in this view of life': N-body simulation models of the Galactic habitable zone
NASA Astrophysics Data System (ADS)
Vukotić, B.; Steinhauser, D.; Martinez-Aviles, G.; Ćirković, M. M.; Micic, M.; Schindler, S.
2016-07-01
We present an isolated Milky-Way-like simulation in the GADGET2 N-body smoothed particle hydrodynamics (SPH) code. The Galactic disc star formation rate (SFR) surface densities and a stellar mass indicative of the Solar neighbourhood are used as thresholds to model the distribution of stellar mass in life-friendly environments. SFR and stellar component density are calculated by averaging the GADGET2 particle properties on a 2D grid mapped on the Galactic plane. The peak values for possibly habitable stellar mass surface density move from 10 to 15 kpc cylindrical galactocentric distance in a 10-Gyr simulated time span. At 10 Gyr, the simulation results imply the following. Stellar particles that have spent almost all of their lifetime in habitable-friendly conditions typically reside at ˜16 kpc from the Galactic Centre and are ˜3 Gyr old. Stellar particles that have spent ≥90 per cent of their 4-5 Gyr long lifetime in habitable-friendly conditions are also predominantly found in the outskirts of the Galactic disc. Fewer than 1 per cent of these particles can be found at a typical Solar system galactocentric distance of 8-10 kpc. Our results imply that the evolution of an isolated spiral galaxy is likely to result in galactic civilizations emerging at the outskirts of the galactic disc around stellar hosts younger than the Sun.
N-body simulations of terrestrial planet formation under the influence of a hot Jupiter
Ogihara, Masahiro; Kobayashi, Hiroshi; Inutsuka, Shu-ichiro E-mail: ogihara@nagoya-u.jp
2014-06-01
We investigate the formation of multiple-planet systems in the presence of a hot Jupiter (HJ) using extended N-body simulations that are performed simultaneously with semianalytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of an HJ. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal mass; the largest planet constitutes more than 50% of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper, we have found a new physical mechanism of induced migration of the HJ, which is called a crowding-out. If the HJ opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the HJ remains. We also discuss angular momentum transfer between the planets and disk.
"Grandeur in this view of life": N-body simulation models of the Galactic habitable zone
NASA Astrophysics Data System (ADS)
Vukotić, B.; Steinhauser, D.; Martinez-Aviles, G.; Ćirković, M. M.; Micic, M.; Schindler, S.
2016-04-01
We present an isolated Milky Way-like simulation in GADGET2 N-body SPH code. The Galactic disk star formation rate (SFR) surface densities and stellar mass indicative of Solar neighbourhood are used as thresholds to model the distribution of stellar mass in life friendly environments. SFR and stellar component density are calculated averaging the GADGET2 particle properties on a 2D grid mapped on the Galactic plane. The peak values for possibly habitable stellar mass surface density move from 10 to 15 kpc cylindrical galactocentric distance in 10 Gyr simulated time span. At 10 Gyr the simulation results imply the following. Stellar particles which have spent almost all of their life time in habitable friendly conditions reside typically at ˜16 kpc from Galactic centre and are ˜3 Gyr old. Stellar particles that have spent ≥90% of their 4 - 5 Gyr long life time in habitable friendly conditions, are also predominantly found in the outskirts of the Galactic disk. Less then 1% of these particles can be found at a typical Solar system galactocentric distance of 8 - 10 kpc. Our results imply that the evolution of an isolated spiral galaxy is likely to result in galactic civilizations emerging at the outskirts of the galactic disk around stellar hosts younger than the Sun.
Direct Large-Scale N-Body Simulations of Planetesimal Dynamics
NASA Astrophysics Data System (ADS)
Richardson, Derek C.; Quinn, Thomas; Stadel, Joachim; Lake, George
2000-01-01
We describe a new direct numerical method for simulating planetesimal dynamics in which N˜10 6 or more bodies can be evolved simultaneously in three spatial dimensions over hundreds of dynamical times. This represents several orders of magnitude improvement in resolution over previous studies. The advance is made possible through modification of a stable and tested cosmological code optimized for massively parallel computers. However, owing to the excellent scalability and portability of the code, modest clusters of workstations can treat problems with N˜10 5 particles in a practical fashion. The code features algorithms for detection and resolution of collisions and takes into account the strong central force field and flattened Keplerian disk geometry of planetesimal systems. We demonstrate the range of problems that can be addressed by presenting simulations that illustrate oligarchic growth of protoplanets, planet formation in the presence of giant planet perturbations, the formation of the jovian moons, and orbital migration via planetesimal scattering. We also describe methods under development for increasing the timescale of the simulations by several orders of magnitude.
Acellular ostrich corneal stroma used as scaffold for construction of tissue-engineered cornea
Liu, Xian-Ning; Zhu, Xiu-Ping; Wu, Jie; Wu, Zheng-Jie; Yin, Yong; Xiao, Xiang-Hua; Su, Xin; Kong, Bin; Pan, Shi-Yin; Yang, Hua; Cheng, Yan; An, Na; Mi, Sheng-Li
2016-01-01
AIM To assess acellular ostrich corneal matrix used as a scaffold to reconstruct a damaged cornea. METHODS A hypertonic saline solution combined with a digestion method was used to decellularize the ostrich cornea. The microstructure of the acellular corneal matrix was observed by transmission electron microscopy (TEM) and hematoxylin and eosin (H&E) staining. The mechanical properties were detected by a rheometer and a tension machine. The acellular corneal matrix was also transplanted into a rabbit cornea and cytokeratin 3 was used to check the immune phenotype. RESULTS The microstructure and mechanical properties of the ostrich cornea were well preserved after the decellularization process. In vitro, the methyl thiazolyl tetrazolium results revealed that extracts of the acellular ostrich corneas (AOCs) had no inhibitory effects on the proliferation of the corneal epithelial or endothelial cells or on the keratocytes. The rabbit lamellar keratoplasty showed that the transplanted AOCs were transparent and completely incorporated into the host cornea while corneal turbidity and graft dissolution occurred in the acellular porcine cornea (APC) transplantation. The phenotype of the reconstructed cornea was similar to a normal rabbit cornea with a high expression of cytokeratin 3 in the superficial epithelial cell layer. CONCLUSION We first used AOCs as scaffolds to reconstruct damaged corneas. Compared with porcine corneas, the anatomical structures of ostrich corneas are closer to those of human corneas. In accordance with the principle that structure determines function, a xenograft lamellar keratoplasty also confirmed that the AOC transplantation generated a superior outcome compared to that of the APC graft. PMID:27158598
Effects of the decellularization method on the local stiffness of acellular lungs.
Melo, Esther; Garreta, Elena; Luque, Tomas; Cortiella, Joaquin; Nichols, Joan; Navajas, Daniel; Farré, Ramon
2014-05-01
Lung bioengineering, a novel approach to obtain organs potentially available for transplantation, is based on decellularizing donor lungs and seeding natural scaffolds with stem cells. Various physicochemical protocols have been used to decellularize lungs, and their performance has been evaluated in terms of efficient decellularization and matrix preservation. No data are available, however, on the effect of different decellularization procedures on the local stiffness of the acellular lung. This information is important since stem cells directly sense the rigidity of the local site they are engrafting to during recellularization, and it has been shown that substrate stiffness modulates cell fate into different phenotypes. The aim of this study was to assess the effects of the decellularization procedure on the inhomogeneous local stiffness of the acellular lung on five different sites: alveolar septa, alveolar junctions, pleura, and vessels' tunica intima and tunica adventitia. Local matrix stiffness was measured by computing Young's modulus with atomic force microscopy after decellularizing the lungs of 36 healthy rats (Sprague-Dawley, male, 250-300 g) with four different protocols with/without perfusion through the lung circulatory system and using two different detergents (sodium dodecyl sulfate [SDS] and 3-[(3-cholamidopropyl) dimethylammonio]-1-propanesulfonate [CHAPS]). The local stiffness of the acellular lung matrix significantly depended on the site within the matrix (p<0.001), ranging from ∼ 15 kPa at the alveolar septum to ∼ 60 kPa at the tunica intima. Acellular lung stiffness (p=0.003) depended significantly, albeit modestly, on the decellularization process. Whereas perfusion did not induce any significant differences in stiffness, the use of CHAPS resulted in a ∼ 35% reduction compared with SDS, the influence of the detergent being more important in the tunica intima. In conclusion, lung matrix stiffness is considerably inhomogeneous, and
Korossis, Sotirios A.; Wilshaw, Stacy-Paul; Jennings, Louise M; Fisher, John; Ingham, Eileen
2014-01-01
Currently available replacement heart valves all have limitations. This study aimed to produce and characterize an acellular, biocompatible porcine pulmonary root conduit for reconstruction of the right ventricular outflow tract e.g., during Ross procedure. A process for the decellularization of porcine pulmonary roots was developed incorporating trypsin treatment of the adventitial surface of the scraped pulmonary artery and sequential treatment with hypotonic Tris buffer (HTB; 10 mM Tris pH 8.0, 0.1% (w/v) EDTA, and 10 KIU aprotinin), 0.1% (w/v) sodium dodecyl sulfate in HTB, two cycles of DNase and RNase, and sterilization with 0.1% (v/v) peracetic acid. Histology confirmed an absence of cells and retention of the gross histoarchitecture. Immunohistochemistry further confirmed cell removal and partial retention of the extracellular matrix, but a loss of collagen type IV. DNA levels were reduced by more than 96% throughout all regions of the acellular tissue and no functional genes were detected using polymerase chain reaction. Total collagen levels were retained but there was a significant loss of glycosaminoglycans following decellularization. The biomechanical, hydrodynamic, and leaflet kinematics properties were minimally affected by the process. Both immunohistochemical labeling and antibody absorption assay confirmed a lack of α-gal epitopes in the acellular porcine pulmonary roots and in vitro biocompatibility studies indicated that acellular leaflets and pulmonary arteries were not cytotoxic. Overall the acellular porcine pulmonary roots have excellent potential for development of a tissue substitute for right ventricular outflow tract reconstruction e.g., during the Ross procedure. PMID:24786313
Large eddy simulation of the flow around bluff body with drag reduction device
NASA Astrophysics Data System (ADS)
Al-Anazi, Khalid Qaied
This thesis focuses on the use of LES to simulate the flow around elliptical bluff body with blunt trailing edge fitted with open base cavity. The main objective of this study is to determine the effects of the cavity on the drag of the body. A secondary but important objective is to demonstrate that LES can provide accurate representation of the flow around this bluff body. Moreover, LES results can complement the available experimental results in order to provide a much better understanding of the flow. The simulations were carried out at a Reynolds number of 2.6×104 based on the height of the body using Spalart-Allmaras RANS model while the LES were performed using Smagorinsky dynamic model. A grid-independence test was conducted using three grids which contain 0.85M, 1.3M and 1.7M cells, respectively. This test shows that the results are grid-independent. The LES results predicted the mean flow field in the near wake with good accuracy as compared to the experimental mean flow field obtained. The base pressure results show that the base pressure coefficient for the base model was around -0.56, which agrees well with the experimental results .By attaching the cavity, the base pressure has increased. The increase in base pressure coefficient was around 44% using 1/3 h cavity and this agrees well with the experimental measurements. The RANS predicted drag coefficient of 0.56 for the base model and 0.471 for the cavity model. This represents a difference of 8% for the base model and 34% for the cavity model when compared with experiment drag coefficients (0.61 for the base model and 0.35 for the cavity model). For the LES, the drag coefficient of the base model was around 0.65 (6.5% difference) and using the cavity, the drag coefficient was reduced to around 0.37 (5.74% difference). Details of the mean velocity components have been compared with experimental data at various locations in the wake region of the flow. Observation on the comparison between LES and
Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations
NASA Astrophysics Data System (ADS)
Kuechler, Erich R.; Giese, Timothy J.; York, Darrin M.
2015-12-01
Accurate modeling of the molecular environment is critical in condensed phase simulations of chemical reactions. Conventional quantum mechanical/molecular mechanical (QM/MM) simulations traditionally model non-electrostatic non-bonded interactions through an empirical Lennard-Jones (LJ) potential which, in violation of intuitive chemical principles, is bereft of any explicit coupling to an atom's local electronic structure. This oversight results in a model whereby short-ranged exchange-repulsion and long-ranged dispersion interactions are invariant to changes in the local atomic charge, leading to accuracy limitations for chemical reactions where significant atomic charge transfer can occur along the reaction coordinate. The present work presents a variational, charge-dependent exchange-repulsion and dispersion model, referred to as the charge-dependent exchange and dispersion (QXD) model, for hybrid QM/MM simulations. Analytic expressions for the energy and gradients are provided, as well as a description of the integration of the model into existing QM/MM frameworks, allowing QXD to replace traditional LJ interactions in simulations of reactive condensed phase systems. After initial validation against QM data, the method is demonstrated by capturing the solvation free energies of a series of small, chlorine-containing compounds that have varying charge on the chlorine atom. The model is further tested on the SN2 attack of a chloride anion on methylchloride. Results suggest that the QXD model, unlike the traditional LJ model, is able to simultaneously obtain accurate solvation free energies for a range of compounds while at the same time closely reproducing the experimental reaction free energy barrier. The QXD interaction model allows explicit coupling of atomic charge with many-body exchange and dispersion interactions that are related to atomic size and provides a more accurate and robust representation of non-electrostatic non-bonded QM/MM interactions.
Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations.
Kuechler, Erich R; Giese, Timothy J; York, Darrin M
2015-12-21
Accurate modeling of the molecular environment is critical in condensed phase simulations of chemical reactions. Conventional quantum mechanical/molecular mechanical (QM/MM) simulations traditionally model non-electrostatic non-bonded interactions through an empirical Lennard-Jones (LJ) potential which, in violation of intuitive chemical principles, is bereft of any explicit coupling to an atom's local electronic structure. This oversight results in a model whereby short-ranged exchange-repulsion and long-ranged dispersion interactions are invariant to changes in the local atomic charge, leading to accuracy limitations for chemical reactions where significant atomic charge transfer can occur along the reaction coordinate. The present work presents a variational, charge-dependent exchange-repulsion and dispersion model, referred to as the charge-dependent exchange and dispersion (QXD) model, for hybrid QM/MM simulations. Analytic expressions for the energy and gradients are provided, as well as a description of the integration of the model into existing QM/MM frameworks, allowing QXD to replace traditional LJ interactions in simulations of reactive condensed phase systems. After initial validation against QM data, the method is demonstrated by capturing the solvation free energies of a series of small, chlorine-containing compounds that have varying charge on the chlorine atom. The model is further tested on the SN2 attack of a chloride anion on methylchloride. Results suggest that the QXD model, unlike the traditional LJ model, is able to simultaneously obtain accurate solvation free energies for a range of compounds while at the same time closely reproducing the experimental reaction free energy barrier. The QXD interaction model allows explicit coupling of atomic charge with many-body exchange and dispersion interactions that are related to atomic size and provides a more accurate and robust representation of non-electrostatic non-bonded QM/MM interactions. PMID
Michikoshi, Shugo; Kokubo, Eiichiro; Inutsuka, Shu-ichiro E-mail: kokubo@th.nao.ac.j
2009-10-01
The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local N-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of nonaxisymmetric structures; the creation of planetesimal seeds; and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to L {sup 3/2} {sub y}, where L{sub y} is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of L{sub x} , where L{sub x} is the size of the computational domain in the radial direction if L{sub x} is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.
Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations
Kuechler, Erich R.; Giese, Timothy J.; York, Darrin M.
2015-12-21
Accurate modeling of the molecular environment is critical in condensed phase simulations of chemical reactions. Conventional quantum mechanical/molecular mechanical (QM/MM) simulations traditionally model non-electrostatic non-bonded interactions through an empirical Lennard-Jones (LJ) potential which, in violation of intuitive chemical principles, is bereft of any explicit coupling to an atom’s local electronic structure. This oversight results in a model whereby short-ranged exchange-repulsion and long-ranged dispersion interactions are invariant to changes in the local atomic charge, leading to accuracy limitations for chemical reactions where significant atomic charge transfer can occur along the reaction coordinate. The present work presents a variational, charge-dependent exchange-repulsion and dispersion model, referred to as the charge-dependent exchange and dispersion (QXD) model, for hybrid QM/MM simulations. Analytic expressions for the energy and gradients are provided, as well as a description of the integration of the model into existing QM/MM frameworks, allowing QXD to replace traditional LJ interactions in simulations of reactive condensed phase systems. After initial validation against QM data, the method is demonstrated by capturing the solvation free energies of a series of small, chlorine-containing compounds that have varying charge on the chlorine atom. The model is further tested on the S{sub N}2 attack of a chloride anion on methylchloride. Results suggest that the QXD model, unlike the traditional LJ model, is able to simultaneously obtain accurate solvation free energies for a range of compounds while at the same time closely reproducing the experimental reaction free energy barrier. The QXD interaction model allows explicit coupling of atomic charge with many-body exchange and dispersion interactions that are related to atomic size and provides a more accurate and robust representation of non-electrostatic non-bonded QM
Digital quantum simulation of many-body non-Markovian dynamics
NASA Astrophysics Data System (ADS)
Sweke, R.; Sanz, M.; Sinayskiy, I.; Petruccione, F.; Solano, E.
2016-08-01
We present an algorithmic method for the digital quantum simulation of many-body locally indivisible non-Markovian open quantum systems. It consists of two parts: first, a Suzuki-Lie-Trotter decomposition of the global system propagator into the product of subsystem propagators, which may not be quantum channels, and second, an algorithmic procedure for the implementation of the subsystem propagators through unitary operations and measurements on a dilated space. By providing rigorous error bounds for the relevant Suzuki-Lie-Trotter decomposition, we are able to analyze the efficiency of the method, and connect it with an appropriate measure of the local indivisibility of the system. In light of our analysis, the proposed method is expected to be experimentally achievable for a variety of interesting cases.
CIRCUMBINARY PLANET FORMATION IN THE KEPLER-16 SYSTEM. I. N-BODY SIMULATIONS
Meschiari, Stefano
2012-06-10
The recently discovered circumbinary planets (Kepler-16 b, Kepler 34-b, Kepler 35-b) represent the first direct evidence of the viability of planet formation in circumbinary orbits. We report on the results of N-body simulations investigating planetesimal accretion in the Kepler-16 b system, focusing on the range of impact velocities under the influence of both stars' gravitational perturbation and friction from a putative protoplanetary disk. Our results show that planet formation might be effectively inhibited for a large range in semimajor axis (1.75 {approx}< a{sub P} {approx}< 4 AU), suggesting that the planetary core must have either migrated from outside 4 AU or formed in situ very close to its current location.
Measured force on elongated bodies in a simulated low-Earth orbit environment
Maldonado, C. A.; Ketsdever, A. D.; Gimelshein, S. F.
2014-12-09
An overview of the development of a magnetically filtered atomic oxygen plasma source and the application of the source to study low-Earth orbit drag on elongated bodies is presented. Plasma diagnostics show that the magnetic filter plasma source produces atomic oxygen ions (O{sup +}) with streaming energies equivalent to the relative orbital environment of approximately 5eV and can supply the appropriate density for LEO simulation. Previous research has demonstrated that momentum transfer between ions and metal surfaces is equivalent to the momentum transfer expected for neutral molecules with similar energy, due to charge exchange occurring prior to momentum transfer. Total drag measurements of aluminum cuboid geometries of varying length to diameter ratios immersed in the extracted plasma plume are presented as a function of streaming ion energy.
REBOUNDx: A library for adding additional effects to N-body simulations
NASA Astrophysics Data System (ADS)
Tamayo, Daniel; Rein, Hanno; Shi, Pengshuai
2016-05-01
Many astrophysical applications involve additional perturbations beyond point-source gravity. We have recently developed REBOUNDx, a library for adding such effects in numerical simulations with the open-source N-body package REBOUND. Various implementations have different numerical properties that in general depend on the underlying integrator employed. In particular, I will discuss adding velocity-dependent/dissipative effects to widely used symplectic integrators, and how one can estimate the introduced numerical errors using the operator-splitting formalism traditionally applied to symplectic integrators. Finally, I will demonstrate how to use the code, and how the Python wrapper we have developed for REBOUND/REBOUNDx makes it easy to interactively leverage powerful analysis, visualization and parallelization libraries.
Measured force on elongated bodies in a simulated low-Earth orbit environment
NASA Astrophysics Data System (ADS)
Maldonado, C. A.; Ketsdever, A. D.; Gimelshein, S. F.
2014-12-01
An overview of the development of a magnetically filtered atomic oxygen plasma source and the application of the source to study low-Earth orbit drag on elongated bodies is presented. Plasma diagnostics show that the magnetic filter plasma source produces atomic oxygen ions (O+) with streaming energies equivalent to the relative orbital environment of approximately 5eV and can supply the appropriate density for LEO simulation. Previous research has demonstrated that momentum transfer between ions and metal surfaces is equivalent to the momentum transfer expected for neutral molecules with similar energy, due to charge exchange occurring prior to momentum transfer. Total drag measurements of aluminum cuboid geometries of varying length to diameter ratios immersed in the extracted plasma plume are presented as a function of streaming ion energy.
Unsteady Navier-Stokes simulation of the canard-wing-body ramp motion
NASA Technical Reports Server (NTRS)
Tu, Eugene L.; Obayashi, Shigeru; Guruswamy, Guru P.
1993-01-01
A time-accurate thin-layer Navier-Stokes simulation of the unsteady flowfield is performed for a typical canard-wing-body configuration undergoing ramp motions. The computations are made at a transonic Mach number of 0.90 and for ramp angles from 0 to 15 degrees. Accuracy is determined by comparisons with steady-state experimental data and with spatial and time-step refinement studies. During the ramp motion, the computational results show improved dynamic lift performance and a strong canard-wing interaction for the canard-on configuration. Formation of the canard leading-edge vortex is inhibited in the early stages of the ramp motion. An analysis performed on the transient flowfield after the ramp motion ends shows that the canard vortex rapidly gains strength and vortex breakdown eventually occurs. These characteristics of the canard vortex have significant influences on wing performance.
The Electrochemical Behavior of TiN/316LSS Material in Simulated Body Fluid Solution.
Thanh, Dinh Thi Mai; Pham, Thi Nam; Huong, Ho Thu; Phuong, Nguyen Thu; Hang, To Thi Xuan; Vy, Uong Van; Hoang, Thai
2015-05-01
We report on the fabrication and the electrochemical behavior of TiN film on the 316L stainless steel (316LSS) material in simulated body fluid (SBF) solution for implant application. The characterization results indicate that the coated TiN is completely crystalline with (111) crystal orientation. Electrochemical results of 316LSS and TiN/316LSS material after 21 days of immersion in SBF show that the durability of the TiN/316LSS is much higher than that of 316LSS, which registers a very low corrosion current density (about tens of nA cm(-2)). The formation of hydroxyapatite on the surface of the TiN/316LSS is also confirmed by SEM, EDX, X-ray and IR spectroscopy. PMID:26505019
Joyce, M.; Marcos, B.
2007-11-15
We apply a recently developed perturbative formalism which describes the evolution under their self-gravity of particles displaced from a perfect lattice to quantify precisely, up to shell crossing, the effects of discreteness in dissipationless cosmological N-body simulations. We give simple expressions, explicitly dependent on the particle density, for the evolution of power in each mode as a function of redshift. For typical starting redshifts the effect of finite particle number is to slow down slightly the growth of power compared to that in the fluid limit (e.g., by about 10% at half the Nyquist frequency), and to induce also dispersion in the growth as a function of direction at a comparable level. In the limit that the initial redshift tends to infinity, at fixed particle density, the evolution in fact diverges from that in the fluid limit (described by the Zeldovich approximation). Contrary to widely held belief, this means that a simulation started at a redshift much higher than the redshift of shell crossing actually gives a worse, rather than a better, result. We also study how these effects are modified when there is a small-scale regularization of the gravitational force. We show that such a smoothing may reduce the anisotropy of the discreteness effects, but it then increases their average effect. This behavior illustrates the fact that the discreteness effects described here are distinct from those usually considered in this context, due to two-body collisions. Indeed the characteristic time for divergence from the collisionless limit is proportional to N{sup 2/3}, rather than N/logN in the latter case.
Sub-discretized surface model with application to contact mechanics in multi-body simulation
Johnson, S; Williams, J
2008-02-28
The mechanics of contact between rough and imperfectly spherical adhesive powder grains are often complicated by a variety of factors, including several which vary over sub-grain length scales. These include several traction factors that vary spatially over the surface of the individual grains, including high energy electron and acceptor sites (electrostatic), hydrophobic and hydrophilic sites (electrostatic and capillary), surface energy (general adhesion), geometry (van der Waals and mechanical), and elasto-plastic deformation (mechanical). For mechanical deformation and reaction, coupled motions, such as twisting with bending and sliding, as well as surface roughness add an asymmetry to the contact force which invalidates assumptions for popular models of contact, such as the Hertzian and its derivatives, for the non-adhesive case, and the JKR and DMT models for adhesive contacts. Though several contact laws have been offered to ameliorate these drawbacks, they are often constrained to particular loading paths (most often normal loading) and are relatively complicated for computational implementation. This paper offers a simple and general computational method for augmenting contact law predictions in multi-body simulations through characterization of the contact surfaces using a hierarchically-defined surface sub-discretization. For the case of adhesive contact between powder grains in low stress regimes, this technique can allow a variety of existing contact laws to be resolved across scales, allowing for moments and torques about the contact area as well as normal and tangential tractions to be resolved. This is especially useful for multi-body simulation applications where the modeler desires statistical distributions and calibration for parameters in contact laws commonly used for resolving near-surface contact mechanics. The approach is verified against analytical results for the case of rough, elastic spheres.
Castro, Nathan J; Tan, Wilhelmina Nanrui; Shen, Charlie; Zhang, Lijie Grace
2016-07-01
Osseous tissue defects caused by trauma present a common clinical problem. Although traditional clinical procedures have been successfully employed, several limitations persist with regards to insufficient donor tissue, disease transmission, and inadequate host-implant integration. Therefore, this work aims to address current limitations regarding inadequate host tissue integration through the use of a novel elastomeric material for three-dimensional (3D) printing biomimetic and bioactive scaffolds. A novel thermoplastic polyurethane-based elastomeric composite filament (Gel-Lay) was used to manufacture porous scaffolds. In an effort to render the scaffolds more bioactive, the flexible scaffolds were subsequently incubated in simulated body fluid at various time points and evaluated for enhanced mechanical properties along with the effects on cell adhesion, proliferation, and 3-week osteogenesis. This work is the first reported use of a novel class of flexible elastomeric materials for the manufacture of 3D printed bioactive scaffold fabrication allowing efficient and effective nucleation of hydroxyapatite (HA) leading to increased nanoscale surface roughness while retaining the bulk geometry of the predesigned structure. Scaffolds with interconnected microfibrous filaments of ∼260 μm were created and nucleated in simulated body fluid that facilitated cell adhesion and spreading after only 24 h in culture. The porous structure further allowed efficient nucleation, exchange of nutrients, and metabolic waste removal during new tissue formation. Through the incorporation of osteoconductive HA, human fetal osteoblast adhesion and differentiation were greatly enhanced thus setting the tone for further exploration of this novel material for biomedical and tissue regenerative applications. PMID:27298115
Modeling flow around bluff bodies and predicting urban dispersion using large eddy simulation.
Tseng, Yu-Heng; Meneveau, Charles; Parlange, Marc B
2006-04-15
Modeling air pollutant transport and dispersion in urban environments is especially challenging due to complex ground topography. In this study, we describe a large eddy simulation (LES) tool including a new dynamic subgrid closure and boundary treatment to model urban dispersion problems. The numerical model is developed, validated, and extended to a realistic urban layout. In such applications fairly coarse grids must be used in which each building can be represented using relatively few grid-points only. By carrying out LES of flow around a square cylinder and of flow over surface-mounted cubes, the coarsest resolution required to resolve the bluff body's cross section while still producing meaningful results is established. Specifically, we perform grid refinement studies showing that at least 6-8 grid points across the bluff body are required for reasonable results. The performance of several subgrid models is also compared. Although effects of the subgrid models on the mean flow are found to be small, dynamic Lagrangian models give a physically more realistic subgrid-scale (SGS) viscosity field. When scale-dependence is taken into consideration, these models lead to more realistic resolved fluctuating velocities and spectra. These results set the minimum grid resolution and subgrid model requirements needed to apply LES in simulations of neutral atmospheric boundary layer flow and scalar transport over a realistic urban geometry. The results also illustrate the advantages of LES over traditional modeling approaches, particularly its ability to take into account the complex boundary details and the unsteady nature of atmospheric boundary layer flow. Thus LES can be used to evaluate probabilities of extreme events (such as probabilities of exceeding threshold pollutant concentrations). Some comments about computer resources required for LES are also included. PMID:16683605
Study of variation in human upper body parameters and motion for use in robotics based simulation.
Lura, Derek J; Carey, Stephanie L; Dubey, Rajiv V
2013-01-01
This paper reviews the variations in human upper body motion of subjects completing activities of daily living. This study was completed to serve as a reference to evaluate the quality of simulated of human motion. In this paper we define the variation in motion as the variation in subjects' parameters (link lengths), joint angles, and hand positions, for a given task. All of these variations are related by forward kinematic equations. Motion data from eight healthy right hand dominant adults performing three activities of daily living (brushing hair, drinking from a cup, and opening a door) were collected using an eight camera Vicon motion analysis system. Subject parameters were calculated using relative positions of functional joint center locations between segments. Joint angles were calculated by Euler angle rotations between body segments. Hand position was defined as the origin of the hand frame relative to the pelvis frame. The variance of recorded human motion was analyzed based on the standard deviations of subject parameters, joint angles, and hand positions. Variances in joint angles were found to be similar in magnitude to root mean squared error of kinematics based motion simulation. To evaluate the relative variance, the forward kinematic solutions of the trials were found after removing subject parameter variance and reducing joint angle variance. The variance in the forward kinematic solution was then compared to the recorded hand position variance. Reductions in subject parameter and joint angle variance produced a proportionally much smaller reduction in the calculated hand position variance. Using the average instead of individual subject parameters had only a small impact on hand position variance. Modifying joint angles to reduce variance had a greater impact on the calculated hand position variance than using average subject parameters, but was still a relatively small change. Future work will focus on using these results to create formalized
A Multirate Variable-timestep Algorithm for N-body Solar System Simulations with Collisions
NASA Astrophysics Data System (ADS)
Sharp, P. W.; Newman, W. I.
2016-03-01
We present and analyze the performance of a new algorithm for performing accurate simulations of the solar system when collisions between massive bodies and test particles are permitted. The orbital motion of all bodies at all times is integrated using a high-order variable-timestep explicit Runge-Kutta Nyström (ERKN) method. The variation in the timestep ensures that the orbital motion of test particles on eccentric orbits or close to the Sun is calculated accurately. The test particles are divided into groups and each group is integrated using a different sequence of timesteps, giving a multirate algorithm. The ERKN method uses a high-order continuous approximation to the position and velocity when checking for collisions across a step. We give a summary of the extensive testing of our algorithm. In our largest simulation—that of the Sun, the planets Earth to Neptune and 100,000 test particles over 100 million years—the relative error in the energy after 100 million years was of the order of 10-11.
Small-body deflection techniques using spacecraft: Techniques in simulating the fate of ejecta
NASA Astrophysics Data System (ADS)
Schwartz, Stephen R.; Yu, Yang; Michel, Patrick; Jutzi, Martin
2016-04-01
We define a set of procedures to numerically study the fate of ejecta produced by the impact of an artificial projectile with the aim of deflecting an asteroid. Here we develop a simplified, idealized model of impact conditions that can be adapted to fit the details of specific deflection-test scenarios, such as what is being proposed for the AIDA project. Ongoing studies based upon the methodology described here can be used to inform observational strategies and safety conditions for an observing spacecraft. To account for ejecta evolution, the numerical strategies we are employing are varied and include a large N-Body component, a smoothed-particle hydrodynamics (SPH) component, and an application of impactor scaling laws. Simulations that use SPH-derived initial conditions show high-speed ejecta escaping at low angles of inclination, and very slowly moving ejecta lofting off the surface at higher inclination angles, some of which reimpacts the small-body surface. We are currently investigating the realism of this and other models' behaviors. Next steps will include the addition of solar perturbations to the model and applying the protocol developed here directly to specific potential mission concepts such as the proposed AIDA scenario.
NASA Astrophysics Data System (ADS)
Kuhl, J. M.; Desjardin, P. E.
2012-01-01
Two-dimensional, fully coupled direct numerical simulations (DNS) are conducted to examine the local energy dynamics of a flexible cantilevered plate in the wake of a two-dimensional circular cylinder. The motion of the cantilevered plate is described using a finite element formulation and a fully compressible, finite volume Navier Stokes solver is used to compute the flow field. A sharp interface level set method is employed in conjunction with a ghost fluid method to describe the immersed boundaries of the bluff body and flexible plate. DNS is first conducted to validate the numerical methodology and compared with previous studies of flexible cantilevered plates and flow over bluff bodies; excellent agreement with previous results is observed. A newly defined power production/loss geometry metric is introduced based on surface curvature and plate velocity. The metric is found to be useful for determining which sections of the plate will produce energy based on curvature and deflection rate. Scatter plots and probability measures are presented showing a high correlation between the direction of energy transfer (i.e., to or from the plate) and the sign of the newly defined curvature-deflection-rate metric. The findings from this study suggest that a simple local geometry/kinematic based metric can be devised to aid in the development and design of flexible wind energy harvesting flutter mills.