Dynamic coupling between fluid flow and vein growth in fractures: a 3D numerical model
NASA Astrophysics Data System (ADS)
Schwarz, J.-O.; Enzmann, F.
2012-04-01
Fluid flow is one of the main mass transport mechanisms in the Earth's crust and abundant mineral vein networks are important indicators for fluid flow and fluid rock interaction. These systems are dynamic and part of the so called RTM processes (reaction-transport-mechanics). Understanding of mineral vein systems requires coupling of these processes. Here we present a conceptional model for dynamic vein growth of syntaxial, posttectonic veins generated by advective fluid flow and show first results of a numerical model for this scenario. Vein generation requires three processes to occur: (i) fracture generation by mechanical stress e.g. hydro-fracturing, (ii) flow of a supersaturated fluid on that fracture and (iii) crystallization of phase(s) on or in the fracture. 3D synthetic fractures are generated with the SynFrac code (Ogilvie, et al. 2006). Subsequently solutions of the Navier-Stokes equation for this fracture are computed by a computational fluid dynamics code called GeoDict (Wiegmann 2007). Transport (advective and diffusive) of chemical species to growth sites in the fracture and vein growth are computed by a self-written MATLAB script. The numerical model discretizes the wall rock and fracture geometry by volumetric pixels (voxels). Based on this representation, the model computes the three basic functions for vein generation: (a) nucleation, (b) fluid flow with transport of chemical species and (c) growth. The following conditions were chosen for these three modules. Nucleation is heterogeneous and occurs instantaneously at the wall rock/fracture interface. Advective and diffusive flow of a supersaturated fluid and related transport of chemical species occurs according to the computed fluid flow field by GeoDict. Concentration of chemical species at the inflow is constant, representing external fluid buffering. Changes/decrease in the concentration of chemical species occurs only due to vein growth. Growth of nuclei is limited either by transport of
User's manual for PELE3D: a computer code for three-dimensional incompressible fluid dynamics
McMaster, W H
1982-05-07
The PELE3D code is a three-dimensional semi-implicit Eulerian hydrodynamics computer program for the solution of incompressible fluid flow coupled to a structure. The fluid and coupling algorithms have been adapted from the previously developed two-dimensional code PELE-IC. The PELE3D code is written in both plane and cylindrical coordinates. The coupling algorithm is general enough to handle a variety of structural shapes. The free surface algorithm is able to accommodate a top surface and several independent bubbles. The code is in a developmental status since all the intended options have not been fully implemented and tested. Development of this code ended in 1980 upon termination of the contract with the Nuclear Regulatory Commission.
Variable Quality Compression of Fluid Dynamical Data Sets Using a 3D DCT Technique
NASA Astrophysics Data System (ADS)
Loddoch, A.; Schmalzl, J.
2005-12-01
In this work we present a data compression scheme that is especially suited for the compression of data sets resulting from computational fluid dynamics (CFD). By adopting the concept of the JPEG compression standard and extending the approach of Schmalzl (Schmalzl, J. Using standard image compression algorithms to store data from computational fluid dynamics. Computers and Geosciences, 29, 10211031, 2003) we employ a three-dimensional discrete cosine transform of the data. The resulting frequency components are rearranged, quantized and finally stored using Huffman-encoding and standard variable length integer codes. The compression ratio and also the introduced loss of accuracy can be adjusted by means of two compression parameters to give the desired compression profile. Using the proposed technique compression ratios of more than 60:1 are possible with an mean error of the compressed data of less than 0.1%.
NASA Astrophysics Data System (ADS)
Stockton, Gregory R.
2011-05-01
Over the last 10 years, very large government, military, and commercial computer and data center operators have spent millions of dollars trying to optimally cool data centers as each rack has begun to consume as much as 10 times more power than just a few years ago. In fact, the maximum amount of data computation in a computer center is becoming limited by the amount of available power, space and cooling capacity at some data centers. Tens of millions of dollars and megawatts of power are being annually spent to keep data centers cool. The cooling and air flows dynamically change away from any predicted 3-D computational fluid dynamic modeling during construction and as time goes by, and the efficiency and effectiveness of the actual cooling rapidly departs even farther from predicted models. By using 3-D infrared (IR) thermal mapping and other techniques to calibrate and refine the computational fluid dynamic modeling and make appropriate corrections and repairs, the required power for data centers can be dramatically reduced which reduces costs and also improves reliability.
NASA Astrophysics Data System (ADS)
Tang, H. S.; Qu, K.; Wu, X. G.
2014-09-01
It is now becoming important to develop our capabilities to simulate coastal ocean flows involved with distinct physical phenomena occurring at a vast range of spatial and temporal scales. This paper presents a hybrid modeling system for such simulation. The system consists of a fully three dimensional (3D) fluid dynamics model and a geophysical fluid dynamics model, which couple with each other in two-way and march in time simultaneously. Particularly, in the hybrid system, the solver for incompressible flow on overset meshes (SIFOM) resolves fully 3D small-scale local flow phenomena, while the unstructured grid finite volume coastal ocean model (FVCOM) captures large-scale background flows. The integration of the two models are realized via domain decomposition implemented with an overset grid method. Numerical experiments on performance of the system in resolving flow patterns and solution convergence rate show that the SIFOM-FVCOM system works as intended, and its solutions compare reasonably with data obtained with measurements and other computational approaches. Its unparalleled capabilities to predict multiphysics and multiscale phenomena with high-fidelity are demonstrated by three typical applications that are beyond the reach of other currently existing models. It is anticipated that the SIFOM-FVCOM system will serve as a new platform to study many emerging coastal ocean problems.
Schmitt, Vivien; Dufresne, Matthieu; Vazquez, Jose; Fischer, Martin; Morin, Antoine
2014-01-01
The aim of this study is to investigate the use of computational fluid dynamics (CFD) to predict the solid separation efficiency of a hydrodynamic separator. The numerical difficulty concerns the discretization of the geometry to simulate both the global behavior and the local phenomena that occur near the screen. In this context, a CFD multiscale approach was used: a global model (at the scale of the device) is used to observe the hydrodynamic behavior within the device; a local model (portion of the screen) is used to determine the local phenomena that occur near the screen. The Eulerian-Lagrangian approach was used to model the particle trajectories in both models. The global model shows the influence of the particles' characteristics on the trapping efficiency. A high density favors the sedimentation. In contrast, particles with small densities (1,040 kg/m(3)) are steered by the hydrodynamic behavior and can potentially be trapped by the separator. The use of the local model allows us to observe the particle trajectories near the screen. A comparison between two types of screens (perforated plate vs expanded metal) highlights the turbulent effects created by the shape of the screen. PMID:24622557
Jungreuthmayer, Christian; Steppert, Petra; Sekot, Gerhard; Zankel, Armin; Reingruber, Herbert; Zanghellini, Jürgen; Jungbauer, Alois
2015-12-18
Polymethacrylate-based monoliths have excellent flow properties. Flow in the wide channel interconnected with narrow channels is theoretically assumed to account for favorable permeability. Monoliths were cut into 898 slices in 50nm distances and visualized by serial block face scanning electron microscopy (SBEM). A 3D structure was reconstructed and used for the calculation of flow profiles within the monolith and for calculation of pressure drop and permeability by computational fluid dynamics (CFD). The calculated and measured permeabilities showed good agreement. Small channels clearly flowed into wide and wide into small channels in a repetitive manner which supported the hypothesis describing the favorable flow properties of these materials. This alternating property is also reflected in the streamline velocity which fluctuated. These findings were corroborated by artificial monoliths which were composed of regular (interconnected) cells where narrow cells followed wide cells. In the real monolith and the artificial monoliths with interconnected flow channels similar velocity fluctuations could be observed. A two phase flow simulation showed a lateral velocity component, which may contribute to the transport of molecules to the monolith wall. Our study showed that the interconnection of small and wide pores is responsible for the excellent pressure flow properties. This study is also a guide for further design of continuous porous materials to achieve good flow properties. PMID:26615711
Saho, Tatsunori; Onishi, Hideo
2015-07-01
In this study, we evaluated hemodynamics using simulated models and determined how cerebral aneurysms develop in simulated and patient-specific models based on medical images. Computational fluid dynamics (CFD) was analyzed by use of OpenFOAM software. Flow velocity, stream line, and wall shear stress (WSS) were evaluated in a simulated model aneurysm with known geometry and in a three-dimensional angiographic model. The ratio of WSS at the aneurysm compared with that at the basilar artery was 1:10 in simulated model aneurysms with a diameter of 10 mm and 1:18 in the angiographic model, indicating similar tendencies. Vortex flow occurred in both model aneurysms, and the WSS decreased in larger model aneurysms. The angiographic model provided accurate CFD information, and the tendencies of simulated and angiographic models were similar. These findings indicate that hemodynamic effects are involved in the development of aneurysms. PMID:25911446
NASA Technical Reports Server (NTRS)
Bartels, Robert E.
2012-01-01
This paper presents the implementation of gust modeling capability in the CFD code FUN3D. The gust capability is verified by computing the response of an airfoil to a sharp edged gust. This result is compared with the theoretical result. The present simulations will be compared with other CFD gust simulations. This paper also serves as a users manual for FUN3D gust analyses using a variety of gust profiles. Finally, the development of an Auto-Regressive Moving-Average (ARMA) reduced order gust model using a gust with a Gaussian profile in the FUN3D code is presented. ARMA simulated results of a sequence of one-minus-cosine gusts is shown to compare well with the same gust profile computed with FUN3D. Proper Orthogonal Decomposition (POD) is combined with the ARMA modeling technique to predict the time varying pressure coefficient increment distribution due to a novel gust profile. The aeroelastic response of a pitch/plunge airfoil to a gust environment is computed with a reduced order model, and compared with a direct simulation of the system in the FUN3D code. The two results are found to agree very well.
NASA Astrophysics Data System (ADS)
Li, Qiang; Yu, Guichang; Liu, Shulian; Zheng, Shuiying
2012-09-01
Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simplified physical model and classic Reynolds equation are always applied. While the application of the general computational fluid dynamics (CFD)-fluid structure interaction (FSI) techniques is more beneficial for analysis of the fluid field in a journal bearing when more detailed solutions are needed. This paper deals with the quasi-coupling calculation of transient fluid dynamics of oil film in journal bearings and rotor dynamics with CFD-FSI techniques. The fluid dynamics of oil film is calculated by applying the so-called "dynamic mesh" technique. A new mesh movement approach is presented while the dynamic mesh models provided by FLUENT are not suitable for the transient oil flow in journal bearings. The proposed mesh movement approach is based on the structured mesh. When the journal moves, the movement distance of every grid in the flow field of bearing can be calculated, and then the update of the volume mesh can be handled automatically by user defined function (UDF). The journal displacement at each time step is obtained by solving the moving equations of the rotor-bearing system under the known oil film force condition. A case study is carried out to calculate the locus of the journal center and pressure distribution of the journal in order to prove the feasibility of this method. The calculating results indicate that the proposed method can predict the transient flow field of a journal bearing in a rotor-bearing system where more realistic models are involved. The presented calculation method provides a basis for studying the nonlinear dynamic behavior of a general rotor-bearing system.
Neidlin, Michael; Steinseifer, Ulrich; Kaufmann, Tim A S
2014-06-01
Neurological complication often occurs during cardiopulmonary bypass (CPB). One of the main causes is hypoperfusion of the cerebral tissue affected by the position of the cannula tip and diminished cerebral autoregulation (CA). Recently, a lumped parameter approach could describe the baroreflex, one of the main mechanisms of cerebral autoregulation, in a computational fluid dynamics (CFD) study of CPB. However, the cerebral blood flow (CBF) was overestimated and the physiological meaning of the variables and their impact on the model was unknown. In this study, we use a 0-D control circuit representation of the Baroreflex mechanism, to assess the parameters with respect to their physiological meaning and their influence on CBF. Afterwards the parameters are transferred to 3D-CFD and the static and dynamic behavior of cerebral autoregulation is investigated. The parameters of the baroreflex mechanism can reproduce normotensive, hypertensive and impaired autoregulation behavior. Further on, the proposed model can mimic the effects of anesthetic agents and other factors controlling dynamic CA. The CFD simulations deliver similar results of static and dynamic CBF as the 0-D control circuit. This study shows the feasibility of a multiscale 0-D/3-D approach to include patient-specific cerebral autoregulation into CFD studies. PMID:24746017
Schwall, James R.; Karim, Naeem U.; Thakkar, Jivan G.; Taylor, Creed; Schulz, Terry; Wright, Richard F.
2006-07-01
The AP1000 is an 1100 MWe advanced nuclear power plant that uses passive safety features to enhance plant safety and to provide significant and measurable improvements in plant simplification, reliability, investment protection and plant costs. The AP1000 received final design approval from the US-NRC in 2004. The AP1000 design is based on the AP600 design that received final design approval in 1999. Wherever possible, the AP1000 plant configuration and layout was kept the same as AP600 to take advantage of the maturity of the design and to minimize new design efforts. As a result, the two-loop configuration was maintained for AP1000, and the containment vessel diameter was kept the same. It was determined that this significant power up-rate was well within the capability of the passive safety features, and that the safety margins for AP1000 were greater than those of operating PWRs. A key feature of the passive core cooling system is the passive residual heat removal heat exchanger (PRHR HX) that provides decay heat removal for postulated LOCA and non-LOCA events. The PRHR HX is a C-tube heat exchanger located in the in-containment refueling water storage tank (IRWST) above the core promoting natural circulation heat removal between the reactor cooling system and the tank. Component testing was performed for the AP600 PRHR HX to determine the heat transfer characteristics and to develop correlations to be used for the AP1000 safety analysis codes. The data from these tests were confirmed by subsequent integral tests at three separate facilities including the ROSA facility in Japan. Owing to the importance of this component, an independent analysis has been performed using the ATHOS-based computational fluid dynamics computer code PRHRCFD. Two separate models of the PRHR HX and IRWST have been developed representing the ROSA test geometry and the AP1000 plant geometry. Confirmation of the ROSA test results were used to validate PRHRCFD, and the AP1000 plant model
Sarrami-Foroushani, Ali; Nasr Esfahany, Mohsen; Nasiraei Moghaddam, Abbas; Saligheh Rad, Hamidreza; Firouznia, Kavous; Shakiba, Madjid; Ghanaati, Hossein; Wilkinson, Iain David; Frangi, Alejandro Federico
2015-01-01
Background: Understanding hemodynamic environment in vessels is important for realizing the mechanisms leading to vascular pathologies. Objectives: Three-dimensional velocity vector field in carotid bifurcation is visualized using TR 3D phase-contrast magnetic resonance imaging (TR 3D PC MRI) and computational fluid dynamics (CFD). This study aimed to present a qualitative and quantitative comparison of the velocity vector field obtained by each technique. Subjects and Methods: MR imaging was performed on a 30-year old male normal subject. TR 3D PC MRI was performed on a 3 T scanner to measure velocity in carotid bifurcation. 3D anatomical model for CFD was created using images obtained from time-of-flight MR angiography. Velocity vector field in carotid bifurcation was predicted using CFD and PC MRI techniques. A statistical analysis was performed to assess the agreement between the two methods. Results: Although the main flow patterns were the same for the both techniques, CFD showed a greater resolution in mapping the secondary and circulating flows. Overall root mean square (RMS) errors for all the corresponding data points in PC MRI and CFD were 14.27% in peak systole and 12.91% in end diastole relative to maximum velocity measured at each cardiac phase. Bland-Altman plots showed a very good agreement between the two techniques. However, this study was not aimed to validate any of methods, instead, the consistency was assessed to accentuate the similarities and differences between Time-resolved PC MRI and CFD. Conclusion: Both techniques provided quantitatively consistent results of in vivo velocity vector fields in right internal carotid artery (RCA). PC MRI represented a good estimation of main flow patterns inside the vasculature, which seems to be acceptable for clinical use. However, limitations of each technique should be considered while interpreting results. PMID:26793288
NASA Astrophysics Data System (ADS)
Meisenheimer, D.; Brueck, C. L.; Wildenschild, D.
2015-12-01
X-ray microtomography imaging of fluid-fluid interfaces in three-dimensional porous media allows for the testing of thermodynamically derived predictions that seek a unique relationship between capillary pressure, fluid saturation, and specific interfacial area (Pc-Sw-Anw). Previous experimental studies sought to test this functional dependence under quasi-equilibrium conditions (assumed static on the imaging time-scale); however, applying predictive models developed under static conditions for dynamic scenarios can lead to substantial flaws in predicted outcomes. Theory and models developed using dynamic data can be verified using fast x-ray microtomography which allows for the unprecedented measurement of developing interfacial areas, curvatures, and trapping behaviors of fluid phases in three-dimensional systems. We will present results of drainage and imbibition experiments of air and water within a mixture of glass beads. The experiments were performed under both quasi-equilibrium and dynamic conditions at the Advanced Photon Source (APS) at Argonne National Laboratory. Fast x-ray microtomography was achieved by utilizing the high brilliance of the x-ray beam at the APS under pink-beam conditions where the white beam is modified with a 4 mm Al absorber and a 0.8 mrad Pt-coated mirror to eliminate low and high-energy photons, respectively. We present a comparison of the results from the quasi-equilibrium and dynamic experiments in an effort to determine if the Pc-Sw-Anw relationship is comparable under either experimental condition and to add to the discussion on whether the Pc-Sw-Anw relationship is unique as hypothesized by existing theory.
TEMPEST/N33.5. Computational Fluid Dynamics Package For Incompressible, 3D, Time Dependent Pro
Trent, Dr.D.S.; Eyler, Dr.L.L.
1991-04-01
TEMPESTN33.5 provides numerical solutions to general incompressible flow problems with coupled heat transfer in fluids and solids. Turbulence is created with a k-e model and gas, liquid or solid constituents may be included with the bulk flow. Problems may be modeled in Cartesian or cylindrical coordinates. Limitations include incompressible flow, Boussinesq approximation, and passive constituents. No direct steady state solution is available; steady state is obtained as the limit of a transient.
Rigano, S; Ferrazzi, E; Boito, S; Pennati, G; Padoan, A; Galan, H
2010-01-01
The primary aim of this pilot study was to study uterine artery (UtA) blood flow volume in uneventful human pregnancies delivered at term, at mid and late gestation by means of 3D and bi-dimensional ultrasound imaging with angio-Doppler combined with fluid-dynamic modeling. Secondary aims were to correlate flow volume to placental site and to UtA Pulsatility Index (PI). Women with singleton, low-risk pregnancies were examined at mid and late gestation. The structure and course of the uterine artery (UtA) was studied in each patient by means of 3D-angio-Doppler and included vessel diameter D, blood flow velocity and PI (measured along the UtA). Fetal weight estimation and placental insertion site were assessed by ultrasound. A robust fluid-dynamic modeling was applied to calculate absolute flow and flow per unit fetal weight. Mean UtA diameter and blood flow velocity increased significantly (p < 0.0001) from mid-gestation to late gestation from 2.6 mm and 67.5 cm/s, to 3.0 mm and 85.3 cm/s, respectively, yielding an increasing absolute flow troughout gestation. h coefficient, derived by fluid-dynamic modeling to calculate mean velocity, increased significantly from 0.52 at mid-gestation to 0.57 at late gestation. UtA blood flow volume ml/min/kg-fetal weight was significantly higher at mid-gestation than at late gestation (535 ml/min/kg vs 193 ml/min/kg; p < 0.0001). In cases with strictly lateral placentas the ipsilateral UtA accommodates at mid and late gestation 63% and 67% of the total UtA flow. In central placentas UtA flow was evenly distributed between the two vessels. An inverse correlation was observed between PI and blood flow volume ml/min/kg (Pearson's coefficient r = -0.54). Our work confirms the technological and methodological limitations in the measurement of uterine artery blood flow. However, Doppler measurements supported by three-dimensional angio imaging of the uterine vessel, high resolution imaging and diameter measurement, and a robust
3D Dynamic Echocardiography with a Digitizer
NASA Astrophysics Data System (ADS)
Oshiro, Osamu; Matani, Ayumu; Chihara, Kunihiro
1998-05-01
In this paper,a three-dimensional (3D) dynamic ultrasound (US) imaging system,where a US brightness-mode (B-mode) imagetriggered with an R-wave of electrocardiogram (ECG)was obtained with an ultrasound diagnostic deviceand the location and orientation of the US probewere simultaneously measured with a 3D digitizer, is described.The obtained B-mode imagewas then projected onto a virtual 3D spacewith the proposed interpolation algorithm using a Gaussian operator.Furthermore, a 3D image was presented on a cathode ray tube (CRT)and stored in virtual reality modeling language (VRML).We performed an experimentto reconstruct a 3D heart image in systole using this system.The experimental results indicatethat the system enables the visualization ofthe 3D and internal structure of a heart viewed from any angleand has potential for use in dynamic imaging,intraoperative ultrasonography and tele-medicine.
Static & Dynamic Response of 3D Solids
Energy Science and Technology Software Center (ESTSC)
1996-07-15
NIKE3D is a large deformations 3D finite element code used to obtain the resulting displacements and stresses from multi-body static and dynamic structural thermo-mechanics problems with sliding interfaces. Many nonlinear and temperature dependent constitutive models are available.
3D dynamic roadmapping for abdominal catheterizations.
Bender, Frederik; Groher, Martin; Khamene, Ali; Wein, Wolfgang; Heibel, Tim Hauke; Navab, Nassir
2008-01-01
Despite rapid advances in interventional imaging, the navigation of a guide wire through abdominal vasculature remains, not only for novice radiologists, a difficult task. Since this navigation is mostly based on 2D fluoroscopic image sequences from one view, the process is slowed down significantly due to missing depth information and patient motion. We propose a novel approach for 3D dynamic roadmapping in deformable regions by predicting the location of the guide wire tip in a 3D vessel model from the tip's 2D location, respiratory motion analysis, and view geometry. In a first step, the method compensates for the apparent respiratory motion in 2D space before backprojecting the 2D guide wire tip into three dimensional space, using a given projection matrix. To countervail the error connected to the projection parameters and the motion compensation, as well as the ambiguity caused by vessel deformation, we establish a statistical framework, which computes a reliable estimate of the guide wire tip location within the 3D vessel model. With this 2D-to-3D transfer, the navigation can be performed from arbitrary viewing angles, disconnected from the static perspective view of the fluoroscopic sequence. Tests on a realistic breathing phantom and on synthetic data with a known ground truth clearly reveal the superiority of our approach compared to naive methods for 3D roadmapping. The concepts and information presented in this paper are based on research and are not commercially available. PMID:18982662
INCORPORATING DYNAMIC 3D SIMULATION INTO PRA
Steven R Prescott; Curtis Smith
2011-07-01
provide superior results and insights. We also couple the state model with the dynamic 3D simulation analysis representing events (such as flooding) to determine which (if any) components fail. Not only does the simulation take into account any failed items from the state model, but any failures caused by the simulation are incorporated back into the state model and factored into the overall results. Using this method we incorporate accurate 3D simulation results, eliminate static-based PRA issues, and have time ordered failure information.
An annotation system for 3D fluid flow visualization
NASA Technical Reports Server (NTRS)
Loughlin, Maria M.; Hughes, John F.
1995-01-01
Annotation is a key activity of data analysis. However, current systems for data analysis focus almost exclusively on visualization. We propose a system which integrates annotations into a visualization system. Annotations are embedded in 3D data space, using the Post-it metaphor. This embedding allows contextual-based information storage and retrieval, and facilitates information sharing in collaborative environments. We provide a traditional database filter and a Magic Lens filter to create specialized views of the data. The system has been customized for fluid flow applications, with features which allow users to store parameters of visualization tools and sketch 3D volumes.
Fast 3D fluid registration of brain magnetic resonance images
NASA Astrophysics Data System (ADS)
Leporé, Natasha; Chou, Yi-Yu; Lopez, Oscar L.; Aizenstein, Howard J.; Becker, James T.; Toga, Arthur W.; Thompson, Paul M.
2008-03-01
Fluid registration is widely used in medical imaging to track anatomical changes, to correct image distortions, and to integrate multi-modality data. Fluid mappings guarantee that the template image deforms smoothly into the target, without tearing or folding, even when large deformations are required for accurate matching. Here we implemented an intensity-based fluid registration algorithm, accelerated by using a filter designed by Bro-Nielsen and Gramkow. We validated the algorithm on 2D and 3D geometric phantoms using the mean square difference between the final registered image and target as a measure of the accuracy of the registration. In tests on phantom images with different levels of overlap, varying amounts of Gaussian noise, and different intensity gradients, the fluid method outperformed a more commonly used elastic registration method, both in terms of accuracy and in avoiding topological errors during deformation. We also studied the effect of varying the viscosity coefficients in the viscous fluid equation, to optimize registration accuracy. Finally, we applied the fluid registration algorithm to a dataset of 2D binary corpus callosum images and 3D volumetric brain MRIs from 14 healthy individuals to assess its accuracy and robustness.
3D imaging using projected dynamic fringes
NASA Astrophysics Data System (ADS)
Shaw, Michael M.; Atkinson, John T.; Harvey, David M.; Hobson, Clifford A.; Lalor, Michael J.
1994-12-01
An instrument capable of highly accurate, non-contact range measurement has been developed, which is based upon the principle of projected rotating fringes. More usually known as dynamic fringe projection, it is this technique which is exploited in the dynamic automated range transducer (DART). The intensity waveform seen at the target and sensed by the detector, contains all the information required to accurately determine the fringe order. This, in turn, allows the range to be evaluated by the substitution of the fringe order into a simple algebraic expression. Various techniques for the analysis of the received intensity signals from the surface of the target have been investigated. The accuracy to which the range can be determined ultimately depends upon the accuracy to which the fringe order can be evaluated from the received intensity waveform. It is extremely important to be able to closely determine the fractional fringe order value, to achieve any meaningful results. This paper describes a number of techniques which have been used to analyze the intensity waveform, and critically appraises their suitability in terms of accuracy and required speed of operation. This work also examines the development of this instrument for three-dimensional measurements based on single or two beam systems. Using CCD array detectors, a 3-D range map of the object's surface may be produced.
Simulating coronal condensation dynamics in 3D
NASA Astrophysics Data System (ADS)
Moschou, S. P.; Keppens, R.; Xia, C.; Fang, X.
2015-12-01
We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.
Complex flow dynamics around 3D microbot prototypes.
Martínez-Aranda, Sergio; Galindo-Rosales, Francisco J; Campo-Deaño, Laura
2016-02-28
A new experimental setup for the study of the complex flow dynamics around 3D microbot prototypes in a straight microchannel has been developed and assessed. The ultimate aim of this work is focused on the analysis of the morphology of different microbot prototypes to get a better insight into their efficiency when they swim through the main conduits of the human circulatory system. The setup consists of a fused silica straight microchannel with a 3D microbot prototype fastened in the center of the channel cross-section by an extremely thin support. Four different prototypes were considered: a cube, a sphere and two ellipsoids with aspect ratios of 1 : 2 and 1 : 4, respectively. Flow visualization and micro-particle image velocimetry (μPIV) measurements were performed using Newtonian and viscoelastic blood analogue fluids. An efficiency parameter, ℑ, to discriminate the prototypes in terms of flow disturbance has been proposed. PMID:26790959
3-D Dynamic Behavior of Generalized Polar Wind
NASA Astrophysics Data System (ADS)
Barakat, A. R.; Schunk, R. W.; Demars, H. G.
2003-12-01
The dynamic behavior of the high-latitude plasma during a representative geomagnetic storm is investigated using a 3-D macroscopic particle-in-cell (mac-PIC) model. In this study, we simulate the behavior of a large number ( ˜100 to 1000) of plasma-filled geomagnetic flux tubes. Each flux tube extends from 1200 km to several Earth radii, includes ˜106 simulation particles, and is followed for ˜12 hours. The lower boundary conditions of the model are provided by a 3-D fluid-like model that extends down to 100 km. Several physical mechanisms are included such as wave-particle interactions, ion-ion collisions, low-altitude ion energization, and magnetospheric particles. The computing-intensive nature of the model requires the utilization of parallel programming techniques. We use a cluster of five nodes, with two (1.6 GHz) processors each, that is available at Utah State University, with the intention of transferring the code to a bigger facility in the future. A 3-D picture is assembled from the temporal evolution of the individual flux tubes by keeping track of their locations. This 3-D picture facilitates comparison with observations, such as radar and satellite measurements. The model and its preliminary results are presented.
Phase unwrapping in the dynamic 3D measurement
NASA Astrophysics Data System (ADS)
Su, Xianyu; Zhang, Qican
2010-04-01
In the dynamic 3D shape measurement phase distribution has 3D character, in which phase changes along x and y directions in space and also along t direction in time. 3D phase unwrapping plays a very important role in the dynamic 3D shape measurement. In the dynamic 3D shape measurement methods based on the structured illumination, Fourier transformation profilometry (FTP) is particularly fit for dynamic 3D measurement, because of only one fringe pattern needed and full field analysis. In this paper some 3D phase unwrapping techniques for dynamic 3D shape measurement mainly in our Lab. are presented and reviewed. The basic methods and algorithm design are introduced. The basic methods include direct 3D phase unwrapping, 3D diamond phase unwrapping, 3D phase unwrapping based on reliability ordering, 3D phase unwrapping based on marked fringe tracing. The advantage of the phase unwrapping based on reliability ordering is that the path of phase unwrapping is always along the direction from the pixel with higher reliability parameter value to the pixel with low reliability parameter value. Therefore, in the worse case the error is limited, if there is any, to local minimum areas.
ICEd-ALE Treatment of 3-D Fluid Flow.
Energy Science and Technology Software Center (ESTSC)
1999-09-13
Version: 00 SALE3D calculates three-dimensional fluid flow at all speeds, from the incompressible limit to highly supersonic. An implicit treatment of the pressure calculation similar to that in the Implicit Continuous-fluid Eulerian (ICE) technique provides this flow speed flexibility. In addition, the computing mesh may move with the fluid in a typical Lagrangian fashion, be held in an Eulerian manner, or move in some arbitrarily specified way to provide a continuous rezoning capability. This latitudemore » results from use of an Arbitrary Lagrangian-Eulerian (ALE) treatment of the mesh. The partial differential equations solved are the Navier-Stokes equations and the mass and internal energy equations. The fluid pressure is determined from an equation of state and supplemented with an artificial viscous pressure for the computation of shock waves. The computing mesh consists of a three-dimensional network of arbitrarily shaped, six-sided deformable cells, and a variety of user-selectable boundary conditions are provided in the program.« less
3D-dynamic representation of DNA sequences.
Wąż, Piotr; Bielińska-Wąż, Dorota
2014-03-01
A new 3D graphical representation of DNA sequences is introduced. This representation is called 3D-dynamic representation. It is a generalization of the 2D-dynamic dynamic representation. The sequences are represented by sets of "material points" in the 3D space. The resulting 3D-dynamic graphs are treated as rigid bodies. The descriptors characterizing the graphs are analogous to the ones used in the classical dynamics. The classification diagrams derived from this representation are presented and discussed. Due to the third dimension, "the history of the graph" can be recognized graphically because the 3D-dynamic graph does not overlap with itself. Specific parts of the graphs correspond to specific parts of the sequence. This feature is essential for graphical comparisons of the sequences. Numerically, both 2D and 3D approaches are of high quality. In particular, a difference in a single base between two sequences can be identified and correctly described (one can identify which base) by both 2D and 3D methods. PMID:24567158
DYNAMIC 3D QSAR TECHNIQUES: APPLICATIONS IN TOXICOLOGY
Two dynamic techniques recently developed to account for conformational flexibility of chemicals in 3D QSARs are presented. In addition to the impact of conformational flexibility of chemicals in 3D QSAR models, the applicability of various molecular descriptors is discussed. The...
Dynamic 3D echocardiography in virtual reality
van den Bosch, Annemien E; Koning, Anton HJ; Meijboom, Folkert J; McGhie, Jackie S; Simoons, Maarten L; van der Spek, Peter J; Bogers, Ad JJC
2005-01-01
Background This pilot study was performed to evaluate whether virtual reality is applicable for three-dimensional echocardiography and if three-dimensional echocardiographic 'holograms' have the potential to become a clinically useful tool. Methods Three-dimensional echocardiographic data sets from 2 normal subjects and from 4 patients with a mitral valve pathological condition were included in the study. The three-dimensional data sets were acquired with the Philips Sonos 7500 echo-system and transferred to the BARCO (Barco N.V., Kortrijk, Belgium) I-space. Ten independent observers assessed the 6 three-dimensional data sets with and without mitral valve pathology. After 10 minutes' instruction in the I-Space, all of the observers could use the virtual pointer that is necessary to create cut planes in the hologram. Results The 10 independent observers correctly assessed the normal and pathological mitral valve in the holograms (analysis time approximately 10 minutes). Conclusion this report shows that dynamic holographic imaging of three-dimensional echocardiographic data is feasible. However, the applicability and use-fullness of this technology in clinical practice is still limited. PMID:16375768
Dynamic 3D Visualization of Vocal Tract Shaping During Speech
Zhu, Yinghua; Kim, Yoon-Chul; Proctor, Michael I.; Narayanan, Shrikanth S.; Nayak, Krishna S.
2014-01-01
Noninvasive imaging is widely used in speech research as a means to investigate the shaping and dynamics of the vocal tract during speech production. 3D dynamic MRI would be a major advance, as it would provide 3D dynamic visualization of the entire vocal tract. We present a novel method for the creation of 3D dynamic movies of vocal tract shaping based on the acquisition of 2D dynamic data from parallel slices and temporal alignment of the image sequences using audio information. Multiple sagittal 2D real-time movies with synchronized audio recordings are acquired for English vowel-consonant-vowel stimuli /ala/, /aɹa/, /asa/ and /aʃa/. Audio data are aligned using mel-frequency cepstral coefficients (MFCC) extracted from windowed intervals of the speech signal. Sagittal image sequences acquired from all slices are then aligned using dynamic time warping (DTW). The aligned image sequences enable dynamic 3D visualization by creating synthesized movies of the moving airway in the coronal planes, visualizing desired tissue surfaces and tube-shaped vocal tract airway after manual segmentation of targeted articulators and smoothing. The resulting volumes allow for dynamic 3D visualization of salient aspects of lingual articulation, including the formation of tongue grooves and sublingual cavities, with a temporal resolution of 78 ms. PMID:23204279
A 3-D perspective of dynamic behaviour of heterogeneous solids
NASA Astrophysics Data System (ADS)
Lu, Yong; Zhou, Rongxin
2015-09-01
The dynamic behaviour of concrete-like materials under high strain rates has been a subject of continuous scrutiny over the years. A prevailing explanation attributes much of the dynamic increase of strength, especially under compression, to the macroscopic inertia confinement. Studies conducted by the authors' group using meso-scale computational models suggest that the heterogeneity of the material composition, in particular the involvement of the aggregates, also plays a sensible part in the process of damage evolution and the increase of the bulk strength under high strain rates, and a detailed investigation into this effect would benefit if a realistic representation of the heterogeneity in 3D can be achieved. This paper presents some recent progress in the development of a 3-D meso-scale computational model incorporating randomly-shaped 3-D aggregate particles, including the general validation of the model, and application in the simulation of the dynamic response of concrete under high strain rate compression.
An Evaluative Review of Simulated Dynamic Smart 3d Objects
NASA Astrophysics Data System (ADS)
Romeijn, H.; Sheth, F.; Pettit, C. J.
2012-07-01
Three-dimensional (3D) modelling of plants can be an asset for creating agricultural based visualisation products. The continuum of 3D plants models ranges from static to dynamic objects, also known as smart 3D objects. There is an increasing requirement for smarter simulated 3D objects that are attributed mathematically and/or from biological inputs. A systematic approach to plant simulation offers significant advantages to applications in agricultural research, particularly in simulating plant behaviour and the influences of external environmental factors. This approach of 3D plant object visualisation is primarily evident from the visualisation of plants using photographed billboarded images, to more advanced procedural models that come closer to simulating realistic virtual plants. However, few programs model physical reactions of plants to external factors and even fewer are able to grow plants based on mathematical and/or biological parameters. In this paper, we undertake an evaluation of plant-based object simulation programs currently available, with a focus upon the components and techniques involved in producing these objects. Through an analytical review process we consider the strengths and weaknesses of several program packages, the features and use of these programs and the possible opportunities in deploying these for creating smart 3D plant-based objects to support agricultural research and natural resource management. In creating smart 3D objects the model needs to be informed by both plant physiology and phenology. Expert knowledge will frame the parameters and procedures that will attribute the object and allow the simulation of dynamic virtual plants. Ultimately, biologically smart 3D virtual plants that react to changes within an environment could be an effective medium to visually represent landscapes and communicate land management scenarios and practices to planners and decision-makers.
On the Dynamic Programming Approach for the 3D Navier-Stokes Equations
Manca, Luigi
2008-06-15
The dynamic programming approach for the control of a 3D flow governed by the stochastic Navier-Stokes equations for incompressible fluid in a bounded domain is studied. By a compactness argument, existence of solutions for the associated Hamilton-Jacobi-Bellman equation is proved. Finally, existence of an optimal control through the feedback formula and of an optimal state is discussed.
3D Fluid-Structure Modeling of a Monofin
NASA Astrophysics Data System (ADS)
Monier, L.; Razafimahery, F.; Rakotomanana, L.
2010-10-01
The purpose of this paper is to develop a numerical modelisation for the behaviour of a monofin. We have developped a fluid struture model simulating the movement of a fin in a swimming pool. We first present the geometry and the equations and then proceed to different numerical experiments in order to validate the model.
The Vibrational Dynamics of 3D HOCl Above Dissociation
NASA Astrophysics Data System (ADS)
Lin, Yi-Der; Reichl, Linda; Jung, Christof
2015-03-01
We have analyzed the vibrational dynamics of HOCl above dissociation using a 3D energy surface which governs the vibrational dynamics of HOCl above dissociation. The dynamics is dominated by an invariant manifold which is transversally unstable for small spacing between Cl and HO complex, and stable for large spacing. Above dissociation, the InM separates two mirror image periodic orbits, embedded in a large chaotic sea, that can hold a large number of quantum states. These periodic orbits have the capability of forming significant quasibound states of the molecule above dissociation. Welch Foundation.
3D two-fluid simulations of turbulence in LAPD
NASA Astrophysics Data System (ADS)
Fisher, Dustin M.
The Large Plasma Device (LAPD) is modeled using a modified version of the 3D Global Braginskii Solver code (GBS) for a nominal Helium plasma. The unbiased low-flow regime is explored in simulations where there is an intrinsic E x B rotation of the plasma. In the simulations this rotation is caused primarily by sheath effects with the Reynolds stress and J x B torque due to a cross-field Pederson conductivity having little effect. Explicit biasing simulations are also explored for the first time where the intrinsic rotation of the plasma is modified through boundary conditions that mimic the biasable limiter used in LAPD. Comparisons to experimental measurements in the unbiased case show strong qualitative agreement with the data, particularly the radial dependence of the density fluctuations, cross-correlation lengths, radial flux dependence outside of the cathode edge, and camera imagery. Kelvin Helmholtz (KH) turbulence at relatively large scales is the dominant driver of cross-field transport in these simulations with smaller-scale drift waves and sheath modes playing a secondary role. Plasma holes and blobs arising from KH vortices are consistent with the scale sizes and overall appearance of those in LAPD camera images. The addition of ion-neutral collisions in the unbiased simulations at previously theorized values reduces the radial particle flux due to a modest stabilizing contribution of the collisions on the KH-modes driving the turbulent transport. In the biased runs the ion-neutral collisions have a much smaller effect due to the modification of the potential from sheath terms. In biasing the plasma to increase the intrinsic rotation, simulations show the emergence of a nonlinearly saturated coherent mode of order m = 6. In addition, the plasma inside of the cathode edge becomes quiescent due to the strong influence of the wall bias in setting up the equilibrium plasma potential. Biasing in the direction opposite to the intrinsic flow reduces the
3-D structure and dynamics of microtubule self-organization
NASA Astrophysics Data System (ADS)
Wang, Jing; Ou-Yang, H. Daniel
2008-03-01
Laser scanning confocal microscopy was used to study the dynamics of 3D assemblies spontaneously formed in microtubule (MT) solutions. Microtubule solutions prepared by mixing and incubating tubulin in the presence of GTP and Oregon Green conjugated taxol in PM buffer were placed in long, sub-millimeter thin glass cells by the capillary action. Within 24 hours, starting with a uniform distribution, microtubules were found to be gradually separated into a few large ``buckled'' bundles along the long direction, and in the middle plane, of the sample cell. A well-defined wavelength of the buckling sinusoids was around 510 μm. The cross section of these round bundles was approximately 40 μm in diameter and the lengths were several centimeters. Detailed analysis of the 3-D image within the bundles revealed that each bundle seemed to consist of loosely packed MTs. It appeared that MTs were phase separated resulting from attractive interactions between charged MT fibers. The ``buckling'' behavior could be the result of geometrical constraints of the repulsive cell walls and the repulsive interaction between bundles. Detailed 3-D observations of the dynamic evolution of MT assembly could provide insight to the mechanisms of cellular MT organization and phase separation of charged colloidal rods.
Fluid and cell behaviors along a 3D printed alginate/gelatin/fibrin channel.
Xu, Yufan; Wang, Xiaohong
2015-08-01
Three-dimensional (3D) cell manipulation is available with the integration of microfluidic technology and rapid prototyping techniques. High-Fidelity (Hi-Fi) constructs hold enormous therapeutic potential for organ manufacturing and regenerative medicine. In the present paper we introduced a quasi-three-dimensional (Q3D) model with parallel biocompatible alginate/gelatin/fibrin hurdles. The behaviors of fluids and cells along the microfluidic channels with various widths were studied. Cells inside the newly designed microfluidic channels attached and grew well. Morphological changes of adipose-derived stem cells (ADSCs) in both two-dimensional (2D) and 3D milieu were found on the printed constructs. Endothelialization occurred with the co-cultures of ADSCs and hepatocytes. This study provides insights into the interactions among fluids, cells and biomaterials, the behaviors of fluids and cells along the microfluidic channels, and the applications of Q3D techniques. PMID:25727058
Correlative Microscopy for 3D Structural Analysis of Dynamic Interactions
Jun, Sangmi; Zhao, Gongpu; Ning, Jiying; Gibson, Gregory A.; Watkins, Simon C.; Zhang, Peijun
2013-01-01
Cryo-electron tomography (cryoET) allows 3D visualization of cellular structures at molecular resolution in a close-to-physiological state1. However, direct visualization of individual viral complexes in their host cellular environment with cryoET is challenging2, due to the infrequent and dynamic nature of viral entry, particularly in the case of HIV-1. While time-lapse live-cell imaging has yielded a great deal of information about many aspects of the life cycle of HIV-13-7, the resolution afforded by live-cell microscopy is limited (~ 200 nm). Our work was aimed at developing a correlation method that permits direct visualization of early events of HIV-1 infection by combining live-cell fluorescent light microscopy, cryo-fluorescent microscopy, and cryoET. In this manner, live-cell and cryo-fluorescent signals can be used to accurately guide the sampling in cryoET. Furthermore, structural information obtained from cryoET can be complemented with the dynamic functional data gained through live-cell imaging of fluorescent labeled target. In this video article, we provide detailed methods and protocols for structural investigation of HIV-1 and host-cell interactions using 3D correlative high-speed live-cell imaging and high-resolution cryoET structural analysis. HeLa cells infected with HIV-1 particles were characterized first by confocal live-cell microscopy, and the region containing the same viral particle was then analyzed by cryo-electron tomography for 3D structural details. The correlation between two sets of imaging data, optical imaging and electron imaging, was achieved using a home-built cryo-fluorescence light microscopy stage. The approach detailed here will be valuable, not only for study of virus-host cell interactions, but also for broader applications in cell biology, such as cell signaling, membrane receptor trafficking, and many other dynamic cellular processes. PMID:23852318
3D precision surface measurement by dynamic structured light
NASA Astrophysics Data System (ADS)
Franke, Ernest A.; Magee, Michael J.; Mitchell, Joseph N.; Rigney, Michael P.
2004-02-01
This paper describes a 3-D imaging technique developed as an internal research project at Southwest Research Institute. The technique is based on an extension of structured light methods in which a projected pattern of parallel lines is rotated over the surface to be measured. A sequence of images is captured and the surface elevation at any location can then be determined from measurements of the temporal pattern, at any point, without considering any other points on the surface. The paper describes techniques for system calibration and surface measurement based on the method of projected quadric shells. Algorithms were developed for image and signal analysis and computer programs were written to calibrate the system and to calculate 3-D coordinates of points on a measured surface. A prototype of the Dynamic Structured Light (DSL) 3-D imaging system was assembled and typical parts were measured. The design procedure was verified and used to implement several different configurations with different measurement volumes and measurement accuracy. A small-parts measurement accuracy of 32 micrometers (.0012") RMS was verified by measuring the surface of a precision-machined plane. Large aircraft control surfaces were measured with a prototype setup that provided .02" depth resolution over a 4" by 8" field of view. Measurement times are typically less than three minutes for 300,000 points. A patent application has been filed.
3D hybrid simulations with gyrokinetic particle ions and fluid electrons
Belova, E.V.; Park, W.; Fu, G.Y.; Strauss, H.R.; Sugiyama, L.E.
1998-12-31
The previous hybrid MHD/particle model (MH3D-K code) represented energetic ions as gyrokinetic (or drift-kinetic) particles coupled to MHD equations using the pressure or current coupling scheme. A small energetic to bulk ion density ratio was assumed, n{sub h}/n{sub b} {much_lt} 1, allowing the neglect of the energetic ion perpendicular inertia in the momentum equation and the use of MHD Ohm`s law E = {minus}v{sub b} {times} B. A generalization of this model in which all ions are treated as gyrokinetic/drift-kinetic particles and fluid description is used for the electron dynamics is considered in this paper.
GPU-Based Visualization of 3D Fluid Interfaces using Level Set Methods
NASA Astrophysics Data System (ADS)
Kadlec, B. J.
2009-12-01
We model a simple 3D fluid-interface problem using the level set method and visualize the interface as a dynamic surface. Level set methods allow implicit handling of complex topologies deformed by evolutions where sharp changes and cusps are present without destroying the representation. We present a highly optimized visualization and computation algorithm that is implemented in CUDA to run on the NVIDIA GeForce 295 GTX. CUDA is a general purpose parallel computing architecture that allows the NVIDIA GPU to be treated like a data parallel supercomputer in order to solve many computational problems in a fraction of the time required on a CPU. CUDA is compared to the new OpenCL™ (Open Computing Language), which is designed to run on heterogeneous computing environments but does not take advantage of low-level features in NVIDIA hardware that provide significant speedups. Therefore, our technique is implemented using CUDA and results are compared to a single CPU implementation to show the benefits of using the GPU and CUDA for visualizing fluid-interface problems. We solve a 1024^3 problem and experience significant speedup using the NVIDIA GeForce 295 GTX. Implementation details for mapping the problem to the GPU architecture are described as well as discussion on porting the technique to heterogeneous devices (AMD, Intel, IBM) using OpenCL. The results present a new interactive system for computing and visualizing the evolution of fluid interface problems on the GPU.
Dynamical Systems Analysis of Fully 3D Ocean Features
NASA Astrophysics Data System (ADS)
Pratt, L. J.
2011-12-01
Dynamical systems analysis of transport and stirring processes has been developed most thoroughly for 2D flow fields. The calculation of manifolds, turnstile lobes, transport barriers, etc. based on observations of the ocean is most often conducted near the sea surface, whereas analyses at depth, usually carried out with model output, is normally confined to constant-z surfaces. At the meoscale and larger, ocean flows are quasi 2D, but smaller scale (submesoscale) motions, including mixed layer phenomena with significant vertical velocity, may be predominantly 3D. The zoology of hyperbolic trajectories becomes richer in such cases and their attendant manifolds are much more difficult to calculate. I will describe some of the basic geometrical features and corresponding Lagrangian Coherent Features expected to arise in upper ocean fronts, eddies, and Langmuir circulations. Traditional GFD models such as the rotating can flow may capture the important generic features. The dynamical systems approach is most helpful when these features are coherent and persistent and the implications and difficulties for this requirement in fully 3D flows will also be discussed.
Dynamic deformable models for 3D MRI heart segmentation
NASA Astrophysics Data System (ADS)
Zhukov, Leonid; Bao, Zhaosheng; Gusikov, Igor; Wood, John; Breen, David E.
2002-05-01
Automated or semiautomated segmentation of medical images decreases interstudy variation, observer bias, and postprocessing time as well as providing clincally-relevant quantitative data. In this paper we present a new dynamic deformable modeling approach to 3D segmentation. It utilizes recently developed dynamic remeshing techniques and curvature estimation methods to produce high-quality meshes. The approach has been implemented in an interactive environment that allows a user to specify an initial model and identify key features in the data. These features act as hard constraints that the model must not pass through as it deforms. We have employed the method to perform semi-automatic segmentation of heart structures from cine MRI data.
NASA Astrophysics Data System (ADS)
Zhang, Zhi-Qian; Liu, G. R.; Khoo, Boo Cheong
2013-02-01
A three-dimensional immersed smoothed finite element method (3D IS-FEM) using four-node tetrahedral element is proposed to solve 3D fluid-structure interaction (FSI) problems. The 3D IS-FEM is able to determine accurately the physical deformation of the nonlinear solids placed within the incompressible viscous fluid governed by Navier-Stokes equations. The method employs the semi-implicit characteristic-based split scheme to solve the fluid flows and smoothed finite element methods to calculate the transient dynamics responses of the nonlinear solids based on explicit time integration. To impose the FSI conditions, a novel, effective and sufficiently general technique via simple linear interpolation is presented based on Lagrangian fictitious fluid meshes coinciding with the moving and deforming solid meshes. In the comparisons to the referenced works including experiments, it is clear that the proposed 3D IS-FEM ensures stability of the scheme with the second order spatial convergence property; and the IS-FEM is fairly independent of a wide range of mesh size ratio.
NASA Astrophysics Data System (ADS)
Vucinic, Dean; Deen, Danny; Oanta, Emil; Batarilo, Zvonimir; Lacor, Chris
This paper focuses on visualization and manipulation of graphical content in distributed network environments. The developed graphical middleware and 3D desktop prototypes were specialized for situational awareness. This research was done in the LArge Scale COllaborative decision support Technology (LASCOT) project, which explored and combined software technologies to support human-centred decision support system for crisis management (earthquake, tsunami, flooding, airplane or oil-tanker incidents, chemical, radio-active or other pollutants spreading, etc.). The performed state-of-the-art review did not identify any publicly available large scale distributed application of this kind. Existing proprietary solutions rely on the conventional technologies and 2D representations. Our challenge was to apply the "latest" available technologies, such Java3D, X3D and SOAP, compatible with average computer graphics hardware. The selected technologies are integrated and we demonstrate: the flow of data, which originates from heterogeneous data sources; interoperability across different operating systems and 3D visual representations to enhance the end-users interactions.
High-Resolution 3D Seismic Imaging of Fluid Flow Anomalies in the Southwest Barents Sea
NASA Astrophysics Data System (ADS)
Planke, S.; Eriksen, F. N.; Eriksen, O. K.; Assad, M.; Stokke, H. H.
2014-12-01
Fluid flow features imaged as gas flares in the water column, pockmarks and mud volcanoes on the seabed, and high-amplitude cross-cutting reflections and bright spots in the sub-surface are abundant in the SW Barents Sea offshore northern Norway. This region is covered by extensive conventional 2D and 3D deep penetration seismic reflection data and multibeam bathymetry. High-resolution 3D P-Cable seismic data have been acquired in the SW Barents Sea over the past few years to image the uppermost ca. 500 m of the sub-surface. The P-Cable system consist of 12 to 16 short streamers (25 m) that are towed on a cross-cable perpendicular to the vessel's steaming direction. This configuration allows for acquisition of seismic data with high trace density, typically with 6 m in-line separation. The vertical resolution is a good as 1-2 m using conventional site survey air gun configurations. The sedimentary succession in the SW Barents Sea consists of upper Paleozoic evaporites overlaid by Mesozoic and Cenozoic clastic sediments. There are several organic-rich intervals in the sequence, including Paleozoic coals and Triassic and Jurassic marine source rocks that are locally in the oil or gas maturation windows. Glacial erosion has locally removed kilometer thick Cenozoic and Mesozoic successions, leaving the Mesozoic and Paleozoic strata in shallow sub-surface. The new high-resolution 3D surveys have targeted shallow fluid anomalies in the subsurface. These are imaged as high-amplitude reflections in fault blocks and structural highs, locally cross-cutting well-defined Mesozoic reflections. Commonly, disturbed reflections are present in overlying sequences, or high-amplitude reflections are imaged in the glacial overburden sediments. Locally, hundreds of pockmarks are imaged by the seafloor reflection. The deep cross-cutting reflections are interpreted as hydrocarbon accumulations that locally migrate towards the surface. The fluids are stored in shallow gas pockets or
3D Dynamic Earthquake Fracture Simulation (Test Case)
NASA Astrophysics Data System (ADS)
Korkusuz Öztürk, Yasemin; Meral Özel, Nurcan; Ando, Ryosuke
2016-04-01
A 3D dynamic earthquake fracture simulation is being developed for the fault structures which are non-planar to understand heterogeneous stress states in the Marmara Sea. Locating in a seismic gap, a large earthquake is expected in the center of the Sea of Marmara. Concerning the fact that more than 14 million inhabitants of İstanbul, located very closely to the Marmara Sea, the importance of the analysis of the Central Marmara Sea is extremely high. A few 3D dynamic earthquake fracture studies have been already done in the Sea of Marmara for pure right lateral strike-slip stress regimes (Oglesby and Mai, 2012; Aochi and Ulrich, 2015). In this study, a 3D dynamic earthquake fracture model with heterogeneous stress patches from the TPV5, a SCEC code validation case, is adapted. In this test model, the fault and the ground surfaces are gridded by a scalene triangulation technique using GMSH program. For a grid size changing between 0.616 km and 1.050 km the number of elements for the fault surface is 1984 and for the ground surface is 1216. When these results are compared with Kaneko's results for TPV5 from SPECFEM3D, reliable findings could be observed for the first 6.5 seconds (stations on the fault) although a stability problem is encountered after this time threshold. To solve this problem grid sizes are made smaller, so the number of elements increase 7986 for the fault surface and 4867 for the ground surface. On the other hand, computational problems arise in that case, since the computation time is directly proportional to the number of total elements and the required memory also increases with the square of that. Therefore, it is expected that this method can be adapted for less coarse grid cases, regarding the main difficulty coming from the necessity of an effective supercomputer and run time limitations. The main objective of this research is to obtain 3D dynamic earthquake rupture scenarios, concerning not only planar and non-planar faults but also
High-resolution 3D seismic data characterize fluid flow systems in the SW Barents Sea
NASA Astrophysics Data System (ADS)
Bünz, Stefan; Mienert, Jürgen; Rajan, Anupama
2010-05-01
The flow of fluids through marine sediments is one of the most dominant and pervasive processes in continental margins. These processes control the evolution of a sedimentary basin and its seafloor environment, and have implications for hydrocarbon exploration and seabed ecosystems. Many seep sites at the seafloor are associated with large but complex faunal communities that have received significant attention in recent years. However, there is a need for a better understanding of the driving mechanism of fluid flow in various geological settings, the accumulation of fluids in the subsurface and their focused flow through conduits and/or faults to the seabed. The Barents Sea is a large hydrocarbon-prone basin of the Norwegian Arctic region. A significant portion of the hydrocarbons has leaked or migrated into the shallow subsurface and is now trapped in gas-hydrate and shallow-gas reservoirs. Furthermore, there are few places in the Barents Sea, where methane gas is leaking from the seafloor into the oceanosphere. Accumulations of free gas in the shallow subsurface are considered a geohazard. They constitute a risk for safe drilling operations and they may pose a threat to global climate if the seal that is trapping them is breached. P-Cable 3D high-resolution seismic data from the Ringvassøya Fault Complex and the Polheim Sub-Platform provide new and detailed insight into fluid flow controls and accumulation mechanisms. The data shows a wide variety of fluid flow features, mostly in the form of pockmarks, bright spots, wipe-out zones or vertical zones of disturbed reflectivity. Fluids migrate by both diapiric mechanism and channelized along sedimentary layers. Glacigenic sediments generally form a strong boundary for fluid flow in the very shallow section. However, we can recognize pockmarks not only at the seafloor but also at one subsurface layer approximately 50 m below sea floor indicating a former venting period in the SW Barents Sea. At few locations high
Modeling tree crown dynamics with 3D partial differential equations.
Beyer, Robert; Letort, Véronique; Cournède, Paul-Henry
2014-01-01
We characterize a tree's spatial foliage distribution by the local leaf area density. Considering this spatially continuous variable allows to describe the spatiotemporal evolution of the tree crown by means of 3D partial differential equations. These offer a framework to rigorously take locally and adaptively acting effects into account, notably the growth toward light. Biomass production through photosynthesis and the allocation to foliage and wood are readily included in this model framework. The system of equations stands out due to its inherent dynamic property of self-organization and spontaneous adaptation, generating complex behavior from even only a few parameters. The density-based approach yields spatially structured tree crowns without relying on detailed geometry. We present the methodological fundamentals of such a modeling approach and discuss further prospects and applications. PMID:25101095
Atomic-Resolution 3D Electron Microscopy with Dynamic Diffraction
O'Keefe, Michael A.; Downing, Kenneth H.; Wenk, Hans-Rudolf; Meisheng, Hu
2005-02-15
Achievement of atomic-resolution electron-beam tomography will allow determination of the three-dimensional structure of nanoparticles (and other suitable specimens) at atomic resolution. Three-dimensional reconstructions will yield ''section'' images that resolve atoms overlapped in normal electron microscope images (projections), resolving lighter atoms such as oxygen in the presence of heavier atoms, and atoms that lie on non-lattice sites such as those in non-periodic defect structures. Lower-resolution electron microscope tomography has been used to produce reconstructed 3D images of nanoparticles [1] but extension to atomic resolution is considered not to be straightforward. Accurate three-dimensional reconstruction from two-dimensional projections generally requires that intensity in the series of 2-D images be a monotonic function of the specimen structure (often specimen density, but in our case atomic potential). This condition is not satisfied in electron microscopy when specimens with strong periodicity are tilted close to zone-axis orientation and produce ''anomalous'' image contrast because of strong dynamic diffraction components. Atomic-resolution reconstructions from tilt series containing zone-axis images (with their contrast enhanced by strong dynamical scattering) can be distorted when the stronger zone-axis images overwhelm images obtained in other ''random'' orientations in which atoms do not line up in neat columns. The first demonstrations of 3-D reconstruction to atomic resolution used five zone-axis images from test specimens of staurolite consisting of a mix of light and heavy atoms [2,3]. Initial resolution was to the 1.6{angstrom} Scherzer limit of a JEOL-ARM1000. Later experiments used focal-series reconstruction from 5 to 10 images to produce staurolite images from the ARM1000 with resolution extended beyond the Scherzer limit to 1.38{angstrom} [4,5]. To obtain a representation of the three-dimensional structure, images were obtained
3D-printed soft microrobot for swimming in biological fluids.
Qiu, Tian; Palagi, Stefano; Fischer, Peer
2015-08-01
Microscopic artificial swimmers hold the potential to enable novel non-invasive medical procedures. In order to ease their translation towards real biomedical applications, simpler designs as well as cheaper yet more reliable materials and fabrication processes should be adopted, provided that the functionality of the microrobots can be kept. A simple single-hinge design could already enable micro-swimming in non-Newtonian fluids, which most bodily fluids are. Here, we address the fabrication of such single-hinge microrobots with a 3D-printed soft material. Firstly, a finite element model is developed to investigate the deformability of the 3D-printed microstructure under typical values of the actuating magnetic fields. Then the microstructures are fabricated by direct 3D-printing of a soft material and their swimming performances are evaluated. The speeds achieved with the 3D-printed microrobots are comparable to those obtained in previous work with complex fabrication procedures, thus showing great promise for 3D-printed microrobots to be operated in biological fluids. PMID:26737396
Image enhancement and segmentation of fluid-filled structures in 3D ultrasound images
NASA Astrophysics Data System (ADS)
Chalana, Vikram; Dudycha, Stephen; McMorrow, Gerald
2003-05-01
Segmentation of fluid-filled structures, such as the urinary bladder, from three-dimensional ultrasound images is necessary for measuring their volume. This paper describes a system for image enhancement, segmentation and volume measurement of fluid-filled structures on 3D ultrasound images. The system was applied for the measurement of urinary bladder volume. Results show an average error of less than 10% in the estimation of the total bladder volume.
A parallel algorithm for 3D dislocation dynamics
NASA Astrophysics Data System (ADS)
Wang, Zhiqiang; Ghoniem, Nasr; Swaminarayan, Sriram; LeSar, Richard
2006-12-01
Dislocation dynamics (DD), a discrete dynamic simulation method in which dislocations are the fundamental entities, is a powerful tool for investigation of plasticity, deformation and fracture of materials at the micron length scale. However, severe computational difficulties arising from complex, long-range interactions between these curvilinear line defects limit the application of DD in the study of large-scale plastic deformation. We present here the development of a parallel algorithm for accelerated computer simulations of DD. By representing dislocations as a 3D set of dislocation particles, we show here that the problem of an interacting ensemble of dislocations can be converted to a problem of a particle ensemble, interacting with a long-range force field. A grid using binary space partitioning is constructed to keep track of node connectivity across domains. We demonstrate the computational efficiency of the parallel micro-plasticity code and discuss how O(N) methods map naturally onto the parallel data structure. Finally, we present results from applications of the parallel code to deformation in single crystal fcc metals.
Quaini, A.; Canic, S.; Glowinski, R.; Igo, S.; Hartley, C.J.; Zoghbi, W.; Little, S.
2011-01-01
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
Quaini, A; Canic, S; Glowinski, R; Igo, S; Hartley, C J; Zoghbi, W; Little, S
2012-01-10
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals
NASA Astrophysics Data System (ADS)
Athavale, M. M.; Przekwas, A. J.; Hendricks, R. C.
1994-01-01
Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.
3-D dynamic rupture simulations by a finite volume method
NASA Astrophysics Data System (ADS)
Benjemaa, M.; Glinsky-Olivier, N.; Cruz-Atienza, V. M.; Virieux, J.
2009-07-01
Dynamic rupture of a 3-D spontaneous crack of arbitrary shape is investigated using a finite volume (FV) approach. The full domain is decomposed in tetrahedra whereas the surface, on which the rupture takes place, is discretized with triangles that are faces of tetrahedra. First of all, the elastodynamic equations are described into a pseudo-conservative form for an easy application of the FV discretization. Explicit boundary conditions are given using criteria based on the conservation of discrete energy through the crack surface. Using a stress-threshold criterion, these conditions specify fluxes through those triangles that have suffered rupture. On these broken surfaces, stress follows a linear slip-weakening law, although other friction laws can be implemented. For The Problem Version 3 of the dynamic-rupture code verification exercise conducted by the SCEC/USGS, numerical solutions on a planar fault exhibit a very high convergence rate and are in good agreement with the reference one provided by a finite difference (FD) technique. For a non-planar fault of parabolic shape, numerical solutions agree satisfactorily well with those obtained with a semi-analytical boundary integral method in terms of shear stress amplitudes, stopping phases arrival times and stress overshoots. Differences between solutions are attributed to the low-order interpolation of the FV approach, whose results are particularly sensitive to the mesh regularity (structured/unstructured). We expect this method, which is well adapted for multiprocessor parallel computing, to be competitive with others for solving large scale dynamic ruptures scenarios of seismic sources in the near future.
A support-operator method for 3-D rupture dynamics
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2009-06-01
We present a numerical method to simulate spontaneous shear crack propagation within a heterogeneous, 3-D, viscoelastic medium. Wave motions are computed on a logically rectangular hexahedral mesh, using the generalized finite-difference method of Support Operators (SOM). This approach enables modelling of non-planar surfaces and non-planar fault ruptures. Our implementation, the Support Operator Rupture Dynamics (SORD) code, is highly scalable, enabling large-scale, multiprocessors calculations. The fault surface is modelled by coupled double nodes, where rupture occurs as dictated by the local stress conditions and a frictional failure law. The method successfully performs test problems developed for the Southern California Earthquake Center (SCEC)/U.S. Geological Survey (USGS) dynamic earthquake rupture code validation exercise, showing good agreement with semi-analytical boundary integral method results. We undertake further dynamic rupture tests to quantify numerical errors introduced by shear deformations to the hexahedral mesh. We generate a family of meshes distorted by simple shearing, in the along-strike direction, up to a maximum of 73°. For SCEC/USGS validation problem number 3, grid-induced errors increase with mesh shear angle, with the logarithm of error approximately proportional to angle over the range tested. At 73°, rms misfits are about 10 per cent for peak slip rate, and 0.5 per cent for both rupture time and total slip, indicating that the method (which, up to now, we have applied mainly to near-vertical strike-slip faulting) is also capable of handling geometries appropriate to low-angle surface-rupturing thrust earthquakes. Additionally, we demonstrate non-planar rupture effects, by modifying the test geometry to include, respectively, cylindrical curvature and sharp kinks.
Fluids escape in subduction zones: new constraints from 3-D microtomography data
NASA Astrophysics Data System (ADS)
Le Roux, V.; Gaetani, G. A.; Slaugenwhite, J.; Miller, K.
2013-12-01
Large amounts of H2O are carried into trenches via subduction of the sediments, basaltic crust and uppermost mantle that make up the oceanic lithosphere. A major question is how much of this subducted H2O is released into the overlying mantle wedge, promoting melting, and how much is carried deeper into the mantle. This depends, at least in part, on whether H2O is able to form an interconnected network among the mineral grains that make up the rock down to very low fluid fractions. In order to achieve connectivity and allow the fluid phase to escape, a minimum amount of fluid (critical porosity) is required when dihedral angles are more than 60 degrees. We investigated the distribution of seawater in simplified sediment analogs (i.e. quartz for siliceous sediments; calcite for carbonate sediments), in natural clays (kaolinite and montmorillonite) and in bulk eclogite. Experiments were performed in a piston-cylinder apparatus at 2 GPa and 650°C. Fluid fractions ranged from ~10% to ~1% to determine the porosity at which connectivity of the seawater network is lost for each rock type. We used synchrotron X-ray microtomographic techniques (at Argonne National Laboratory, IL) to obtain 3-D images of the pore space network in order to constrain the grain scale distribution of fluids in a subducted slab. This nondestructive 3-D imaging technique has a spatial resolution of 0.7 μm and provides quantitative information on geometrical parameters of fluid topology, such as porosity, dihedral angle distribution, fluid channel sizes and connectivity. The geometrical parameters were extracted using the VSG Avizo software. This study lays the groundwork for determining the 3-D grain scale distribution of fluids in a range of subducted lithologies. Results from this study provide important new insights into the amount of fluid that can be transported into the deep mantle by subduction.
Computational Fluid Dynamics Library
Energy Science and Technology Software Center (ESTSC)
2005-03-04
CFDLib05 is the Los Alamos Computational Fluid Dynamics LIBrary. This is a collection of hydrocodes using a common data structure and a common numerical method, for problems ranging from single-field, incompressible flow, to multi-species, multi-field, compressible flow. The data structure is multi-block, with a so-called structured grid in each block. The numerical method is a Finite-Volume scheme employing a state vector that is fully cell-centered. This means that the integral form of the conservation lawsmore » is solved on the physical domain that is represented by a mesh of control volumes. The typical control volume is an arbitrary quadrilateral in 2D and an arbitrary hexahedron in 3D. The Finite-Volume scheme is for time-unsteady flow and remains well coupled by means of time and space centered fluxes; if a steady state solution is required, the problem is integrated forward in time until the user is satisfied that the state is stationary.« less
Dynamics of free subduction from 3-D boundary element modeling
NASA Astrophysics Data System (ADS)
Li, Zhong-Hai; Ribe, Neil M.
2012-06-01
In order better to understand the physical mechanisms underlying free subduction, we perform three-dimensional boundary-element numerical simulations of a dense fluid sheet with thickness h and viscosity η2 sinking in an `ambient mantle' with viscosity η1. The mantle layer is bounded above by a traction-free surface, and is either (1) infinitely deep or (2) underlain by a rigid boundary at a finite depth H + d, similar to the typical geometry used in laboratory experiments. Instantaneous solutions in configuration (1) show that the sheet's dimensionless `stiffness' S determines whether the slab's sinking speed is controlled by the viscosity of the ambient mantle (S < 1) or the viscosity of the sheet itself (S > 10). Time-dependent solutions with tracers in configuration (2) demonstrate a partial return flow around the leading edge of a retreating slab and return flow around its sides. The extra `edge drag' exerted by the flow around the sides causes transverse deformation of the slab, and makes the sinking speed of a 3-D slab up to 40% less than that of a 2-D slab. A systematic investigation of the slab's interaction with the bottom boundary as a function of η2/η1 and H/h delineates a rich regime diagram of different subduction modes (trench retreating, slab folding, trench advancing) and reveals a new `advancing-folding' mode in which slab folding is preceded by advancing trench motion. The solutions demonstrate that mode selection is controlled by the dip of the leading edge of the slab at the time when it first encounters the bottom boundary.
Computational Analysis of the Transonic Dynamics Tunnel Using FUN3D
NASA Technical Reports Server (NTRS)
Chwalowski, Pawel; Quon, Eliot; Brynildsen, Scott E.
2016-01-01
This paper presents results from an exploratory two-year effort of applying Computational Fluid Dynamics (CFD) to analyze the empty-tunnel flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The TDT is a continuous-flow, closed circuit, 16- x 16-foot slotted-test-section wind tunnel, with capabilities to use air or heavy gas as a working fluid. In this study, experimental data acquired in the empty tunnel using the R-134a test medium was used to calibrate the computational data. The experimental calibration data includes wall pressures, boundary-layer profiles, and the tunnel centerline Mach number profiles. Subsonic and supersonic flow regimes were considered, focusing on Mach 0.5, 0.7 and Mach 1.1 in the TDT test section. This study discusses the computational domain, boundary conditions, and initial conditions selected and the resulting steady-state analyses using NASA's FUN3D CFD software.
RV functional imaging: 3-D echo-derived dynamic geometry and flow field simulations.
Pasipoularides, Ares D; Shu, Ming; Womack, Michael S; Shah, Ashish; Von Ramm, Olaf; Glower, Donald D
2003-01-01
We describe a novel functional imaging approach for quantitative analysis of right ventricular (RV) blood flow patterns in specific experimental animals (or humans) using real-time, three-dimensional (3-D) echocardiography (RT3D). The method is independent of the digital imaging modality used. It comprises three parts. First, a semiautomated segmentation aided by intraluminal contrast medium locates the RV endocardial surface. Second, a geometric scheme for dynamic RV chamber reconstruction applies a time interpolation procedure to the RT3D data to quantify wall geometry and motion at 400 Hz. A volumetric prism method validated the dynamic geometric reconstruction against simultaneous sonomicrometric canine measurements. Finally, the RV endocardial border motion information is used for mesh generation on a computational fluid dynamics solver to simulate development of the early RV diastolic inflow field. Boundary conditions (tessellated endocardial surface nodal velocities) for the solver are directly derived from the endocardial geometry and motion information. The new functional imaging approach may yield important kinematic information on the distribution of instantaneous velocities in the RV diastolic flow field of specific normal or diseased hearts. PMID:12388220
NASA Astrophysics Data System (ADS)
Druzgalski, Clara; Mani, Ali
2014-11-01
We have investigated the transport dynamics of an electrokinetic instability that occurs when ions are driven from bulk fluids to ion-selective membranes due to externally applied electric fields. This phenomenon is relevant to a wide range of electrochemical applications including electrodialysis for fresh water production. Using data from our 3D DNS, we show how electroconvective instability, arising from concentration polarization, results in a chaotic flow that significantly alters the net ion transport rate across the membrane surface. The 3D DNS results, which fully resolve the spatiotemporal scales including the electric double layers, enable visualization of instantaneous snapshots of current density directly on the membrane surface, as well as analysis of transport statistics such as concentration variance and fluctuating advective fluxes. Furthermore, we present a full spectral analysis revealing broadband spectra in both concentration and flow fields and deduce the key parameter controlling the range of contributing scales.
A 3D GCL compatible cell-centered Lagrangian scheme for solving gas dynamics equations
NASA Astrophysics Data System (ADS)
Georges, Gabriel; Breil, Jérôme; Maire, Pierre-Henri
2016-01-01
Solving the gas dynamics equations under the Lagrangian formalism enables to simulate complex flows with strong shock waves. This formulation is well suited to the simulation of multi-material compressible fluid flows such as those encountered in the domain of High Energy Density Physics (HEDP). These types of flows are characterized by complex 3D structures such as hydrodynamic instabilities (Richtmyer-Meshkov, Rayleigh-Taylor, etc.). Recently, the 3D extension of different Lagrangian schemes has been proposed and appears to be challenging. More precisely, the definition of the cell geometry in the 3D space through the treatment of its non-planar faces and the limiting of a reconstructed field in 3D in the case of a second-order extension are of great interest. This paper proposes two new methods to solve these problems. A systematic and symmetric geometrical decomposition of polyhedral cells is presented. This method enables to define a discrete divergence operator leading to the respect of the Geometric Conservation Law (GCL). Moreover, a multi-dimensional minmod limiter is proposed. This new limiter constructs, from nodal gradients, a cell gradient which enables to ensure the monotonicity of the numerical solution even in presence of strong discontinuity. These new ingredients are employed into a cell-centered Lagrangian scheme. Robustness and accuracy are assessed against various representative test cases.
Conservation of Fluid Mass and Energy by RELAP5-3D during a SBLOCA
Cliff B. Davis
2009-08-01
Mass and energy balances were performed to check the accuracy of RELAP5-3D’s solution during a loss-of-coolant accident initiated by a small break in a typical pressurized water reactor. Mass and energy balances were performed for the combined liquid and gas phases and the gas phase by itself. The analysis showed that RELAP5-3D adequately conserved mass and energy for the combined fluid and the gas phase.
Navier-Stokes Neutral and Plasma Fluid Modelling in 3D
Riemann, J; Borchardt, M; Schneider, R; Mutzke, A; Rognlien, T; Umansky, M
2004-05-17
The 3D finite volume transport code BoRiS is applied to a system of coupled plasma and neutral fluid equations in a slab. Demonstrating easy implementation of new equations, a new parallel BoRiS version is tested on three different models for the neutral fluid - diffusive, parallel Navier-Stokes and full Navier-Stokes - and the results are compared to each other. Typical effects like density enhancement by ionization of recycled neutrals in front of a target plate can be seen and differences are linked to the neutral models in use.
Dynamics of tokamak plasma surface current in 3D ideal MHD model
NASA Astrophysics Data System (ADS)
Galkin, Sergei A.; Svidzinski, V. A.; Zakharov, L. E.
2013-10-01
Interest in the surface current which can arise on perturbed sharp plasma vacuum interface in tokamaks was recently generated by a few papers (see and references therein). In dangerous disruption events with plasma-touching-wall scenarios, the surface current can be shared with the wall leading to the strong, damaging forces acting on the wall A relatively simple analytic definition of δ-function surface current proportional to a jump of tangential component of magnetic field nevertheless leads to a complex computational problem on the moving plasma-vacuum interface, requiring the incorporation of non-linear 3D plasma dynamics even in one-fluid ideal MHD. The Disruption Simulation Code (DSC), which had recently been developed in a fully 3D toroidal geometry with adaptation to the moving plasma boundary, is an appropriate tool for accurate self-consistent δfunction surface current calculation. Progress on the DSC-3D development will be presented. Self-consistent surface current calculation under non-linear dynamics of low m kink mode and VDE will be discussed. Work is supported by the US DOE SBIR grant #DE-SC0004487.
North Cascadia heat flux and fluid flow from gas hydrates: Modeling 3-D topographic effects
NASA Astrophysics Data System (ADS)
Li, Hong-lin; He, Tao; Spence, George D.
2014-01-01
The bottom-simulating reflector (BSR) of gas hydrate is well imaged from two perpendicular seismic grids in the region of a large carbonate mound, informally called Cucumber Ridge off Vancouver Island. We use a new method to calculate 3-D heat flow map from the BSR depths, in which we incorporate 3-D topographic corrections after calibrated by the drilling results from nearby (Integrated) Ocean Drilling Program Site 889 and Site U1327. We then estimate the associated fluid flow by relating it to the topographically corrected heat flux anomalies. In the midslope region, a heat flux anomaly of 1 mW/m2 can be associated with an approximate focused fluid flow rate of 0.09 mm/yr. Around Cucumber Ridge, high rates of focused fluid flow were observed at steep slopes with values more than double the average regional diffusive fluid discharge rate of 0.56 mm/yr. As well, in some areas of relatively flat seafloor, the focused fluid flow rates still exceeded 0.5 mm/yr. On the seismic lines the regions of focused fluid flow were commonly associated with seismic blanking zones above the BSR and sometimes with strong reflectors below the BSR, indicating that the faults/fractures provide high-permeability pathways for fluids to carry methane from BSR depths to the seafloor. These high fluid flow regions cover mostly the western portion of our area with gas hydrate concentration estimations of ~6% based on empirical correlations from Hydrate Ridge in south off Oregon, significantly higher than previously recognized values of ~2.5% in the eastern portion determined from Site U1327.
NASA Technical Reports Server (NTRS)
1989-01-01
An overview of computational fluid dynamics (CFD) activities at the Langley Research Center is given. The role of supercomputers in CFD research, algorithm development, multigrid approaches to computational fluid flows, aerodynamics computer programs, computational grid generation, turbulence research, and studies of rarefied gas flows are among the topics that are briefly surveyed.
ERIC Educational Resources Information Center
Bird, R. Byron
1980-01-01
Problems in polymer fluid dynamics are described, including development of constitutive equations, rheometry, kinetic theory, flow visualization, heat transfer studies, flows with phase change, two-phase flow, polymer unit operations, and drag reduction. (JN)
Dynamical Lorentz symmetry breaking in 3D and charge fractionalization
Charneski, B.; Gomes, M.; Silva, A. J. da; Mariz, T.; Nascimento, J. R.
2009-03-15
We analyze the breaking of Lorentz invariance in a 3D model of fermion fields self-coupled through four-fermion interactions. The low-energy limit of the theory contains various submodels which are similar to those used in the study of graphene or in the description of irrational charge fractionalization.
3D numerical simulations of vesicle and inextensible capsule dynamics
NASA Astrophysics Data System (ADS)
Farutin, Alexander; Biben, Thierry; Misbah, Chaouqi
2014-10-01
Vesicles are locally-inextensible fluid membranes, capsules are endowed with in-plane shear elasticity mimicking the cytoskeleton of red blood cells (RBCs), but are extensible, while RBCs are inextensible. We use boundary integral (BI) methods based on the Green function techniques to model and solve numerically their dynamics. We regularize the single layer integral by subtraction of exact identities for the terms involving the normal and the tangential components of the force. The stability and precision of BI calculation is enhanced by taking advantage of additional quadrature nodes located in vertices of an auxiliary mesh, constructed by a standard refinement procedure from the main mesh. We extend the partition of unity technique to boundary integral calculation on triangular meshes. The proposed algorithm offers the same treatment of near-singular integration regardless whether the source and the target points belong to the same surface or not. Bending forces are calculated by using expressions derived from differential geometry. Membrane incompressibility is handled by using two penalization parameters per suspended entity: one for deviation of the global area from prescribed value and another for the sum of squares of local strains defined on each vertex. Extensible or inextensible capsules, a model of RBC, are studied by storing the position in the reference configuration for each vertex. The elastic force is then calculated by direct variation of the elastic energy. Various nonequilibrium physical examples on vesicles and capsules will be presented and the convergence and precision tests highlighted. Overall, a good convergence is observed with numerical error inversely proportional to the number of vertices used for surface discretization, the highest order of convergence allowed by piece-wise linear interpolation of the surface.
Computational fluid dynamic control
NASA Technical Reports Server (NTRS)
Hartley, Tom T.; Deabreu-Garcia, Alex
1989-01-01
A general technique is presented for modeling fluid, or gas, dynamic systems specifically for the development of control systems. The numerical methods which are generally used in computational fluid dynamics are borrowed to create either continuous-time or discrete-time models of the particular fluid system. The resulting equations can be either left in a nonlinear form, or easily linearized about an operating point. As there are typically very many states in these systems, the usual linear model reduction methods can be used on them to allow a low-order controller to be designed. A simple example is given which typifies many internal flow control problems. The resulting control is termed computational fluid dynamic control.
DREAM3D simulations of inner-belt dynamics
Cunningham, Gregory Scott
2015-05-26
A 1973 paper by Lyons and Thorne explains the two-belt structure for electrons in the inner magnetosphere as a balance between inward radial diffusion and loss to the atmosphere, where the loss to the atmosphere is enabled by pitch-angle scattering from Coulomb and wave-particle interactions. In the 1973 paper, equilibrium solutions to a decoupled set of 1D radial diffusion equations, one for each value of the first invariant of motion, μ, were computed to produce the equilibrium two-belt structure. Each 1D radial diffusion equation incorporated an L-and μ-dependent `lifetime' due to the Coulomb and wave-particle interactions. This decoupling of the problem is appropriate under the assumption that radial diffusion is slow in comparison to pitch-angle scattering. However, for some values of μ and L the lifetime associated with pitch-angle scattering is comparable to the timescale associated with radial diffusion, suggesting that the true equilibrium solutions might reflect `coupled modes' involving pitch-angle scattering and radial diffusion and thus requiring a 3D diffusion model. In the work we show here, we have computed the equilibrium solutions using our 3D diffusion model, DREAM3D, that allows for such coupling. We find that the 3D equilibrium solutions are quite similar to the solutions shown in the 1973 paper when we use the same physical models for radial diffusion and pitch-angle scattering from hiss. However, we show that the equilibrium solutions are quite sensitive to various aspects of the physics model employed in the 1973 paper that can be improved, suggesting that additional work needs to be done to understand the two-belt structure.
Salinity effects on cracking morphology and dynamics in 3-D desiccating clays
NASA Astrophysics Data System (ADS)
DeCarlo, Keita F.; Shokri, Nima
2014-04-01
Saline conditions induce not only chemical but physical changes in swelling clays, and have a significant influence on the crack dynamics and morphology of desiccating clays. In this study, we used X-ray microtomography to experimentally investigate the effects of sodium chloride on the morphology and dynamics of desiccation cracks in three-dimensional mixtures of sand-bentonite slurry under varying rheological conditions. Rectangular glass containers were packed with slurries of different salt concentrations, with the top boundary exposed to air for evaporation. The growth and propagation of the cracking network that subsequently formed was visualized in 3-D at multiple intervals. The characterization of cracking and branching behavior shows a high extent of localized surficial crack networks at low salinity, with a transition to less extensive but more centralized crack networks with increased salinity. The observed behavior was described in the context of the physicochemical properties of the montmorillonite clay, where shifts from an "entangled" (large platelet spacing, small pore structure) to a "stacked" (small platelet spacing, open pore structure) network influence fluid distribution and thus extent of cracking and branching behavior. This is further corroborated by vertical profiles of water distribution, which shows localized desiccation fronts that shift to uniform desaturation with increasing salt concentration. Our results provide new insights regarding the formation, dynamics, and patterns of desiccation cracks formed during evaporation from 3-D saline clay structures, which will be useful in hydrological applications including water management, land surface evaporation, and subsurface contaminant transport.
Fluid flow pathways study from the 3D seismic data offshore southwestern Taiwan
NASA Astrophysics Data System (ADS)
Chen, L.; Chi, W. C.; Chiang, H. T.; Lin, S.
2014-12-01
3D seismic reflection data provide detailed information on the physical properties of the crust, which can be used for hydrocarbon exploration. Recently, scientists from Taiwan and Germany are collaborating on a project to use a portable 3D seismic system, called P-Cable, to study gas hydrates offshore southwest Taiwan. We have collected 3 cubes, covering the active and passive margins. At these three sites, there is a wide-spread bottom-simulating reflector (BSR). We use the BSR to study the shallow thermal structures of these prospect sites, and use the temperature field information to study fluid migration patterns. We have also done in-situ heat flow measurements, and found similar results, showing focused fluid flow migrations in some pathways. Some of the high temperature fields also correlate with gas chimneys found through seismic attribute analyses. Preliminary results show that there might be active fluid migration above the BSR in the gas hydrate stability zone. In September and October of 2014, we will collect additional P-Cable datasets to be incorporated into this study. Such results will be used to evaluate some proposed sites for future drilling programs.
NASA Astrophysics Data System (ADS)
Krebs, Isabel; Jardin, Stephen C.; Igochine, Valentin; Guenter, Sibylle; Hoelzl, Matthias; ASDEX Upgrade Team
2014-10-01
We study sawtooth reconnection in ASDEX Upgrade tokamak plasmas by means of 3D non-linear two-fluid MHD simulations in toroidal geometry using the high-order finite element code M3D-C1. Parameters and equilibrium of the simulations are based on typical sawtoothing ASDEX Upgrade discharges. The simulation results are compared to features of the experimental observations such as the sawtooth crash time and frequency, the evolution of the safety factor profile and the 3D evolution of the temperature. 2D ECE imaging measurements during sawtooth crashes in ASDEX Upgrade indicate that the heat is transported out of the core through a narrow poloidally localized region. We investigate if incomplete sawtooth reconnection can be seen in the simulations which is suggested by soft X-ray tomography measurements in ASDEX Upgrade showing that an (m = 1, n = 1) perturbation is typically observed to survive the sawtooth crash and approximately maintain its radial position.
Relativistic fluid dynamics. Proceedings.
NASA Astrophysics Data System (ADS)
Anile, A. M.; Choquet-Bruhat, Y.
Contents: 1. Covariant theory of conductivity in ideal fluid or solid media (B. Carter). 2. Hamiltonian techniques for relativistic fluid dynamics and stability theory (D. D. Holm). 3. Covariant fluid mechanics and thermodynamics: an introduction (W. Israel). 4. Relativistic plasmas (H. Weitzner). 5. An improved relativistic warm plasma model (A. M. Anile, S. Pennisi). 6. Relativistic extended thermodynamics II (I. Müller). 7. Relativistic extended thermodynamics: general assumptions and mathematical procedure (T. Ruggeri). 8. Relativistic hydrodynamics and heavy ion reactions (D. Strottman). 9. Some problems in relativistic hydrodynamics (C. G. van Weert).
3D topographic correction of the BSR heat flow and detection of focused fluid flow
NASA Astrophysics Data System (ADS)
He, Tao; Li, Hong-Lin; Zou, Chang-Chun
2014-06-01
The bottom-simulating reflector (BSR) is a seismic indicator of the bottom of a gas hydrate stability zone. Its depth can be used to calculate the seafloor surface heat flow. The calculated BSR heat flow variations include disturbances from two important factors: (1) seafloor topography, which focuses the heat flow over regions of concave topography and defocuses it over regions of convex topography, and (2) the focused warm fluid flow within the accretionary prism coming from depths deeper than BSR. The focused fluid flow can be detected if the contribution of the topography to the BSR heat flow is removed. However, the analytical equation cannot solve the topographic effect at complex seafloor regions. We prove that 3D finite element method can model the topographic effect on the regional background heat flow with high accuracy, which can then be used to correct the topographic effect and obtain the BSR heat flow under the condition of perfectly flat topography. By comparing the corrected BSR heat flow with the regional background heat flow, focused fluid flow regions can be detected that are originally too small and cannot be detected using present-day equipment. This method was successfully applied to the midslope region of northern Cascadia subducting margin. The results suggest that the Cucumber Ridge and its neighboring area are positive heat flow anomalies, about 10%-20% higher than the background heat flow after 3D topographic correction. Moreover, the seismic imaging associated the positive heat flow anomaly areas with seabed fracture-cavity systems. This suggests flow of warm gas-carrying fluids along these high-permeability pathways, which could result in higher gas hydrate concentrations.
How effective can optical-CT 3D dosimetry be without refractive fluid matching?
NASA Astrophysics Data System (ADS)
Rankine, L.; Oldham
2013-06-01
Achieving accurate optical CT 3D dosimetry without the use of viscous refractive index (RI) matching fluids would greatly increase convenience. Software has been developed to simulate optical CT 3D dosimetry for a range of scanning configurations including parallel-beam, point and converging light sources. For each configuration the efficacy of 3 refractive media were investigated: air, water, and a fluid closely matched to Presage (RI = 1.00, 1.33 and 1.49 respectively). The results revealed that the useable radius of the dosimeter (i.e. where data was within 2% of truth) reduced to 68% for water-matching, and 31% for dry-scanning in air. Point source incident ray geometry produced slightly more favourable results, although variation between the three geometries was relatively small. The required detector size however, increased by a factor six for dry-scanning, introducing cost penalties. For applications where dose information is not required in the periphery, some dry and low-viscous matching configurations may be feasible.
3D Global Two-Fluid Simulations of Turbulence in LAPD
NASA Astrophysics Data System (ADS)
Fisher, Dustin; Rogers, Barrett; Ricci, Paolo
2012-10-01
3D global two-fluid simulations are presented in an ongoing effort to identify and understand the physics of instabilities that arise in the Large Plasma Device (LAPD) at UCLA's Basic Science Facility. The LAPD, with its wide range of tunable parameters and device configurations, is ideally suited for studying space and laboratory plasmas. Moreover, the highly detailed and reproducible measurements of the LAPD lend themselves amicably to comparisons with simulations. Ongoing modeling is done using a modified version of the Global Braginskii Solver (GBS) [1] that models the plasma from source to edge region in a fully 3D two-fluid code. The reduced Braginskii equations are solved on a field-aligned grid using a finite difference method and 4th order Runge-Kutta time stepping and are parallelized on Dartmouth's Discovery cluster. Recent progress has been made to account for the thermionic cathode emission of fast electrons at the source, the axial dependence of the plasma source, and it is now possible to vary the potential on the front and side walls. Preliminary results, seen from the density and temperature profiles, show that the low frequency Kelvin Helmholtz instability still dominates the turbulence in the device.[4pt] [1] B. Rogers and P. Ricci. Phys. Rev. Lett. 104:225002, 2010
Adhikari, Deepak; Gemmell, Brad J; Hallberg, Michael P; Longmire, Ellen K; Buskey, Edward J
2015-11-01
We describe an automated, volumetric particle image velocimetry (PIV) and tracking method that measures time-resolved, 3D zooplankton trajectories and surrounding volumetric fluid velocity fields simultaneously and non-intrusively. The method is demonstrated for groups of copepods flowing past a wall-mounted cylinder. We show that copepods execute escape responses when subjected to a strain rate threshold upstream of a cylinder, but the same threshold range elicits no escape responses in the turbulent wake downstream. The method was also used to document the instantaneous slip velocity of zooplankton and the resulting differences in trajectory between zooplankton and non-inertial fluid particles in the unsteady wake flow, showing the method's capability to quantify drift for both passive and motile organisms in turbulent environments. Applications of the method extend to any group of organisms interacting with the surrounding fluid environment, where organism location, larger-scale eddies and smaller-scale fluid deformation rates can all be tracked and analyzed. PMID:26486364
Code verification for unsteady 3-D fluid-solid interaction problems
NASA Astrophysics Data System (ADS)
Yu, Kintak Raymond; Étienne, Stéphane; Hay, Alexander; Pelletier, Dominique
2015-12-01
This paper describes a procedure to synthesize Manufactured Solutions for Code Verification of an important class of Fluid-Structure Interaction (FSI) problems whose behaviors can be modeled as rigid body vibrations in incompressible fluids. We refer this class of FSI problems as Fluid-Solid Interaction problems, which can be found in many practical engineering applications. The methodology can be utilized to develop Manufactured Solutions for both 2-D and 3-D cases. We demonstrate the procedure with our numerical code. We present details of the formulation and methodology. We also provide the reasonings behind our proposed approach. Results from grid and time step refinement studies confirm the verification of our solver and demonstrate the versatility of the simple synthesis procedure. In addition, the results also demonstrate that the modified decoupled approach to verify flow problems with high-order time-stepping schemes can be employed equally well to verify code for multi-physics problems (here, those of the Fluid-Solid Interaction) when the numerical discretization is based on the Method of Lines.
The vibrational dynamics of 3D HOCl above dissociation
Lin, Yi-Der; Reichl, L. E.; Jung, Christof
2015-03-28
We explore the classical vibrational dynamics of the HOCl molecule for energies above the dissociation energy of the molecule. Above dissociation, we find that the classical dynamics is dominated by an invariant manifold which appears to stabilize two periodic orbits at energies significantly above the dissociation energy. These stable periodic orbits can hold a large number of quantum states and likely can support a significant quasibound state of the molecule, well above the dissociation energy. The classical dynamics and the lifetime of quantum states on the invariant manifold are determined.
The vibrational dynamics of 3D HOCl above dissociation
NASA Astrophysics Data System (ADS)
Lin, Yi-Der; Reichl, L. E.; Jung, Christof
2015-03-01
We explore the classical vibrational dynamics of the HOCl molecule for energies above the dissociation energy of the molecule. Above dissociation, we find that the classical dynamics is dominated by an invariant manifold which appears to stabilize two periodic orbits at energies significantly above the dissociation energy. These stable periodic orbits can hold a large number of quantum states and likely can support a significant quasibound state of the molecule, well above the dissociation energy. The classical dynamics and the lifetime of quantum states on the invariant manifold are determined.
NASA Astrophysics Data System (ADS)
Chung, T. J.
2002-03-01
Computational fluid dynamics (CFD) techniques are used to study and solve complex fluid flow and heat transfer problems. This comprehensive text ranges from elementary concepts for the beginner to state-of-the-art CFD for the practitioner. It discusses and illustrates the basic principles of finite difference (FD), finite element (FE), and finite volume (FV) methods, with step-by-step hand calculations. Chapters go on to examine structured and unstructured grids, adaptive methods, computing techniques, and parallel processing. Finally, the author describes a variety of practical applications to problems in turbulence, reacting flows and combustion, acoustics, combined mode radiative heat transfer, multiphase flows, electromagnetic fields, and relativistic astrophysical flows. Students and practitioners--particularly in mechanical, aerospace, chemical, and civil engineering--will use this authoritative text to learn about and apply numerical techniques to the solution of fluid dynamics problems.
Guidelines in the experimental validation of a 3D heat and fluid flow model of keyhole laser welding
NASA Astrophysics Data System (ADS)
Courtois, Mickael; Carin, Muriel; Le Masson, Philippe; Gaied, Sadok; Balabane, Mikhaël
2016-04-01
During the past few years, numerous sophisticated models have been proposed to predict in a self-consistent way the dynamics of the keyhole, together with the melt pool and vapor jet. However, these models are only partially compared to experimental data, so the reliability of these models is questionable. The present paper aims to propose a more complete experimental set-up in order to validate the most relevant results calculated by these models. A complete heat transfer and fluid flow three-dimensional (3D) model is first proposed in order to describe laser welding in keyhole regimes. The interface is tracked with a level set method and fluid flows are calculated in liquid and gas. The mechanisms of recoil pressure and keyhole creation are highlighted in a fusion line configuration chosen as a reference. Moreover, a complete validation of the model is proposed with guidelines on the variables to observe. Numerous comparisons with dedicated experiments (thermocouples, pyrometry, high-speed camera) are proposed to estimate the validity of the model. In addition to traditional geometric measurements, the main variables calculated, temperatures, and velocities in the melt pool are at the center of this work. The goal is to propose a reference validation for complex 3D models proposed over the last few years.
Environmental Fluid Dynamics Code
The Environmental Fluid Dynamics Code (EFDC)is a state-of-the-art hydrodynamic model that can be used to simulate aquatic systems in one, two, and three dimensions. It has evolved over the past two decades to become one of the most widely used and technically defensible hydrodyn...
First 3-D simulations of meteor plasma dynamics and turbulence
NASA Astrophysics Data System (ADS)
Oppenheim, Meers M.; Dimant, Yakov S.
2015-02-01
Millions of small but detectable meteors hit the Earth's atmosphere every second, creating trails of hot plasma that turbulently diffuse into the background atmosphere. For over 60 years, radars have detected meteor plasmas and used these signals to infer characteristics of the meteoroid population and upper atmosphere, but, despite the importance of meteor radar measurements, the complex processes by which these plasmas evolve have never been thoroughly explained or modeled. In this paper, we present the first fully 3-D simulations of meteor evolution, showing meteor plasmas developing instabilities, becoming turbulent, and inhomogeneously diffusing into the background ionosphere. These instabilities explain the characteristics and strength of many radar observations, in particular the high-resolution nonspecular echoes made by large radars. The simulations reveal how meteors create strong electric fields that dig out deep plasma channels along the Earth's magnetic fields. They also allow researchers to explore the impacts of the intense winds and wind shears, commonly found at these altitudes, on meteor plasma evolution. This study will allow the development of more sophisticated models of meteor radar signals, enabling the extraction of detailed information about the properties of meteoroid particles and the atmosphere.
Quantifying cellular interaction dynamics in 3-D fluorescence microscopy data
Klauschen, Frederick; Ishii, Masaru; Qi, Hai; Bajénoff, Marc; Egen, Jackson G.; Germain, Ronald N.; Meier-Schellersheim, Martin
2012-01-01
The wealth of information available from advanced fluorescence imaging techniques used to analyze biological processes with high spatial and temporal resolution calls for high-throughput image analysis methods. Here, we describe a fully automated approach to analyzing cellular interaction behavior in 3-D fluorescence microscopy images. As example application we present the analysis of drug-induced and S1P1-knock-out-related changes in bone-osteoclast interactions. Moreover, we apply our approach to images showing the spatial association of dendritic cells with the fibroblastic reticular cell network within lymph nodes and to microscopy data about T-B lymphocyte synapse formation. Such analyses that yield important information about the molecular mechanisms determining cellular interaction behavior would be very difficult to perform with approaches that rely on manual/semi-automated analyses. This protocol integrates adaptive threshold segmentation, object detection, adaptive color channel merging and neighborhood analysis and permits rapid, standardized, quantitative analysis and comparison of the relevant features in large data sets. PMID:19696749
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
NASA Astrophysics Data System (ADS)
Ogilvie, Gordon I.
2016-06-01
These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is `frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.
NASA Astrophysics Data System (ADS)
Ogilvie, Gordon I.
2016-06-01
> These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is `frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.
3D visualization of deformation structures and potential fluid pathways at the Grimsel Test Site
NASA Astrophysics Data System (ADS)
Schneeberger, Raphael; Kober, Florian; Berger, Alfons; Spillmann, Thomas; Herwegh, Marco
2015-04-01
Knowledge on the ability of fluids to infiltrate subsurface rocks is of major importance for underground constructions, geothermal or radioactive waste disposal projects. In this study, we focus on the characterization of water infiltration pathways, their 3D geometries and origins. Based on surface and subsurface mapping in combination with drill core data, we developed by the use of MoveTM (Midland Valley Exploration Ltd.) a 3D structural model of the Grimsel Test Site (GTS). GTS is an underground laboratory operated by NAGRA, the Swiss organisation responsible for the management of nuclear waste. It is located within a suite of post-Variscan magmatic bodies comprising former granitic and granodioritic melts, which are dissected by mafic and aplitic dikes. During Alpine orogeny, the suite was tectonically overprinted within two stages of ductile deformation (Wehrens et al., in prep.) followed by brittle overprint of some of the shear zones during the retrograde exhumation history. It is this brittle deformation, which controls today's water infiltration network. However, the associated fractures, cataclasites and fault gouges are controlled themselves by aforementioned pre-existing mechanical discontinuities, whose origin ranges back as far as to the magmatic stage. For example, two sets of vertically oriented mafic dikes (E-W and NW-SE striking) and compositional heterogeneities induced by magmatic segregation processes in the plutonic host rocks served as nucleation sites for Alpine strain localization. Subsequently, NE-SW, E-W and NW-SE striking ductile shear zones were formed, in combination with high temperature fracturing while dissecting the host rocks in a complex 3D pattern (Wehrens et al, in prep.). Whether the ductile shear zones have been subjected to brittle reactivation and can serve as infiltration pathways or not, depends strongly on their orientations with respect to the principal stress field. Especially where deformation structures intersect
Nanoscale 3D distribution of low melt and fluid fractions in mantle rocks
NASA Astrophysics Data System (ADS)
Gardes, Emmanuel; Morales, Luiz; Heinrich, Wilhelm; Sifre, David; Hashim, Leila; Gaillard, Fabrice; Katharina, Marquardt
2016-04-01
The presence of melts or fluids in the intergranular medium of rocks strongly influences their bulk physico-chemical properties (e.g. mass transport and chemical reactivity, electrical conductivity, seismic wave velocity, etc). Actually, the effects can be so large that only small melt or fluid fractions must sometimes be involved for explaining mantle geophysical discontinuities and anomalies. The investigation of the distribution of such small fractions in the intergranular medium of mantle rocks is therefore crucial for relating them to bulk and large scale properties. However, it involves submicrometric structures which are hardly characterizable using conventional techniques. Here we present how the FIB-SEM-STEM microscope can be used to produce 3D imaging at unequalled resolution. We show that low melt and fluid fractions can form films as thin as 20 nm at olivine grain boundaries, and that they can modify the physico-chemical properties of mantle rocks by orders of magnitude. The fine relationships between films at grain boundaries, tubules at triple junctions and pockets at grain corners can be explored, and appear to be complex and to differ from usual visions.
Introducing a New 3D Dynamical Model for Barred Galaxies
NASA Astrophysics Data System (ADS)
Jung, Christof; Zotos, Euaggelos E.
2015-11-01
The regular or chaotic dynamics of an analytical realistic three dimensional model composed of a spherically symmetric central nucleus, a bar and a flat disk is investigated. For describing the properties of the bar, we introduce a new simple dynamical model and we explore the influence on the character of orbits of all the involved parameters of it, such as the mass and the scale length of the bar, the major semi-axis and the angular velocity of the bar, as well as the energy. Regions of phase space with ordered and chaotic motion are identified in dependence on these parameters and for breaking the rotational symmetry. First, we study in detail the dynamics in the invariant plane z = pz = 0 using the Poincaré map as a basic tool and then study the full three-dimensional case using the Smaller Alignment index method as principal tool for distinguishing between order and chaos. We also present strong evidence obtained through the numerical simulations that our new bar model can realistically describe the formation and the evolution of the observed twin spiral structure in barred galaxies.
Fluid force and static symmetry breaking modes of 3D bluff bodies.
NASA Astrophysics Data System (ADS)
Cadot, Olivier; Evrard, Antoine; DFA Team
2015-11-01
A cavity at the base of the squareback Ahmed model at Re =6.106 is able to reduce the base suction by 18% and the drag coefficient by 9%, while the flow at the separation remains unaffected. Instantaneous pressure measurements at the body base, fluid force measurements and wake velocity measurements are investigated varying the cavity depth from 0 to 35% of the base height. Due to the reflectional symmetry of the rectangular base, there are two Reflectional Symmetry Breaking (RSB) mirror modes present in the natural wake that switch from one to the other randomly in accordance with the recent findings of Grandemange et al. (2013). It is shown that these modes exhibit an energetic 3D static vortex system close to the base of the body. A sufficiently deep cavity is able to stabilize the wake toward a symmetry preserved wake, thus suppressing the RSB modes and leading to a weaker elliptical toric recirculation. The stabilization can be modeled with a Langevin equation. The plausible mechanism for drag reduction with the base cavity is based on the interaction of the static 3D vortex system of the RSB modes with the base and their suppression by stabilization. There are some strong evidences that this mechanism may be generalized to axisymmetric bodies with base cavity.
Influence of Young's moduli in 3D fluid-structure coupled models of the human cochlea
NASA Astrophysics Data System (ADS)
Böhnke, Frank; Semmelbauer, Sebastian; Marquardt, Torsten
2015-12-01
The acoustic wave propagation in the human cochlea was studied using a tapered box-model with linear assumptions respective to all mechanical parameters. The discretisation and evaluation is conducted by a commercial finite element package (ANSYS). The main difference to former models of the cochlea was the representation of the basilar membrane by a 3D elastic solid. The Young's moduli of this solid were modified to study their influence on the travelling wave. The lymph in the scala vestibuli and scala tympani was represented by a viscous and nearly incompressible fluid finite element approach. Our results show the maximum displacement for f = 2kHz at half of the length of the cochlea in accordance with former experiments. For low frequencies f <200 Hz nearly zero phase shifts were found, whereas for f =1 kHz it reaches values up to -12 cycles depending on the degree of orthotropy.
Kent, J C; Eaton, A R
1982-03-01
A new technique has been developed for studies of fluid motion within the cylinder of a reciprocating piston engine during the air induction process. Helium-filled bubbles, serving as neutrally buoyant flow tracer particles, enter the cylinder along with the inducted air charge. The bubble motion is recorded by stereo cine photography through the transparent cylinder of a specially designed research engine. Quantitative data on the 3-D velocity field generated during induction is obtained from frame-to-frame analysis of the stereo images, taking into account refraction of the rays due to the transparent cylinder. Other applications for which this technique appears suitable include measurements of velocity fields within intake ports and flow-field dynamics within intake manifolds of multicylinder engines. PMID:20372559
Computational fluid dynamics research
NASA Technical Reports Server (NTRS)
Chandra, Suresh; Jones, Kenneth; Hassan, Hassan; Mcrae, David Scott
1992-01-01
The focus of research in the computational fluid dynamics (CFD) area is two fold: (1) to develop new approaches for turbulence modeling so that high speed compressible flows can be studied for applications to entry and re-entry flows; and (2) to perform research to improve CFD algorithm accuracy and efficiency for high speed flows. Research activities, faculty and student participation, publications, and financial information are outlined.
Perspectives in Fluid Dynamics
NASA Astrophysics Data System (ADS)
Batchelor, G. K.; Moffatt, H. K.; Worster, M. G.
2002-12-01
With applications ranging from modelling the environment to automotive design and physiology to astrophysics, conventional textbooks cannot hope to give students much information on what topics in fluid dynamics are currently being researched, or how to choose between them. This book rectifies matters. It consists of eleven chapters that introduce and review different branches of the subject for graduate-level courses, or for specialists seeking introductions to other areas. Hb ISBN (2001): 0-521-78061-6
Cotranslational processing mechanisms: towards a dynamic 3D model.
Giglione, Carmela; Fieulaine, Sonia; Meinnel, Thierry
2009-08-01
Recent major advances have been made in understanding how cotranslational events are achieved in the course of protein biosynthesis. Specifically, several studies have shed light into the dynamic process of how nascent chains emerging from the ribosome are supported by protein biogenesis factors to ensure both processing and folding mechanisms. To take into account the awareness that coordination is needed, a new 'concerted model' recently proposed simultaneous action of both processes on the ribosome. In the model, any emerging nascent chain is first encountered by the chaperone trigger factor (TF), which forms an open cradle underneath the ribosomal exit tunnel. This cradle serves as a passive router that channels the nascent chains to the first cotranslational event, the proteolysis event performed by the N-terminal methionine excision machinery. Although fascinating, this model clearly raises more questions than it answers. Does the data used to develop this model stand up to scrutiny and, if not, what are the alternative mechanisms that the data suggest? PMID:19647435
3D video-based deformation measurement of the pelvis bone under dynamic cyclic loading
2011-01-01
Background Dynamic three-dimensional (3D) deformation of the pelvic bones is a crucial factor in the successful design and longevity of complex orthopaedic oncological implants. The current solutions are often not very promising for the patient; thus it would be interesting to measure the dynamic 3D-deformation of the whole pelvic bone in order to get a more realistic dataset for a better implant design. Therefore we hypothesis if it would be possible to combine a material testing machine with a 3D video motion capturing system, used in clinical gait analysis, to measure the sub millimetre deformation of a whole pelvis specimen. Method A pelvis specimen was placed in a standing position on a material testing machine. Passive reflective markers, traceable by the 3D video motion capturing system, were fixed to the bony surface of the pelvis specimen. While applying a dynamic sinusoidal load the 3D-movement of the markers was recorded by the cameras and afterwards the 3D-deformation of the pelvis specimen was computed. The accuracy of the 3D-movement of the markers was verified with 3D-displacement curve with a step function using a manual driven 3D micro-motion-stage. Results The resulting accuracy of the measurement system depended on the number of cameras tracking a marker. The noise level for a marker seen by two cameras was during the stationary phase of the calibration procedure ± 0.036 mm, and ± 0.022 mm if tracked by 6 cameras. The detectable 3D-movement performed by the 3D-micro-motion-stage was smaller than the noise level of the 3D-video motion capturing system. Therefore the limiting factor of the setup was the noise level, which resulted in a measurement accuracy for the dynamic test setup of ± 0.036 mm. Conclusion This 3D test setup opens new possibilities in dynamic testing of wide range materials, like anatomical specimens, biomaterials, and its combinations. The resulting 3D-deformation dataset can be used for a better estimation of material
Optical-CT 3D Dosimetry Using Fresnel Lenses with Minimal Refractive-Index Matching Fluid.
Bache, Steven; Malcolm, Javian; Adamovics, John; Oldham, Mark
2016-01-01
Telecentric optical computed tomography (optical-CT) is a state-of-the-art method for visualizing and quantifying 3-dimensional dose distributions in radiochromic dosimeters. In this work a prototype telecentric system (DFOS-Duke Fresnel Optical-CT Scanner) is evaluated which incorporates two substantial design changes: the use of Fresnel lenses (reducing lens costs from $10-30K t0 $1-3K) and the use of a 'solid tank' (which reduces noise, and the volume of refractively matched fluid from 1 ltr to 10 cc). The efficacy of DFOS was evaluated by direct comparison against commissioned scanners in our lab. Measured dose distributions from all systems were compared against the predicted dose distributions from a commissioned treatment planning system (TPS). Three treatment plans were investigated including a simple four-field box treatment, a multiple small field delivery, and a complex IMRT treatment. Dosimeters were imaged within 2 h post irradiation, using consistent scanning techniques (360 projections acquired at 1 degree intervals, reconstruction at 2mm). DFOS efficacy was evaluated through inspection of dose line-profiles, and 2D and 3D dose and gamma maps. DFOS/TPS gamma pass rates with 3%/3mm dose difference/distance-to-agreement criteria ranged from 89.3% to 92.2%, compared to from 95.6% to 99.0% obtained with the commissioned system. The 3D gamma pass rate between the commissioned system and DFOS was 98.2%. The typical noise rates in DFOS reconstructions were up to 3%, compared to under 2% for the commissioned system. In conclusion, while the introduction of a solid tank proved advantageous with regards to cost and convenience, further work is required to improve the image quality and dose reconstruction accuracy of the new DFOS optical-CT system. PMID:27019460
Optical-CT 3D Dosimetry Using Fresnel Lenses with Minimal Refractive-Index Matching Fluid
Bache, Steven; Malcolm, Javian; Adamovics, John; Oldham, Mark
2016-01-01
Telecentric optical computed tomography (optical-CT) is a state-of-the-art method for visualizing and quantifying 3-dimensional dose distributions in radiochromic dosimeters. In this work a prototype telecentric system (DFOS—Duke Fresnel Optical-CT Scanner) is evaluated which incorporates two substantial design changes: the use of Fresnel lenses (reducing lens costs from $10-30K t0 $1-3K) and the use of a ‘solid tank’ (which reduces noise, and the volume of refractively matched fluid from 1ltr to 10cc). The efficacy of DFOS was evaluated by direct comparison against commissioned scanners in our lab. Measured dose distributions from all systems were compared against the predicted dose distributions from a commissioned treatment planning system (TPS). Three treatment plans were investigated including a simple four-field box treatment, a multiple small field delivery, and a complex IMRT treatment. Dosimeters were imaged within 2h post irradiation, using consistent scanning techniques (360 projections acquired at 1 degree intervals, reconstruction at 2mm). DFOS efficacy was evaluated through inspection of dose line-profiles, and 2D and 3D dose and gamma maps. DFOS/TPS gamma pass rates with 3%/3mm dose difference/distance-to-agreement criteria ranged from 89.3% to 92.2%, compared to from 95.6% to 99.0% obtained with the commissioned system. The 3D gamma pass rate between the commissioned system and DFOS was 98.2%. The typical noise rates in DFOS reconstructions were up to 3%, compared to under 2% for the commissioned system. In conclusion, while the introduction of a solid tank proved advantageous with regards to cost and convenience, further work is required to improve the image quality and dose reconstruction accuracy of the new DFOS optical-CT system. PMID:27019460
NASA Technical Reports Server (NTRS)
Fowlis, W. W.
1981-01-01
Systematic scaling or dimensional analysis reveals that certain scales of geophysical fluid flows (such as stellar, ocean, and planetary atmosphere circulations) can be accurately modeled in the laboratory using a procedure which differs from conventional engineering modeling. Rather than building a model to obtain numbers for a specific design problem, the relative effects of the significant forces are systematically varied in an attempt to deepen understanding of the effects of these forces. Topics covered include: (1) modeling a large-scale planetary atmospheric flow in a rotating cylindrical annulus; (2) achieving a radial dielectric body force; (3) spherical geophysical fluid dynamics experiments for Spacelab flights; (4) measuring flow and temperature; and (5) the possible effect of rotational or precessional disturbances on the flow in the rotating spherical containers.
3D Global Braginskii Simulations of Plasma Dynamics and Turbulence in LAPD
NASA Astrophysics Data System (ADS)
Fisher, Dustin; Rogers, Barrett
2013-10-01
3D global two-fluid simulations are presented in an ongoing effort to identify and understand the plasma dynamics in the Large Plasma Device (LAPD) at UCLA's Basic Science Facility. Modeling is done using a modified version of the Global Braginskii Solver (GBS) that models the plasma from source to edge region on a field-aligned grid using a finite difference method and 4th order Runge-Kutta time stepping. Progress has been made to account for the thermionic cathode emission of fast electrons at the source, the axial dependence of the plasma source, and biasing the front and side walls. Along with trying to understand the effect sheath's and neutrals have in setting the plasma potential, work is being done to model the biasable limiter recently used by colleagues at UCLA to better understand flow shear and particle transport in the LAPD. Comparisons of the zero bias case are presented along with analysis of the growth and dynamics of turbulent structures (such as drift waves) seen in the simulations. Supported through CICART under the auspices of the DOE's EPSCoR Grant No. DE-FG02-10ER46372.
Doyle, Andrew D; Carvajal, Nicole; Jin, Albert; Matsumoto, Kazue; Yamada, Kenneth M
2015-01-01
The physical properties of two-dimensional (2D) extracellular matrices (ECMs) modulate cell adhesion dynamics and motility, but little is known about the roles of local microenvironmental differences in three-dimensional (3D) ECMs. Here we generate 3D collagen gels of varying matrix microarchitectures to characterize their regulation of 3D adhesion dynamics and cell migration. ECMs containing bundled fibrils demonstrate enhanced local adhesion-scale stiffness and increased adhesion stability through balanced ECM/adhesion coupling, whereas highly pliable reticular matrices promote adhesion retraction. 3D adhesion dynamics are locally regulated by ECM rigidity together with integrin/ECM association and myosin II contractility. Unlike 2D migration, abrogating contractility stalls 3D migration regardless of ECM pore size. We find force is not required for clustering of activated integrins on 3D native collagen fibrils. We propose that efficient 3D migration requires local balancing of contractility with ECM stiffness to stabilize adhesions, which facilitates the detachment of activated integrins from ECM fibrils. PMID:26548801
Doyle, Andrew D.; Carvajal, Nicole; Jin, Albert; Matsumoto, Kazue; Yamada, Kenneth M.
2015-01-01
The physical properties of two-dimensional (2D) extracellular matrices (ECMs) modulate cell adhesion dynamics and motility, but little is known about the roles of local microenvironmental differences in three-dimensional (3D) ECMs. Here we generate 3D collagen gels of varying matrix microarchitectures to characterize their regulation of 3D adhesion dynamics and cell migration. ECMs containing bundled fibrils demonstrate enhanced local adhesion-scale stiffness and increased adhesion stability through balanced ECM/adhesion coupling, whereas highly pliable reticular matrices promote adhesion retraction. 3D adhesion dynamics are locally regulated by ECM rigidity together with integrin/ECM association and myosin II contractility. Unlike 2D migration, abrogating contractility stalls 3D migration regardless of ECM pore size. We find force is not required for clustering of activated integrins on 3D native collagen fibrils. We propose that efficient 3D migration requires local balancing of contractility with ECM stiffness to stabilize adhesions, which facilitates the detachment of activated integrins from ECM fibrils. PMID:26548801
NASA Astrophysics Data System (ADS)
Doyle, Andrew D.; Carvajal, Nicole; Jin, Albert; Matsumoto, Kazue; Yamada, Kenneth M.
2015-11-01
The physical properties of two-dimensional (2D) extracellular matrices (ECMs) modulate cell adhesion dynamics and motility, but little is known about the roles of local microenvironmental differences in three-dimensional (3D) ECMs. Here we generate 3D collagen gels of varying matrix microarchitectures to characterize their regulation of 3D adhesion dynamics and cell migration. ECMs containing bundled fibrils demonstrate enhanced local adhesion-scale stiffness and increased adhesion stability through balanced ECM/adhesion coupling, whereas highly pliable reticular matrices promote adhesion retraction. 3D adhesion dynamics are locally regulated by ECM rigidity together with integrin/ECM association and myosin II contractility. Unlike 2D migration, abrogating contractility stalls 3D migration regardless of ECM pore size. We find force is not required for clustering of activated integrins on 3D native collagen fibrils. We propose that efficient 3D migration requires local balancing of contractility with ECM stiffness to stabilize adhesions, which facilitates the detachment of activated integrins from ECM fibrils.
PFLOW: A 3-D Numerical Modeling Tool for Calculating Fluid-Pressure Diffusion from Coulomb Strain
NASA Astrophysics Data System (ADS)
Wolf, L. W.; Lee, M.; Meir, A.; Dyer, G.; Ma, K.; Chan, C.
2009-12-01
A new 3D time-dependent pore-pressure diffusion model PFLOW is developed to investigate the response of pore fluids to the crustal deformation generated by strong earthquakes in heterogeneous geologic media. Given crustal strain generated by changes in Coulomb stress, this MATLAB-based code uses Skempton's coefficient to calculate resulting changes fluid pressure. Pore-pressure diffusion can be tracked over time in a user-defined model space with user-prescribed Neumann or Dirchilet boundary conditions and with spatially variable values of permeability. PFLOW employs linear or quadratic finite elements for spatial discretization and first order or second order, explicit or implicit finite difference discretization in time. PFLOW is easily interfaced with output from deformation modeling programs such as Coulomb (Toda et al., 2007) or 3D-DEF (Gomberg and Ellis, 1994). The code is useful for investigating to first-order the evolution of pore pressure changes induced by changes in Coulomb stress and their possible relation to water-level changes in wells or changes in stream discharge. It can also be used for student research and classroom instruction. As an example application, we calculate the coseismic pore pressure changes and diffusion induced by volumetric strain associated with the 1999 Chi-Chi earthquake (Mw = 7.6) in Taiwan. The Chi-Chi earthquake provides an unique opportunity to investigate the spatial and time-dependent poroelastic response of near-field rocks and sediments because there exist extensive observational data of water-level changes and crustal deformation. The integrated model allows us to explore whether changes in Coulomb stress can adequately explain hydrologic anomalies observed in areas such as Taiwan’s western foothills and the Choshui River alluvial plain. To calculate coseismic strain, we use the carefully calibrated finite fault-rupture model of Ma et al. (2005) and the deformation modeling code Coulomb 3.1 (Toda et al., 2007
Dynamic visual image modeling for 3D synthetic scenes in agricultural engineering
NASA Astrophysics Data System (ADS)
Gao, Li; Yan, Juntao; Li, Xiaobo; Ji, Yatai; Li, Xin
The dynamic visual image modeling for 3D synthetic scenes by using dynamic multichannel binocular visual image based on the mobile self-organizing network. Technologies of 3D modeling synthetic scenes have been widely used in kinds of industries. The main purpose of this paper is to use multiple networks of dynamic visual monitors and sensors to observe an unattended area, to use the advantages of mobile network in rural areas for improving existing mobile network information service further and providing personalized information services. The goal of displaying is to provide perfect representation of synthetic scenes. Using low-power dynamic visual monitors and temperature/humidity sensor or GPS installed in the node equipment, monitoring data will be sent at scheduled time. Then through the mobile self-organizing network, 3D model is rebuilt by synthesizing the returned images. On this basis, we formalize a novel algorithm for multichannel binocular visual 3D images based on fast 3D modeling. Taking advantage of these low prices mobile, mobile self-organizing networks can get a large number of video from where is not suitable for human observation or unable to reach, and accurately synthetic 3D scene. This application will play a great role in promoting its application in agriculture.
NASA Astrophysics Data System (ADS)
Yang, Xiaofeng; Ghafourian, Pegah; Sharma, Puneet; Salman, Khalil; Martin, Diego; Fei, Baowei
2012-02-01
We have applied image analysis methods in the assessment of human kidney perfusion based on 3D dynamic contrast-enhanced (DCE) MRI data. This approach consists of 3D non-rigid image registration of the kidneys and fuzzy C-mean classification of kidney tissues. The proposed registration method reduced motion artifacts in the dynamic images and improved the analysis of kidney compartments (cortex, medulla, and cavities). The dynamic intensity curves show the successive transition of the contrast agent through kidney compartments. The proposed method for motion correction and kidney compartment classification may be used to improve the validity and usefulness of further model-based pharmacokinetic analysis of kidney function.
Parallel contact detection algorithm for transient solid dynamics simulations using PRONTO3D
Attaway, S.W.; Hendrickson, B.A.; Plimpton, S.J.
1996-09-01
An efficient, scalable, parallel algorithm for treating material surface contacts in solid mechanics finite element programs has been implemented in a modular way for MIMD parallel computers. The serial contact detection algorithm that was developed previously for the transient dynamics finite element code PRONTO3D has been extended for use in parallel computation by devising a dynamic (adaptive) processor load balancing scheme.
Qian Chen
2008-08-18
The generation, motion, and interaction of dislocations play key roles during the plastic deformation process of crystalline solids. 3D Dislocation Dynamics has been employed as a mesoscale simulation algorithm to investigate the collective and cooperative behavior of dislocations. Most current research on 3D Dislocation Dynamics is based on the solutions available in the framework of classical isotropic elasticity. However, due to some degree of elastic anisotropy in almost all crystalline solids, it is very necessary to extend 3D Dislocation Dynamics into anisotropic elasticity. In this study, first, the details of efficient and accurate incorporation of the fully anisotropic elasticity into 3D discrete Dislocation Dynamics by numerically evaluating the derivatives of Green's functions are described. Then the intrinsic properties of perfect dislocations, including their stability, their core properties and disassociation characteristics, in newly discovered rare earth-based intermetallics and in conventional intermetallics are investigated, within the framework of fully anisotropic elasticity supplemented with the atomistic information obtained from the ab initio calculations. Moreover, the evolution and interaction of dislocations in these intermetallics as well as the role of solute segregation are presented by utilizing fully anisotropic 3D dislocation dynamics. The results from this work clearly indicate the role and the importance of elastic anisotropy on the evolution of dislocation microstructures, the overall ductility and the hardening behavior in these systems.
NASA Astrophysics Data System (ADS)
Barnhart, Donald H.; Chan, Victor S. S.; Halliwell, Neil A.; Coupland, Jeremy M.
1999-10-01
This paper introduces a new approach to 3D displacement and velocity measurements that unifies the disciplines of holographic interferometry and holographic particle image velocimetry (HPIV). Equally applicable to fluid and solid mechanics, the overall system enables quantitative displacement measurements between two holographically recorded events from either particle or surface scattering sites, working with both pulsed and continuous-wave laser systems. The resulting measurements exhibit an accuracy corresponding to interferometric system, but with a dynamic range found with PIV systems. Most importantly, this paper introduces the novel use of an optical fiber to specify the measurement points, remove optical aberrations of windows, and eliminate directional ambiguity. An optical fiber is used to probe the recorded holographic image space at each 3D measurement point in order to extract the 3D displacement vectors. This fiber system also employs a novel optical image shifting method to eliminate the problem of directional ambiguity. In addition, the reported system uses 3D complex optical correlation rather than 2D real digital correlation. It is therefore a simple matter to directly obtain 3D displacement and velocity measurements at precisely known 3D locations in the object space. By correlating both the amplitude and phase information in the holographic image, this system can measure spatial distributions of displacements even when the presence of severe aberrations preclude the detection of sharp images.
NASA Astrophysics Data System (ADS)
Jackson, Derek; Cooper, Andrew; Green, Andrew; Beyers, Meiring; Wiles, Errol; Benallack, Keegan
2016-04-01
Un-vegetated dune fields provide excellent opportunities to examine airflow dynamics over various types and scales of dune landforms. The three dimensional surface over which lower boundary layers travel, help adjust surface airflow and consequently the aeolian response of the dunes themselves. The use of computational fluid dynamic (CFD) modelling in recent studies now enables investigation of the 3D behaviour of airflow over complex terrain, providing new insights into heterogeneous surface flow and aeolian response of dune surfaces on a large (dunefield) scale. Using a largely un-vegetated coastal dune field site at Mpekweni, Eastern Cape, South Africa, a detailed (0.1m gridded) terrestrial laser scanning survey was conducted to create a high resolution topographical surface. Using local wind flow measurements and local met station records as input, CFD modelling was performed for a number of scenarios involving variable direction and magnitude to examine surface flow patterns across multiple dune forms. Near surface acceleration, expansion and separation of airflow inducing convergence and divergence (steering) of flow velocity streamlines are investigated. Flow acceleration over dune crests/brink lines is a key parameter in driving dune migration and slip face dynamics. Dune aspect ratio (height to length) is also important in determining the degree of crestal flow acceleration, with an increase in flow associated with increasing aspect ratios. Variations in dune height appear to be the most important parameter in driving general flow acceleration. The results from the study provide new insights into dune migration behaviour at this site as well as surface flow behaviour across multiple dune configurations and length scales within un-vegetated dune fields.
3D Dynamics of Freshwater Lenses in the Near-Surface Layer of the Tropical Ocean
NASA Astrophysics Data System (ADS)
Soloviev, Alexander; Dean, Cayla
2015-04-01
Convective rains in the Intertropical Convergence Zone (ITCZ) produce lenses of freshened water on the ocean surface. These lenses are localized in space and typically involve both salinity and temperature anomalies. Due to significant density anomalies, strong pressure gradients develop, which result in lateral spreading of freshwater lenses in a form resembling gravity currents. Gravity currents inherently involve three-dimensional dynamics. As a type of organized structure, gravity currents in the upper layer of the ocean may also interact with, and be shaped by, the ambient oceanic environment and atmospheric conditions. Among the important factors are the background stratification, wind stress, wind/wave mixing and spatially coherent organized motions in the near-surface layer of the ocean. Under certain conditions, a resonant interaction between a propagating freshwater lens and internal waves in the underlying pycnocline (e.g., barrier layer) may develop, whereas interaction with wind stress may produce an asymmetry in the freshwater lens and associated mixing. These two types of interactions working in concert may explain the series of sharp frontal interfaces, which have been observed in association with freshwater lenses during TOGA COARE. In this work, we have conducted a series of numerical experiments using computational fluid dynamics tools. These numerical simulations were designed to elucidate the relationship between vertical mixing and horizontal advection of salinity under various environmental conditions and potential impact on the Aquarius and SMOS satellite image formation. Available near-surface data from field experiments served as a guidance for numerical simulations. The results of this study indicate that 3D dynamics of freshwater lenses are essential within a certain range of wind/wave conditions and the freshwater influx in the surface layer of the ocean.
Costa, Pedro F; Hutmacher, Dietmar W; Theodoropoulos, Christina; Gomes, Manuela E; Reis, Rui L; Vaquette, Cédryck
2015-04-22
The ability to test large arrays of cell and biomaterial combinations in 3D environments is still rather limited in the context of tissue engineering and regenerative medicine. This limitation can be generally addressed by employing highly automated and reproducible methodologies. This study reports on the development of a highly versatile and upscalable method based on additive manufacturing for the fabrication of arrays of scaffolds, which are enclosed into individualized perfusion chambers. Devices containing eight scaffolds and their corresponding bioreactor chambers are simultaneously fabricated utilizing a dual extrusion additive manufacturing system. To demonstrate the versatility of the concept, the scaffolds, while enclosed into the device, are subsequently surface-coated with a biomimetic calcium phosphate layer by perfusion with simulated body fluid solution. 96 scaffolds are simultaneously seeded and cultured with human osteoblasts under highly controlled bidirectional perfusion dynamic conditions over 4 weeks. Both coated and noncoated resulting scaffolds show homogeneous cell distribution and high cell viability throughout the 4 weeks culture period and CaP-coated scaffolds result in a significantly increased cell number. The methodology developed in this work exemplifies the applicability of additive manufacturing as a tool for further automation of studies in the field of tissue engineering and regenerative medicine. PMID:25721231
Magnetic levitation fluid dynamics
NASA Astrophysics Data System (ADS)
Bojarevics, V.; Pericleous, K.
2001-06-01
This work is concerned with the accurate computation of flow in a rapidly deforming liquid metal droplet, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope, which is initially spherical. The relative positional change between the liquid metal surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. A spectral technique is used to solve this problem, which is assumed axisymmetric. The computed results are compared against a physical experiment and "ideal sphere" analytic solutions. A comparison between the "magnetic pressure" approximation and the full electromagnetic force solutions, shows fundamental differences; the full electromagnetic force solution is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation. Figs 9, Refs 13.
Fluid dynamics of heart development.
Santhanakrishnan, Arvind; Miller, Laura A
2011-09-01
The morphology, muscle mechanics, fluid dynamics, conduction properties, and molecular biology of the developing embryonic heart have received much attention in recent years due to the importance of both fluid and elastic forces in shaping the heart as well as the striking relationship between the heart's evolution and development. Although few studies have directly addressed the connection between fluid dynamics and heart development, a number of studies suggest that fluids may play a key role in morphogenic signaling. For example, fluid shear stress may trigger biochemical cascades within the endothelial cells of the developing heart that regulate chamber and valve morphogenesis. Myocardial activity generates forces on the intracardiac blood, creating pressure gradients across the cardiac wall. These pressures may also serve as epigenetic signals. In this article, the fluid dynamics of the early stages of heart development is reviewed. The relevant work in cardiac morphology, muscle mechanics, regulatory networks, and electrophysiology is also reviewed in the context of intracardial fluid dynamics. PMID:21327946
First application of the 3D-MHB on dynamic compressive behavior of UHPC
NASA Astrophysics Data System (ADS)
Cadoni, Ezio; Dotta, Matteo; Forni, Daniele; Riganti, Gianmario; Albertini, Carlo
2015-09-01
In order to study the dynamic behaviour of material in confined conditions a new machine was conceived and called 3D-Modified Hopkinson Bar (3D-MHB). It is a Modified Hopkinson Bar apparatus designed to apply dynamic loading in materials having a tri-axial stress state. It consists of a pulse generator system (with pre-tensioned bar and brittle joint), 1 input bar, and 5 output bars. The first results obtained on Ultra High Performance Concrete in compression with three different mono-axial compression states are presented. The results show how the pre-stress states minimize the boundary condition and a more uniform response is obtained.
Parallel load balancing strategy for Volume-of-Fluid methods on 3-D unstructured meshes
NASA Astrophysics Data System (ADS)
Jofre, Lluís; Borrell, Ricard; Lehmkuhl, Oriol; Oliva, Assensi
2015-02-01
Volume-of-Fluid (VOF) is one of the methods of choice to reproduce the interface motion in the simulation of multi-fluid flows. One of its main strengths is its accuracy in capturing sharp interface geometries, although requiring for it a number of geometric calculations. Under these circumstances, achieving parallel performance on current supercomputers is a must. The main obstacle for the parallelization is that the computing costs are concentrated only in the discrete elements that lie on the interface between fluids. Consequently, if the interface is not homogeneously distributed throughout the domain, standard domain decomposition (DD) strategies lead to imbalanced workload distributions. In this paper, we present a new parallelization strategy for general unstructured VOF solvers, based on a dynamic load balancing process complementary to the underlying DD. Its parallel efficiency has been analyzed and compared to the DD one using up to 1024 CPU-cores on an Intel SandyBridge based supercomputer. The results obtained on the solution of several artificially generated test cases show a speedup of up to ∼12× with respect to the standard DD, depending on the interface size, the initial distribution and the number of parallel processes engaged. Moreover, the new parallelization strategy presented is of general purpose, therefore, it could be used to parallelize any VOF solver without requiring changes on the coupled flow solver. Finally, note that although designed for the VOF method, our approach could be easily adapted to other interface-capturing methods, such as the Level-Set, which may present similar workload imbalances.
3D dynamic simulation of crack propagation in extracorporeal shock wave lithotripsy
NASA Astrophysics Data System (ADS)
Wijerathne, M. L. L.; Hori, Muneo; Sakaguchi, Hide; Oguni, Kenji
2010-06-01
Some experimental observations of Shock Wave Lithotripsy(SWL), which include 3D dynamic crack propagation, are simulated with the aim of reproducing fragmentation of kidney stones with SWL. Extracorporeal shock wave lithotripsy (ESWL) is the fragmentation of kidney stones by focusing an ultrasonic pressure pulse onto the stones. 3D models with fine discretization are used to accurately capture the high amplitude shear shock waves. For solving the resulting large scale dynamic crack propagation problem, PDS-FEM is used; it provides numerically efficient failure treatments. With a distributed memory parallel code of PDS-FEM, experimentally observed 3D photoelastic images of transient stress waves and crack patterns in cylindrical samples are successfully reproduced. The numerical crack patterns are in good agreement with the experimental ones, quantitatively. The results shows that the high amplitude shear waves induced in solid, by the lithotriptor generated shock wave, play a dominant role in stone fragmentation.
Stellar Astrophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Thompson, Michael J.; Christensen-Dalsgaard, Jørgen
2003-05-01
Preface; 1. A selective overview Jørgen Christensen-Dalsgaard and Michael J. Thompson; Part I. Stellar Convection and Oscillations: 2. On the diversity of stellar pulsations Wojciech A. Dziembowski; 3. Acoustic radiation and mode excitation by turbulent convection Günter Houdek; 4. Understanding roAp stars Margarida S. Cunha; 5. Waves in the magnetised solar atmosphere Colin S. Rosenthal; Part II. Stellar Rotation and Magnetic Fields: 6. Stellar rotation: a historical survey Leon Mestel; 7. The oscillations of rapidly rotating stars Michel Rieutord; 8. Solar tachocline dynamics: eddy viscosity, anti-friction, or something in between? Michael E. McIntyre; 9. Dynamics of the solar tachocline Pascale Garaud; 10. Dynamo processes: the interaction of turbulence and magnetic fields Michael Proctor; 11. Dynamos in planets Chris Jones; Part III. Physics and Structure of Stellar Interiors: 12. Solar constraints on the equation of state Werner Däppen; 13. 3He transport and the solar neutrino problem Chris Jordinson; 14. Mixing in stellar radiation zones Jean-Paul Zahn; 15. Element settling and rotation-induced mixing in slowly rotating stars Sylvie Vauclair; Part IV. Helio- and Asteroseismology: 16. Solar structure and the neutrino problem Hiromoto Shibahashi; 17. Helioseismic data analysis Jesper Schou; 18. Seismology of solar rotation Takashi Sekii; 19. Telechronohelioseismology Alexander Kosovichev; Part V. Large-Scale Numerical Experiments: 20. Bridges between helioseismology and models of convection zone dynamics Juri Toomre; 21. Numerical simulations of the solar convection zone Julian R. Elliott; 22. Modelling solar and stellar magnetoconvection Nigel Weiss; 23. Nonlinear magnetoconvection in the presence of a strong oblique field Keith Julien, Edgar Knobloch and Steven M. Tobias; 24. Simulations of astrophysical fluids Marcus Brüggen; Part VI. Dynamics: 25. A magic electromagnetic field Donald Lynden-Bell; 26. Continuum equations for stellar dynamics Edward A
Stellar Astrophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Thompson, Michael J.; Christensen-Dalsgaard, Jørgen
2008-02-01
Preface; 1. A selective overview Jørgen Christensen-Dalsgaard and Michael J. Thompson; Part I. Stellar Convection and Oscillations: 2. On the diversity of stellar pulsations Wojciech A. Dziembowski; 3. Acoustic radiation and mode excitation by turbulent convection Günter Houdek; 4. Understanding roAp stars Margarida S. Cunha; 5. Waves in the magnetised solar atmosphere Colin S. Rosenthal; Part II. Stellar Rotation and Magnetic Fields: 6. Stellar rotation: a historical survey Leon Mestel; 7. The oscillations of rapidly rotating stars Michel Rieutord; 8. Solar tachocline dynamics: eddy viscosity, anti-friction, or something in between? Michael E. McIntyre; 9. Dynamics of the solar tachocline Pascale Garaud; 10. Dynamo processes: the interaction of turbulence and magnetic fields Michael Proctor; 11. Dynamos in planets Chris Jones; Part III. Physics and Structure of Stellar Interiors: 12. Solar constraints on the equation of state Werner Däppen; 13. 3He transport and the solar neutrino problem Chris Jordinson; 14. Mixing in stellar radiation zones Jean-Paul Zahn; 15. Element settling and rotation-induced mixing in slowly rotating stars Sylvie Vauclair; Part IV. Helio- and Asteroseismology: 16. Solar structure and the neutrino problem Hiromoto Shibahashi; 17. Helioseismic data analysis Jesper Schou; 18. Seismology of solar rotation Takashi Sekii; 19. Telechronohelioseismology Alexander Kosovichev; Part V. Large-Scale Numerical Experiments: 20. Bridges between helioseismology and models of convection zone dynamics Juri Toomre; 21. Numerical simulations of the solar convection zone Julian R. Elliott; 22. Modelling solar and stellar magnetoconvection Nigel Weiss; 23. Nonlinear magnetoconvection in the presence of a strong oblique field Keith Julien, Edgar Knobloch and Steven M. Tobias; 24. Simulations of astrophysical fluids Marcus Brüggen; Part VI. Dynamics: 25. A magic electromagnetic field Donald Lynden-Bell; 26. Continuum equations for stellar dynamics Edward A
Trans3D: a free tool for dynamical visualization of EEG activity transmission in the brain.
Blinowski, Grzegorz; Kamiński, Maciej; Wawer, Dariusz
2014-08-01
The problem of functional connectivity in the brain is in the focus of attention nowadays, since it is crucial for understanding information processing in the brain. A large repertoire of measures of connectivity have been devised, some of them being capable of estimating time-varying directed connectivity. Hence, there is a need for a dedicated software tool for visualizing the propagation of electrical activity in the brain. To this aim, the Trans3D application was developed. It is an open access tool based on widely available libraries and supporting both Windows XP/Vista/7(™), Linux and Mac environments. Trans3D can create animations of activity propagation between electrodes/sensors, which can be placed by the user on the scalp/cortex of a 3D model of the head. Various interactive graphic functions for manipulating and visualizing components of the 3D model and input data are available. An application of the Trans3D tool has helped to elucidate the dynamics of the phenomena of information processing in motor and cognitive tasks, which otherwise would have been very difficult to observe. Trans3D is available at: http://www.eeg.pl/. PMID:24967953
NASA Technical Reports Server (NTRS)
Bancroft, Gregory N.; Sikavitsas, Vassilios I.; van den Dolder, Juliette; Sheffield, Tiffany L.; Ambrose, Catherine G.; Jansen, John A.; Mikos, Antonios G.; McIntire, L. V. (Principal Investigator)
2002-01-01
Bone is a complex highly structured mechanically active 3D tissue composed of cellular and matrix elements. The true biological environment of a bone cell is thus derived from a dynamic interaction between responsively active cells experiencing mechanical forces and a continuously changing 3D matrix architecture. To investigate this phenomenon in vitro, marrow stromal osteoblasts were cultured on 3D scaffolds under flow perfusion with different rates of flow for an extended period to permit osteoblast differentiation and significant matrix production and mineralization. With all flow conditions, mineralized matrix production was dramatically increased over statically cultured constructs with the total calcium content of the cultured scaffolds increasing with increasing flow rate. Flow perfusion induced de novo tissue modeling with the formation of pore-like structures in the scaffolds and enhanced the distribution of cells and matrix throughout the scaffolds. These results represent reporting of the long-term effects of fluid flow on primary differentiating osteoblasts and indicate that fluid flow has far-reaching effects on osteoblast differentiation and phenotypic expression in vitro. Flow perfusion culture permits the generation and study of a 3D, actively modeled, mineralized matrix and can therefore be a valuable tool for both bone biology and tissue engineering.
Computational 3D fluid-structure interaction for the aortic valve
NASA Astrophysics Data System (ADS)
Luo, Haoxiang; Chen, Ye; Sun, Wei
2015-11-01
Three-dimensional fluid-structure interaction (FSI) involving large deformations of flexible bodies is common in biological systems. A typical example is the heart valves. Accurate and efficient numerical approaches for modeling such systems are still lacking. In this work, we report a successful case of combining an immersed-boundary flow solver with a nonlinear finite-element solid-dynamics solver, both in-house programs, specifically for three-dimensional simulations. Based on the Cartesian grid, the viscous incompressible flow solver can handle boundaries of large displacements with simple mesh generation. The solid-dynamics solver has separate subroutines for analyzing general three-dimensional bodies and thin-walled structures composed of frames, membranes, and plates. Both geometric nonlinearity associated with large displacements and material nonlinearity associated with large strains are incorporated in the solver. The FSI is achieved through a strong coupling and partitioned approach. We have performed several benchmarking cases to validate the FSI solver. Application to the native aortic valve will be demonstrated. Supported by the NSF grant (CBET-1066962).
ERIC Educational Resources Information Center
Berney, Sandra; Bétrancourt, Mireille; Molinari, Gaëlle; Hoyek, Nady
2015-01-01
The emergence of dynamic visualizations of three-dimensional (3D) models in anatomy curricula may be an adequate solution for spatial difficulties encountered with traditional static learning, as they provide direct visualization of change throughout the viewpoints. However, little research has explored the interplay between learning material…
Effect of Ductile Agents on the Dynamic Behavior of SiC3D Network Composites
NASA Astrophysics Data System (ADS)
Zhu, Jingbo; Wang, Yangwei; Wang, Fuchi; Fan, Qunbo
2016-07-01
Co-continuous SiC ceramic composites using pure aluminum, epoxy, and polyurethane (PU) as ductile agents were developed. The dynamic mechanical behavior and failure mechanisms were investigated experimentally using the split Hopkinson pressure bar (SHPB) method and computationally by finite element (FE) simulations. The results show that the SiC3D/Al composite has the best overall performance in comparison with SiC3D/epoxy and SiC3D/PU composites. FE simulations are generally consistent with experimental data. These simulations provide valuable help in predicting mechanical strength and in interpreting the experimental results and failure mechanisms. They may be combined with micrographs for fracture characterizations of the composites. We found that interactions between the SiC phase and ductile agents under dynamic compression in the SHPB method are complex, and that interfacial condition is an important parameter that determines the mechanical response of SiC3D composites with a characteristic interlocking structure during dynamic compression. However, the effect of the mechanical properties of ductile agents on dynamic behavior of the composites is a second consideration in the production of the composites.
3D Flow Visualization Using Texture Advection
NASA Technical Reports Server (NTRS)
Kao, David; Zhang, Bing; Kim, Kwansik; Pang, Alex; Moran, Pat (Technical Monitor)
2001-01-01
Texture advection is an effective tool for animating and investigating 2D flows. In this paper, we discuss how this technique can be extended to 3D flows. In particular, we examine the use of 3D and 4D textures on 3D synthetic and computational fluid dynamics flow fields.
Sketch on dynamic gesture tracking and analysis exploiting vision-based 3D interface
NASA Astrophysics Data System (ADS)
Woo, Woontack; Kim, Namgyu; Wong, Karen; Tadenuma, Makoto
2000-12-01
In this paper, we propose a vision-based 3D interface exploiting invisible 3D boxes, arranged in the personal space (i.e. reachable space by the body without traveling), which allows robust yet simple dynamic gesture tracking and analysis, without exploiting complicated sensor-based motion tracking systems. Vision-based gesture tracking and analysis is still a challenging problem, even though we have witnessed rapid advances in computer vision over the last few decades. The proposed framework consists of three main parts, i.e. (1) object segmentation without bluescreen and 3D box initialization with depth information, (2) movement tracking by observing how the body passes through the 3D boxes in the personal space and (3) movement feature extraction based on Laban's Effort theory and movement analysis by mapping features to meaningful symbols using time-delay neural networks. Obviously, exploiting depth information using multiview images improves the performance of gesture analysis by reducing the errors introduced by simple 2D interfaces In addition, the proposed box-based 3D interface lessens the difficulties in both tracking movement in 3D space and in extracting low-level features of the movement. Furthermore, the time-delay neural networks lessens the difficulties in movement analysis by training. Due to its simplicity and robustness, the framework will provide interactive systems, such as ATR I-cubed Tangible Music System or ATR Interactive Dance system, with improved quality of the 3D interface. The proposed simple framework also can be extended to other applications requiring dynamic gesture tracking and analysis on the fly.
Computational fluid dynamic applications
Chang, S.-L.; Lottes, S. A.; Zhou, C. Q.
2000-04-03
The rapid advancement of computational capability including speed and memory size has prompted the wide use of computational fluid dynamics (CFD) codes to simulate complex flow systems. CFD simulations are used to study the operating problems encountered in system, to evaluate the impacts of operation/design parameters on the performance of a system, and to investigate novel design concepts. CFD codes are generally developed based on the conservation laws of mass, momentum, and energy that govern the characteristics of a flow. The governing equations are simplified and discretized for a selected computational grid system. Numerical methods are selected to simplify and calculate approximate flow properties. For turbulent, reacting, and multiphase flow systems the complex processes relating to these aspects of the flow, i.e., turbulent diffusion, combustion kinetics, interfacial drag and heat and mass transfer, etc., are described in mathematical models, based on a combination of fundamental physics and empirical data, that are incorporated into the code. CFD simulation has been applied to a large variety of practical and industrial scale flow systems.
Tunesi, Marta; Fusco, Federica; Fiordaliso, Fabio; Corbelli, Alessandro; Biella, Gloria; Raimondi, Manuela T
2016-01-01
Frontotemporal lobar degeneration (FTLD) is a severe neurodegenerative disorder that is diagnosed with increasing frequency in clinical setting. Currently, no therapy is available and in addition the molecular basis of the disease are far from being elucidated. Consequently, it is of pivotal importance to develop reliable and cost-effective in vitro models for basic research purposes and drug screening. To this respect, recent results in the field of Alzheimer's disease have suggested that a tridimensional (3D) environment is an added value to better model key pathologic features of the disease. Here, we have tried to add complexity to the 3D cell culturing concept by using a microfluidic bioreactor, where cells are cultured under a continuous flow of medium, thus mimicking the interstitial fluid movement that actually perfuses the body tissues, including the brain. We have implemented this model using a neuronal-like cell line (SH-SY5Y), a widely exploited cell model for neurodegenerative disorders that shows some basic features relevant for FTLD modeling, such as the release of the FTLD-related protein progranulin (PRGN) in specific vesicles (exosomes). We have efficiently seeded the cells on 3D scaffolds, optimized a disease-relevant oxidative stress experiment (by targeting mitochondrial function that is one of the possible FTLD-involved pathological mechanisms) and evaluated cell metabolic activity in dynamic culture in comparison to static conditions, finding that SH-SY5Y cells cultured in 3D scaffold are susceptible to the oxidative damage triggered by a mitochondrial-targeting toxin (6-OHDA) and that the same cells cultured in dynamic conditions kept their basic capacity to secrete PRGN in exosomes once recovered from the bioreactor and plated in standard 2D conditions. We think that a further improvement of our microfluidic system may help in providing a full device where assessing basic FTLD-related features (including PRGN dynamic secretion) that may be
Dynamic scattering theory for dark-field electron holography of 3D strain fields.
Lubk, Axel; Javon, Elsa; Cherkashin, Nikolay; Reboh, Shay; Gatel, Christophe; Hÿtch, Martin
2014-01-01
Dark-field electron holography maps strain in crystal lattices into reconstructed phases over large fields of view. Here we investigate the details of the lattice strain-reconstructed phase relationship by applying dynamic scattering theory both analytically and numerically. We develop efficient analytic linear projection rules for 3D strain fields, facilitating a straight-forward calculation of reconstructed phases from 3D strained materials. They are used in the following to quantify the influence of various experimental parameters like strain magnitude, specimen thickness, excitation error and surface relaxation. PMID:24012934
Atmospheric and Oceanic Fluid Dynamics
NASA Astrophysics Data System (ADS)
Vallis, Geoffrey K.
2006-11-01
Fluid dynamics is fundamental to our understanding of the atmosphere and oceans. Although many of the same principles of fluid dynamics apply to both the atmosphere and oceans, textbooks tend to concentrate on the atmosphere, the ocean, or the theory of geophysical fluid dynamics (GFD). This textbook provides a comprehensive unified treatment of atmospheric and oceanic fluid dynamics. The book introduces the fundamentals of geophysical fluid dynamics, including rotation and stratification, vorticity and potential vorticity, and scaling and approximations. It discusses baroclinic and barotropic instabilities, wave-mean flow interactions and turbulence, and the general circulation of the atmosphere and ocean. Student problems and exercises are included at the end of each chapter. Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation will be an invaluable graduate textbook on advanced courses in GFD, meteorology, atmospheric science and oceanography, and an excellent review volume for researchers. Additional resources are available at www.cambridge.org/9780521849692. Includes end of chapter review questions to aid understanding Unified and comprehensive treatment of both atmospheric and oceanic fluid dynamics Covers many modern topics and provides up to date knowledge
Oh, Heeseok; Lee, Sanghoon; Bovik, Alan Conrad
2016-02-01
The human visual system perceives 3D depth following sensing via its binocular optical system, a series of massively parallel processing units, and a feedback system that controls the mechanical dynamics of eye movements and the crystalline lens. The process of accommodation (focusing of the crystalline lens) and binocular vergence is controlled simultaneously and symbiotically via cross-coupled communication between the two critical depth computation modalities. The output responses of these two subsystems, which are induced by oculomotor control, are used in the computation of a clear and stable cyclopean 3D image from the input stimuli. These subsystems operate in smooth synchronicity when one is viewing the natural world; however, conflicting responses can occur when viewing stereoscopic 3D (S3D) content on fixed displays, causing physiological discomfort. If such occurrences could be predicted, then they might also be avoided (by modifying the acquisition process) or ameliorated (by changing the relative scene depth). Toward this end, we have developed a dynamic accommodation and vergence interaction (DAVI) model that successfully predicts visual discomfort on S3D images. The DAVI model is based on the phasic and reflex responses of the fast fusional vergence mechanism. Quantitative models of accommodation and vergence mismatches are used to conduct visual discomfort prediction. Other 3D perceptual elements are included in the proposed method, including sharpness limits imposed by the depth of focus and fusion limits implied by Panum's fusional area. The DAVI predictor is created by training a support vector machine on features derived from the proposed model and on recorded subjective assessment results. The experimental results are shown to produce accurate predictions of experienced visual discomfort. PMID:26672036
Analysis of the formation and evolution of vortex rings in non Newtonian fluids using 3D PTV
NASA Astrophysics Data System (ADS)
Bajpayee, Abhishek; Techet, Alexandra
2013-11-01
Formation and evolution of vortex rings have been studied for a long time but mostly only in Newtonian fluids. However, many fluids in nature and in the industry such as blood, crude oil, etc., exhibit non Newtonian characteristics. Palacios-Morales and Zenit recently studied the formation of vortex rings in shear thinning liquids for the first time using 2D PIV and compared experimental findings with theoretical predictions. The authors recently demonstrated the applicability of Light Field (LF) imaging to conduct 3D Particle Tracking Velocimetry (PTV) to study densely seeded flow fields and their evolution over time using synthetic data. LF based 3D PTV is now used to quantitatively study vortex rings created in Glycerin based on multiple parameters and the results are compared with previous findings. ONR (Grant #N00014-12-1-0787, Dr. Steven Russell), Naval Engineering Education Center.
Computational Fluid Dynamics Technology for Hypersonic Applications
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.
2003-01-01
Several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented from code validation and code benchmarking efforts to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified. Highlights of diverse efforts to address these challenges are then discussed. One such effort to re-engineer and synthesize the existing analysis capability in LAURA, VULCAN, and FUN3D will provide context for these discussions. The critical (and evolving) role of agile software engineering practice in the capability enhancement process is also noted.
SAGE 2D and 3D Simulations of the Explosive Venting of Supercritical Fluids Through Porous Media
NASA Astrophysics Data System (ADS)
Weaver, R.; Gisler, G.; Svensen, H.; Mazzini, A.
2008-12-01
Magmatic intrusive events in large igneous provinces heat sedimentary country rock leading to the eventual release of volatiles. This has been proposed as a contributor to climate change and other environmental impacts. By means of numerical simulations, we examine ways in which these volatiles can be released explosively from depth. Gases and fluids cooked out of country rock by metamorphic heating may be confined for a time by impermeable clays or other barriers, developing high pressures and supercritical fluids. If confinement is suddenly breached (by an earthquake for example) in such a way that the fluid has access to porous sediments, a violent eruption of a non-magmatic mixture of fluid and sediment may result. Surface manifestations of these events could be hydrothermal vent complexes, kimberlite pipes, pockmarks, or mud volcanoes. These are widespread on Earth, especially in large igneous provinces, as in the Karoo Basin of South Africa, the North Sea off the Norwegian margin, and the Siberian Traps. We have performed 2D and 3D simulations with the Sage hydrocode (from Los Alamos and Science Applications International) of supercritical venting in a variety of geometries and configurations. The simulations show several different patterns of propagation and fracturing in porous or otherwise weakened overburden, dependent on depth, source conditions (fluid availability, temperature, and pressure), and manner of confinement breach. Results will be given for a variety of 2D and 3D simulations of these events exploring the release of volatiles into the atmosphere.
Santee, G.E. Jr.; Chang, F.H.; Mortensen, G.A.; Brockett, G.F.; Gross, M.B.; Belytschko, T.B.
1982-11-01
This report, the third in a series of reports for RP-1065, describes the final step in the stepwise approach for developing the three-dimensional, nonlinear, fluid-structure interaction methodology to assess the hydroloads on a large PWR during the subcooled portions of a hypothetical LOCA. The final step in the methodology implements enhancements and special modifications to the STEALTH 3D computer program and the WHAMSE 3D computer program. After describing the enhancements, the individual and the coupled computer programs are assessed by comparing calculational results with either analytical solutions or with experimental data. The coupled 3D STEALTH/WHAMSE computer program is then applied to the simulation of HDR Test V31.1 to further assess the program and to investigate the role that fluid-structure interaction plays in the hydrodynamic loading of reactor internals during subcooled blowdown.
Predictions of bubbly flows in vertical pipes using two-fluid models in CFDS-FLOW3D code
Banas, A.O.; Carver, M.B.; Unrau, D.
1995-09-01
This paper reports the results of a preliminary study exploring the performance of two sets of two-fluid closure relationships applied to the simulation of turbulent air-water bubbly upflows through vertical pipes. Predictions obtained with the default CFDS-FLOW3D model for dispersed flows were compared with the predictions of a new model (based on the work of Lee), and with the experimental data of Liu. The new model, implemented in the CFDS-FLOW3D code, included additional source terms in the {open_quotes}standard{close_quotes} {kappa}-{epsilon} transport equations for the liquid phase, as well as modified model coefficients and wall functions. All simulations were carried out in a 2-D axisymmetric format, collapsing the general multifluid framework of CFDS-FLOW3D to the two-fluid (air-water) case. The newly implemented model consistently improved predictions of radial-velocity profiles of both phases, but failed to accurately reproduce the experimental phase-distribution data. This shortcoming was traced to the neglect of anisotropic effects in the modelling of liquid-phase turbulence. In this sense, the present investigation should be considered as the first step toward the ultimate goal of developing a theoretically sound and universal CFD-type two-fluid model for bubbly flows in channels.
Progress in the Peeling-Ballooning Model of ELMs: Numerical Studies of 3D Nonlinear ELM Dynamics
Snyder, P B; Wilson, H R; Xu, X Q
2004-12-13
Nonlinear simulations with the 3D electromagnetic two-fluid BOUT code are employed to study the dynamics of edge localized modes (ELMs) driven by intermediate wavelength peeling-ballooning modes. It is found that the early behavior of the modes is similar to expectations from linear, ideal peeling-ballooning mode theory, with the modes growing linearly at a fraction of the Alfven frequency. In the non-linear phase, the modes grow explosively, forming a number of extended filaments which propagate rapidly from the outer closed flux region into the open flux region toward the outer wall. Similarities to non-linear linear ballooning theory, as well as additional complexities are observed. Comparison to observations reveals a number of similarities. Implications of the simulations and proposals for the dynamics of the full ELM crash are discussed.
PROGRESS IN THE PEELING-BALLOONING MODEL OF ELMS: NUMERICAL STUDIES OF 3D NONLINEAR ELM DYNAMICS
SNYDER,P.B; WILSON,H.R; XU,X.Q
2004-11-01
Nonlinear simulations with the 3D electromagnetic two-fluid BOUT code are employed to study the dynamics of edge localized modes (ELMs) driven by intermediate wavelength peeling-ballooning modes. It is found that the early behavior of the modes is similar to expectations from linear, ideal peeling-ballooning mode theory, with the modes growing linearly at a fraction of the Alfven frequency. In the nonlinear phase, the modes grow explosively, forming a number of extended filaments which propagate rapidly from the outer closed flux region into the open flux region toward the outboard wall. Similarities to non-linear ballooning theory, as well as additional complexities are observed. Comparison to observations reveals a number of similarities. Implications of the simulations and proposals for the dynamics of the full ELM crash are discussed.
Fluid Dynamics with Free Surfaces
Energy Science and Technology Software Center (ESTSC)
1992-02-01
RIPPLE is a two-dimensional, transient, free surface incompressible fluid dynamics program. It allows multiple free surfaces with surface tension and wall adhesion forces and has a partial cell treatment which allows curved boundaries and interior obstacles.
The Variety of Fluid Dynamics.
ERIC Educational Resources Information Center
Barnes, Francis; And Others
1980-01-01
Discusses three research topics which are concerned with eminently practical problems and deal at the same time with fundamental fluid dynamical problems. These research topics come from the general areas of chemical and biological engineering, geophysics, and pure mathematics. (HM)
3D PLLA/ibuprofen composite scaffolds obtained by a supercritical fluids assisted process.
Cardea, S; Baldino, L; Scognamiglio, M; Reverchon, E
2014-04-01
The emerging next generation of engineered tissues is based on the development of loaded scaffolds containing bioactive molecules in order to control the cellular function or to interact on the surrounding tissues. Indeed, implantation of engineered biomaterials might cause local inflammation because of the host's immune response; thereby, the use of anti-inflammatory agents, whether steroidal or nonsteroidal is required. One of the most important stages of tissue engineering is the design and the generation of a porous 3D structure, with high porosity, high interconnectivity and homogenous morphology. Various techniques have been reported in the literature for the fabrication of biodegradable scaffolds, but they suffer several limitations. In this study, for the first time, the possibility of generating 3D polymeric scaffolds loaded with an active compound by supercritical freeze extraction process is evaluated; this innovative process combines the advantages of the thermally induced phase separation process and of the supercritical carbon dioxide drying. Poly-L-lactid acid/ibuprofen composite scaffolds characterized by a 3D geometry, micrometric cellular structures and wrinkled pores walls have been obtained; moreover, homogeneous drug distribution and controlled release of the active principle have been assured. PMID:24366467
Dynamic 3D micropatterned cell co-cultures within photocurable and chemically degradable hydrogels
Sugiura, Shinji; Cha, Jae Min; Yanagawa, Fumiki; Zorlutuna, Pinar; Bae, Hojae; Khademhosseini, Ali
2014-01-01
In this paper we report on the development of dynamically controlled 3D micropatterned cellular co-cultures within photocurable and chemically degradable hydrogels. Specifically, we generated dynamic co-cultures of micropatterned murine embryonic stem (mES) cells with human hepatocellular carcinoma (HepG2) cells within 3D hydrogels. HepG2 cells were used due to their ability to direct the differentiation of mES cells through secreted paracrine factors. To generate dynamic co-cultures, mES cells were first encapsulated within micropatterned photocurable poly(ethylene glycol) (PEG) hydrogels. These micropatterned cell-laden PEG hydrogels were subsequently surrounded by calcium alginate (Ca-Alg) hydrogels containing HepG2 cells. After 4 days, the co-culture step was halted by exposing the system to sodium citrate solution, which removed the alginate gels and the encapsulated HepG2 cells. The encapsulated mES cells were then maintained in the resulting cultures for 16 days and cardiac differentiation was analyzed. We observed that the mES cells that were exposed to HepG2 cells in the co-cultures, generated cells with higher expression of cardiac genes and proteins as well as increased spontaneous beating. Due to its ability to control the 3D microenvironment of cells in a spatially and temporally regulated manner the method presented in this study is useful for a range of cell culture applications related to tissue engineering and regenerative medicine. PMID:24170301
Dynamical D4-D8 and D3-D7 branes in supergravity
Binetruy, Pierre; Sasaki, Misao; Uzawa, Kunihito
2009-07-15
We present a class of dynamical solutions for intersecting D4-D8 and D3-D7 brane systems in ten-dimensional type IIA and IIB supergravity. We discuss if these solutions can be recovered in lower-dimensional effective theories for the warped compactification of a general p-brane system. It is found that an effective p+1-dimensional description is not possible in general due to the entanglement of the transverse coordinates and the p+1-dimensional coordinates in the metric components. For the D4-D8 brane system, the dynamical solutions reduces to a static warped AdS{sub 6}xS{sup 4} geometry in a certain spacetime region. For the D3-D7 brane system, we find a dynamical solution whose metric form is similar to that of a D3-brane solution. The main difference is the existence of a nontrivial dilaton configuration in the D3-D7 solution. Then we discuss cosmology of these solutions. We find that they behave like a Kasner-type cosmological solution at {tau}{yields}{infinity}, while it reduces to a warped static solution at {tau}{yields}0, where {tau} is the cosmic time.
Dynamical D4-D8 and D3-D7 branes in supergravity
NASA Astrophysics Data System (ADS)
Binetruy, Pierre; Sasaki, Misao; Uzawa, Kunihito
2009-07-01
We present a class of dynamical solutions for intersecting D4-D8 and D3-D7 brane systems in ten-dimensional type IIA and IIB supergravity. We discuss if these solutions can be recovered in lower-dimensional effective theories for the warped compactification of a general p-brane system. It is found that an effective p+1-dimensional description is not possible in general due to the entanglement of the transverse coordinates and the p+1-dimensional coordinates in the metric components. For the D4-D8 brane system, the dynamical solutions reduces to a static warped AdS6×S4 geometry in a certain spacetime region. For the D3-D7 brane system, we find a dynamical solution whose metric form is similar to that of a D3-brane solution. The main difference is the existence of a nontrivial dilaton configuration in the D3-D7 solution. Then we discuss cosmology of these solutions. We find that they behave like a Kasner-type cosmological solution at τ→∞, while it reduces to a warped static solution at τ→0, where τ is the cosmic time.
PRONTO3D users` instructions: A transient dynamic code for nonlinear structural analysis
Attaway, S.W.; Mello, F.J.; Heinstein, M.W.; Swegle, J.W.; Ratner, J.A.; Zadoks, R.I.
1998-06-01
This report provides an updated set of users` instructions for PRONTO3D. PRONTO3D is a three-dimensional, transient, solid dynamics code for analyzing large deformations of highly nonlinear materials subjected to extremely high strain rates. This Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion. Eight-node, uniform strain, hexahedral elements and four-node, quadrilateral, uniform strain shells are used in the finite element formulation. An adaptive time step control algorithm is used to improve stability and performance in plasticity problems. Hourglass distortions can be eliminated without disturbing the finite element solution using either the Flanagan-Belytschko hourglass control scheme or an assumed strain hourglass control scheme. All constitutive models in PRONTO3D are cast in an unrotated configuration defined using the rotation determined from the polar decomposition of the deformation gradient. A robust contact algorithm allows for the impact and interaction of deforming contact surfaces of quite general geometry. The Smooth Particle Hydrodynamics method has been embedded into PRONTO3D using the contact algorithm to couple it with the finite element method.
Can symmetry transitions of complex fields enable 3-d control of fluid vorticity?
Martin, James E.; Solis, Kyle Jameson
2015-08-01
Methods of inducing vigorous noncontact fluid flow are important to technologies involving heat and mass transfer and fluid mixing, since they eliminate the need for moving parts, pipes and seals, all of which compromise system reliability. Unfortunately, traditional noncontact flow methods are few, and have limitations of their own. We have discovered two classes of fields that can induce fluid vorticity without requiring either gravity or a thermal gradient. The first class we call Symmetry-Breaking Rational Fields. These are triaxial fields comprised of three orthogonal components, two ac and one dc. The second class is Rational Triad Fields, which differ in that all three components are alternating. In this report we quantify the induced vorticity for a wide variety of fields and consider symmetry transitions between these field types. These transitions give rise to orbiting vorticity vectors, a technology for non-contact, non-stationary fluid mixing.
Dynamic WIFI-Based Indoor Positioning in 3D Virtual World
NASA Astrophysics Data System (ADS)
Chan, S.; Sohn, G.; Wang, L.; Lee, W.
2013-11-01
A web-based system based on the 3DTown project was proposed using Google Earth plug-in that brings information from indoor positioning devices and real-time sensors into an integrated 3D indoor and outdoor virtual world to visualize the dynamics of urban life within the 3D context of a city. We addressed limitation of the 3DTown project with particular emphasis on video surveillance camera used for indoor tracking purposes. The proposed solution was to utilize wireless local area network (WLAN) WiFi as a replacement technology for localizing objects of interest due to the wide spread availability and large coverage area of WiFi in indoor building spaces. Indoor positioning was performed using WiFi without modifying existing building infrastructure or introducing additional access points (AP)s. A hybrid probabilistic approach was used for indoor positioning based on previously recorded WiFi fingerprint database in the Petrie Science and Engineering building at York University. In addition, we have developed a 3D building modeling module that allows for efficient reconstruction of outdoor building models to be integrated with indoor building models; a sensor module for receiving, distributing, and visualizing real-time sensor data; and a web-based visualization module for users to explore the dynamic urban life in a virtual world. In order to solve the problems in the implementation of the proposed system, we introduce approaches for integration of indoor building models with indoor positioning data, as well as real-time sensor information and visualization on the web-based system. In this paper we report the preliminary results of our prototype system, demonstrating the system's capability for implementing a dynamic 3D indoor and outdoor virtual world that is composed of discrete modules connected through pre-determined communication protocols.
NASA Technical Reports Server (NTRS)
Gayman, W. H.
1974-01-01
Test method and apparatus determine fluid effective mass and damping in frequency range where effective mass may be considered as total mass less sum of slosh masses. Apparatus is designed so test tank and its mounting yoke are supported from structural test wall by series of flexures.
Berney, Sandra; Bétrancourt, Mireille; Molinari, Gaëlle; Hoyek, Nady
2015-01-01
The emergence of dynamic visualizations of three-dimensional (3D) models in anatomy curricula may be an adequate solution for spatial difficulties encountered with traditional static learning, as they provide direct visualization of change throughout the viewpoints. However, little research has explored the interplay between learning material presentation formats, spatial abilities, and anatomical tasks. First, to understand the cognitive challenges a novice learner would be faced with when first exposed to 3D anatomical content, a six-step cognitive task analysis was developed. Following this, an experimental study was conducted to explore how presentation formats (dynamic vs. static visualizations) support learning of functional anatomy, and affect subsequent anatomical tasks derived from the cognitive task analysis. A second aim was to investigate the interplay between spatial abilities (spatial visualization and spatial relation) and presentation formats when the functional anatomy of a 3D scapula and the associated shoulder flexion movement are learned. Findings showed no main effect of the presentation formats on performances, but revealed the predictive influence of spatial visualization and spatial relation abilities on performance. However, an interesting interaction between presentation formats and spatial relation ability for a specific anatomical task was found. This result highlighted the influence of presentation formats when spatial abilities are involved as well as the differentiated influence of spatial abilities on anatomical tasks. PMID:25689057
Temporal dynamics and developmental memory of 3D chromatin architecture at Hox gene loci
Noordermeer, Daan; Leleu, Marion; Schorderet, Patrick; Joye, Elisabeth; Chabaud, Fabienne; Duboule, Denis
2014-01-01
Hox genes are essential regulators of embryonic development. Their step-wise transcriptional activation follows their genomic topology and the various states of activation are subsequently memorized into domains of progressively overlapping gene products. We have analyzed the 3D chromatin organization of Hox clusters during their early activation in vivo, using high-resolution circular chromosome conformation capture. Initially, Hox clusters are organized as single chromatin compartments containing all genes and bivalent chromatin marks. Transcriptional activation is associated with a dynamic bi-modal 3D organization, whereby the genes switch autonomously from an inactive to an active compartment. These local 3D dynamics occur within a framework of constitutive interactions within the surrounding Topological Associated Domains, indicating that this regulation process is mostly cluster intrinsic. The step-wise progression in time is fixed at various body levels and thus can account for the chromatin architectures previously described at a later stage for different anterior to posterior levels. DOI: http://dx.doi.org/10.7554/eLife.02557.001 PMID:24843030
Using articulated scene models for dynamic 3d scene analysis in vista spaces
NASA Astrophysics Data System (ADS)
Beuter, Niklas; Swadzba, Agnes; Kummert, Franz; Wachsmuth, Sven
2010-09-01
In this paper we describe an efficient but detailed new approach to analyze complex dynamic scenes directly in 3D. The arising information is important for mobile robots to solve tasks in the area of household robotics. In our work a mobile robot builds an articulated scene model by observing the environment in the visual field or rather in the so-called vista space. The articulated scene model consists of essential knowledge about the static background, about autonomously moving entities like humans or robots and finally, in contrast to existing approaches, information about articulated parts. These parts describe movable objects like chairs, doors or other tangible entities, which could be moved by an agent. The combination of the static scene, the self-moving entities and the movable objects in one articulated scene model enhances the calculation of each single part. The reconstruction process for parts of the static scene benefits from removal of the dynamic parts and in turn, the moving parts can be extracted more easily through the knowledge about the background. In our experiments we show, that the system delivers simultaneously an accurate static background model, moving persons and movable objects. This information of the articulated scene model enables a mobile robot to detect and keep track of interaction partners, to navigate safely through the environment and finally, to strengthen the interaction with the user through the knowledge about the 3D articulated objects and 3D scene analysis. [Figure not available: see fulltext.
NASA Astrophysics Data System (ADS)
Jung, Jaewoon; Kobayashi, Chigusa; Imamura, Toshiyuki; Sugita, Yuji
2016-03-01
Three-dimensional Fast Fourier Transform (3D FFT) plays an important role in a wide variety of computer simulations and data analyses, including molecular dynamics (MD) simulations. In this study, we develop hybrid (MPI+OpenMP) parallelization schemes of 3D FFT based on two new volumetric decompositions, mainly for the particle mesh Ewald (PME) calculation in MD simulations. In one scheme, (1d_Alltoall), five all-to-all communications in one dimension are carried out, and in the other, (2d_Alltoall), one two-dimensional all-to-all communication is combined with two all-to-all communications in one dimension. 2d_Alltoall is similar to the conventional volumetric decomposition scheme. We performed benchmark tests of 3D FFT for the systems with different grid sizes using a large number of processors on the K computer in RIKEN AICS. The two schemes show comparable performances, and are better than existing 3D FFTs. The performances of 1d_Alltoall and 2d_Alltoall depend on the supercomputer network system and number of processors in each dimension. There is enough leeway for users to optimize performance for their conditions. In the PME method, short-range real-space interactions as well as long-range reciprocal-space interactions are calculated. Our volumetric decomposition schemes are particularly useful when used in conjunction with the recently developed midpoint cell method for short-range interactions, due to the same decompositions of real and reciprocal spaces. The 1d_Alltoall scheme of 3D FFT takes 4.7 ms to simulate one MD cycle for a virus system containing more than 1 million atoms using 32,768 cores on the K computer.
Härmä, Ville; Schukov, Hannu-Pekka; Happonen, Antti; Ahonen, Ilmari; Virtanen, Johannes; Siitari, Harri; Åkerfelt, Malin; Lötjönen, Jyrki; Nees, Matthias
2014-01-01
Glandular epithelial cells differentiate into complex multicellular or acinar structures, when embedded in three-dimensional (3D) extracellular matrix. The spectrum of different multicellular morphologies formed in 3D is a sensitive indicator for the differentiation potential of normal, non-transformed cells compared to different stages of malignant progression. In addition, single cells or cell aggregates may actively invade the matrix, utilizing epithelial, mesenchymal or mixed modes of motility. Dynamic phenotypic changes involved in 3D tumor cell invasion are sensitive to specific small-molecule inhibitors that target the actin cytoskeleton. We have used a panel of inhibitors to demonstrate the power of automated image analysis as a phenotypic or morphometric readout in cell-based assays. We introduce a streamlined stand-alone software solution that supports large-scale high-content screens, based on complex and organotypic cultures. AMIDA (Automated Morphometric Image Data Analysis) allows quantitative measurements of large numbers of images and structures, with a multitude of different spheroid shapes, sizes, and textures. AMIDA supports an automated workflow, and can be combined with quality control and statistical tools for data interpretation and visualization. We have used a representative panel of 12 prostate and breast cancer lines that display a broad spectrum of different spheroid morphologies and modes of invasion, challenged by a library of 19 direct or indirect modulators of the actin cytoskeleton which induce systematic changes in spheroid morphology and differentiation versus invasion. These results were independently validated by 2D proliferation, apoptosis and cell motility assays. We identified three drugs that primarily attenuated the invasion and formation of invasive processes in 3D, without affecting proliferation or apoptosis. Two of these compounds block Rac signalling, one affects cellular cAMP/cGMP accumulation. Our approach supports
Dynamic 3D imaging based on acousto-optic heterodyne fringe interferometry.
Guan, Yingjian; Yin, Yongkai; Li, Ameng; Liu, Xiaoli; Peng, Xiang
2014-06-15
An acoustic-optics heterodyne fringe interferometry coupled with a three-camera system is developed for dynamic 3D imaging. In this system, first-order beams with a slight frequency difference diffracted from two acousto-optic deflectors (AODs) form a beat intensity fringe pattern. Setting the frequency of the trigger signal for the CCD cameras into four times the beat frequency, four-step phase-shifting fringe patterns can be obtained, and the wrapped phase map (WPM) can be calculated. Under the epipolar constraint among three cameras, the homologous points can be determined unambiguously with the assistant of a WPM; thus the 3D shape can be reconstructed while skipping the phase unwrapping step. Experimental results are presented to validate this approach. PMID:24978566
User's manuals for DYNA3D and DYNAP: nonlinear dynamic analysis of solids in three dimensions
Hallquist, J.O.
1981-07-01
This report provides a user's manual for DYNA3D, an explicit three-dimensional finite element code for analyzing the large deformation dynamic response of inelastic solids. A contact-impact algorithm permits gaps and sliding along material interfaces. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning without the need of transition regions. Spatial discretization is achieved by the use of 8-node solid elements, and the equations-of-motion are integrated by the central difference method. Post-processors for DYNA3D include GRAPE for plotting deformed shapes and stress contours and DYNAP for plotting time histories. A user's manual for DYNAP is also provided in this report.
Radial electric field 3D modeling for wire arrays driving dynamic hohlraums on Z.
Mock, Raymond Cecil
2007-06-01
The anode-cathode structure of the Z-machine wire array results in a higher negative radial electric field (Er) on the wires near the cathode relative to the anode. The magnitude of this field has been shown to anti-correlate with the axial radiation top/bottom symmetry in the DH (Dynamic Hohlraum). Using 3D modeling, the structure of this field is revealed for different wire-array configurations and for progressive mechanical alterations, providing insight for minimizing the negative Er on the wire array in the anode-to-cathode region of the DH. Also, the 3D model is compared to Sasorov's approximation, which describes Er at the surface of the wire in terms of wire-array parameters.
Blob Dynamics in 3D BOUT Simulations of Tokamak Edge Turbulence
Russell, D; D'Ippolito, D; Myra, J; Nevins, W; Xu, X
2004-08-23
Propagating filaments of enhanced plasma density, or blobs, observed in 3D numerical simulations of a diverted, neutral-fueled tokamak are studied. Fluctuations of vorticity, electrical potential {phi}, temperature T{sub e} and current density J{sub {parallel}} associated with the blobs have a dipole structure perpendicular to the magnetic field and propagate radially with large E {center_dot} B drift velocities (> 1 km/s). The simulation results are consistent with a 3D blob dynamics model that incorporates increased parallel plasma resistivity (from neutral cooling of the X-point region), blob disconnection from the divertor sheath, X-point closure of the current loops, and collisional physics to sustain the {phi}, T{sub e}, J{sub {parallel}} dipoles.
Nonlinear dynamics of Airy-vortex 3D wave packets: emission of vortex light waves.
Driben, Rodislav; Meier, Torsten
2014-10-01
The dynamics of 3D Airy-vortex wave packets is studied under the action of strong self-focusing Kerr nonlinearity. Emissions of nonlinear 3D waves out of the main wave packets with the topological charges were demonstrated. Because of the conservation of the total angular momentum, charges of the emitted waves are equal to those carried by the parental light structure. The rapid collapse imposes a severe limitation on the propagation of multidimensional waves in Kerr media. However, the structure of the Airy beam carrier allows the coupling of light from the leading, most intense peak into neighboring peaks and consequently strongly postpones the collapse. The dependence of the critical input amplitude for the appearance of a fast collapse on the beam width is studied for wave packets with zero and nonzero topological charges. Wave packets carrying angular momentum are found to be much more resistant to the rapid collapse. PMID:25360922
Hand Gesture Spotting Based on 3D Dynamic Features Using Hidden Markov Models
NASA Astrophysics Data System (ADS)
Elmezain, Mahmoud; Al-Hamadi, Ayoub; Michaelis, Bernd
In this paper, we propose an automatic system that handles hand gesture spotting and recognition simultaneously in stereo color image sequences without any time delay based on Hidden Markov Models (HMMs). Color and 3D depth map are used to segment hand regions. The hand trajectory will determine in further step using Mean-shift algorithm and Kalman filter to generate 3D dynamic features. Furthermore, k-means clustering algorithm is employed for the HMMs codewords. To spot meaningful gestures accurately, a non-gesture model is proposed, which provides confidence limit for the calculated likelihood by other gesture models. The confidence measures are used as an adaptive threshold for spotting meaningful gestures. Experimental results show that the proposed system can successfully recognize isolated gestures with 98.33% and meaningful gestures with 94.35% reliability for numbers (0-9).
An exact solution for the 3D MHD stagnation-point flow of a micropolar fluid
NASA Astrophysics Data System (ADS)
Borrelli, A.; Giantesio, G.; Patria, M. C.
2015-01-01
The influence of a non-uniform external magnetic field on the steady three dimensional stagnation-point flow of a micropolar fluid over a rigid uncharged dielectric at rest is studied. The total magnetic field is parallel to the velocity at infinity. It is proved that this flow is possible only in the axisymmetric case. The governing nonlinear partial differential equations are reduced to a system of ordinary differential equations by a similarity transformation, before being solved numerically. The effects of the governing parameters on the fluid flow and on the magnetic field are illustrated graphically and discussed.
3D kinematic and dynamic analysis of the front crawl tumble turn in elite male swimmers.
Puel, F; Morlier, J; Avalos, M; Mesnard, M; Cid, M; Hellard, P
2012-02-01
The aim of this study was to identify kinematic and dynamic variables related to the best tumble turn times (3mRTT, the turn time from 3-m in to 3-m out, independent variable) in ten elite male swimmers using a three-dimensional (3D) underwater analysis protocol and the Lasso (least absolute shrinkage and selection operator) as statistical method. For each swimmer, the best-time turn was analyzed with five stationary and synchronized underwater cameras. The 3D reconstruction was performed using the Direct Linear Transformation algorithm. An underwater piezoelectric 3D force platform completed the set-up to compute dynamic variables. Data were smoothed by the Savitzky-Golay filtering method. Three variables were considered relevant in the best Lasso model (3mRTT=2.58-0.425 RD+0.204 VPe+0.0046 TD): the head-wall distance where rotation starts (RD), the horizontal speed at the force peak (VPe), and the 3D length of the path covered during the turn (TD). Furthermore, bivariate analysis showed that upper body (CUBei) and lower limb extension indexes at first contact (CLLei) were also linked to the turn time (r=-0.65 and p<0.05 for both variables). Thus the best turn times were associated with a long RD, slower VPe and reduced TD. By an early transverse rotation, male elite swimmers reach the wall with a slightly flexed posture that results in fast extension. These swimmers opt for a movement that is oriented forward and they focus on reducing the distance covered. PMID:22176710
Gierke, Sarah; Wittmann, Torsten
2012-01-01
SUMMARY Background Epithelial remodeling, in which apical-basal polarized cells switch to a migratory phenotype, plays a central role in development and disease of multicellular organisms. Although dynamic microtubules (MTs) are required for directed migration on flat surfaces, how MT dynamics are controlled or contribute to epithelial remodeling in a more physiological three-dimensional (3D) environment is not understood. We use confocal live cell imaging to analyze MT function and dynamics during 3D epithelial morphogenesis and remodeling of polarized Madin-Darby canine kidney (MDCK) epithelial cells that undergo partial epithelial-to-mesenchymal transition (EMT) in response to hepatocyte growth factor (HGF). Results We find that HGF treatment increases MT growth rate before morphological changes are evident, and that large numbers of MTs grow into HGF-induced cell extensions independent of centrosome reorientation. Using lentivirus-mediated shRNA, we demonstrate that EB1, an adaptor protein that mediates recruitment of numerous other +TIP proteins to growing MT plus ends, is required for this HGF-induced MT reorganization. We further show that protrusion and adhesion dynamics are disorganized, and that vesicular trafficking to the tip of HGF-induced cell extensions is disrupted in EB1-depleted cells. Conclusions We conclude that EB1-mediated interactions with growing MTs are important to coordinate cell shape changes and directed migration into the surrounding extracellular matrix during epithelial remodeling in a physiological 3D environment. In contrast, EB1 is not required for the establishment or maintenance of apical-basal cell polarity, suggesting different functions of +TIPs and MTs in different types of cell polarity. PMID:22483942
Non-linear dynamics of viscous bilayers subjected to an electric field: 3D phase field simulations
NASA Astrophysics Data System (ADS)
Dritselis, Christos; Karapetsas, George; Bontozoglou, Vasilis
2014-11-01
The scope of this work is to investigate the non-linear dynamics of the electro-hydrodynamic instability of a bilayer of immiscible liquids. We consider the case of two viscous films which is separated from the top electrode by air. We assume that the liquids are perfect dielectrics and consider the case of both flat and patterned electrodes. We develop a computational model using the diffuse interface method and carry out 3D numerical simulations fully accounting for the flow and electric field in all phases. We perform a parametric study and investigate the influence of the electric properties of fluids, applied voltage and various geometrical characteristics of the mask. The authors acknowledge the support by the General Secretariat of Research and Technology of Greece under the action ``Supporting Postdoctoral Researchers'' (Grant Number PE8/906), co-funded by the European Social Fund and National Resources.
Prodanović, M; Lindquist, W B; Seright, R S
2006-06-01
Using oil-wet polyethylene core models, we present the development of robust throat finding techniques for the extraction, from X-ray microtomographic images, of a pore network description of porous media having porosity up to 50%. Measurements of volume, surface area, shape factor, and principal diameters are extracted for pores and area, shape factor and principal diameters for throats. We also present results on the partitioning of wetting and non-wetting phases in the pore space at fixed volume increments of the injected fluid during a complete cycle of drainage and imbibition. We compare these results with fixed fractional flow injection, where wetting and non-wetting phase are simultaneously injected at fixed volume ratio. Finally we demonstrate the ability to differentiate three fluid phases (oil, water, air) in the pore space. PMID:16364351
On the 3D Steady Flow of a Second Grade Fluid Past an Obstacle
NASA Astrophysics Data System (ADS)
Konieczny, Paweł; Kreml, Ondřej
2012-06-01
We study steady flow of a second grade fluid past an obstacle in three space dimensions. We prove existence of solution in weighted Lebesgue spaces with anisotropic weights and thus existence of the wake region behind the obstacle. We use properties of the fundamental Oseen tensor together with results achieved in Koch (Quad Mat 15:59-122, 2004) and properties of solutions to steady transport equation to get up to arbitrarily small ɛ the same decay as the Oseen fundamental solution.
Noninvasive Imaging of 3D Dynamics in Thickly Fluorescent Specimens Beyond the Diffraction Limit
Gao, Liang; Shao, Lin; Higgins, Christopher D.; Poulton, John S.; Peifer, Mark; Davidson, Michael W.; Wu, Xufeng; Goldstein, Bob; Betzig, Eric
2013-01-01
SUMMARY Optical imaging of the dynamics of living specimens involves tradeoffs between spatial resolution, temporal resolution, and phototoxicity, made more difficult in three-dimensions. Here, however, we report that rapid 3D dynamics can be studied beyond the diffraction limit in thick or densely fluorescent living specimens over many time points by combining ultra-thin planar illumination produced by scanned Bessel beams with superresolution structured illumination microscopy. We demonstrate in vivo karyotyping of chromosomes during mitosis and identify different dynamics for the actin cytoskeleton at the dorsal and ventral surfaces of fibroblasts. Compared to spinning disk confocal microscopy, we demonstrate substantially reduced photodamage when imaging rapid morphological changes in D. discoideum cells, as well as improved contrast and resolution at depth within developing C. elegans embryos. Bessel beam structured plane illumination thus promises new insights into complex biological phenomena that require 4D subcellular spatiotemporal detail in either a single or multicellular context. PMID:23217717
Quantification of Diaphragm Mechanics in Pompe Disease Using Dynamic 3D MRI
Mogalle, Katja; Perez-Rovira, Adria; Ciet, Pierluigi; Wens, Stephan C. A.; van Doorn, Pieter A.; Tiddens, Harm A. W. M.; van der Ploeg, Ans T.; de Bruijne, Marleen
2016-01-01
Background Diaphragm weakness is the main reason for respiratory dysfunction in patients with Pompe disease, a progressive metabolic myopathy affecting respiratory and limb-girdle muscles. Since respiratory failure is the major cause of death among adult patients, early identification of respiratory muscle involvement is necessary to initiate treatment in time and possibly prevent irreversible damage. In this paper we investigate the suitability of dynamic MR imaging in combination with state-of-the-art image analysis methods to assess respiratory muscle weakness. Methods The proposed methodology relies on image registration and lung surface extraction to quantify lung kinematics during breathing. This allows for the extraction of geometry and motion features of the lung that characterize the independent contribution of the diaphragm and the thoracic muscles to the respiratory cycle. Results Results in 16 3D+t MRI scans (10 Pompe patients and 6 controls) of a slow expiratory maneuver show that kinematic analysis from dynamic 3D images reveals important additional information about diaphragm mechanics and respiratory muscle involvement when compared to conventional pulmonary function tests. Pompe patients with severely reduced pulmonary function showed severe diaphragm weakness presented by minimal motion of the diaphragm. In patients with moderately reduced pulmonary function, cranial displacement of posterior diaphragm parts was reduced and the diaphragm dome was oriented more horizontally at full inspiration compared to healthy controls. Conclusion Dynamic 3D MRI provides data for analyzing the contribution of both diaphragm and thoracic muscles independently. The proposed image analysis method has the potential to detect less severe diaphragm weakness and could thus be used to determine the optimal start of treatment in adult patients with Pompe disease in prospect of increased treatment response. PMID:27391236
NASA Astrophysics Data System (ADS)
Makhfudz, Imam
2016-04-01
Axion electrodynamics, first proposed in the context of particle physics, manifests itself in condensed matter physics in the topological field theory description of 3 d topological insulators and gives rise to magnetoelectric effect, where applying magnetic (electric) field B (E ) induces polarization (magnetization) p (m ) . We use linear response theory to study the associated topological current using the Fu-Kane-Mele model of 3 d topological insulators in the presence of time-dependent uniform weak magnetic field. By computing the dynamical current susceptibility χij jpjp(ω ) , we discover from its static limit an `order parameter' of the topological phase transition between weak topological (or ordinary) insulator and strong topological insulator, found to be continuous. The χij jpjp(ω ) shows a sign-changing singularity at a critical frequency with suppressed strength in the topological insulating state. Our results can be verified in current noise experiment on 3 d TI candidate materials for the detection of such topological phase transition.
Semi-automatic segmentation for 3D motion analysis of the tongue with dynamic MRI.
Lee, Junghoon; Woo, Jonghye; Xing, Fangxu; Murano, Emi Z; Stone, Maureen; Prince, Jerry L
2014-12-01
Dynamic MRI has been widely used to track the motion of the tongue and measure its internal deformation during speech and swallowing. Accurate segmentation of the tongue is a prerequisite step to define the target boundary and constrain the tracking to tissue points within the tongue. Segmentation of 2D slices or 3D volumes is challenging because of the large number of slices and time frames involved in the segmentation, as well as the incorporation of numerous local deformations that occur throughout the tongue during motion. In this paper, we propose a semi-automatic approach to segment 3D dynamic MRI of the tongue. The algorithm steps include seeding a few slices at one time frame, propagating seeds to the same slices at different time frames using deformable registration, and random walker segmentation based on these seed positions. This method was validated on the tongue of five normal subjects carrying out the same speech task with multi-slice 2D dynamic cine-MR images obtained at three orthogonal orientations and 26 time frames. The resulting semi-automatic segmentations of a total of 130 volumes showed an average dice similarity coefficient (DSC) score of 0.92 with less segmented volume variability between time frames than in manual segmentations. PMID:25155697
Geophysical Fluid Dynamics Outreach Films
NASA Astrophysics Data System (ADS)
Aurnou, J. M.; Schwarz, J. W.; Noguez, G.
2012-12-01
Here we will present high definition films of laboratory experiments demonstrating basic fluid motions similar to those occurring in atmospheres and oceans. In these experiments, we use water to simulate the fluid dynamics of both the liquid (oceans) and gaseous (atmospheric) envelopes. To simulate the spinning of the earth, we carry out the experiments on a rotating table. For each experiment, we begin by looking at our system first without the effects of rotation. Then, we include rotation to see how the behavior of the fluid changes due to the Coriolis accelerations. Our hope is that by viewing these experiments one will develop a sense for how fluids behave both in rotating and non-rotating systems. By noting the differences between the experiments, it should then be possible to establish a basis to think about large-scale fluid motions that exist in Earth's oceans and atmospheres as well as on planets other than Earth.Plan view image of vortices in a rotating tank of fluid. Movies of such flows make accessible the often difficult to comprehend fluid dynamical processes that occur in planetary atmospheres and oceans.
An easy implementation of displacement calculations in 3D discrete dislocation dynamics codes
NASA Astrophysics Data System (ADS)
Fivel, Marc; Depres, Christophe
2014-10-01
Barnett's coordinate-free expression of the displacement field of a triangular loop in an isotropic media is revisited in a view to be implemented in 3D discrete dislocation dynamics codes. A general meshing procedure solving the problems of non-planar loops is presented. The method is user-friendly and can be used in numerical simulations since it gives the contribution of each dislocation segment to the global displacement field without defining the connectivity of closed loops. Easy to implement in parallel calculations, this method is successfully applied to large-scale simulations.
Parallel 3-D particle-in-cell modelling of charged ultrarelativistic beam dynamics
NASA Astrophysics Data System (ADS)
Boronina, Marina A.; Vshivkov, Vitaly A.
2015-12-01
> ) in supercolliders. We use the 3-D set of Maxwell's equations for the electromagnetic fields, and the Vlasov equation for the distribution function of the beam particles. The model incorporates automatically the longitudinal effects, which can play a significant role in the cases of super-high densities. We present numerical results for the dynamics of two focused ultrarelativistic beams with a size ratio 10:1:100. The results demonstrate high efficiency of the proposed computational methods and algorithms, which are applicable to a variety of problems in relativistic plasma physics.
Time resolved 3D momentum imaging of ultrafast dynamics by coherent VUV-XUV radiation.
Sturm, F P; Wright, T W; Ray, D; Zalyubovskaya, I; Shivaram, N; Slaughter, D S; Ranitovic, P; Belkacem, A; Weber, Th
2016-06-01
We present a new experimental setup for measuring ultrafast nuclear and electron dynamics of molecules after photo-excitation and ionization. We combine a high flux femtosecond vacuum ultraviolet (VUV) and extreme ultraviolet (XUV) source with an internally cold molecular beam and a 3D momentum imaging particle spectrometer to measure electrons and ions in coincidence. We describe a variety of tools developed to perform pump-probe studies in the VUV-XUV spectrum and to modify and characterize the photon beam. First benchmark experiments are presented to demonstrate the capabilities of the system. PMID:27370429
Simulation studies of defect textures and dynamics in 3-d cholesteric droplets
NASA Astrophysics Data System (ADS)
Gimenez-Pinto, Vianney; Lu, Shin-Ying; Selinger, Jonathan; Selinger, Robin
2010-03-01
We model defect texture evolution in droplets of cholesteric liquid crystals by solving for the dynamics of the nematic director field. In order to accommodate defects in the simulated texture, we use a finite difference formulation that is explicitly independent of sign reversal of the director at any position in the sample. Textures are visualized using either the Berreman 4x4 matrix method or by mapping free energy density. We study both planar and focal conic cholesteric textures in 3-d spherical and cylindrical droplets, with the goal to optimize device geometries for bistable display applications.
Simulation studies of dynamics and defect textures in 3-d cholesteric droplets
NASA Astrophysics Data System (ADS)
Gimenez-Pinto, Vianney; Lu, Shin-Ying; Selinger, Jonathan; Selinger, Robin
2010-04-01
We model defect texture evolution in droplets of cholesteric liquid crystals by solving for the dynamics of the nematic director field. In order to accommodate defects in the simulated texture, we use a finite difference formulation that is explicitly independent of sign reversal of the director at any position in the sample. Textures are visualized using either the Berreman 4x4 matrix method or by mapping free energy density. We study both planar and focal conic cholesteric textures in 3-d spherical and cylindrical droplets, with the goal to optimize device geometries for bistable display applications.
Dynamic Characteristics of a Model and Prototype for 3D-RC Structure
NASA Astrophysics Data System (ADS)
Moniuddin, Md. Khaja; Vasanthalakshmi, G.; Chethan, K.; Babu, R. Ramesh
2016-06-01
Infill walls provide durable and economical partitions that have relatively excellent thermal and sound insulation with high fire resistance. Monolithic infilled walls are provided within RC structures without being analyzed as a combination of concrete and brick elements, although in reality they act as a single unit during earthquakes. The performance of such structures during earthquakes has proved to be superior in comparison to bare frames in terms of stiffness, strength and energy dissipation. To know the dynamic characteristics of monolithic infill wall panels and masonry infill, modal, response spectrum and time history analyses have been carried out on a model and prototype of a 3D RC structure for a comparative study.
Time resolved 3D momentum imaging of ultrafast dynamics by coherent VUV-XUV radiation
NASA Astrophysics Data System (ADS)
Sturm, F. P.; Wright, T. W.; Ray, D.; Zalyubovskaya, I.; Shivaram, N.; Slaughter, D. S.; Ranitovic, P.; Belkacem, A.; Weber, Th.
2016-06-01
We present a new experimental setup for measuring ultrafast nuclear and electron dynamics of molecules after photo-excitation and ionization. We combine a high flux femtosecond vacuum ultraviolet (VUV) and extreme ultraviolet (XUV) source with an internally cold molecular beam and a 3D momentum imaging particle spectrometer to measure electrons and ions in coincidence. We describe a variety of tools developed to perform pump-probe studies in the VUV-XUV spectrum and to modify and characterize the photon beam. First benchmark experiments are presented to demonstrate the capabilities of the system.
Tunesi, Marta; Fusco, Federica; Fiordaliso, Fabio; Corbelli, Alessandro; Biella, Gloria; Raimondi, Manuela T.
2016-01-01
Frontotemporal lobar degeneration (FTLD) is a severe neurodegenerative disorder that is diagnosed with increasing frequency in clinical setting. Currently, no therapy is available and in addition the molecular basis of the disease are far from being elucidated. Consequently, it is of pivotal importance to develop reliable and cost-effective in vitro models for basic research purposes and drug screening. To this respect, recent results in the field of Alzheimer’s disease have suggested that a tridimensional (3D) environment is an added value to better model key pathologic features of the disease. Here, we have tried to add complexity to the 3D cell culturing concept by using a microfluidic bioreactor, where cells are cultured under a continuous flow of medium, thus mimicking the interstitial fluid movement that actually perfuses the body tissues, including the brain. We have implemented this model using a neuronal-like cell line (SH-SY5Y), a widely exploited cell model for neurodegenerative disorders that shows some basic features relevant for FTLD modeling, such as the release of the FTLD-related protein progranulin (PRGN) in specific vesicles (exosomes). We have efficiently seeded the cells on 3D scaffolds, optimized a disease-relevant oxidative stress experiment (by targeting mitochondrial function that is one of the possible FTLD-involved pathological mechanisms) and evaluated cell metabolic activity in dynamic culture in comparison to static conditions, finding that SH-SY5Y cells cultured in 3D scaffold are susceptible to the oxidative damage triggered by a mitochondrial-targeting toxin (6-OHDA) and that the same cells cultured in dynamic conditions kept their basic capacity to secrete PRGN in exosomes once recovered from the bioreactor and plated in standard 2D conditions. We think that a further improvement of our microfluidic system may help in providing a full device where assessing basic FTLD-related features (including PRGN dynamic secretion) that may
Lattice Boltzmann simulation of dynamics of plunge and pitch of 3D flexible wing
NASA Astrophysics Data System (ADS)
Qi, Dewei; Shyy, Wei
2008-11-01
The method of lattice Boltzmann (LB) simulation has been used to simulate fluid structures and motion of a flexible insect wing in a 3D space. In the method, a beam has been discretized into a chain of rigid segments. Each segment is connected through ball and socket joints at its ends. One segment may be bent and twisted with its neighboring segment. A constraint force is applied to each joint to ensure the solid structure moving as a whole flexible elastic body.We have demonstrated that the LB method is suitable for modeling of aerodynamics of insects flight at low Reynolds numbers. First, a simulation of plunging and pitching of a rigid wing is performed at Re=75 in a 2D space and the results of lift forces and flow structures are in excellent agreement with the previous results. Second, plunging and pitching of a flexible wing in span-wise direction is simulated at Re=136 in a 3D space. We found that when twisting elasticity is large enough the twisting angle could be controlled at a level of smaller than 0.2 degree. It is shown that as bending and twisting elasticity is large enough, the motion of flexible wing approaches that of a rigid membrane wing. The simulation results show that the optimization of flexibility in span-wise direction will benefit thrust and an intermediate level is favorable. The results are consistent with experimental finding.
Internal waves patterns in the wake of a 3D body towed in a two-layer fluid
NASA Astrophysics Data System (ADS)
Lacaze, Laurent; Mercier, Matthieu; Thual, Olivier; Paci, Alexandre
2014-11-01
Stratified flows over obstacles are important features in meteorology and oceanography. The characterization of these flows is crucial in order to propose models of geophysical processes such as mixing and ocean circulation or orographic drag in the atmosphere. For some specific stratification profiles, the energy of internal waves generated by the obstacle can be trapped at a given depth, at the base of the oceanic mixing layer or at the top of the atmospheric boundary layer for instance. This scenario can be modelled by a two-layer stratified fluid for which gravity waves spread at the interface between the two layers. The work presented here focuses on a two-layer flow over a 3D obstacle, or equivalently, an obstacle towed in a fluid at rest. Experiments performed both in the large-scale flume of CNRM-GAME Toulouse (METEO-FRANCE & CNRS) and in a smaller tank apparatus, are presented with a specific attention on the measurement of the 3D wave patterns. A non-hydrostatic linear analysis is used to describe the observed wave patterns. The experiments highlight the strong influence of the Froude number on the generated waves. More specifically, we investigate the nature of the wake angle obtained from the wave pattern, and discuss a transition from Kelvin to Mach angle.
Insights from 3D numerical simulations on the dynamics of the India-Asia collision zone
NASA Astrophysics Data System (ADS)
Pusok, A. E.; Kaus, B.; Popov, A.
2013-12-01
The dynamics of the India-Asia collision zone remains one of the most remarkable topics of the current research interest: the transition from subduction to collision and uplift, followed by the rise of the abnormally thick Tibetan plateau, and the deformation at its Eastern and Western syntaxes, are processes still not fully understood. Models that have addressed this topic include wholescale underthrusting of Indian lithospheric mantle under Tibet, distributed homogeneous shortening or the thin-sheet model, slip-line field model for lateral extrusion or lower crustal flow models for the exhumation of the Himalayan units and lateral spreading of the Tibetan plateau. Of these, the thin-sheet model has successfully illustrated some of the basic physics of continental collision and has the advantage of a 3D model being reduced to 2D, but one of its major shortcomings is that it cannot simultaneously represent channel flow and gravitational collapse of the mantle lithosphere, since these mechanisms require the lithosphere to interact with the underlying mantle, or to have a vertically non-homogeneous rheology. As a consequence, 3D models are emerging as powerful tools to understand the dynamics of coupled systems. However, because of yet recent developments and various complexities, the current 3D models simulating the dynamics of continent collision zones have relied on certain explicit assumptions, such as replacing part of the asthenosphere with various types of boundary conditions that mimic the effect of mantle flow, in order to focus on the lithospheric/crustal deformation. Here, we employ the parallel 3D code LaMEM (Lithosphere and Mantle Evolution Model), with a finite difference staggered grid solver, which is capable of simulating lithospheric deformation while simultaneously taking mantle flow and a free surface into account. We present qualitative results on lithospheric and upper-mantle scale simulations in which the Indian lithosphere is subducted and
Rapid 3D dynamic arterial spin labeling with a sparse model-based image reconstruction.
Zhao, Li; Fielden, Samuel W; Feng, Xue; Wintermark, Max; Mugler, John P; Meyer, Craig H
2015-11-01
Dynamic arterial spin labeling (ASL) MRI measures the perfusion bolus at multiple observation times and yields accurate estimates of cerebral blood flow in the presence of variations in arterial transit time. ASL has intrinsically low signal-to-noise ratio (SNR) and is sensitive to motion, so that extensive signal averaging is typically required, leading to long scan times for dynamic ASL. The goal of this study was to develop an accelerated dynamic ASL method with improved SNR and robustness to motion using a model-based image reconstruction that exploits the inherent sparsity of dynamic ASL data. The first component of this method is a single-shot 3D turbo spin echo spiral pulse sequence accelerated using a combination of parallel imaging and compressed sensing. This pulse sequence was then incorporated into a dynamic pseudo continuous ASL acquisition acquired at multiple observation times, and the resulting images were jointly reconstructed enforcing a model of potential perfusion time courses. Performance of the technique was verified using a numerical phantom and it was validated on normal volunteers on a 3-Tesla scanner. In simulation, a spatial sparsity constraint improved SNR and reduced estimation errors. Combined with a model-based sparsity constraint, the proposed method further improved SNR, reduced estimation error and suppressed motion artifacts. Experimentally, the proposed method resulted in significant improvements, with scan times as short as 20s per time point. These results suggest that the model-based image reconstruction enables rapid dynamic ASL with improved accuracy and robustness. PMID:26169322
Incompressible Viscous Fluid Dynamics
Energy Science and Technology Software Center (ESTSC)
1992-02-13
NACHOS2 is a finite element program designed for the analysis of two-dimensional, incompressible viscous fluid flow problems. The basic flows considered may be isothermal, nonisothermal, or may involve other physical processes, such as mass transport. Both steady and transient flows may be analyzed. The class of problems treated are those described by the two-dimensional (plane or axisymmetric) incompressible form of the Navier-Stokes equations. An energy transport equation is included in the formulation for problems inmore » which heat transfer effects are important. Two auxiliary transport equations can be added to describe other physical processes,e.g. mass transfer, chemical reactions. Among the specific types of flow problems treated are: isothermal flow; forced, free, or mixed convection; conjugate heat transfer; flow in saturated porous media with or without heat transfer; and inelastic, non-Newtonian flows with or without heat transfer. Other problem classes are possible depending on the specific definitions applied to the auxiliary transport equations.« less
NASA Astrophysics Data System (ADS)
Hu, Guangli; Liu, Junwei; Graham-Brady, Lori; Ramesh, K. T.
2015-05-01
We present a validated fully 3D mechanism-based micromechanical constitutive model for brittle solids under dynamic multiaxial loading conditions. Flaw statistics are explicitly incorporated through a defect density, and evolving flaw distributions in both orientation and size. Interactions among cracks are modeled by means of a crack-matrix-effective-medium approach. A tensorial damage parameter is defined based upon the crack length and orientation development under local effective stress fields. At low confining stresses, the wing-cracking mechanism dominates, leading to the degradation of the modulus and peak strength of the material, whereas at high enough confining stresses, the cracking mechanism is completely shut-down and dislocation mechanisms become dominant. The model handles general multiaxial stress states, accounts for evolving internal variables in the form of evolving flaw size and orientation distributions, includes evolving anisotropic damage and irreversible damage strains in a thermodynamically consistent fashion, incorporates rate-dependence through the micromechanics, and includes dynamic bulking based on independent experimental data. Simulation results are discussed and compared with experimental results on one specific structural ceramic, aluminum nitride. We demonstrate that this 3D constitutive model is capable of capturing the general constitutive response of structural ceramics.
Grounding line dynamics inferred from a 3D full-Stokes model solving the contact problem
NASA Astrophysics Data System (ADS)
Favier, Lionel; Gagliardini, Olivier; Durand, Gael; Zwinger, Thomas
2010-05-01
The mass balance of marine ice-sheets, such as the West Antarctic Ice Sheet, is mostly controlled by their grounding line dynamics. Most numerical models simulating marine ice-sheets involve simplifications and do not include all the stress gradients. First results obtained with a 3D full-Stokes model for the grounded ice-sheet / floating ice-shelf transition, using the finite-element code Elmer/Ice, are presented. The initial geometry, which takes into account a dome and a calving front, has been laterally extruded from a previously investigated 2D flowline geometry. The grounding line migration is computed by solving the contact problem between the ice and the rigid downward sloping bedrock, where a non linear friction law is applied in the two horizontal directions. The evolutions of the sea-air and sea-ice interfaces are determined by the solution of a local transport equation. The consistency between the 3D model and the analogous results of the flowline model is shown by comparing the results in the basic extruded case, with no normal flux through lateral boundaries. Thereafter, spatially non uniform perturbations are introduced, to simulate the grounding line dynamics under fully three-dimensional perturbations.
NASA Astrophysics Data System (ADS)
Petersen, C.; Buenz, S.; Hustoft, S.; Mienert, J.
2007-12-01
High-resolution seismic data were acquired using the 3D seismic P-Cable system of the University of Tromsoe to investigate how the fluid flow penetrates gas hydrate systems of the Vestnesa Ridge. The ridge represents a current-controlled sediment drift on the continental margin offshore western Svalbard. The survey area is located at the northwestern part of the Vestnesa Ridge and centered at the ridge crest that resembles an anticline in a water depth of 1250-1320 m. The seafloor morphology at the crest is characterized by an abundance of pockmarks with a diameter between 50-500 m indicating recent fluid flow activity. Since the area is within the gas hydrate stability zone (GHSZ), it is an ideal site to understand where and how fluids escape through a hydrated sediment drift. 35 reflection seismic profiles with a spacing of about 40-60 m were shot resulting in a seismic cube covering an area of approximately 22 km2. In addition, regional single channel streamer (SCS) seismic lines were acquired across the ridge perpendicular to the crest to connect the 3D area with the regional structural setting. The seismic data provide images of the subsurface to about 500 ms TWT (two-way time) below the seafloor (bsf), where gas accumulations cause acoustic attenuations that hinder deeper acoustic signal penetration. The well-stratified sediments exhibit a bottom simulating reflector (BSR) at about 200 ms TWT bsf at the base of the GHSZ. The BSR is difficult to identify due to the stratification, but it is accompanied by the onset of an ubiquitous band of strong reflectivity indicating free gas accumulation zones beneath the GHSZ. Fluid flow activity is evident from a link between gas accumulations (bright spots), gas wipeouts and disturbed reflectivity in the seismic data. These features are observed not only beneath the pockmark structures, but also in the sediment without seafloor expressions of fluid venting. The fluid source might be related to deep tectonic processes at
3D microscale laser dynamic forming: Multiscale modeling and experimental validation
Gao Huang; Cheng, Gary J.
2011-05-15
Microscale laser dynamic forming ({mu}LDF) shows great potential in fabricating robust and high-aspect-ratio metallic microcomponents. Experiments revealed that strain rate and sample size play important roles in determining the dynamic plasticity and final results of {mu}LDF. To further understand these effects, a multiscale modeling methodology is adopted to characterize the microscale dynamic plasticity considering the evolutions of nano-to-submicron dislocations avalanches under shock loading. In this methodology, 3D discrete dislocation dynamics simulations are implemented to derive the yield strength and the initial strain hardening dependence on size and strain rate. It is observed that there exist three dynamic stages during deformation process. The initial strain hardening rate in Stage II increases with strain rate. The mechanical threshold stress model, intrinsically equipped with strain-rate-dependent flow stress and initial hardening, is chosen and modified to incorporate size effect quantitatively. This scale-dependent model, implemented in abaqus/explicit, provides deformation depths and thickness variations in good agreement with experimental results in {mu}LDF.
Ito, Takashi; Pilat, Marcin L; Suzuki, Reiji; Arita, Takaya
2016-01-01
Recent studies have reported that population dynamics and evolutionary dynamics, occurring at different time scales, can be affected by each other. Our purpose is to explore the interaction between population and evolutionary dynamics using an artificial life approach based on a 3D physically simulated environment in the context of predator-prey and morphology-behavior coevolution. The morphologies and behaviors of virtual prey creatures are evolved using a genetic algorithm based on the predation interactions between predators and prey. Both population sizes are also changed, depending on the fitness. We observe two types of cyclic behaviors, corresponding to short-term and long-term dynamics. The former can be interpreted as a simple population dynamics of Lotka-Volterra type. It is shown that the latter cycle is based on the interaction between the changes in the prey strategy against predators and the long-term change in both population sizes, resulting partly from a tradeoff between their defensive success and the cost of defense. PMID:26934093
Progress in geophysical fluid dynamics
NASA Astrophysics Data System (ADS)
Robinson, Allan R.
Geophysical fluid dynamics deals with the motions and physics of the atmosphere, oceans and interior of the earth and other planets: the winds, the swirls, the currents that occur on myriads of scales from millimeter to climatological. Explanations of natural phenomena, basic processes and abstractions are sought. The rotation of the earth, the buoyancy of its fluids and the tendency towards large-scale turbulence characterize these flows. But geophysical fluid dynamics is importantly a part of modern fluid dynamics which is contributing to the development of nonlinear mechanics generally. Some general insights are emerging for nonlinear systems which must be regarded as partly deterministic and partly random or which are complex and aperiodic. Contributions from geophysical fluid dynamics come from its methodology, from the experience of examples, and from the perspective provided by its unique scale. Contributions have been made to turbulent, chaotic and coherently structured nonlinear process research. Turbulent vortices larger than man himself naturally invite detailed investigation and deterministic physical studies. Examples are storms in the atmosphere and large ring vortices spun off by the Gulf Stream current in mid-ocean. The statistics of these events determine critical aspects of the general circulations. Fluid dynamicists generally now know that it is often relevant or necessary to study local dynamical processes of typical eddies even though only the average properties of the flow are of interest; progress in understanding the turbulent boundary layer in pipes involves the study of millimeter-scale vortices. Weather-related studies were seminal to the construction of the new scientific field of chaos. Coherent vortices abound of which the Great Red Spot of Jupiter is a spectacular example. Geophysical fluid dynamicists have been among forefront researchers in exploiting the steadily increasing speed and capacity of modern computers. Supercomputers
Three-fluid, 3D MHD solar wind modeling with turbulence transport and eddy viscosity
NASA Astrophysics Data System (ADS)
Usmanov, A. V.; Goldstein, M. L.; Matthaeus, W. H.
2014-12-01
We present results from a three-fluid, fully three-dimensional MHD solar wind model that includes turbulence transport, eddy viscosity, turbulent resistivity, and turbulent heating. The solar wind plasma is described as a co-moving system of three species: the solar wind protons, electrons, and interstellar pickup protons. Separate energy equations are employed for each species. We obtain numerical solutions of Reynolds-averaged solar wind equations coupled with turbulence transport equations in the region from 0.3 to 100 AU. The integrated system of equations includes the effects of electron heat conduction, Coulomb collisions, photoionization of interstellar hydrogen atoms and their charge exchange with the solar wind protons, turbulence energy generation by pickup protons, and turbulent heating of solar wind protons and electrons. Using either a dipole approximation for the solar magnetic field or synoptic solar magnetograms from the Wilcox Solar Observatory for assigning boundary conditions at the coronal base, we apply the model to study the global structure of the solar wind and its three-dimensional properties, including turbulence parameters, throughout the heliosphere. The model results are compared with observations on WIND, Ulysses and Voyager 2 spacecraft. This work is partially supported by LWS and Heliophysics Grand Challenges programs.
Tian, Fang-Bao; Dai, Hu; Luo, Haoxiang; Doyle, James F; Rousseau, Bernard
2014-02-01
Three-dimensional fluid-structure interaction (FSI) involving large deformations of flexible bodies is common in biological systems, but accurate and efficient numerical approaches for modeling such systems are still scarce. In this work, we report a successful case of combining an existing immersed-boundary flow solver with a nonlinear finite-element solid-mechanics solver specifically for three-dimensional FSI simulations. This method represents a significant enhancement from the similar methods that are previously available. Based on the Cartesian grid, the viscous incompressible flow solver can handle boundaries of large displacements with simple mesh generation. The solid-mechanics solver has separate subroutines for analyzing general three-dimensional bodies and thin-walled structures composed of frames, membranes, and plates. Both geometric nonlinearity associated with large displacements and material nonlinearity associated with large strains are incorporated in the solver. The FSI is achieved through a strong coupling and partitioned approach. We perform several validation cases, and the results may be used to expand the currently limited database of FSI benchmark study. Finally, we demonstrate the versatility of the present method by applying it to the aerodynamics of elastic wings of insects and the flow-induced vocal fold vibration. PMID:24415796
NASA Technical Reports Server (NTRS)
Walatka, Pamela P.; Buning, Pieter G.; Pierce, Larry; Elson, Patricia A.
1990-01-01
PLOT3D is a computer graphics program designed to visualize the grids and solutions of computational fluid dynamics. Seventy-four functions are available. Versions are available for many systems. PLOT3D can handle multiple grids with a million or more grid points, and can produce varieties of model renderings, such as wireframe or flat shaded. Output from PLOT3D can be used in animation programs. The first part of this manual is a tutorial that takes the reader, keystroke by keystroke, through a PLOT3D session. The second part of the manual contains reference chapters, including the helpfile, data file formats, advice on changing PLOT3D, and sample command files.
Electromagnetics and Fluid Dynamics
NASA Technical Reports Server (NTRS)
Gaitonde, Datta
1998-01-01
Previous efforts focused on developing tools for design of low observables were sustained. The final product was the maturation of a high-order accurate finite-volume based code to solve Maxwell's equations. One of the primary achievements was the development and implementation of efficient filtering techniques which enhance the robustness of high-order and optimized schemes without significant adverse impact on accuracy. This has eliminated the stability barrier which restrains the common use of high-order schemes for conservative wave propagation phenomena on curvilinear meshes. A study was performed of crossing shock interactions under conditions of increasing interaction strength and asymmetry. In the first category, the observed computed topological bifurcations were correlated with the formation of various lines of coalescence and divergence evident in experimental and computed surf-ace oil maps. ne flow structure arising from asymmetric interactions was investigated with particular emphasis on: 1) vorticity dynamics, 2) shock-structure and 3) sidewall vortex loading. Several efforts of the prior year were successfully published in archival journals. The high-order algorithms developed for CEM have been implemented into the FDL3DI CFD code are presently undergoing extensive testing. Preliminary results are highly encouraging.
3D Dynamics of the Near-Surface Layer of the Ocean in the Presence of Freshwater Influx
NASA Astrophysics Data System (ADS)
Dean, C.; Soloviev, A.
2015-12-01
Freshwater inflow due to convective rains or river runoff produces lenses of freshened water in the near surface layer of the ocean. These lenses are localized in space and typically involve both salinity and temperature anomalies. Due to significant density anomalies, strong pressure gradients develop, which result in lateral spreading of freshwater lenses in a form resembling gravity currents. Gravity currents inherently involve three-dimensional dynamics. The gravity current head can include the Kelvin-Helmholtz billows with vertical density inversions. In this work, we have conducted a series of numerical experiments using computational fluid dynamics tools. These numerical simulations were designed to elucidate the relationship between vertical mixing and horizontal advection of salinity under various environmental conditions and potential impact on the pollution transport including oil spills. The near-surface data from the field experiments in the Gulf of Mexico during the SCOPE experiment were available for validation of numerical simulations. In particular, we observed a freshwater layer within a few-meter depth range and, in some cases, a density inversion at the edge of the freshwater lens, which is consistent with the results of numerical simulations. In conclusion, we discuss applicability of these results to the interpretation of Aquarius and SMOS sea surface salinity satellite measurements. The results of this study indicate that 3D dynamics of the near-surface layer of the ocean are essential in the presence of freshwater inflow.
Fluid dynamic bowtie attenuators
NASA Astrophysics Data System (ADS)
Szczykutowicz, Timothy P.; Hermus, James
2015-03-01
Fluence field modulated CT allows for improvements in image quality and dose reduction. To date, only 1-D modulators have been proposed, the extension to 2-D modulation is difficult with solid-metal attenuation-based modulators. This work proposes to use liquids and gas to attenuate the x-ray beam which can be arrayed allowing for 2-D fluence modulation. The thickness of liquid and the pressure for a given path length of gas were determined that provided the same attenuation as 30 cm of soft tissue at 80, 100, 120, and 140 kV. Gaseous Xenon and liquid Iodine, Zinc Chloride, and Cerium Chloride were studied. Additionally, we performed some proof-of-concept experiments in which (1) a single cell of liquid was connected to a reservoir which allowed the liquid thickness to be modulated and (2) a 96 cell array was constructed in which the liquid thickness in each cell was adjusted manually. Liquid thickness varied as a function of kV and chemical composition, with Zinc Chloride allowing for the smallest thickness; 1.8, 2.25, 3, and 3.6 cm compensated for 30 cm of soft tissue at 80, 100, 120, and 140 kV respectively. The 96 cell Iodine attenuator allowed for a reduction in both dynamic range to the detector and scatter to primary ratio. Successful modulation of a single cell was performed at 0, 90, and 130 degrees using a simple piston/actuator. The thickness of liquids and the Xenon gas pressure seem logistically implementable within the constraints of CBCT and diagnostic CT systems.
NASA Astrophysics Data System (ADS)
Inogamov, Nail A.; Zhakhovsky, Vasily V.
2016-02-01
There are many important applications in which the ultrashort diffraction-limited and therefore tightly focused laser pulses irradiates metal films mounted on dielectric substrate. Here we present the detailed picture of laser peeling and 3D structure formation of the thin (relative to a depth of a heat affected zone in the bulk targets) gold films on glass substrate. The underlying physics of such diffraction-limited laser peeling was not well understood previously. Our approach is based on a physical model which takes into consideration the new calculations of the two-temperature (2T) equation of state (2T EoS) and the two-temperature transport coefficients together with the coupling parameter between electron and ion subsystems. The usage of the 2T EoS and the kinetic coefficients is required because absorption of an ultrashort pulse with duration of 10-1000 fs excites electron subsystem of metal and transfers substance into the 2T state with hot electrons (typical electron temperatures 1-3 eV) and much colder ions. It is shown that formation of submicrometer-sized 3D structures is a result of the electron-ion energy transfer, melting, and delamination of film from substrate under combined action of electron and ion pressures, capillary deceleration of the delaminated liquid metal or semiconductor, and ultrafast freezing of molten material. We found that the freezing is going in non-equilibrium regime with strongly overcooled liquid phase. In this case the Stefan approximation is non-applicable because the solidification front speed is limited by the diffusion rate of atoms in the molten material. To solve the problem we have developed the 2T Lagrangian code including all this reach physics in. We also used the high-performance combined Monte- Carlo and molecular dynamics code for simulation of surface 3D nanostructuring at later times after completion of electron-ion relaxation.
Yang, Yang; Pei, Zhiqiang; Li, Zhen; Wei, Yen; Ji, Yan
2016-02-24
Making dynamic three-dimensional (3D) structures capable of reversible shape changes or locomotion purely out of dry polymers is very difficult. Meanwhile, no previous dynamic 3D structures can be remade into new configurations while being resilient to mechanical damages and low temperature. Here, we show that light-activated transesterification in carbon nanotube dispersed liquid crystalline vitrimers enables flexible design and easy building of dynamic 3D structures out of flat films upon irradiation of light without screws, glues, or molds. Shining light also enables dynamic 3D structures to be quickly modified on demand, restored from distortion, repaired if broken, in situ healed when microcrack appears, assembled for more sophisticated structures, reconfigured, and recycled after use. Furthermore, the fabrication, reconfiguration, actuation, reparation, and assembly as well as healing can be performed even at extremely low temperatures (e.g., -130 °C). PMID:26840838
Computational fluid dynamics uses in fluid dynamics/aerodynamics education
NASA Technical Reports Server (NTRS)
Holst, Terry L.
1994-01-01
The field of computational fluid dynamics (CFD) has advanced to the point where it can now be used for the purpose of fluid dynamics physics education. Because of the tremendous wealth of information available from numerical simulation, certain fundamental concepts can be efficiently communicated using an interactive graphical interrogation of the appropriate numerical simulation data base. In other situations, a large amount of aerodynamic information can be communicated to the student by interactive use of simple CFD tools on a workstation or even in a personal computer environment. The emphasis in this presentation is to discuss ideas for how this process might be implemented. Specific examples, taken from previous publications, will be used to highlight the presentation.
Observing molecular dynamics with time-resolved 3D momentum imaging
NASA Astrophysics Data System (ADS)
Sturm, F. P.; Wright, T.; Bocharova, I.; Ray, D.; Shivaram, N.; Cryan, J.; Belkacem, A.; Weber, T.; Dörner, R.
2014-05-01
Photo-excitation and ionization trigger rich dynamics in molecular systems which play a key role in many important processes in nature such as vision, photosynthesis or photoprotection. Observing those reactions in real-time without significantly disturbing the molecules by a strong electric field has been a great challenge. Recent experiments using Time-of-Flight and Velocity Map Imaging techniques have revealed important information on the dynamics of small molecular systems upon photo-excitation. We have developed an apparatus for time-resolved momentum imaging of electrons and ions in all three spatial dimensions that employs two-color femtosecond laser pulses in the vacuum and extreme ultraviolet (VUV, XUV) for probing molecular dynamics. Our COLTRIMS style reaction microscope can measure electrons and ions in coincidence and reconstruct the momenta of the reaction fragments in 3D. We use a high power 800 nm laser in a loose focusing geometry gas cell to efficinetly drive High Harmonic Generation. The resulting photon flux is sufficient to perform 2-photon pump-probe experiments using VUV and XUV pulses for both pump and probe. With this setup we investigate non-Born-Oppenheimer dynamics in small molecules such as C2H4 and CO2 on a femtosecond time scale. Supported by Chemical Sciences, Geosciences and Biosciences division of BES/DOE.
McCullough, Dean P.; Gudla, Prabhakar R.; Harris, Bradley S.; Collins, Jason A.; Meaburn, Karen J.; Nakaya, Masa-Aki; Yamaguchi, Terry P.; Misteli, Tom; Lockett, Stephen J.
2009-01-01
Communications between cells in large part drive tissue development and function, as well as disease-related processes such as tumorigenesis. Understanding the mechanistic bases of these processes necessitates quantifying specific molecules in adjacent cells or cell nuclei of intact tissue. However, a major restriction on such analyses is the lack of an efficient method that correctly segments each object (cell or nucleus) from 3-D images of an intact tissue specimen. We report a highly reliable and accurate semi-automatic algorithmic method for segmenting fluorescence-labeled cells or nuclei from 3-D tissue images. Segmentation begins with semi-automatic, 2-D object delineation in a user-selected plane, using dynamic programming (DP) to locate the border with an accumulated intensity per unit length greater that any other possible border around the same object. Then the two surfaces of the object in planes above and below the selected plane are found using an algorithm that combines DP and combinatorial searching. Following segmentation, any perceived errors can be interactively corrected. Segmentation accuracy is not significantly affected by intermittent labeling of object surfaces, diffuse surfaces, or spurious signals away from surfaces. The unique strength of the segmentation method was demonstrated on a variety of biological tissue samples where all cells, including irregularly shaped cells, were accurately segmented based on visual inspection. PMID:18450544
Mutual information as a measure of image quality for 3D dynamic lung imaging with EIT
Crabb, M G; Davidson, J L; Little, R; Wright, P; Morgan, A R; Miller, C A; Naish, J H; Parker, G J M; Kikinis, R; McCann, H; Lionheart, W R B
2014-01-01
We report on a pilot study of dynamic lung electrical impedance tomography (EIT) at the University of Manchester. Low-noise EIT data at 100 frames per second (fps) were obtained from healthy male subjects during controlled breathing, followed by magnetic resonance imaging (MRI) subsequently used for spatial validation of the EIT reconstruction. The torso surface in the MR image and electrode positions obtained using MRI fiducial markers informed the construction of a 3D finite element model extruded along the caudal-distal axis of the subject. Small changes in the boundary that occur during respiration were accounted for by incorporating the sensitivity with respect to boundary shape into a robust temporal difference reconstruction algorithm. EIT and MRI images were co-registered using the open source medical imaging software, 3D Slicer. A quantitative comparison of quality of different EIT reconstructions was achieved through calculation of the mutual information with a lung-segmented MR image. EIT reconstructions using a linear shape correction algorithm reduced boundary image artefacts, yielding better contrast of the lungs, and had 10% greater mutual information compared with a standard linear EIT reconstruction. PMID:24710978
Description of patellar movement by 3D parameters obtained from dynamic CT acquisition
NASA Astrophysics Data System (ADS)
de Sá Rebelo, Marina; Moreno, Ramon Alfredo; Gobbi, Riccardo Gomes; Camanho, Gilberto Luis; de Ávila, Luiz Francisco Rodrigues; Demange, Marco Kawamura; Pecora, Jose Ricardo; Gutierrez, Marco Antonio
2014-03-01
The patellofemoral joint is critical in the biomechanics of the knee. The patellofemoral instability is one condition that generates pain, functional impairment and often requires surgery as part of orthopedic treatment. The analysis of the patellofemoral dynamics has been performed by several medical image modalities. The clinical parameters assessed are mainly based on 2D measurements, such as the patellar tilt angle and the lateral shift among others. Besides, the acquisition protocols are mostly performed with the leg laid static at fixed angles. The use of helical multi slice CT scanner can allow the capture and display of the joint's movement performed actively by the patient. However, the orthopedic applications of this scanner have not yet been standardized or widespread. In this work we present a method to evaluate the biomechanics of the patellofemoral joint during active contraction using multi slice CT images. This approach can greatly improve the analysis of patellar instability by displaying the physiology during muscle contraction. The movement was evaluated by computing its 3D displacements and rotations from different knee angles. The first processing step registered the images in both angles based on the femuŕs position. The transformation matrix of the patella from the images was then calculated, which provided the rotations and translations performed by the patella from its position in the first image to its position in the second image. Analysis of these parameters for all frames provided real 3D information about the patellar displacement.
Visualization of unsteady computational fluid dynamics
NASA Technical Reports Server (NTRS)
Haimes, Robert
1994-01-01
A brief summary of the computer environment used for calculating three dimensional unsteady Computational Fluid Dynamic (CFD) results is presented. This environment requires a super computer as well as massively parallel processors (MPP's) and clusters of workstations acting as a single MPP (by concurrently working on the same task) provide the required computational bandwidth for CFD calculations of transient problems. The cluster of reduced instruction set computers (RISC) is a recent advent based on the low cost and high performance that workstation vendors provide. The cluster, with the proper software can act as a multiple instruction/multiple data (MIMD) machine. A new set of software tools is being designed specifically to address visualizing 3D unsteady CFD results in these environments. Three user's manuals for the parallel version of Visual3, pV3, revision 1.00 make up the bulk of this report.
Visualization of unsteady computational fluid dynamics
NASA Astrophysics Data System (ADS)
Haimes, Robert
1994-11-01
A brief summary of the computer environment used for calculating three dimensional unsteady Computational Fluid Dynamic (CFD) results is presented. This environment requires a super computer as well as massively parallel processors (MPP's) and clusters of workstations acting as a single MPP (by concurrently working on the same task) provide the required computational bandwidth for CFD calculations of transient problems. The cluster of reduced instruction set computers (RISC) is a recent advent based on the low cost and high performance that workstation vendors provide. The cluster, with the proper software can act as a multiple instruction/multiple data (MIMD) machine. A new set of software tools is being designed specifically to address visualizing 3D unsteady CFD results in these environments. Three user's manuals for the parallel version of Visual3, pV3, revision 1.00 make up the bulk of this report.
Remote Visualization and Remote Collaboration On Computational Fluid Dynamics
NASA Technical Reports Server (NTRS)
Watson, Val; Lasinski, T. A. (Technical Monitor)
1995-01-01
A new technology has been developed for remote visualization that provides remote, 3D, high resolution, dynamic, interactive viewing of scientific data (such as fluid dynamics simulations or measurements). Based on this technology, some World Wide Web sites on the Internet are providing fluid dynamics data for educational or testing purposes. This technology is also being used for remote collaboration in joint university, industry, and NASA projects in computational fluid dynamics and wind tunnel testing. Previously, remote visualization of dynamic data was done using video format (transmitting pixel information) such as video conferencing or MPEG movies on the Internet. The concept for this new technology is to send the raw data (e.g., grids, vectors, and scalars) along with viewing scripts over the Internet and have the pixels generated by a visualization tool running on the viewer's local workstation. The visualization tool that is currently used is FAST (Flow Analysis Software Toolkit).
3D Plasma Clusters: Analysis of dynamical evolution and individual particle interaction
Antonova, T.; Thomas, H. M.; Morfill, G. E.; Annaratone, B. M.
2008-09-07
3D plasma clusters (up to 100 particles) have been built inside small (32 mm{sup 3}) plasma volume in gravity. It has been estimated that the external confinement has a negligible influence on the processes inside the clusters. At such conditions the analysis of dynamical evolution and individual particle interactions have shown that the binary interaction among particles in addition to the repelling Coulomb force exhibits also an attractive part. The tendency of the systems to approach the state with minimum energy by rearranging particles inside has been detected. The measured 63 particles' cluster vibrations are in close agreement with vibrations of a drop with surface tension. This indicates that even a 63 particle cluster already exhibits properties normally associated with the cooperative regime.
NASA Astrophysics Data System (ADS)
Hu, Shengsun; Guo, Chaobo; Wang, Dongpo; Wang, Zhijiang
2016-07-01
The nonuniform distributions of the residual stress were simulated by a 3D finite element model to analyze the elastic-plastic dynamic ultrasonic impact treatment (UIT) process of multiple impacts on the 2024 aluminum alloy. The evolution of the stress during the impact process was discussed. The successive impacts during the UIT process improve the uniformity of the plastic deformation and decrease the maximum compressive residual stress beneath the former impact indentations. The influences of different controlled parameters, including the initial impact velocity, pin diameter, pin tip, device moving, and offset distances, on the residual stress distributions were analyzed. The influences of the controlled parameters on the residual stress distributions are apparent in the offset direction due to the different surface coverage in different directions. The influences can be used to understand the UIT process and to obtain the desired residual stress by optimizing the controlled parameters.
Self-Consistent 3D Modeling of Electron Cloud Dynamics and Beam Response
Furman, Miguel; Furman, M.A.; Celata, C.M.; Kireeff-Covo, M.; Sonnad, K.G.; Vay, J.-L.; Venturini, M.; Cohen, R.; Friedman, A.; Grote, D.; Molvik, A.; Stoltz, P.
2007-04-02
We present recent advances in the modeling of beam electron-cloud dynamics, including surface effects such as secondary electron emission, gas desorption, etc, and volumetric effects such as ionization of residual gas and charge-exchange reactions. Simulations for the HCX facility with the code WARP/POSINST will be described and their validity demonstrated by benchmarks against measurements. The code models a wide range of physical processes and uses a number of novel techniques, including a large-timestep electron mover that smoothly interpolates between direct orbit calculation and guiding-center drift equations, and a new computational technique, based on a Lorentz transformation to a moving frame, that allows the cost of a fully 3D simulation to be reduced to that of a quasi-static approximation.
Dynamic Implicit 3D Adaptive Mesh Refinement for Non-Equilibrium Radiation Diffusion
Philip, Bobby; Wang, Zhen; Berrill, Mark A; Rodriguez Rodriguez, Manuel; Pernice, Michael
2014-01-01
The time dependent non-equilibrium radiation diffusion equations are important for solving the transport of energy through radiation in optically thick regimes and find applications in several fields including astrophysics and inertial confinement fusion. The associated initial boundary value problems that are encountered exhibit a wide range of scales in space and time and are extremely challenging to solve. To efficiently and accurately simulate these systems we describe our research on combining techniques that will also find use more broadly for long term time integration of nonlinear multiphysics systems: implicit time integration for efficient long term time integration of stiff multiphysics systems, local control theory based step size control to minimize the required global number of time steps while controlling accuracy, dynamic 3D adaptive mesh refinement (AMR) to minimize memory and computational costs, Jacobian Free Newton Krylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent linear solver convergence.
The Quantum Dynamics of a Dilute Gas in a 3D BCC Optical Lattice
NASA Astrophysics Data System (ADS)
Reichl, Linda; Boretz, Yingyue
2015-03-01
The classical and quantum dynamics of a dilute gas of rubidium atoms, in a 3D body-centered cubic optical lattice, is studied for a range of polarizations of the laser beams forming the lattice. The relative polarization of the lasers determines the the structure of the potential energy seen by the rubidium atoms. If three pairs of in-phase mutually perpendicular laser beams, with the same wavelength, form the lattice, only a limited range of possible couplings can be realized in the lab. We have determined the band structure of the BCC optical lattice for all theoretically possible couplings, and find that the band structure for lattices realizable in the lab, differs significantly from that expected for a BCC crystal. As coupling is increased, the lattice becomes increasingly chaotic and it becomes possible to produce band structure that has qualitative similarity to a BCC. Welch Foundation
He, W.; Anderson, R.N.
1998-08-25
A method is disclosed for inverting 3-D seismic reflection data obtained from seismic surveys to derive impedance models for a subsurface region, and for inversion of multiple 3-D seismic surveys (i.e., 4-D seismic surveys) of the same subsurface volume, separated in time to allow for dynamic fluid migration, such that small scale structure and regions of fluid and dynamic fluid flow within the subsurface volume being studied can be identified. The method allows for the mapping and quantification of available hydrocarbons within a reservoir and is thus useful for hydrocarbon prospecting and reservoir management. An iterative seismic inversion scheme constrained by actual well log data which uses a time/depth dependent seismic source function is employed to derive impedance models from 3-D and 4-D seismic datasets. The impedance values can be region grown to better isolate the low impedance hydrocarbon bearing regions. Impedance data derived from multiple 3-D seismic surveys of the same volume can be compared to identify regions of dynamic evolution and bypassed pay. Effective Oil Saturation or net oil thickness can also be derived from the impedance data and used for quantitative assessment of prospective drilling targets and reservoir management. 20 figs.
He, Wei; Anderson, Roger N.
1998-01-01
A method is disclosed for inverting 3-D seismic reflection data obtained from seismic surveys to derive impedance models for a subsurface region, and for inversion of multiple 3-D seismic surveys (i.e., 4-D seismic surveys) of the same subsurface volume, separated in time to allow for dynamic fluid migration, such that small scale structure and regions of fluid and dynamic fluid flow within the subsurface volume being studied can be identified. The method allows for the mapping and quantification of available hydrocarbons within a reservoir and is thus useful for hydrocarbon prospecting and reservoir management. An iterative seismic inversion scheme constrained by actual well log data which uses a time/depth dependent seismic source function is employed to derive impedance models from 3-D and 4-D seismic datasets. The impedance values can be region grown to better isolate the low impedance hydrocarbon bearing regions. Impedance data derived from multiple 3-D seismic surveys of the same volume can be compared to identify regions of dynamic evolution and bypassed pay. Effective Oil Saturation or net oil thickness can also be derived from the impedance data and used for quantitative assessment of prospective drilling targets and reservoir management.
3D Reconstruction of Human Laryngeal Dynamics Based on Endoscopic High-Speed Recordings.
Semmler, Marion; Kniesburges, Stefan; Birk, Veronika; Ziethe, Anke; Patel, Rita; Dollinger, Michael
2016-07-01
Standard laryngoscopic imaging techniques provide only limited two-dimensional insights into the vocal fold vibrations not taking the vertical component into account. However, previous experiments have shown a significant vertical component in the vibration of the vocal folds. We present a 3D reconstruction of the entire superior vocal fold surface from 2D high-speed videoendoscopy via stereo triangulation. In a typical camera-laser set-up the structured laser light pattern is projected on the vocal folds and captured at 4000 fps. The measuring device is suitable for in vivo application since the external dimensions of the miniaturized set-up barely exceed the size of a standard rigid laryngoscope. We provide a conservative estimate on the resulting resolution based on the hardware components and point out the possibilities and limitations of the miniaturized camera-laser set-up. In addition to the 3D vocal fold surface, we extended previous approaches with a G2-continuous model of the vocal fold edge. The clinical applicability was successfully established by the reconstruction of visual data acquired from 2D in vivo high-speed recordings of a female and a male subject. We present extracted dynamic parameters like maximum amplitude and velocity in the vertical direction. The additional vertical component reveals deeper insights into the vibratory dynamics of the vocal folds by means of a non-invasive method. The successful miniaturization allows for in vivo application giving access to the most realistic model available and hence enables a comprehensive understanding of the human phonation process. PMID:26829782
Visualization of Computational Fluid Dynamics
NASA Technical Reports Server (NTRS)
Gerald-Yamasaki, Michael; Hultquist, Jeff; Bryson, Steve; Kenwright, David; Lane, David; Walatka, Pamela; Clucas, Jean; Watson, Velvin; Lasinski, T. A. (Technical Monitor)
1995-01-01
Scientific visualization serves the dual purpose of exploration and exposition of the results of numerical simulations of fluid flow. Along with the basic visualization process which transforms source data into images, there are four additional components to a complete visualization system: Source Data Processing, User Interface and Control, Presentation, and Information Management. The requirements imposed by the desired mode of operation (i.e. real-time, interactive, or batch) and the source data have their effect on each of these visualization system components. The special requirements imposed by the wide variety and size of the source data provided by the numerical simulation of fluid flow presents an enormous challenge to the visualization system designer. We describe the visualization system components including specific visualization techniques and how the mode of operation and source data requirements effect the construction of computational fluid dynamics visualization systems.
NASA Astrophysics Data System (ADS)
Kluesner, J. W.; Silver, E. A.; Gibson, J. C.; Bangs, N. L.; McIntosh, K.; von Huene, R.; Orange, D.; Ranero, C. R.
2012-12-01
Offshore southern Costa Rica we have identified 161 potential fluid seepage sites on the shelf and slope regions within an 11 x 55 km strip where no fluid indicators had been reported previously using conventional deep-water mutlibeam bathymetry (100 m grid cell size) and deep towed side scan sonar. Evidence includes large and small pockmarks, mounds, ridges, and slope failure features with localized anomalous high-amplitude backscatter strength. The majority of seepage indicators are associated with shallow sub-bottom reversed polarity bright spots and flat spots imaged within the CRISP 3D seismic grid. Data were collected ~50 km west of Osa Peninsula, Costa Rica onboard the R/V Marcus G. Langseth during the spring of 2011. We obtained EM122 multibeam data using fixed, closely spaced receiver beams and 9-10 times swath overlap, which greatly improved the signal-to-noise ratio and sounding density and allowed for very small grid and mosaic cell sizes (2-10 m). A gas plume in the water column, seen on a 3.5 kHz profile, is located along a fault trace and above surface and subsurface seep indicators. Fluid indicators on the outer shelf occur largely on a dense array of faults, some of which cut through the reflective basement. Seismic flat spots commonly underlie axes of large anticlines on the shelf and slope. Pockmarks are also located at the foot of mid-slope canyons, very near to the upper end of the BSR. These pockmarks appear to be associated with canyon abandonment and folded beds that channel fluids upward, causing hydrate instability. Our findings suggest that significant amounts of methane are venting into ocean and potentially into the atmosphere across the heavily deformed shelf and slope of Costa Rica.
3D Dynamics of Magnetopause Reconnection Using Hall-MHD Global Simulations
NASA Astrophysics Data System (ADS)
Maynard, K.; Germaschewski, K.; Raeder, J.; Bhattacharjee, A.
2011-12-01
Magnetic reconnection at Earth's magnetopause and in the magnetotail is of crucial importance for the dynamics of the global magnetosphere and space weather. Even though the plasma conditions in the magnetosphere are largely in the collisionless regime, most of the existing research using global computational models employ single-fluid magnetohydrodynamics (MHD) with artificial resistivity. Studies of reconnection in simplified, two-dimensional geometries have established that two-fluid and kinetic effects can dramatically alter dynamics and reconnection rates when compared with single-fluid models. These enhanced models also introduce particular signatures, for example a quadrupolar out-of-plane magnetic field component that has already been observed in space by satellite measurements. However, results from simplified geometries cannot be translated directly to the dynamics of three-dimensional magnetospheric reconnection. For instance, magnetic flux originating from the solar wind and arriving at the magnetopause can either reconnect or be advected around the magnetosphere. In this study, we use a new version of the OpenGGCM code that incorporates the Hall term in a Generalized Ohm's Law to study magnetopause reconnection under synthetic solar wind conditions and investigate how reconnection rates and dynamics of flux transfer events depend on the strength of the Hall term. The OpenGGCM, a global model of Earth's magnetosphere, has recently been ported to exploit modern computing architectures like the Cell processor and SIMD capabilities of conventional processors using an automatic code generator. These enhancements provide us with the performance needed to include the computationally expensive Hall physics.
Transforming 2d Cadastral Data Into a Dynamic Smart 3d Model
NASA Astrophysics Data System (ADS)
Tsiliakou, E.; Labropoulos, T.; Dimopoulou, E.
2013-08-01
3D property registration has become an imperative need in order to optimally reflect all complex cases of the multilayer reality of property rights and restrictions, revealing their vertical component. This paper refers to the potentials and multiple applications of 3D cadastral systems and explores the current state-of-the art, especially the available software with which 3D visualization can be achieved. Within this context, the Hellenic Cadastre's current state is investigated, in particular its data modeling frame. Presenting the methodologies and specifications addressing the registration of 3D properties, the operating cadastral system's shortcomings and merits are pointed out. Nonetheless, current technological advances as well as the availability of sophisticated software packages (proprietary or open source) call for 3D modeling. In order to register and visualize the complex reality in 3D, Esri's CityEngine modeling software has been used, which is specialized in the generation of 3D urban environments, transforming 2D GIS Data into Smart 3D City Models. The application of the 3D model concerns the Campus of the National Technical University of Athens, in which a complex ownership status is established along with approved special zoning regulations. The 3D model was built using different parameters based on input data, derived from cadastral and urban planning datasets, as well as legal documents and architectural plans. The process resulted in a final 3D model, optimally describing the cadastral situation and built environment and proved to be a good practice example of 3D visualization.
Fundamentals of Geophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
McWilliams, James C.
2006-07-01
Earth's atmosphere and oceans exhibit complex patterns of fluid motion over a vast range of space and time scales. These patterns combine to establish the climate in response to solar radiation that is inhomogeneously absorbed by the materials comprising air, water, and land. Spontaneous, energetic variability arises from instabilities in the planetary-scale circulations, appearing in many different forms such as waves, jets, vortices, boundary layers, and turbulence. Geophysical fluid dynamics (GFD) is the science of all these types of fluid motion. This textbook is a concise and accessible introduction to GFD for intermediate to advanced students of the physics, chemistry, and/or biology of Earth's fluid environment. The book was developed from the author's many years of teaching a first-year graduate course at the University of California, Los Angeles. Readers are expected to be familiar with physics and mathematics at the level of general dynamics (mechanics) and partial differential equations. Covers the essential GFD required for atmospheric science and oceanography courses Mathematically rigorous, concise coverage of basic theory and applications to both oceans and atmospheres Author is a world expert; this book is based on the course he has taught for many years Exercises are included, with solutions available to instructors from solutions@cambridge.org
Visualization of unsteady computational fluid dynamics
NASA Technical Reports Server (NTRS)
Haimes, Robert
1995-01-01
The current computing environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array) provide the required computation bandwidth for CFD calculations of transient problems. Work is in progress on a set of software tools designed specifically to address visualizing 3D unsteady CFD results in these super-computer-like environments. The visualization is concurrently executed with the CFD solver. The parallel version of Visual3, pV3 required splitting up the unsteady visualization task to allow execution across a network of workstation(s) and compute servers. In this computing model, the network is almost always the bottleneck so much of the effort involved techniques to reduce the size of the data transferred between machines.
Visualization of unsteady computational fluid dynamics
NASA Astrophysics Data System (ADS)
Haimes, Robert
1995-10-01
The current computing environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array) provide the required computation bandwidth for CFD calculations of transient problems. Work is in progress on a set of software tools designed specifically to address visualizing 3D unsteady CFD results in these super-computer-like environments. The visualization is concurrently executed with the CFD solver. The parallel version of Visual3, pV3 required splitting up the unsteady visualization task to allow execution across a network of workstation(s) and compute servers. In this computing model, the network is almost always the bottleneck so much of the effort involved techniques to reduce the size of the data transferred between machines.
Yu, Yingjie; Zhang, Qi; Buscaglia, Jonathan; Chang, Chung-Chueh; Liu, Ying; Yang, Zhenhua; Guo, Yichen; Wang, Yantian; Levon, Kalle; Rafailovich, Miriam
2016-07-21
In this study, a sensitive, yet robust, biosensing system with real-time electrochemical readout was developed. The biosensor system was applied to the detection of carcinoembryonic antigen (CEA), which is a common marker for many cancers such as pancreatic, breast, and colon cancer. Real time detection of CEA during a medical procedure can be used to make critical decisions regarding further surgical intervention. CEA was templated on gold surface (RMS roughness ∼3-4 nm) coated with a hydrophilic self-assembled monolayer (SAM) on the working electrode of an open circuit potentiometric network. The subsequent removal of template CEA makes the biosensor capable of CEA detection based on its specific structure and conformation. The molecular imprinting (MI) biosensor was further calibrated using the potentiometric responses in solutions with known CEA concentrations and a detection limit of 0.5 ng ml(-1) was achieved. Potentiometric sensing was then applied to pancreatic cyst fluid samples obtained from 18 patients when the cyst fluid was also evaluated using ELISA in a certified pathology laboratory. Excellent agreement was obtained between the quantitation of CEA obtained by both the ELISA and MI biosensor detection for CEA. A 3-D MI model, using the natural rms roughness of PVD gold layers, is presented to explain the high degree of sensitivity and linearity observed in those experiments. PMID:27193921
Using 3D dynamic models to reproduce X-ray properties of colliding wind binaries
NASA Astrophysics Data System (ADS)
Russell, Christopher Michael Post
Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The two massive stars contained in these systems have powerful radiatively driven stellar winds, and the conversion of their kinetic energy to heat (up to 108 K) at the wind-wind collision region generates hard thermal X-rays (up to 10 keV). Rich data sets exist of several multi-year-period systems, as well as key observations of shorter period systems, and detailed models are required to disentangle the phase-locked emission and absorption processes in these systems. To interpret these X-ray light curves and spectra, this dissertation models the wind-wind interaction of CWBs using 3D smoothed particle hydrodynamics (SPH), and solves the 3D formal solution of radiative transfer to synthesize the model X-ray properties, allowing direct comparison with the colliding-wind X-ray spectra observed by, e.g., RXTE and XMM. The multi-year-period, highly eccentric CWBs we examine are eta Carinae and WR140. For the commonly inferred primary mass loss rate of ˜10 -3 Msun/yr, eta Carinae's 3D model reproduces quite well the 2-10 keV RXTE light curve, hardness ratio, and dynamic spectra in absolute units. This agreement includes the ˜3 month X-ray minimum associated with the 1998.0 and 2003.5 periastron passages, which we find to occur as the primary wind encroaches into the secondary wind's acceleration region. This modeling provides further evidence that the observer is mainly viewing the system through the secondary's shock cone, and suggests that periastron occurs ~1 month after the onset of the X-ray minimum. The model RXTE observables of WR140 match the data well in absolute units, although the decrease in model X-rays around periastron is less than observed. There is very good agreement between the observed XMM spectrum taken on the rise before periastron and the model. We also model two short-period CWBs, HD150136, which has a wind-star collision, and delta Orionis A, the closest eclipsing
Agour, Mostafa; Falldorf, Claas; Bergmann, Ralf B
2016-06-27
We present a new method for the generation of a dynamic wave field with high space bandwidth product (SBP). The dynamic wave field is generated from several wave fields diffracted by a display which comprises multiple spatial light modulators (SLMs) each having a comparably low SBP. In contrast to similar approaches in stereoscopy, we describe how the independently generated wave fields can be coherently superposed. A major benefit of the scheme is that the display system may be extended to provide an even larger display. A compact experimental configuration which is composed of four phase-only SLMs to realize the coherent combination of independent wave fields is presented. Effects of important technical parameters of the display system on the wave field generated across the observation plane are investigated. These effects include, e.g., the tilt of the individual SLM and the gap between the active areas of multiple SLMs. As an example of application, holographic reconstruction of a 3D object with parallax effects is demonstrated. PMID:27410593
A Real-time, 3D Musculoskeletal Model for Dynamic Simulation of Arm Movements
Chadwick, Edward K.; Blana, Dimitra; van den Bogert, Antonie J.; Kirsch, Robert F.
2010-01-01
Neuroprostheses can be used to restore movement of the upper limb in individuals with high-level spinal cord injury. Development and evaluation of command and control schemes for such devices typically requires real-time, “patient-in-the-loop” experimentation. A real-time, three-dimensional, musculoskeletal model of the upper limb has been developed for use in a simulation environment to allow such testing to be carried out non-invasively. The model provides real-time feedback of human arm dynamics that can be displayed to the user in a virtual reality environment. The model has a three degree-of-freedom gleno-humeral joint as well as elbow flexion/extension and pronation/supination, and contains 22 muscles of the shoulder and elbow divided into multiple elements. The model is able to run in real time on modest desktop hardware and demonstrates that a large-scale, 3D model can be made to run in real time. This is a prerequisite for a real-time, whole arm model that will form part of a dynamic arm simulator for use in the development, testing and user training of neural prosthesis systems. PMID:19272926
A digital holography set-up for 3D vortex flow dynamics
NASA Astrophysics Data System (ADS)
Lebon, Benoît; Perret, Gaële; Coëtmellec, Sébastien; Godard, Gilles; Gréhan, Gérard; Lebrun, Denis; Brossard, Jérôme
2016-06-01
In the present paper, a digital in-line holography (DIH) set-up, with a converging beam, is used to take three-dimensional (3D) velocity measurements of vortices. The vortices are formed periodically at the edges of a submerged horizontal plate submitted to regular waves. They take the form of vortex filaments that extend from side to side of the channel. They undergo strongly three-dimensional instability mechanisms that remain very complicated to characterize experimentally. The experiments are performed in a 10 × 0.3 × 0.3 m3 wave flume. The DIH set-up is performed using a modulated laser diode emitting at the wavelength of 640 nm and a lensless CCD camera. The beam crosses the channel side to side. To reveal the flow dynamics, 30-μm hydrogen bubbles are generated at the edge of the plate to serve as tracers. Their locations are recorded on the holograms multiple times to access the dynamics of the flow. This method leads to an accuracy in the order of 100 μm on the axial location. Those measurements have been validated with stereo-PIV measurements. A very good agreement is found on time-averaged velocity fields between the two techniques.
NASA Astrophysics Data System (ADS)
Javidi, Bahram; Yeom, Seokwon; Moon, Inkyu; Daneshpanah, Mehdi
2006-05-01
In this paper, we present an overview of three-dimensional (3D) optical imaging techniques for real-time automated sensing, visualization, and recognition of dynamic biological microorganisms. Real time sensing and 3D reconstruction of the dynamic biological microscopic objects can be performed by single-exposure on-line (SEOL) digital holographic microscopy. A coherent 3D microscope-based interferometer is constructed to record digital holograms of dynamic micro biological events. Complex amplitude 3D images of the biological microorganisms are computationally reconstructed at different depths by digital signal processing. Bayesian segmentation algorithms are applied to identify regions of interest for further processing. A number of pattern recognition approaches are addressed to identify and recognize the microorganisms. One uses 3D morphology of the microorganisms by analyzing 3D geometrical shapes which is composed of magnitude and phase. Segmentation, feature extraction, graph matching, feature selection, and training and decision rules are used to recognize the biological microorganisms. In a different approach, 3D technique is used that are tolerant to the varying shapes of the non-rigid biological microorganisms. After segmentation, a number of sampling patches are arbitrarily extracted from the complex amplitudes of the reconstructed 3D biological microorganism. These patches are processed using a number of cost functions and statistical inference theory for the equality of means and equality of variances between the sampling segments. Also, we discuss the possibility of employing computational integral imaging for 3D sensing, visualization, and recognition of biological microorganisms illuminated under incoherent light. Experimental results with several biological microorganisms are presented to illustrate detection, segmentation, and identification of micro biological events.
Chisholm, Kelsey; Miles, Devin; Rankine, Leith; Oldham, Mark
2015-01-01
Purpose: In optical-CT, the use of a refractively matched polyurethane solid-tank in place of a fluid bath has the potential to greatly increase practical convenience, reduce cost, and possibly improve the efficacy of flood corrections. This work investigates the feasibility of solid-tank optical-CT imaging for 3D dosimetry through computer simulation. Methods: A matlab ray-tracing simulation platform, ScanSim, was used to model a parallel-source telecentric optical-CT imaging system through a polyurethane solid-tank containing a central cylindrical hollow into which PRESAGE radiochromic dosimeters can be placed. A small amount of fluid fills the 1–5 mm gap between the dosimeter and the walls of the tank. The use of the solid-tank reduces the required amount of fluid by approximately 97%. To characterize the efficacy of solid-tank, optical-CT scanning simulations investigated sensitivity to refractive index (RI) mismatches between dosimeter, solid-tank, and fluid, for a variety of dosimeter (RI = 1.5–1.47) and fluid (RI = 1.55–1.0) combinations. Efficacy was evaluated through the usable radius (ru) metric, defined as the fraction of the radius of the dosimeter where measured dose is predicted to be within 2% of the ground truth entered into the simulation. Additional simulations examined the effect of increasing gap size (1–5 mm) between the dosimeter and solid-tank well. The effects of changing the lens tolerance (0.5°–5.0°) were also investigated. Results: As the RI mismatch between the dosimeter and solid-tank increased from 0 to 0.02, the usable radius decreased from 97.6% to 50.2%. The optimal fluid RI decreased nonlinearly from 1.5 to 1.34 as the mismatch increased and was up to 9% lower than the tank. Media mismatches between the dosimeter and solid-tank also exacerbate the effects of changing the gap size, with no easily quantifiable relationship with usable radius. Generally, the optimal fluid RI value increases as gap size increases and is
Chisholm, Kelsey; Miles, Devin; Rankine, Leith; Oldham, Mark
2015-05-15
Purpose: In optical-CT, the use of a refractively matched polyurethane solid-tank in place of a fluid bath has the potential to greatly increase practical convenience, reduce cost, and possibly improve the efficacy of flood corrections. This work investigates the feasibility of solid-tank optical-CT imaging for 3D dosimetry through computer simulation. Methods: A MATLAB ray-tracing simulation platform, ScanSim, was used to model a parallel-source telecentric optical-CT imaging system through a polyurethane solid-tank containing a central cylindrical hollow into which PRESAGE radiochromic dosimeters can be placed. A small amount of fluid fills the 1–5 mm gap between the dosimeter and the walls of the tank. The use of the solid-tank reduces the required amount of fluid by approximately 97%. To characterize the efficacy of solid-tank, optical-CT scanning simulations investigated sensitivity to refractive index (RI) mismatches between dosimeter, solid-tank, and fluid, for a variety of dosimeter (RI = 1.5–1.47) and fluid (RI = 1.55–1.0) combinations. Efficacy was evaluated through the usable radius (r{sub u}) metric, defined as the fraction of the radius of the dosimeter where measured dose is predicted to be within 2% of the ground truth entered into the simulation. Additional simulations examined the effect of increasing gap size (1–5 mm) between the dosimeter and solid-tank well. The effects of changing the lens tolerance (0.5°–5.0°) were also investigated. Results: As the RI mismatch between the dosimeter and solid-tank increased from 0 to 0.02, the usable radius decreased from 97.6% to 50.2%. The optimal fluid RI decreased nonlinearly from 1.5 to 1.34 as the mismatch increased and was up to 9% lower than the tank. Media mismatches between the dosimeter and solid-tank also exacerbate the effects of changing the gap size, with no easily quantifiable relationship with usable radius. Generally, the optimal fluid RI value increases as gap size increases and
Lectures on Geophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Samelson, Roger M.
The fluid kaleidoscope of the Earth's ocean and atmosphere churns and sparkles with jets, gyres, eddies, waves, streams, and cyclones. These vast circulations, essential elements of the physical environment that support human life, are given a special character by the Earth's rotation and by their confinement to a shallow surficial layer, thin relative to the solid Earth in roughly the same proportion as an apple skin is to an apple. Geophysical fluid dynamics exploits this special character to develop a unified theoretical approach to the physics of the ocean and atmosphere.With Lectures on Geophysical Fluid Dynamics, Rick Salmon has added an insightful and provocative volume to the handful of authoritative texts currently available on the subject. The book is intended for first-year graduate students, but advanced students and researchers also will find it useful. It is divided into seven chapters, the first four of these adapted from course lectures. The book is well written and presents a fresh and stimulating perspective that complements existing texts. It would serve equally well either as the main text for a core graduate curriculum or as a supplementary resource for students and teachers seeking new approaches to both classical and contemporary problems. A lively set of footnotes contains many references to very recent work. The printing is attractive, the binding is of high quality, and typographical errors are few.
NASA Astrophysics Data System (ADS)
Campbell, Nicholas
2014-11-01
The success of dielectric barrier discharge (DBD) plasma actuators as flow control devices in transducing electrical energy directly into near instantaneous fluid motion has been limited due to momentum loss near the wall. To increase the feasibility of these devices, they have been used to drive a channel flow, creating a jet under quiescent conditions. Electrostatic Fluid Accelerators (EFA) have also been shown to drive internal gas flows. The present work draws on the success of the DBD driven plasma channels, while exploring a new electrode configuration that stems from EFA designs, in order to actuate more of the bulk fluid. Major parameters, applied voltage and operating frequency as well as electrode gap and choice of electrode (material, shape, size); were experimentally investigated using Particle Image Velocimetry to obtain time averaged, 2D velocity fields. Results indicate significant variation of performance with these parameters and suggest that in comparison to surface DBD actuators an order magnitude improvement in efficiency is possible. Furthermore, the qualitative aspect of an electro-fluid dynamic jet shows greater versatility in application for use as both boundary layer flow control and driving internal gas flows.
NASA Astrophysics Data System (ADS)
Duru, Kenneth; Dunham, Eric M.
2016-01-01
Dynamic propagation of shear ruptures on a frictional interface in an elastic solid is a useful idealization of natural earthquakes. The conditions relating discontinuities in particle velocities across fault zones and tractions acting on the fault are often expressed as nonlinear friction laws. The corresponding initial boundary value problems are both numerically and computationally challenging. In addition, seismic waves generated by earthquake ruptures must be propagated for many wavelengths away from the fault. Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods. We present a high order accurate finite difference method for: a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration; b) dynamic propagation of earthquake ruptures along nonplanar faults; and c) accurate propagation of seismic waves in heterogeneous media with free surface topography. We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts (SBP) finite difference operators in space. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. The finite difference stencils used in this paper are sixth order accurate in the interior and third order accurate close to the boundaries. However, the method is applicable to any spatial operator with a diagonal norm satisfying the SBP property. Time stepping is performed with a 4th order accurate explicit low storage Runge-Kutta scheme, thus yielding a globally fourth order accurate method in both space and time. We show numerical simulations on band limited self-similar fractal faults revealing the complexity of rupture dynamics
NASA Astrophysics Data System (ADS)
Duru, K.; Dunham, E. M.; Bydlon, S. A.; Radhakrishnan, H.
2014-12-01
Dynamic propagation of shear ruptures on a frictional interface is a useful idealization of a natural earthquake.The conditions relating slip rate and fault shear strength are often expressed as nonlinear friction laws.The corresponding initial boundary value problems are both numerically and computationally challenging.In addition, seismic waves generated by earthquake ruptures must be propagated, far away from fault zones, to seismic stations and remote areas.Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods.We present a numerical method for:a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration;b) dynamic propagation of earthquake ruptures along rough faults; c) accurate propagation of seismic waves in heterogeneous media with free surface topography.We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts finite differences in space. The finite difference stencils are 6th order accurate in the interior and 3rd order accurate close to the boundaries. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. Time stepping is performed with a 4th order accurate explicit low storage Runge-Kutta scheme. We have performed extensive numerical experiments using a slip-weakening friction law on non-planar faults, including recent SCEC benchmark problems. We also show simulations on fractal faults revealing the complexity of rupture dynamics on rough faults. We are presently extending our method to rate-and-state friction laws and off-fault plasticity.
NASA Astrophysics Data System (ADS)
Bartzke, Gerhard; Rogers, Benedict D.; Fourtakas, Georgios; Mokos, Athanasios; Huhn, Katrin
2016-04-01
The processes that cause the creation of a variety of sediment morphological features, e.g. laminated beds, ripples, or dunes, are based on the initial motion of individual sediment grains. However, with experimental techniques it is difficult to measure the flow characteristics, i.e., the velocity of the pore water flow in sediments, at a sufficient resolution and in a non-intrusive way. As a result, the role of fluid infiltration at the surface and in the interior affecting the initiation of motion of a sediment bed is not yet fully understood. Consequently, there is a strong need for numerical models, since these are capable of quantifying fluid driven sediment transport processes of complex sediment beds composed of irregular shapes. The numerical method Smoothed Particle Hydrodynamics (SPH) satisfies this need. As a meshless and Lagrangian technique, SPH is ideally suited to simulating flows in sediment beds composed of various grain shapes, but also flow around single grains at a high temporal and spatial resolution. The solver chosen is DualSPHysics (www.dual.sphysics.org) since this is validated for a range of flow conditions. For the present investigation a 3-D numerical flume model was generated using SPH with a length of 4.0 cm, a width of 0.05 cm and a height of 0.2 cm where mobile sediment particles were deposited in a recess. An experimental setup was designed to test sediment configurations composed of irregular grain shapes (grain diameter, D50=1000 μm). Each bed consisted of 3500 mobile objects. After the bed generation process, the entire domain was flooded with 18 million fluid particles. To drive the flow, an oscillating motion perpendicular to the bed was applied to the fluid, reaching a peak value of 0.3 cm/s, simulating 4 seconds of real time. The model results showed that flow speeds decreased logarithmically from the top of the domain towards the surface of the beds, indicating a fully developed boundary layer. Analysis of the fluid
NASA Astrophysics Data System (ADS)
Dégi, Júlia; Török, Kálmán; Németh, Bianca; Rhede, Dieter; Takács, Ágnes; Habler, Gerlinde; Abart, Rainer
2016-04-01
Complex kelyphitic rims after garnet were studied in a lower crustal mafic granulite xenolith from the Bakony-Balaton Highland Volcanic Field, W-Hungary. The xenolith is dominated by a garnet granulite mineral assemblage equilibrated at 900 ° C, 1.4 GPa within the overthickened orogenic root of the Alps. Garnet breakdown was initiated during the extension of the Pannonian Basin and remained continuous until the xenoliths reached the surface. This resulted in the formation of various microstructural domains within the kelyphitic rims which were distinguished to three main types: fine-grained symplectites, recrystallized symplectites and cyrstallized melts. Fine-grained symplectites are the products of isochemical breakdown of garnet to pure anorthite, Al-orthopyroxene and hercynitic spinel. Nanoscale topography built up by curved chains of humps on the garnet surface showing regular spatial distribution is observed at the reaction front in 3D reconstructions. These patterns follow the contours of anorthites nucleating at the reaction front. This suggests that diffusion-controlled material transfer in solid state together with surface energy minimization determines symplectite microstructure. The latter leads to continuous isochemical coarsening getting further from the reaction front. Slight increase of Na and Ti-content in fine-grained symplectites is associated with sudden changes in 3D microstructure reflecting the effect of aquaeous fluids infiltrating to the reaction rim. A similar aquaeous fluid may have induced the formation of nearly isochemical melting and in situ recrystallization of the symplectites to form recrystallized symplectites. Some of these remained pristine, while some experienced ductile deformation and re-melting due to a reaction with an external melt in the lower crust. Well-crystallized melt pockets consisting of complexly zoned pyroxene, spinel and plagioclase grains were formed around 1000 ° C in this stage. Following this, interaction
Toward a 3D dynamic model of a faulty duplex ball bearing
NASA Astrophysics Data System (ADS)
Kogan, Gideon; Klein, Renata; Kushnirsky, Alex; Bortman, Jacob
2015-03-01
Bearings are vital components for safe and proper operation of machinery. Increasing efficiency of bearing diagnostics usually requires training of health and usage monitoring systems via expensive and time-consuming ground calibration tests. The main goal of this research, therefore, is to improve bearing dynamics modeling tools in order to reduce the time and budget needed to implement the health and usage monitoring approach. The proposed three-dimensional ball bearing dynamic model is based on the classic dynamic and kinematic equations. Interactions between the bodies are simulated using non-linear springs combined with dampers described by Hertz-type contact relation. The force friction is simulated using the hyperbolic-tangent function. The model allows simulation of a wide range of mechanical faults. It is validated by comparison to known bearing behavior and to experimental results. The model results are verified by demonstrating numerical convergence. The model results for the two cases of single and duplex angular ball bearings with axial deformation in the outer ring are presented. The qualitative investigation provides insight into bearing dynamics, the sensitivity study generalizes the qualitative findings for similar cases, and the comparison to the test results validates model reliability. The article demonstrates the variety of the cases that the 3D bearing model can simulate and the findings to which it may lead. The research allowed the identification of new patterns generated by single and duplex bearings with axially deformed outer race. It also enlightened the difference between single and duplex bearing manifestation. In the current research the dynamic model enabled better understanding of the physical behavior of the faulted bearings. Therefore, it is expected that the modeling approach has the potential to simplify and improve the development process of diagnostic algorithms. • A deformed outer race of a single axially loaded bearing is
Fluid dynamics and vibration of tube banks in fluid flow
Zukauskas, A.; Ulinskas, R.; Katinas, V.
1988-01-01
This work presents results derived in fluid dynamics, hydraulic drag and flow-induced vibrations within transverse and yawed tube banks. The studies encompass banks of smooth, rough and finned tubes at Reynolds numbers from 1 to 2x10/sup 6/. Highlighted in the text are fluid dynamic parameters of tube banks measured at inter-tube spaces and tube surfaces.
NASA Astrophysics Data System (ADS)
Barnhart, Donald H.; Chan, Victor S. S.; Halliwell, Neil A.; Coupland, Jeremy M.
1999-10-01
This paper reports the use of a new holographic measurement system in the study of 3D surface displacements. Although equally applicable to fluid and solid mechanics, the aim of this report is to demonstrate the system's use in quantitative surface displacement measurements with a classical cantilever experiment, using a continuous-wave diode-pumped YAG laser system. The reported results exhibit an accuracy corresponding to other interferometric systems, but with a much larger displacement range. The measurement system employs a novel optical image shifting method to eliminate the problem of directional ambiguity. In addition, the reported system uses 3D complex correlation rather than 2D real correlation, thereby offering a direct method for measuring 3D displacement in 3D space. FInally, with the novel use of an optical fiber to probe the recorded holographic image space, it is found to be a simple matter to directly obtain 3D displacement measurements at precisely known surface locations.
NASA Technical Reports Server (NTRS)
Folta, David; Bauer, Frank H. (Technical Monitor)
2001-01-01
The autonomous formation flying control algorithm developed by the Goddard Space Flight Center (GSFC) for the New Millennium Program (NMP) Earth Observing-1 (EO-1) mission is investigated for applicability to libration point orbit formations. In the EO-1 formation-flying algorithm, control is accomplished via linearization about a reference transfer orbit with a state transition matrix (STM) computed from state inputs. The effect of libration point orbit dynamics on this algorithm architecture is explored via computation of STMs using the flight proven code, a monodromy matrix developed from a N-body model of a libration orbit, and a standard STM developed from the gravitational and coriolis effects as measured at the libration point. A comparison of formation flying Delta-Vs calculated from these methods is made to a standard linear quadratic regulator (LQR) method. The universal 3-D approach is optimal in the sense that it can be accommodated as an open-loop or closed-loop control using only state information.