Coupled 1D-3D hydrodynamic modelling, with application to the Pearl River Delta
NASA Astrophysics Data System (ADS)
Twigt, Daniel J.; de Goede, Erik D.; Zijl, Firmijn; Schwanenberg, Dirk; Chiu, Alex Y. W.
2009-12-01
Within the hydrodynamic modelling community, it is common practice to apply different modelling systems for coastal waters and river systems. Whereas for coastal waters 3D finite difference or finite element grids are commonly used, river systems are generally modelled using 1D networks. Each of these systems is tailored towards specific applications. Three-dimensional coastal water models are designed to model the horizontal and vertical variability in coastal waters and are less well suited for representing the complex geometry and cross-sectional areas of river networks. On the other hand, 1D river network models are designed to accurately represent complex river network geometries and complex structures like weirs, barrages and dams. A disadvantage, however, is that they are unable to resolve complex spatial flow variability. In real life, however, coastal oceans and rivers interact. In deltaic estuaries, both tidal intrusion of seawater into the upstream river network and river discharge into open waters play a role. This is frequently approached by modelling the systems independently, with off-line coupling of the lateral boundary forcing. This implies that the river and the coastal model run sequentially, providing lateral discharge (1D) and water level (3D) forcing to each other without the possibility of direct feedback or interaction between these processes. An additional disadvantage is that due to the time aggregation usually applied to exchanged quantities, mass conservation is difficult to ensure. In this paper, we propose an approach that couples a 3D hydrodynamic modelling system for coastal waters (Delft3D) with a 1D modelling system for river hydraulics (SOBEK) online. This implies that contrary to off-line coupling, the hydrodynamic quantities are exchanged between the 1D and 3D domains during runtime to resolve the real-time exchange and interaction between the coastal waters and river network. This allows for accurate and mass conserving
Coupling Magnetic Fields and ALE Hydrodynamics for 3D Simulations of MFCG's
White, D; Rieben, R; Wallin, B
2006-09-20
We review the development of a full 3D multiphysics code for the simulation of explosively driven Magnetic Flux Compression Generators (MFCG) and related pulse power devices. In a typical MFCG the device is seeded with an initial electric current and the device is then detonated. The detonation compresses the magnetic field and amplifies the current. This is a multiphysics problem in that detonation kinetics, electromagnetic diffusion and induction, material deformation, and thermal effects are all important. This is a tightly coupled problem in that the different physical quantities have comparable spatial and temporal variation, and hence should be solved simultaneously on the same computational mesh.
NASA Astrophysics Data System (ADS)
Granato, Gian Luigi; Ragone-Figueroa, Cinthia; Domínguez-Tenreiro, Rosa; Obreja, Aura; Borgani, Stefano; De Lucia, Gabriella; Murante, Giuseppe
2015-06-01
We compute and study the infrared and sub-mm properties of high-redshift (z ≳ 1) simulated clusters and protoclusters. The results of a large set of hydrodynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z = 0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and spectral energy distributions (SEDs) of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far-IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation (SF) rather than accretion on to super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode.
3D Hydrodynamic Simulations of Relativistic Jets
NASA Astrophysics Data System (ADS)
Hughes, P. A.; Miller, M. A.; Duncan, G. C.; Swift, C. M.
1998-12-01
We present the results of validation runs and the first extragalactic jet simulations performed with a 3D relativistic numerical hydrodynamic code employing a solver of the RHLLE type and using adaptive mesh refinement (AMR; Duncan & Hughes, 1994, Ap. J., 436, L119). Test problems include the shock tube, blast wave and spherical shock reflection (implosion). Trials with the code show that as a consequence of AMR it is viable to perform exploratory runs on workstation class machines (with no more than 128Mb of memory) prior to production runs. In the former case we achieve a resolution not much less than that normally regarded as the minimum needed to capture the essential physics of a problem, which means that such runs can provide valuable guidance allowing the optimum use of supercomputer resources. We present initial results from a program to explore the 3D stability properties of flows previously studied using a 2D axisymmetric code, and our first attempt to explore the structure and morphology of a relativistic jet encountering an ambient density gradient that mimics an ambient inhomogeneity or cloud.
Explicit 3-D Hydrodynamic FEM Program
Energy Science and Technology Software Center (ESTSC)
2000-11-07
DYNA3D is a nonlinear explicit finite element code for analyzing 3-D structures and solid continuum. The code is vectorized and available on several computer platforms. The element library includes continuum, shell, beam, truss and spring/damper elements to allow maximum flexibility in modeling physical problems. Many materials are available to represent a wide range of material behavior, including elasticity, plasticity, composites, thermal effects and rate dependence. In addition, DYNA3D has a sophisticated contact interface capability, includingmore » frictional sliding, single surface contact and automatic contact generation.« less
Explicit 3-D Hydrodynamic FEM Program
2000-11-07
DYNA3D is a nonlinear explicit finite element code for analyzing 3-D structures and solid continuum. The code is vectorized and available on several computer platforms. The element library includes continuum, shell, beam, truss and spring/damper elements to allow maximum flexibility in modeling physical problems. Many materials are available to represent a wide range of material behavior, including elasticity, plasticity, composites, thermal effects and rate dependence. In addition, DYNA3D has a sophisticated contact interface capability, including frictional sliding, single surface contact and automatic contact generation.
FARGO3D: Hydrodynamics/magnetohydrodynamics code
NASA Astrophysics Data System (ADS)
Benítez Llambay, Pablo; Masset, Frédéric
2015-09-01
A successor of FARGO (ascl:1102.017), FARGO3D is a versatile HD/MHD code that runs on clusters of CPUs or GPUs, with special emphasis on protoplanetary disks. FARGO3D offers Cartesian, cylindrical or spherical geometry; 1-, 2- or 3-dimensional calculations; and orbital advection (aka FARGO) for HD and MHD calculations. As in FARGO, a simple Runge-Kutta N-body solver may be used to describe the orbital evolution of embedded point-like objects. There is no need to know CUDA; users can develop new functions in C and have them translated to CUDA automatically to run on GPUs.
3D hydrodynamic focusing microfluidics for emerging sensing technologies.
Daniele, Michael A; Boyd, Darryl A; Mott, David R; Ligler, Frances S
2015-05-15
While the physics behind laminar flows has been studied for 200 years, understanding of how to use parallel flows to augment the capabilities of microfluidic systems has been a subject of study primarily over the last decade. The use of one flow to focus another within a microfluidic channel has graduated from a two-dimensional to a three-dimensional process and the design principles are only now becoming established. This review explores the underlying principles for hydrodynamic focusing in three dimensions (3D) using miscible fluids and the application of these principles for creation of biosensors, separation of cells and particles for sample manipulation, and fabrication of materials that could be used for biosensors. Where sufficient information is available, the practicality of devices implementing fluid flows directed in 3D is evaluated and the advantages and limitations of 3D hydrodynamic focusing for the particular application are highlighted. PMID:25041926
3-D HYDRODYNAMIC MODELING IN A GEOSPATIAL FRAMEWORK
Bollinger, J; Alfred Garrett, A; Larry Koffman, L; David Hayes, D
2006-08-24
3-D hydrodynamic models are used by the Savannah River National Laboratory (SRNL) to simulate the transport of thermal and radionuclide discharges in coastal estuary systems. Development of such models requires accurate bathymetry, coastline, and boundary condition data in conjunction with the ability to rapidly discretize model domains and interpolate the required geospatial data onto the domain. To facilitate rapid and accurate hydrodynamic model development, SRNL has developed a pre- and post-processor application in a geospatial framework to automate the creation of models using existing data. This automated capability allows development of very detailed models to maximize exploitation of available surface water radionuclide sample data and thermal imagery.
Regularity criterion for the 3D Hall-magneto-hydrodynamics
NASA Astrophysics Data System (ADS)
Dai, Mimi
2016-07-01
This paper studies the regularity problem for the 3D incompressible resistive viscous Hall-magneto-hydrodynamic (Hall-MHD) system. The Kolmogorov 41 phenomenological theory of turbulence [14] predicts that there exists a critical wavenumber above which the high frequency part is dominated by the dissipation term in the fluid equation. Inspired by this idea, we apply an approach of splitting the wavenumber combined with an estimate of the energy flux to obtain a new regularity criterion. The regularity condition presented here is weaker than conditions in the existing criteria (Prodi-Serrin type criteria) for the 3D Hall-MHD system.
3D hydrodynamical and radiative transfer modeling of η Carinae's colliding winds
NASA Astrophysics Data System (ADS)
Madura, T. I.; Clementel, N.; Gull, T. R.; Kruip, C. J. H.; Paardekooper, J.-P.; Icke, V.
We present results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system η Carinae. We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We investigate several computational domain sizes and Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing. While we initially focus on η Car, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty `pinwheel' (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where various observed time-variable emission and absorption lines form in these unique objects.
3D Hydrodynamical and Radiative Transfer Modeling of Eta Carinae's Colliding Winds
NASA Astrophysics Data System (ADS)
Madura, Thomas Ignatius; Clementel, Nicola; Gull, Theodore R.; Kruip, Chael J. H.; Paardekooper, Jan-Pieter; Icke, Vincent
2015-08-01
We present the results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system Eta Carinae (Clementel, Madura, et al. 2014, MNRAS, 443, 2475 and Clementel, Madura, et al. 2015, MNRAS, 447, 2445). We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to 3D smoothed particle hydrodynamics simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium in 3D. We investigate several computational domain sizes and Luminous Blue Variable primary-star mass-loss rates. We show how the SimpleX simulations can be used to generate synthetic spectral data cubes for comparison to data obtained with the Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph as part of a multi-cycle program to map changes in Eta Carinae's spatially extended interacting wind structures across one binary cycle. Comparison of the HST observations to the SimpleX models can help lead to more accurate constraints on the orbital, stellar, and wind parameters of the Eta Carinae system, such as the LBV primary's mass-loss rate and the companion star's temperature and luminosity. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing (Madura et al. 2015, arXiv:1503.00716). While we initially focus specifically on Eta Carinae, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty ‘pinwheel’ (WR 112, WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where dust can form and survive in these unique objects.
3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock
NASA Astrophysics Data System (ADS)
Mohamed, S.; Mackey, J.; Langer, N.
2013-05-01
Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few × 10-3 M⊙), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.
Photospheric Emission from Collapsar Jets in 3D Relativistic Hydrodynamics
NASA Astrophysics Data System (ADS)
Ito, Hirotaka; Matsumoto, Jin; Nagataki, Shigehiro; Warren, Donald C.; Barkov, Maxim V.
2015-12-01
We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show that non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra even though only thermal photons are injected initially and the effect of non-thermal particles is not considered.
NASA Astrophysics Data System (ADS)
Martin, Laurent; Emma, Gouze
2010-05-01
. Since 2008, these three buoys have been also equipped with nitrate, chlorophyll and oxygen probes. Thanks to these measurements, a 3D hydrodynamic model (TELEMAC©) has been developped and validated to better qualify and quantify the relationships between the salinity of the lagoon, the fresh water inputs (from the powerplant and from the rivers), the water exchanges through the Caronte channel and the wind mixing. This model is currently used by the electricity producer to manage fresh water discharges complying with salinity indicators fixed by the European Court of Justice. Then, a biogeochemical model (DelWAQ©) coupled with the hydrodynamic model has been developped to understand the ecosystem functioning and to assess the hydroelectric powerplant implication in the eutrophication of the lagoon. Simulations reproduce quite well: 1/ the seasonal variations of nutrients, 2/ biogeochemical processes, 3/ anoxia events in connection with stratification periods at deep stations and 4/ are able to calculate nutrients budgets over a year. The results show that high primary production rates are based on high dynamical mineralization processes. The allochtonous nutrients sources are not sufficient to feed the phytoplanktonic demand (less than 1%). This models coupling is the only way to compile the physical and biogeochemical variables and processes. It's a tool aiming at a better assessment of the high complexity of the lagoon nutrients cycles. It will help us to understand the powerplant implication for the eutrophication with respect to the anthropised rivers. Moreover we would be able to test various managing scenarii (e.g. drop of nutrients loads) and to suggest new rehabilitation strategies.
3D Kinematics and Hydrodynamic Analysis of Freely Swimming Cetacean
NASA Astrophysics Data System (ADS)
Ren, Yan; Sheinberg, Dustin; Liu, Geng; Dong, Haibo; Fish, Frank; Javed, Joveria
2015-11-01
It's widely thought that flexibility and the ability to control flexibility are crucial elements in determining the performance of animal swimming. However, there is a lack of quantification of both span-wise and chord-wise deformation of Cetacean's flukes and associated hydrodynamic performance during actively swimming. To fill this gap, we examined the motion and flexure of both dolphin fluke and orca fluke in steady swimming using a combined experimental and computational approach. It is found that the fluke surface morphing can effectively modulate the flow structures and influence the propulsive performance. Findings from this work are fundamental for understanding key kinematic features of effective Cetacean propulsors, and for quantifying the hydrodynamic force production that naturally occurs during different types of swimming. This work is supported by ONR MURI N00014-14-1-0533 and NSF CBET-1313217.
Using the symmetries in 3D Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Cruz-Pérez, Juan P.; González, José A.
2011-10-01
We apply the cartoon SPH technique to solve numerically the Newtonian Euler's equations in scenarios with spherical symmetry with the standard Smoothed Particle Hydrodynamics method, using the equations written in cartesian coordinates. This implementation allow us to increase the resolution of the simulations in order to obtain accurate results. We test it using the Shock tube and the isothermal collapse of a perfect fluid in spherical symmetry.
NASA Astrophysics Data System (ADS)
Moreno Navas, Juan; Telfer, Trevor C.; Ross, Lindsay G.
2011-04-01
Hydrographic conditions, and particularly current speeds, have a strong influence on the management of fish cage culture. These hydrodynamic conditions can be used to predict particle movement within the water column and the results used to optimise environmental conditions for effective site selection, setting of environmental quality standards, waste dispersion, and potential disease transfer. To this end, a 3D hydrodynamic model, MOHID, has been coupled to a particle tracking model to study the effects of mean current speed, quiescent water periods and bulk water circulation in Mulroy Bay, Co. Donegal Ireland, an Irish fjard (shallow fjordic system) important to the aquaculture industry. A Lagangrian method simulated the instantaneous release of "particles" emulating discharge from finfish cages to show the behaviour of waste in terms of water circulation and water exchange. The 3D spatial models were used to identify areas of mixed and stratified water using a version of the Simpson-Hunter criteria, and to use this in conjunction with models of current flow for appropriate site selection for salmon aquaculture. The modelled outcomes for stratification were in good agreement with the direct measurements of water column stratification based on observed density profiles. Calculations of the Simpson-Hunter tidal parameter indicated that most of Mulroy Bay was potentially stratified with a well mixed region over the shallow channels where the water is faster flowing. The fjard was characterised by areas of both very low and high mean current speeds, with some areas having long periods of quiescent water. The residual current and the particle tracking animations created through the models revealed an anticlockwise eddy that may influence waste dispersion and potential for disease transfer, among salmon cages and which ensures that the retention time of waste substances from cages is extended. The hydrodynamic model results were incorporated into the ArcView TM GIS
3D hydrodynamic interactions lead to divergences in 2D diffusion.
Bleibel, Johannes; Domínguez, Alvaro; Oettel, Martin
2015-05-20
We investigate the influence of 3D hydrodynamic interactions on confined colloidal suspensions, where only the colloids are restricted to one or two dimensions. In the absence of static interactions among the colloids, i.e., an ideal gas of colloidal particles with a finite hydrodynamic radius, we find a divergent collective diffusion coefficient. The origin of the divergence is traced back to the dimensional mismatch of 3D hydrodynamic interactions and the colloidal particles moving only in 1D or 2D. Our results from theory are confirmed by Stokesian dynamics simulations and supported by light scattering observational data for particles at a fluid interface. PMID:25923320
3D hydrodynamic interactions lead to divergences in 2D diffusion
NASA Astrophysics Data System (ADS)
Bleibel, Johannes; Domínguez, Alvaro; Oettel, Martin
2015-05-01
We investigate the influence of 3D hydrodynamic interactions on confined colloidal suspensions, where only the colloids are restricted to one or two dimensions. In the absence of static interactions among the colloids, i.e., an ideal gas of colloidal particles with a finite hydrodynamic radius, we find a divergent collective diffusion coefficient. The origin of the divergence is traced back to the dimensional mismatch of 3D hydrodynamic interactions and the colloidal particles moving only in 1D or 2D. Our results from theory are confirmed by Stokesian dynamics simulations and supported by light scattering observational data for particles at a fluid interface.
Ruh, Dominic; Tränkle, Benjamin; Rohrbach, Alexander
2011-10-24
Multi-dimensional, correlated particle tracking is a key technology to reveal dynamic processes in living and synthetic soft matter systems. In this paper we present a new method for tracking micron-sized beads in parallel and in all three dimensions - faster and more precise than existing techniques. Using an acousto-optic deflector and two quadrant-photo-diodes, we can track numerous optically trapped beads at up to tens of kHz with a precision of a few nanometers by back-focal plane interferometry. By time-multiplexing the laser focus, we can calibrate individually all traps and all tracking signals in a few seconds and in 3D. We show 3D histograms and calibration constants for nine beads in a quadratic arrangement, although trapping and tracking is easily possible for more beads also in arbitrary 2D arrangements. As an application, we investigate the hydrodynamic coupling and diffusion anomalies of spheres trapped in a 3 × 3 arrangement. PMID:22109012
The hydrodynamic part of the 3D CEMBS model for the Baltic Sea
NASA Astrophysics Data System (ADS)
Dzierzbicka-Glowacka, Lidia; Jakacki, Jaromir; Janecki, Maciej; Nowicki, Artur
2013-04-01
The paper presents a hydrodynamic part of the coupled ice-ocean model 3D CEMBS designed for the Baltic Sea. It is based on the Community Earth System Model (CESM from the National Center for Atmospheric Research). It was adopted for the Baltic Sea as a coupled sea-ice model. It consists of the Community Ice Code (CICE model, version 4.0) and the Parallel Ocean Program (version 2.1). The models are linked through the coupler (CPL7) based on the Model Coupling Toolkit library. The ocean model has 21 vertical levels and horizontal grid of 600x640 cells. Horizontal resolution is approximately 2km. It is forced by atmospheric fields from European Centre for Medium-Range Weather Forecasts and in operational mode from 48-hour atmospheric forecasts provided by the UM model from the Interdisciplinary Centre for Mathematical and Computational Modelling of Warsaw University (ICM). The study was financially supported by the Polish State Committee of Scientific Research (grants: No N N305 111636, N N306 353239). The partial support for this study was also provided by the project Satellite Monitoring of the Baltic Sea Environment - SatBaltyk founded by European Union through European Regional Development Fund contract no. POIG 01.01.02-22-011/09. Calculations were carried out at the Academy Computer Centre in Gdansk.
A coupled RELAPS-3D/CFD methodology with a proof-of-principle calculation
Aumiller, D.L.; Tomlinson, E.T.; Bauer, R.C.
2000-01-01
The RELAP5-3D computer code was modified to make the explicit coupling capability in the code fully functional. As a test of the modified code, a coupled RELAP5/RELAP5 analysis of the Edwards-O'Brien blowdown problem was performed which showed no significant deviations from the standard RELAP5-3D predictions. In addition, a multiphase Computational Fluid Dynamics (CFD) code was modified to permit explicit coupling to RELAP5-3D. Several calculations were performed with this code. The first analysis used the experimental pressure history from a point just upstream of the break as a boundary condition. This analysis showed that a multiphase CFD code could calculate the thermodynamic and hydrodynamic conditions during a rapid blowdown transient. Finally, a coupled RELAP5/CFD analysis was performed. The results are presented in this paper.
Exact Relativistic Ideal Hydrodynamical Solutions in (1+3)D with Longitudinal and Transverse Flows
Liao, Jinfeng; Koch, Volker
2009-05-20
A new method for solving relativistic ideal hydrodynamics in (1+3)D is developed. Longitudinal and transverse radial flows are explicitly embedded into the ansatz for velocity field and the hydrodynamic equations are reduced to a single equation for the transverse velocity field only, which is analytically more tractable as compared with the full hydrodynamic equations. As an application we use the method to find analytically all possible solutions whose transverse velocity fields have power dependence on proper time and transverse radius. Possible application to the Relativistic Heavy Ion Collisions and possible generalizations of the method are discussed.
Recent Hydrodynamics Improvements to the RELAP5-3D Code
Richard A. Riemke; Cliff B. Davis; Richard.R. Schultz
2009-07-01
The hydrodynamics section of the RELAP5-3D computer program has been recently improved. Changes were made as follows: (1) improved turbine model, (2) spray model for the pressurizer model, (3) feedwater heater model, (4) radiological transport model, (5) improved pump model, and (6) compressor model.
Monolithic cell counter based on 3D hydrodynamic focusing in microfluidic channels
NASA Astrophysics Data System (ADS)
Paiè, Petra; Bragheri, Francesca; Osellame, Roberto
2014-03-01
Hydrodynamic focusing is a powerful technique frequently used in microfluidics that presents a wide range of applications since it allows focusing the sample flowing in the device to a narrow region in the center of the microchannel. In fact thanks to the laminarity of the fluxes in microchannels it is possible to confine the sample solution with a low flow rate by using a sheath flow with a higher flow rate. This in turn allows the flowing of one sample element at a time in the detection region, thus enabling analysis on single particles. Femtosecond laser micromachining is ideally suited to fabricate device integrating full hydrodynamic focusing functionalities thanks to the intrinsic 3D nature of this technique, especially if compared to expensive and complicated lithographic multi-step fabrication processes. Furthermore, because of the possibility to fabricate optical waveguides with the same technology, it is possible to obtain compact optofluidic devices to perform optical analysis of the sample even at the single cell level, as is the case for optical cell stretchers and sorters. In this work we show the fabrication and the fluidic characterization of extremely compact devices having only two inlets for 2D (both in vertical and horizontal planes) as well as full 3D symmetric hydrodynamic focusing. In addition we prove one of the possible application of the hydrodynamic focusing module, by fabricating and validating (both with polystyrene beads and erythrocytes) a monolithic cell counter obtained by integrating optical waveguides in the 3D hydrodynamic focusing device.
Oxygen spectral line synthesis: 3D non-LTE with CO5BOLD hydrodynamical model atmospheres.
NASA Astrophysics Data System (ADS)
Prakapavičius, D.; Steffen, M.; Kučinskas, A.; Ludwig, H.-G.; Freytag, B.; Caffau, E.; Cayrel, R.
In this work we present first results of our current project aimed at combining the 3D hydrodynamical stellar atmosphere approach with non-LTE (NLTE) spectral line synthesis for a number of key chemical species. We carried out a full 3D-NLTE spectrum synthesis of the oxygen IR 777 nm triplet, using a modified and improved version of our NLTE3D package to calculate departure coefficients for the atomic levels of oxygen in a CO5BOLD 3D hydrodynamical solar model atmosphere. Spectral line synthesis was subsequently performed with the Linfor3D code. In agreement with previous studies, we find that the lines of the oxygen triplet produce deeper cores under NLTE conditions, due to the diminished line source function in the line forming region. This means that the solar oxygen IR 777 nm lines should be stronger in NLTE, leading to negative 3D NLTE-LTE abundance corrections. Qualitatively this result would support previous claims for a relatively low solar oxygen abundance. Finally, we outline several further steps that need to be taken in order to improve the physical realism and numerical accuracy of our current 3D-NLTE calculations.
Wind forcing of upland lake hydrodynamics: implementation and validation of a 3D numerical model
NASA Astrophysics Data System (ADS)
Morales, L.; French, J.; Burningham, H.; Evans, C.; Battarbee, R.
2010-12-01
Upland lakes act as important archives of environmental change, yet inferences based on the analysis of sediment cores are frequently compromised by an incomplete understanding of the hydrodynamic processes controlling the distribution and completeness of lake sediment sequences and their linkages to wider environmental factors. Many upland lakes are characterized by complex vertical and horizontal circulation patterns induced by the action of wind on the water surface. Wind forcing is important not only for the resuspension of bottom sediments in shallow marginal areas, but may also control the broader distribution of sediment accumulation. The work presented here represents the first stage of a project aimed at elucidating the linkages between wind forcing and the distribution of bottom sediments in upland lakes and the extent to which simple 'sediment focusing' models provide an adequate basis for predicting optimal locations for the acquisition of core samples for palaeolimnological analysis. As a first step, a 3D numerical hydrodynamic model is implemented for Llyn Conwy, a small oligotrophic upland lake in North Wales, UK. This utilises the community ocean model, FVCOM, that solves the Navier-Stokes equations in 3D on an unstructured triangular mesh using the finite volume method. A new graphical user interface has been developed for FVCOM to facilitate pre- and post-processing of lake modelling problems. At Llyn Conwy, the model is forced using local meteorological data and validated against vertical temperature profiles recorded by a long-term buoy deployment and short-term observations of vertical current structure measured using an upward-looking acoustic doppler profiler and surface circulation obtained from GPS drifters. Challenges in the application of FVCOM to a small lake include the design of a mesh that ensures numerical stability whilst resolving a complex bathymetry, and the need for careful treatment of model 'spin-up'. Once calibrated, the
Description of a parallel, 3D, finite element, hydrodynamics-diffusion code
Milovich, J L; Prasad, M K; Shestakov, A I
1999-04-11
We describe a parallel, 3D, unstructured grid finite element, hydrodynamic diffusion code for inertial confinement fusion (ICF) applications and the ancillary software used to run it. The code system is divided into two entities, a controller and a stand-alone physics code. The code system may reside on different computers; the controller on the user's workstation and the physics code on a supercomputer. The physics code is composed of separate hydrodynamic, equation-of-state, laser energy deposition, heat conduction, and radiation transport packages and is parallelized for distributed memory architectures. For parallelization, a SPMD model is adopted; the domain is decomposed into a disjoint collection of subdomains, one per processing element (PE). The PEs communicate using MPI. The code is used to simulate the hydrodynamic implosion of a spherical bubble.
Sub-micrometer-precision, three-dimensional (3D) hydrodynamic focusing via “microfluidic drifting”
Nawaz, Ahmad Ahsan; Zhang, Xiangjun; Mao, Xiaole; Rufo, Joseph; Lin, Sz-Chin Steven; Guo, Feng; Zhao, Yanhui; Lapsley, Michael; Li, Peng; McCoy, J. Philip; Levine, Stewart J.; Huang, Tony Jun
2014-01-01
In this article, we demonstrate single-layered, “microfluidic drifting” based three-dimensional (3D) hydrodynamic focusing devices with particle/cell focal positioning approaching submicron precision along both lateral and vertical directions. By systematically optimizing channel geometries and sample/sheath flow rates, a series of “microfluidic drifting” based 3D hydrodynamic focusing devices with different curvature angles are designed and fabricated. Their performances are then evaluated by confocal microscopy, fast camera imaging, and side-view imaging techniques. Using a device with a curvature angle of 180°, we have achieved a standard deviation of ±0.45 µm in particle focal position and a coefficient of variation (CV) of 2.37% in flow cytometric measurements. To the best of our knowledge, this is the best CV that has been achieved by a microfluidic flow cytometry device. Moreover, the device showed the capability to distinguish 8 peaks when subjected to a stringent 8-peak rainbow calibration test, signifying the ability to perform sensitive, accurate tests similar to commercial flow cytometers. We have further tested and validated our device by detection of HEK-293 cells. With its advantages in simple fabrication (i.e., single-layered device), precise 3D hydrodynamic focusing (i.e., submicrometer precision along both lateral and vertical directions), and high detection resolution (i.e., low CV), our method could serve as an important basis for high-performance, mass-producible microfluidic flow cytometry. PMID:24287742
Sub-micrometer-precision, three-dimensional (3D) hydrodynamic focusing via "microfluidic drifting".
Nawaz, Ahmad Ahsan; Zhang, Xiangjun; Mao, Xiaole; Rufo, Joseph; Lin, Sz-Chin Steven; Guo, Feng; Zhao, Yanhui; Lapsley, Michael; Li, Peng; McCoy, J Philip; Levine, Stewart J; Huang, Tony Jun
2014-01-21
In this article, we demonstrate single-layered, "microfluidic drifting" based three-dimensional (3D) hydrodynamic focusing devices with particle/cell focal positioning approaching submicron precision along both lateral and vertical directions. By systematically optimizing channel geometries and sample/sheath flow rates, a series of "microfluidic drifting" based 3D hydrodynamic focusing devices with different curvature angles are designed and fabricated. Their performances are then evaluated using confocal microscopy, fast camera imaging, and side-view imaging techniques. Using a device with a curvature angle of 180°, we have achieved a standard deviation of ±0.45 μm in particle focal position and a coefficient of variation (CV) of 2.37% in flow cytometric measurements. To the best of our knowledge, this is the best CV that has been achieved using a microfluidic flow cytometry device. Moreover, the device showed the capability to distinguish 8 peaks when subjected to a stringent 8-peak rainbow calibration test, signifying the ability to perform sensitive, accurate tests similar to commercial flow cytometers. We have further tested and validated our device by detection of HEK-293 cells. With its advantages in simple fabrication (i.e., single-layered device), precise 3D hydrodynamic focusing (i.e., submicrometer precision along both lateral and vertical directions), and high detection resolution (i.e., low CV), our method could serve as an important basis for high-performance, mass-producible microfluidic flow cytometry. PMID:24287742
How do hydrodynamic instabilities affect 3D transport in geophysical vortices?
NASA Astrophysics Data System (ADS)
Wang, Peng; Özgökmen, Tamay M.
2015-03-01
Three-dimensional (3D) transport within geophysical vortices (e.g. ocean eddies) is important in understanding processes at a variety of scales, ranging from plankton production to climate variability. 3D transport can be affected by hydrodynamic instabilities of geophysical vortices; however, how the instabilities affecting 3D transport is not clear. Focusing on barotropic, inertial and 3D instabilities, we investigate the joint impacts of instabilities on 3D transport by using analytical methods and direct numerical simulations. We discover for the first time that material can be exchanged through 3D pathways which link a family of vortices generated by the instabilities in a single, initially unstable vortex. We also show that instabilities can increase the magnitude of vertical velocity, mixing rate and vertical material exchange. Besides, we find that instabilities can cause the kinetic energy wavenumber spectrum to have a power-law regime different than the classic regimes of k - 5 / 3 and k-3, and propose a new energy spectrum to interpret the non-classic regime.
Hydrodynamic growth experiments with the 3-D, “native-roughness” modulations on NIF
NASA Astrophysics Data System (ADS)
Smalyuk, V. A.; Weber, S. V.; Casey, D.; Clark, D. S.; Coppari, F.; Field, J. E.; Haan, S. W.; Hammel, B. A.; Hamza, A.; Hsing, W.; Landen, O.; Nikroo, A.; Robey, H. F.; Weber, C. R.
2016-05-01
Hydrodynamic instability growth experiments with threedimensional (3-D) surface-roughness modulations were performed on plastic (CH) shell spherical implosions at the National Ignition Facility (NIF). The initial capsule outer-surface roughness was similar to the standard specifications (“native roughness”) used in a majority of implosions on NIF. At a convergence ratio of ∼3, the measured tent modulations were close to those predicted by 3-D simulations (within ∼15-20%), while measured 3-D, broadband modulations were ∼3-4 times larger than those simulated based on the growth of the known imposed initial surface modulations. One of the hypotheses to explain the results is based on the increased instability amplitudes due to modulations of the oxygen content in the bulk of the capsule. These new experiments results have prompted looking for ways to reduce UV light exposure during target fabrication.
How Do Hydrodynamic Instabilities Affect 3D Transport in Geophysical Vortices?
NASA Astrophysics Data System (ADS)
Wang, P.; Ozgokmen, T. M.
2014-12-01
Understanding three-dimensional (3D) transport in ocean eddies is important for processes at a variety of scales, ranging from plankton production to climate variability. It is well known that geophysical vortices are subject to various hydrodynamic instabilities. Yet the influence of these instabilities on 3D material transport in vortex systems is not well investigated. Focusing on barotropic, inertial and 3D instabilities, we analyze these instabilities with normal-mode method, and reproduce their characteristics via highly-resolved numerical simulations using a spectral element Navier-Stokes solver. By comparing the simulation results of stable and unstable vortices, we investigate the joint impacts of instabilities on 3D transport through three major aspects: (i) energy transfer, (ii) overturning transport of the secondary circulation, and (iii) rates of vertical exchange and mixing. It is found that instabilities can enhance local nonlinear interactions and cause the kinetic energy wavenumber spectrum to have slopes between the conventional -5/3 and -3 at inertial ranges. The cascade of a new quantity is proposed to explain these non-conventional slopes. One of our main results is the discovery of material exchange between the central vortex and satellite vortices through 3D pathways, called funnels. These funnels modify the concept of elliptic regions that can trap material when confined to 2D dynamics. Thus, we show that a family of vortices, created by the hydrodynamic instabilities of the initially unstable vortex, can still continue to operate in unity in order to complete the 3D transport in these systems. We also show that flow instabilities can double the magnitude of vertical velocity, increase the rate of vertical exchange by an order of magnitude and enhance mixing rate more than 100%.
Numerical 3D Hydrodynamics Study of Gravitational Instabilities in a Circumbinary Disk
NASA Astrophysics Data System (ADS)
Desai, Karna Mahadev; Steiman-Cameron, Thomas Y.; Michael, Scott; Cai, Kai; Durisen, Richard H.
2016-01-01
We present a 3D hydrodynamical study of gravitational instabilities (GIs) in a circumbinary protoplanetary disk around a Solar mass star and a brown dwarf companion (0.02 M⊙). GIs can play an important, and at times dominant, role in driving the structural evolution of protoplanetary disks. The reported simulations were performed employing CHYMERA, a radiative 3D hydrodynamics code developed by the Indiana University Hydrodynamics Group. The simulations include disk self-gravity and radiative cooling governed by realistic dust opacities. We examine the role of GIs in modulating the thermodynamic state of the disks, and determine the strengths of GI-induced density waves, non-axisymmetric density structures, radial mass transport, and gravitational torques. The principal goal of this study is to determine how the presence of the companion affects the nature and strength of GIs. Results are compared with a parallel simulation of a protoplanetary disk without the presence of the brown dwarf binary companion. We detect no fragmentation in either disk. A persistent vortex forms in the inner region of both disks. The vortex seems to be stabilized by the presence of the binary companion.
Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Zhang, Mingyu
2010-09-01
The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503-1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710-1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.
Interchain coupling and 3D modeling of trans-polyacetylene
Bronold, F.; Saxena, A.; Bishop, A.R.
1992-01-01
In spite of the success of the SSH model for trans-polyacetylene in interpreting many experimental results (e.g. optical and magnetic properties) there remain some aspects of the real material which are outside the scope of the simple 1D model. Especially ordering phenomena of doped and undoped trans-polyacetylene as well as transport properties (e.g. electronic and thermal conductivity) are beyond a 1D description. There are many attempts to construct a transport theory for this novel class of materials using solitons or polaxons as the basic ingredients. But so far it is not yet clear whether these typical 1D excitations still exist in crystalline transpolyacetylene. Therefore, to clarify the role which intrinsic self-localized nonlinear excitations characteristic of 1D models play in the bulk (3D) material, we study the stability of a polaronic excitation against interchain coupling. As a preliminary step we consider first two coupled t-(CH){sub x}-chains where the {pi}-electrons are allowed to hop from one chain to the other. Then we introduce a 3D generalization of the SSH model and study a polaron in a 3D crystalline environment.
Interchain coupling and 3D modeling of trans-polyacetylene
Bronold, F.; Saxena, A.; Bishop, A.R.
1992-09-01
In spite of the success of the SSH model for trans-polyacetylene in interpreting many experimental results (e.g. optical and magnetic properties) there remain some aspects of the real material which are outside the scope of the simple 1D model. Especially ordering phenomena of doped and undoped trans-polyacetylene as well as transport properties (e.g. electronic and thermal conductivity) are beyond a 1D description. There are many attempts to construct a transport theory for this novel class of materials using solitons or polaxons as the basic ingredients. But so far it is not yet clear whether these typical 1D excitations still exist in crystalline transpolyacetylene. Therefore, to clarify the role which intrinsic self-localized nonlinear excitations characteristic of 1D models play in the bulk (3D) material, we study the stability of a polaronic excitation against interchain coupling. As a preliminary step we consider first two coupled t-(CH){sub x}-chains where the {pi}-electrons are allowed to hop from one chain to the other. Then we introduce a 3D generalization of the SSH model and study a polaron in a 3D crystalline environment.
3D RNA and Functional Interactions from Evolutionary Couplings.
Weinreb, Caleb; Riesselman, Adam J; Ingraham, John B; Gross, Torsten; Sander, Chris; Marks, Debora S
2016-05-01
Non-coding RNAs are ubiquitous, but the discovery of new RNA gene sequences far outpaces the research on the structure and functional interactions of these RNA gene sequences. We mine the evolutionary sequence record to derive precise information about the function and structure of RNAs and RNA-protein complexes. As in protein structure prediction, we use maximum entropy global probability models of sequence co-variation to infer evolutionarily constrained nucleotide-nucleotide interactions within RNA molecules and nucleotide-amino acid interactions in RNA-protein complexes. The predicted contacts allow all-atom blinded 3D structure prediction at good accuracy for several known RNA structures and RNA-protein complexes. For unknown structures, we predict contacts in 160 non-coding RNA families. Beyond 3D structure prediction, evolutionary couplings help identify important functional interactions-e.g., at switch points in riboswitches and at a complex nucleation site in HIV. Aided by increasing sequence accumulation, evolutionary coupling analysis can accelerate the discovery of functional interactions and 3D structures involving RNA. PMID:27087444
Influence of a postural change of the swimmer's head in hydrodynamic performances using 3D CFD.
Popa, Catalin Viorel; Arfaoui, Ahlem; Fohanno, Stéphane; Taïar, Redha; Polidori, Guillaume
2014-01-01
This study deals with recent researches undertaken by the authors in the field of hydrodynamics of human swimming. The aim of this numerical study was to investigate the flow around the entire swimmer's body. The results presented in this article focus on the combination of a 3D computational fluid dynamics code and the use of the k-ω turbulence model, in the range of Reynolds numbers representative of a swimming level varying from national to international competition. Emphasis is placed on the influence of a postural change of the swimmer's head in hydrodynamic performances, which is directly related to the reduction of overall drag. These results confirm and complete those, less accurate, of a preliminary 2D study recently published by the authors and allow the authors to optimise the swimmer's head position in underwater swimming. PMID:22587390
Nonlinear dynamic analysis of hydrodynamically-coupled stainless steel structures
Zhao, Y.
1996-12-01
Spent nuclear fuel is usually stored temporarily on the site of nuclear power plants. The spent fuel storage racks are nuclear-safety-related stainless steel structures required to be analyzed for seismic loads. When the storage pool is subjected to three-dimensional (3-D) floor seismic excitations, rack modules, stored fuel bundles, adjacent racks and pool walls, and surrounding water are hydrodynamically coupled. Hydrodynamic coupling (HC) significantly affects the dynamic responses of the racks that are free-standing and submerged in water within the pool. A nonlinear time-history dynamic analysis is usually needed to describe the motion behavior of the racks that are both geometrically nonlinear and material nonlinear in nature. The nonlinearities include the friction resistance between the rack supporting legs and the pool floor, and various potential impacts of fuel-rack, rack-rack, and rack-pool wall. The HC induced should be included in the nonlinear dynamic analysis using the added-hydrodynamic-mass concept based on potential theory per the US Nuclear Regulatory Commission (USNRC) acceptance criteria. To this end, a finite element analysis constitutes a feasible and effective tool. However, most people perform somewhat simplified 1-D, or 2-D, or 3-D single rack and 2-D multiple rack analyses. These analyses are incomplete because a 3-D single rack model behaves quite differently from a 2-D mode. Furthermore, a 3-D whole pool multi-rack model behaves differently than a 3-D single rack model, especially when the strong HC effects are unsymmetrical. In this paper 3-D nonlinear dynamic time-history analyses were performed in a more quantitative manner using sophisticated finite element models developed for a single rack as well as all twelve racks in the whole-pool. Typical response results due to different HC effects are determined and discussed.
Simulations of implosions with a 3D, parallel, unstructured-grid, radiation-hydrodynamics code
Kaiser, T B; Milovich, J L; Prasad, M K; Rathkopf, J; Shestakov, A I
1998-12-28
An unstructured-grid, radiation-hydrodynamics code is used to simulate implosions. Although most of the problems are spherically symmetric, they are run on 3D, unstructured grids in order to test the code's ability to maintain spherical symmetry of the converging waves. Three problems, of increasing complexity, are presented. In the first, a cold, spherical, ideal gas bubble is imploded by an enclosing high pressure source. For the second, we add non-linear heat conduction and drive the implosion with twelve laser beams centered on the vertices of an icosahedron. In the third problem, a NIF capsule is driven with a Planckian radiation source.
Shao, Yan-Lin Faltinsen, Odd M.
2014-10-01
We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.
Delft3D-FLOW on PRACE infrastructures for real life hydrodynamic applications.
NASA Astrophysics Data System (ADS)
Donners, John; Genseberger, Menno; Jagers, Bert; de Goede, Erik; Mourits, Adri
2013-04-01
PRACE, the Partnership for Advanced Computing in Europe, offers access to the largest high-performance computing systems in Europe. PRACE invites and helps industry to increase their innovative potential through the use of the PRACE infrastructure. This poster describes different efforts to assist Deltares with porting the open-source simulation software Delft3D-FLOW to PRACE infrastructures. Analysis of the performance on these infrastructures has been done for real life flow applications. Delft3D-FLOW is a 2D and 3D shallow water solver which calculates non-steady flow and transport phenomena resulting from tidal and meteorological forcing on a curvilinear, boundary fitted grid in Cartesian or spherical coordinates. It also includes a module which sediment transport (both suspended and bed total load) and morphological changes for an arbitrary number of cohesive and non-cohesive fractions. As Delft3D-FLOW has been developed over several decades, with a variety of functionality and over 350k lines of source code, porting to PRACE infrastructures needs some effort. At the moment Delft3D-FLOW uses MPI with domain decomposition in one direction as its parallellisation approach. Because it is hard to identify scaling issues if one immediately starts with a complex case with many features enabled, different cases with increasing complexity have been used to investigate scaling of this parallellisation approach on several PRACE platforms. As a base reference case we started with a schematic high-resolution 2D hydrodynamic model of the river Waal that turned out to be surprisingly well-suited to the highly-parallel PRACE machines. Although Delft3D-FLOW employs a sophisticated build system, several modifications were required to port it to most PRACE systems due to the use of specific, highly-tuned compilers and MPI-libraries. After this we moved to a 3D hydrodynamic model of Rotterdam harbour that includes sections of the rivers Rhine and Meuse and a part of the North
Carafa, R; Marinov, D; Dueri, S; Wollgast, J; Ligthart, J; Canuti, E; Viaroli, P; Zaldívar, J M
2006-10-01
Sacca di Goro is a shallow coastal microtidal lagoon with a surface area of 26 km2, and an average depth of about 1.5m. Fresh water pollutant loads from Po River branches and several drainage canals lead to anthropogenic eutrophication, frequent summer anoxia crises and chemical contamination. Such events not only affect the lagoon ecosystem but also cause serious economic losses, the lagoon being the second largest producer of clams in Italy. The present work aims at using a fate model coupled with COHERENS 3D hydrodynamic model to simulate and to explain the spatial distribution and temporal variations of s-triazines herbicides in the Sacca di Goro lagoon. The simulation results of spatial and temporal dynamic behaviour of atrazine, simazine and terbuthylazine have been compared with experimental data obtained during an annual monitoring programme. PMID:16643962
Non-LTE oxygen line formation in 3D hydrodynamic model stellar atmospheres
NASA Astrophysics Data System (ADS)
Amarsi, A. M.; Asplund, M.; Collet, R.; Leenaarts, J.
2016-02-01
The O I 777 nm lines are among the most commonly used diagnostics for the oxygen abundances in the atmospheres of FGK-type stars. However, they form in conditions that are far from local thermodynamic equilibrium (LTE). We explore the departures from LTE of atomic oxygen, and their impact on O I lines, across the STAGGER-grid of three-dimensional hydrodynamic model atmospheres. For the O I 777 nm triplet, we find significant departures from LTE. These departures are larger in stars with larger effective temperatures, smaller surface gravities, and larger oxygen abundances. We present grids of predicted 3D non-LTE based equivalent widths for the O I 616 nm, [O I] 630 nm, [O I] 636 nm, and O I 777 nm lines, as well as abundance corrections to 1D LTE based results.
Size-Controlled Fabrication of Polyaniline Microfibers Based on 3D Hydrodynamic Focusing Approach.
Yoo, Imsung; Song, Simon; Uh, Kyungchan; Lee, Chan Woo; Kim, Jong-Man
2015-07-01
Owing to the relatively high conductivity and unique redox behavior, polyaniline (PANI) has been one of the most technologically promising conducting polymers. Although various methodologies have been developed, fabrication of PANI microfibers has been a challenging task owing to the poor solubility in most organic solvents. By taking advantage of a microfluidic technology and organic soluble acid labile t-Boc-protected PANI (t-Boc-PANI) as the conducting polymer precursor, fabrication of PANI microfibers in a size-controlled manner is possible. Introduction of a THF solution containing t-Boc-PANI, and dodecylbenzenesulfonic acid (DBSA) as a core flow, and water as a sheath flow into a microfluidic channel with a 3D hydrodynamic focusing effect results in crystallization of the polymer fiber. By changing the flow rate, linear PANI microfibers that range from 16.2 to 39.4 μm in diameter are readily obtained. PMID:25882095
3D hydrodynamic simulations of the Galactic supernova remnant CTB 109
NASA Astrophysics Data System (ADS)
Bolte, J.; Sasaki, M.; Breitschwerdt, D.
2015-10-01
Context. Using detailed 3D hydrodynamic simulations we study the nature of the Galactic supernova remnant (SNR) CTB 109 (G109.1-1.0), which is well known for its semicircular shape and a bright diffuse X-ray emission feature inside the SNR. Aims: Our model has been designed to explain the observed morphology, with a special emphasis on the bright emission feature inside the SNR. Moreover, we determine the age of the remnant and compare our findings with X-ray observations. With CTB 109 we test a new method of detailed numerical simulations of diffuse young objects, using realistic initial conditions derived directly from observations. Methods: We performed numerical 3D simulations with the RAMSES code. The initial density structure has been directly taken from 12CO emission data, adding an additional dense cloud, which, when it is shocked, causes the bright emission feature. Results: From parameter studies we obtained the position (ℓ,b) = (109.1545°,-1.0078°) for an elliptical cloud with ncloud = 25 cm-3 based on the preshock density from Chandra data and a maximum diameter of 4.54 pc, whose encounter with the supernova (SN) shock wave generates the bright X-ray emission inside the SNR. The calculated age of the remnant is about 11 000 yr according to our simulations. In addition, we can also determine the most probable site of the SN explosion. Conclusions: Hydrodynamic simulations can reproduce the morphology and the observed size of the SNR CTB 109 remarkably well. Moreover, the simulations show that it is very plausible that the bright X-ray emission inside the SNR is the result of an elliptical dense cloud shocked by the SN explosion wave. We show that numerical simulations using observational data for an initial model can produce meaningful results.
New Insights on Pulsating White Dwarfs from 3D Radiation-Hydrodynamical Simulations
NASA Astrophysics Data System (ADS)
Tremblay, Pier-Emmanuel; Fontaine, Gilles; Ludwig, Hans-Günter
2015-08-01
We have recently computed a grid of 3D radiation-hydrodynamical simulations for the atmosphere of 70 pure-hydrogen DA white dwarfs in the range 7.0 < log g < 9.0. This includes the full ZZ Ceti instability strip where DA white dwarfs are pulsating, by far the most common type of degenerate pulsators. We have significantly improved the theoretical framework to study these objects by removing the free parameters of 1D convection, which were previously a major modeling hurdle. We will compare our new models with the observed sample of ZZ Ceti stars and highlight the improved derived properties of these objects. In particular, the new spectroscopically determined 3D atmospheric parameters allow for an improved definition of instability strip edges. We have also made new predictions for the size of convection zones, which significantly impact the position where the pulsations are driven, and the region of the HR diagram where white dwarfs are expected to pulsate. Finally, we will present new results from non-adiabatic pulsation calculations.
NASA Astrophysics Data System (ADS)
Lacis, K.; Muiznieks, A.; Ratnieks, G.
2005-06-01
A system of three-dimensional numerical models is described to analyse the melt hydrodynamics in the floating zone crystal growth by the needle-eye technique under a rotating magnetic field for the production of high quality silicon single crystals of large diameters big( 100dots 200 mm big). Since the pancake inductor has only one turn, the high frequency (HF) electromagnetic (EM) field and the distribution of heat sources and EM forces on the melt free surface have distinct asymmetric features. This asymmetry together with the displacement of the crystal and feed rod axis and crystal rotation manifests itself as three dimensional hydrodynamic, thermal and dopant concentration fields in the molten zone and causes variations of resistivity in the grown single crystal, which are known as the so-called rotational striations. Additionally, the rotating magnetic field can be used to influence the melt hydrodynamics and to reduce the flow asymmetry. In the present 3D model system, the shape of the molten zone is obtained from symmetric FZ shape calculations. The asymmetric HF EM field is calculated by the 3D boundary element method. The low-frequency rotating magnetic field and a corresponding force density distribution in the melt are calculated by the 3D finite element method. The obtained asymmetric HF field power distribution on the free melt surface, the corresponding HF EM forces and force density of the rotating magnetic field are used for the coupled calculation of 3D steady-state hydrodynamic and temperature fields in the molten zone on a body fitted structured 3D grid by a commercial program package with a control volume approach. Beside the EM forces, also the buoyancy and Marangoni forces are considered. After HD calculations a corresponding 3D dopant concentration field is calculated and used to derive the variations resistivity in the grown crystal. The capability of the system of models is illustrated by a calculation example of a realistic FZ system
NASA Astrophysics Data System (ADS)
Madura, Thomas; Gull, Theodore R.; Clementel, Nicola; Paardekooper, Jan-Pieter; Kruip, Chael; Corcoran, Michael F.; Hamaguchi, Kenji; Teodoro, Mairan
2015-01-01
We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 MSun), highly eccentric (e ~ 0.9) binary Eta Carinae. Using a consumer-grade 3D printer (Makerbot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics simulations of Eta Carinae's inner (r ~110 AU) wind-wind collision interface at multiple orbital phases. These 3D prints reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (φ ~1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the hot, adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unknown physical features highlight the important role 3D printing can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.
High-resolution modelling of 3D hydrodynamics in coastal archipelagos
NASA Astrophysics Data System (ADS)
Miettunen, Elina; Tuomi, Laura; Ropponen, Janne; Lignell, Risto
2016-04-01
Dynamics of the coastal seas are affected by eutrophication, over-fishing, coastal construction and climate change. To enable the sustainable development of these areas, monitoring and modelling of the state of the sea are needed. The Archipelago Sea, located in the northern part of the semi-enclosed and brackish water Baltic Sea, is one of the most complex coastal areas with over 40 000 small islands and islets. It is also very vulnerable area already heavily stressed with eutrophication. Applicable modelling tools are needed to support the decision making and to provide sufficiently reliable information on the effects of the planned actions on the state of the coastal waters. We used 3D hydrodynamic model COHERENS to model the Archipelago Sea area with high spatial resolution of 0.25 nmi. Boundary conditions for this limited area were provided from coarser resolution, 2 nmi, Baltic Sea grid. In order to evaluate the performance of the high-resolution coastal model implementation a comprehensive measurement dataset was gathered, including hydrographic data from three intensive monitoring stations and several more rarely visited monitoring or research stations. The hydrodynamic model was able to simulate the surface temperature and salinity fields and their seasonal variation with good accuracy in this complex area. The sharp depth gradients typical for this area provided some challenges to the modelling. There was some over mixing and related to too strong vertical currents in the steep slopes of the deeper fault lines. Also the water exchange between the more open sea and coastal areas through narrow channels between the islands is not sufficiently well reproduced with the current resolution, leading to too high bottom temperatures.
Excitation of anti-symmetric coupled spoof SPPs in 3D SIS waveguides based on coupling
NASA Astrophysics Data System (ADS)
Li-li, Tian; Yang, Chen; Jian-long, Liu; Kai, Guo; Ke-ya, Zhou; Yang, Gao; Shu-tian, Liu
2016-07-01
According to the electromagnetic field distributions, there exist two kinds of coupled spoof surface plasmon polaritons (SSPPs), the symmetric and anti-symmetric modes, in the three-dimensional (3D) subwavelength spoof–insulator–spoof (SIS) waveguide. We study the dispersion and excitation of the two kinds of coupled SSPPs supported by the 3D SIS waveguide. The evolution of the dispersion with the thickness and gap width of the waveguide is numerically investigated, and we give a theoretical analysis according to the coupling mechanism. Specially, based on the coupling mechanism, we design a zipper structure, through which the excitation and propagation of the anti-symmetric coupled modes can be realized effectively. Project supported by the National Basic Research Program of China (Grant No. 2013CBA01702) and the National Natural Science Foundation of China (Grant Nos. 61377016, 61575055, 10974039, 61307072, 61308017, and 61405056).
Ameur, Houari; Bouzit, Mohamed
2013-08-16
The 3D hydrodynamics and shear rates distributions within the United States Pharmacopeia Apparatus 2 have been investigated in this paper. With the help of a CFD package, several geometric modifications to the device were evaluated in this study. Specially, we examine the influence of impeller clearance, blade diameter, shape of the vessel base and shape of the lower part of blade. Increasing the impeller clearance was observed to exacerbate the heterogeneity in shear and would likely result in greater variability in dissolution measurements. Use of moderate blade diameter and dished bottom were observed to reduce shear heterogeneity in the regions where tablets are most likely to visit during testing. The comparative analysis shows better reproducibility and accelerated dissolution rates with the modified vessel shape, the dished bottom can enhance mixing near the vessel base when compared with the flat bottom. Increasing length of the lower edge of the paddle was observed to generate high radial pumping and to enlarge the dead zone located at the center of the vessel base. PMID:23680733
3-D hydrodynamic modelling of flood impacts on a building and indoor flooding processes
NASA Astrophysics Data System (ADS)
Gems, Bernhard; Mazzorana, Bruno; Hofer, Thomas; Sturm, Michael; Gabl, Roman; Aufleger, Markus
2016-06-01
Given the current challenges in flood risk management and vulnerability assessment of buildings exposed to flood hazards, this study presents three-dimensional numerical modelling of torrential floods and its interaction with buildings. By means of a case study application, the FLOW-3D software is applied to the lower reach of the Rio Vallarsa torrent in the village of Laives (Italy). A single-family house on the flood plain is therefore considered in detail. It is exposed to a 300-year flood hydrograph. Different building representation scenarios, including an entire impervious building envelope and the assumption of fully permeable doors, light shafts and windows, are analysed. The modelling results give insight into the flooding process of the building's interior, the impacting hydrodynamic forces on the exterior and interior walls, and further, they quantify the impact of the flooding of a building on the flow field on the surrounding flood plain. The presented study contributes to the development of a comprehensive physics-based vulnerability assessment framework. For pure water floods, this study presents the possibilities and limits of advanced numerical modelling techniques within flood risk management and, thereby, the planning of local structural protection measures.
Linking 1D evolutionary to 3D hydrodynamical simulations of massive stars
NASA Astrophysics Data System (ADS)
Cristini, A.; Meakin, C.; Hirschi, R.; Arnett, D.; Georgy, C.; Viallet, M.
2016-03-01
Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The ‘stiffness’ of a convective boundary can be quantified using the bulk Richardson number ({{Ri}}{{B}}), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A ‘stiff’ boundary ({{Ri}}{{B}}˜ {10}4) will suppress CBM, whereas in the opposite case a ‘soft’ boundary ({{Ri}}{{B}}˜ 10) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are ‘stiffer’ than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett (2007 Astrophys. J. 667 448-75). This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in super-asymptotic giant branch stars and also the onset of novae.
Interface requirements to couple thermal-hydraulic codes to 3D neutronic codes
Langenbuch, S.; Austregesilo, H.; Velkov, K.
1997-07-01
The present situation of thermalhydraulics codes and 3D neutronics codes is briefly described and general considerations for coupling of these codes are discussed. Two different basic approaches of coupling are identified and their relative advantages and disadvantages are discussed. The implementation of the coupling for 3D neutronics codes in the system ATHLET is presented. Meanwhile, this interface is used for coupling three different 3D neutronics codes.
Dynamic coupling of three hydrodynamic models
NASA Astrophysics Data System (ADS)
Hartnack, J. N.; Philip, G. T.; Rungoe, M.; Smith, G.; Johann, G.; Larsen, O.; Gregersen, J.; Butts, M. B.
2008-12-01
The need for integrated modelling is evidently present within the field of flood management and flood forecasting. Engineers, modellers and managers are faced with flood problems which transcend the classical hydrodynamic fields of urban, river and coastal flooding. Historically the modeller has been faced with having to select one hydrodynamic model to cover all the aspects of the potentially complex dynamics occurring in a flooding situation. Such a single hydrodynamic model does not cover all dynamics of flood modelling equally well. Thus the ideal choice may in fact be a combination of models. Models combining two numerical/hydrodynamic models are becoming more standard, typically these models combine a 1D river model with a 2D overland flow model or alternatively a 1D sewer/collection system model with a 2D overland solver. In complex coastal/urban areas the flood dynamics may include rivers/streams, collection/storm water systems along with the overland flow. The dynamics within all three areas is of the same time scale and there is feedback in the system across the couplings. These two aspects dictate a fully dynamic three way coupling as opposed to running the models sequentially. It will be shown that the main challenges of the three way coupling are time step issues related to the difference in numerical schemes used in the three model components and numerical instabilities caused by the linking of the model components. MIKE FLOOD combines the models MIKE 11, MIKE 21 and MOUSE into one modelling framework which makes it possible to couple any combination of river, urban and overland flow fully dynamically. The MIKE FLOOD framework will be presented with an overview of the coupling possibilities. The flood modelling concept will be illustrated through real life cases in Australia and in Germany. The real life cases reflect dynamics and interactions across all three model components which are not possible to reproduce using a two-way coupling alone. The
PCB modeling in the Gulf of Lions using a 3D coupled model
NASA Astrophysics Data System (ADS)
Alekseenko, Elena; Thouvenin, Bénédicte; Tixier, Céline; Tronczynski, Jacek; Garreau, Pierre; Verney, Romaric; Carlotti, Francois; Espinasse, Boris; Queguiner, Bernard; Baklouti, Melika
2013-04-01
Polychlorobiphenyls (PCBs) are synthetic chlorinated organic compounds, which were widely used in many industrial materials. These compounds are persistent, bioaccumulable and toxic for living organisms. The riverine and atmospheric fluxes are the major routes of entry for these chemicals into marine ecosystems, where they are now embedded in natural biogeochemical cycles (Lohmann et al. 2007). Because of bioaccumulation and biomagnification processes in food webs, even nowadays, these compounds may attain dangerous concentration levels especially in the top predators including marine mammals. The contamination of marine biota by PCBs in Mediterranean has also become a matter of concern as the concentrations in some species are at levels putting them at risk for significant biological effects. This may pose potential human health risks in commercial edible species (Carpenter 2006). Planktonic populations play a key role in the trophic food webs in marine ecosystems by the mobilisation and transfer of energy and organic matter towards higher trophic levels. This work aims at a better understanding of the role of plankton in the transfer of PCBs to higher trophic levels in the Gulf of Lions (Mediterranean) by coupling of biogeochemical, ecological and hydrodynamical processes. Modeling is a powerful tool for coupling processes of different disciplines and scales. The recent development of 3D hydrodynamic, hydrosedimentary and biogeochemical models in the Mediterranean (André et al, 2005,2009, Ulses et al, 2008, Dufois et al, 2008, Auger et al, 2011), enables feasibility testing of coupling these models with transfer processes of chemical contaminants. The lack of detailed observations in the sea and the significant uncertainty on contaminants inputs prevent from a proper validation of such modeling tests. However, these tools are very useful to assess the influence of fast processes on the transfer of contaminants to bioaccumulative species. Sensitivity analysis
Andréfouët, S; Ouillon, S; Brinkman, R; Falter, J; Douillet, P; Wolk, F; Smith, R; Garen, P; Martinez, E; Laurent, V; Lo, C; Remoissenet, G; Scourzic, B; Gilbert, A; Deleersnijder, E; Steinberg, C; Choukroun, S; Buestel, D
2006-10-01
A workshop organized in French Polynesia in November 2004 allowed reviewing the current methods to model the three-dimensional hydrodynamic circulation in semi-enclosed atoll lagoons for aquaculture applications. Mollusk (e.g. pearl oyster, clam) aquaculture is a major source of income for South Pacific countries such as French Polynesia or Cook Islands. This aquaculture now requires a better understanding of circulation patterns to improve the spatial use of the lagoons, especially to define the best area to set larvae collectors. The pelagic larval duration of the relevant species (<20 days) and the size of the semi-closed lagoons (few hundreds of km2) drive the specifications of the model in terms of the spatial and temporal scale. It is considered that, in contrast with fish, mollusk larvae movements are limited and that their cycle occurs completely in the lagoon, without an oceanic stage. Atolls where aquaculture is productive are generally well-bounded, or semi-closed, without significant large and deep openings to the ocean. Nevertheless part of the lagoon circulation is driven by oceanic water inputs through the rim, ocean swells, tides and winds. Therefore, boundary conditions of the lagoon system are defined by the spatial structure of a very shallow rim (exposition and number of hoas), the deep ocean swell climate, tides and wind regimes. To obtain a realistic 3D numerical model of lagoon circulation with adequate forcing, it is thus necessary to connect in an interdisciplinary way a variety of methods (models, remote sensing and in situ data collection) to accurately represent the different components of the lagoon system and its specific boundary conditions. We review here the current methods and tools used to address these different components for a hypothetical atoll of the Tuamotu Archipelago (French Polynesia), representative of the semi-closed lagoons of the South Pacific Ocean. We hope this paper will serve as a guide for similar studies
Generalized hydrodynamics model for strongly coupled plasmas
NASA Astrophysics Data System (ADS)
Diaw, A.; Murillo, M. S.
2015-07-01
Beginning with the exact equations of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, we obtain the density, momentum, and stress tensor-moment equations. We close the moment equations with two closures, one that guarantees an equilibrium state given by density-functional theory and another that includes collisions in the relaxation of the stress tensor. The introduction of a density functional-theory closure ensures self-consistency in the equation-of-state properties of the plasma (ideal and excess pressure, electric fields, and correlations). The resulting generalized hydrodynamics thus includes all impacts of Coulomb coupling, viscous damping, and the high-frequency (viscoelastic) response. We compare our results with those of several known models, including generalized hydrodynamic theory and models obtained using the Singwi-Tosi-Land-Sjolander approximation and the quasilocalized charge approximation. We find that the viscoelastic response, including both the high-frequency elastic generalization and viscous wave damping, is important for correctly describing ion-acoustic waves. We illustrate this result by considering three very different systems: ultracold plasmas, dusty plasmas, and dense plasmas. The new model is validated by comparing its results with those of the current autocorrelation function obtained from molecular-dynamics simulations of Yukawa plasmas, and the agreement is excellent. Generalizations of this model to mixtures and quantum systems should be straightforward.
3D neutronic/thermal-hydraulic coupled analysis of MYRRHA
Vazquez, M.; Martin-Fuertes, F.
2012-07-01
The current tendency in multiphysics calculations applied to reactor physics is the use of already validated computer codes, coupled by means of an iterative approach. In this paper such an approach is explained concerning neutronics and thermal-hydraulics coupled analysis with MCNPX and COBRA-IV codes using a driver program and file exchange between codes. MCNPX provides the neutronic analysis of heterogeneous nuclear systems, both in critical and subcritical states, while COBRA-IV is a subchannel code that can be used for rod bundles or core thermal-hydraulics analysis. In our model, the MCNP temperature dependence of nuclear data is handled via pseudo-material approach, mixing pre-generated cross section data set to obtain the material with the desired cross section temperature. On the other hand, COBRA-IV has been updated to allow for the simulation of liquid metal cooled reactors. The coupled computational tool can be applied to any geometry and coolant, as it is the case of single fuel assembly, at pin-by-pin level, or full core simulation with the average pin of each fuel-assembly. The coupling tool has been applied to the critical core layout of the SCK-CEN MYRRHA concept, an experimental LBE cooled fast reactor presently in engineering design stage. (authors)
Transitional behavior in hydrodynamically coupled oscillators
NASA Astrophysics Data System (ADS)
Box, S.; Debono, L.; Phillips, D. B.; Simpson, S. H.
2015-02-01
In this article we consider the complete set of synchronized and phase-locked states available to pairs of hydrodynamically coupled colloidal rotors, consisting of spherical beads driven about circular paths in the same, and in opposing senses. Oscillators such as these have previously been used as coarse grained, minimal models of beating cilia. Two mechanisms are known to be important in establishing synchrony. The first involves perturbation of the driving force, and the second involves deformation of the rotor trajectory. We demonstrate that these mechanisms are of similar strength, in the regime of interest, and interact to determine observed behavior. Combining analysis and simulation with experiments performed using holographic optical tweezers, we show how varying the amplitude of the driving force perturbation leads to a transition from synchronized to phase-locked states. Analogies with biological systems are discussed, as are implications for the design of biomimetic devices.
One-layer microfluidic device for hydrodynamic 3D self-flow-focusing operating in low flow speed
NASA Astrophysics Data System (ADS)
Daghighi, Yasaman; Gnyawali, Vaskar; Strohm, Eric M.; Tsai, Scott S. H.; Kolios, Michael C.
2016-03-01
Hydrodynamic 3D flow-focusing techniques in microfluidics are categorized as (a) sheathless techniques which require high flow rates and long channels, resulting in high operating cost and high flow rates which are inappropriate for applications with flow rate limitations, and (b) sheath-flow based techniques which usually require excessive sheath flow rate to achieve hydrodynamic 3D flow-focusing. Many devices based on these principles use complicated fabrication methods to create multi-layer microchannels. We have developed a sheath-flow based microfluidic device that is capable of hydrodynamic 3D self-flow-focusing. In this device the main flow (black ink) in a low speed, and a sheath flow, enter through two inlets and enter a 180 degree curved channel (300 × 300 μm cross-section). Main flow migrates outwards into the sheath-flow due to centrifugal effects and consequently, vertical focusing is achieved at the end of the curved channel. Then, two other sheath flows horizontally confine the main flow to achieve horizontal focusing. Thus, the core flow is three-dimensionally focused at the center of the channel at the downstream. Using centrifugal force for 3D flow-focusing in a single-layer fabricated microchannel has been previously investigated by few groups. However, their demonstrated designs required high flow speed (>1 m/s) which is not suitable for many applications that live biomedical specie are involved. Here, we introduce a new design which is operational in low flow speed (<0.05 m/s) and is suitable for applications involving live cells. This microfluidic device can be used in detecting, counting and isolating cells in many biomedical applications.
NASA Technical Reports Server (NTRS)
Asenov, Asen; Brown, A. R.; Slavcheva, G.; Davies, J. H.
2000-01-01
When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain
Matter coupling in 3D ‘minimal massive gravity’
NASA Astrophysics Data System (ADS)
Arvanitakis, Alex S.; Routh, Alasdair J.; Townsend, Paul K.
2014-12-01
The ‘minimal massive gravity’ model of massive gravity in three spacetime dimensions (which has the same anti-de Sitter (AdS) bulk properties as ‘topologically massive gravity’ but improved boundary properties) is coupled to matter. Consistency requires a particular matter source tensor, which is quadratic in the stress tensor. The consequences are explored for an ideal fluid in the context of asymptotically de Sitter (dS) cosmological solutions, which bounce smoothly from contraction to expansion. Various vacuum solutions are also found, including warped (A)dS, and (for special values of parameters) static black holes and an (A)dS2× {{S}1} vacuum.
NASA Astrophysics Data System (ADS)
Hostache, Renaud; Krein, Andreas; Barrière, Julien
2014-05-01
Coupling the 3D hydro-morphodynamic model Telemac-3D-sisyphe and seismic measurements to estimate bedload transport rates in a small gravel-bed river. Renaud Hostache, Andreas Krein, Julien Barrière During flood events, amounts of river bed material are transported via bedload. This causes problems, like the silting of reservoirs or the disturbance of biological habitats. Some current bedload measuring techniques have limited possibilities for studies in high temporal resolutions. Optical systems are usually not applicable because of high turbidity due to concentrated suspended sediment transported. Sediment traps or bedload samplers yield only summative information on bedload transport with low temporal resolution. An alternative bedload measuring technique is the use of seismological systems installed next to the rivers. The potential advantages are observations in real time and under undisturbed conditions. The study area is a 120 m long reach of River Colpach (21.5 km2), a small gravel bed river in Northern Luxembourg. A combined approach of hydro-climatological observations, hydraulic measurements, sediment sampling, and seismological measurements is used in order to investigate bedload transport phenomena. Information derived from seismic measurements and results from a 3-dimensional hydro-morphodynamic model are exemplarily discussed for a November 2013 flood event. The 3-dimensional hydro-morphodynamic model is based on the Telemac hydroinformatic system. This allows for dynamically coupling a 3D hydrodynamic model (Telemac-3D) and a morphodynamic model (Sisyphe). The coupling is dynamic as these models exchange their information during simulations. This is a main advantage as it allows for taking into account the effects of the morphologic changes of the riverbed on the water hydrodynamic and the bedload processes. The coupled model has been calibrated using time series of gauged water depths and time series of bed material collected sequentially (after
3D conductive coupling for efficient generation of prominent Fano resonances in metamaterials.
Liu, Zhiguang; Liu, Zhe; Li, Jiafang; Li, Wuxia; Li, Junjie; Gu, Changzhi; Li, Zhi-Yuan
2016-01-01
We demonstrate a 3D conductive coupling mechanism for the efficient generation of prominent and robust Fano resonances in 3D metamaterials (MMs) formed by integrating vertical U-shape split-ring resonators (SRRs) or vertical rectangular plates along a planar metallic hole array with extraordinary optical transmission (EOT). In such a configuration, intensified vertical E-field is induced along the metallic holes and naturally excites the electric resonances of the vertical structures, which form non-radiative "dark" modes. These 3D conductive "dark" modes strongly interfere with the "bright" resonance mode of the EOT structure, generating significant Fano resonances with both prominent destructive and constructive interferences. The demonstrated 3D conductive coupling mechanism is highly universal in that both 3D MMs with vertical SRRs and vertical plates exhibit the same prominent Fano resonances despite their dramatic structural difference, which is conceptually different from conventional capacitive and inductive coupling mechanisms that degraded drastically upon small structural deviations. PMID:27296109
Cholla: 3D GPU-based hydrodynamics code for astrophysical simulation
NASA Astrophysics Data System (ADS)
Schneider, Evan E.; Robertson, Brant E.
2016-07-01
Cholla (Computational Hydrodynamics On ParaLLel Architectures) models the Euler equations on a static mesh and evolves the fluid properties of thousands of cells simultaneously using GPUs. It can update over ten million cells per GPU-second while using an exact Riemann solver and PPM reconstruction, allowing computation of astrophysical simulations with physically interesting grid resolutions (>256^3) on a single device; calculations can be extended onto multiple devices with nearly ideal scaling beyond 64 GPUs.
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
Coupling hydrodynamic and wave propagation modeling for waveform modeling of SPE.
NASA Astrophysics Data System (ADS)
Larmat, C. S.; Steedman, D. W.; Rougier, E.; Delorey, A.; Bradley, C. R.
2015-12-01
The goal of the Source Physics Experiment (SPE) is to bring empirical and theoretical advances to the problem of detection and identification of underground nuclear explosions. This paper presents effort to improve knowledge of the processes that affect seismic wave propagation from the hydrodynamic/plastic source region to the elastic/anelastic far field thanks to numerical modeling. The challenge is to couple the prompt processes that take place in the near source region to the ones taking place later in time due to wave propagation in complex 3D geologic environments. In this paper, we report on results of first-principles simulations coupling hydrodynamic simulation codes (Abaqus and CASH), with a 3D full waveform propagation code, SPECFEM3D. Abaqus and CASH model the shocked, hydrodynamic region via equations of state for the explosive, borehole stemming and jointed/weathered granite. LANL has been recently employing a Coupled Euler-Lagrange (CEL) modeling capability. This has allowed the testing of a new phenomenological model for modeling stored shear energy in jointed material. This unique modeling capability has enabled highfidelity modeling of the explosive, the weak grout-filled borehole, as well as the surrounding jointed rock. SPECFEM3D is based on the Spectral Element Method, a direct numerical method for full waveform modeling with mathematical accuracy (e.g. Komatitsch, 1998, 2002) thanks to its use of the weak formulation of the wave equation and of high-order polynomial functions. The coupling interface is a series of grid points of the SEM mesh situated at the edge of the hydrodynamic code domain. Displacement time series at these points are computed from output of CASH or Abaqus (by interpolation if needed) and fed into the time marching scheme of SPECFEM3D. We will present validation tests and waveforms modeled for several SPE tests conducted so far, with a special focus on effect of the local topography.
NASA Astrophysics Data System (ADS)
Shi, Jian; Zhang, Qian
2016-03-01
A uniqueness result of weak solution for the 3D viscous magneto-hydrodynamics equations in {B^1_{infty,infty}} is proved by means of the Fourier localization technique and the losing derivative estimates.
Hydrodynamics Versus Intracellular Coupling in the Synchronization of Eukaryotic Flagella
NASA Astrophysics Data System (ADS)
Quaranta, Greta; Aubin-Tam, Marie-Eve; Tam, Daniel
2015-12-01
The influence of hydrodynamic forces on eukaryotic flagella synchronization is investigated by triggering phase locking between a controlled external flow and the flagella of C. reinhardtii. Hydrodynamic forces required for synchronization are over an order of magnitude larger than hydrodynamic forces experienced in physiological conditions. Our results suggest that synchronization is due instead to coupling through cell internal fibers connecting the flagella. This conclusion is confirmed by observations of the vfl3 mutant, with impaired mechanical connection between the flagella.
PCB modeling in the Gulf of Lions using a 3D coupled model
NASA Astrophysics Data System (ADS)
Alekseenko, Elena; Thouvenin, Bénédicte; Tixier, Céline; Tronczynski, Jacek; Garreau, Pierre; Verney, Romaric; Carlotti, Francois; Espinasse, Boris; Queguiner, Bernard; Baklouti, Melika
2013-04-01
Polychlorobiphenyls (PCBs) are synthetic chlorinated organic compounds, which were widely used in many industrial materials. These compounds are persistent, bioaccumulable and toxic for living organisms. The riverine and atmospheric fluxes are the major routes of entry for these chemicals into marine ecosystems, where they are now embedded in natural biogeochemical cycles (Lohmann et al. 2007). Because of bioaccumulation and biomagnification processes in food webs, even nowadays, these compounds may attain dangerous concentration levels especially in the top predators including marine mammals. The contamination of marine biota by PCBs in Mediterranean has also become a matter of concern as the concentrations in some species are at levels putting them at risk for significant biological effects. This may pose potential human health risks in commercial edible species (Carpenter 2006). Planktonic populations play a key role in the trophic food webs in marine ecosystems by the mobilisation and transfer of energy and organic matter towards higher trophic levels. This work aims at a better understanding of the role of plankton in the transfer of PCBs to higher trophic levels in the Gulf of Lions (Mediterranean) by coupling of biogeochemical, ecological and hydrodynamical processes. Modeling is a powerful tool for coupling processes of different disciplines and scales. The recent development of 3D hydrodynamic, hydrosedimentary and biogeochemical models in the Mediterranean (André et al, 2005,2009, Ulses et al, 2008, Dufois et al, 2008, Auger et al, 2011), enables feasibility testing of coupling these models with transfer processes of chemical contaminants. The lack of detailed observations in the sea and the significant uncertainty on contaminants inputs prevent from a proper validation of such modeling tests. However, these tools are very useful to assess the influence of fast processes on the transfer of contaminants to bioaccumulative species. Sensitivity analysis
3D hydrodynamic lift force model for AREVA fuel assembly in EDF PWRs
Ekomie, S.; Bigot, J.; Dolleans, Ph.; Vallory, J.
2007-07-01
The accurate knowledge of the hydrodynamic lift force acting on a fuel assembly in PWR core is necessary to design the hold-down system of this assembly. This paper presents the model used by AREVA NP and EDF for computing this force. It results from a post-processing of sub-channel thermal-hydraulic codes respectively porous medium approach code THYC (EDF) and sub-channel type code FLICA III-F (AREVA NP). This model is based on the application of the Euler's theorem. Some hypotheses used to simplify the complexity of fuel assembly geometry are supported by CFD calculations. Then the model is compared to some experimental results obtained on a single fuel assembly inserted in the HERMES-T test facility located in CEA - Cadarache. Finally, the model is applied to calculate the lift force for the whole core. Various loading patterns including homogenous and mixed cores have been investigated and compared. (authors)
3D analysis of eddy current loss in the permanent magnet coupling
NASA Astrophysics Data System (ADS)
Zhu, Zina; Meng, Zhuo
2016-07-01
This paper first presents a 3D analytical model for analyzing the radial air-gap magnetic field between the inner and outer magnetic rotors of the permanent magnet couplings by using the Amperian current model. Based on the air-gap field analysis, the eddy current loss in the isolation cover is predicted according to the Maxwell's equations. A 3D finite element analysis model is constructed to analyze the magnetic field spatial distributions and vector eddy currents, and then the simulation results obtained are analyzed and compared with the analytical method. Finally, the current losses of two types of practical magnet couplings are measured in the experiment to compare with the theoretical results. It is concluded that the 3D analytical method of eddy current loss in the magnet coupling is viable and could be used for the eddy current loss prediction of magnet couplings.
3D analysis of eddy current loss in the permanent magnet coupling.
Zhu, Zina; Meng, Zhuo
2016-07-01
This paper first presents a 3D analytical model for analyzing the radial air-gap magnetic field between the inner and outer magnetic rotors of the permanent magnet couplings by using the Amperian current model. Based on the air-gap field analysis, the eddy current loss in the isolation cover is predicted according to the Maxwell's equations. A 3D finite element analysis model is constructed to analyze the magnetic field spatial distributions and vector eddy currents, and then the simulation results obtained are analyzed and compared with the analytical method. Finally, the current losses of two types of practical magnet couplings are measured in the experiment to compare with the theoretical results. It is concluded that the 3D analytical method of eddy current loss in the magnet coupling is viable and could be used for the eddy current loss prediction of magnet couplings. PMID:27475575
Evolution of M82-like starburst winds revisited: 3D radiative cooling hydrodynamical simulations
NASA Astrophysics Data System (ADS)
Melioli, C.; de Gouveia Dal Pino, E. M.; Geraissate, F. G.
2013-04-01
In this study we present three-dimensional radiative cooling hydrodynamical simulations of galactic winds generated particularly in M82-like starburst galaxies. We have considered intermittent winds induced by supernova (SN) explosions within super star clusters randomly distributed (in space and time) in the central region of the galaxy (within a radius of R = 150 pc) and were able to reproduce the observed M82 wind conditions with its complex morphological outflow structure. We have found that the environmental conditions in the disc in the nearly recent past are crucial to determine whether the wind will develop a large-scale rich filamentary structure, as in M82 wind, or not. If a sufficiently large number of super stellar clusters are built up in a starburst mainly over a period of a few million years, then the simulations reproduce the multiphase gas observed in M82-like winds, i.e. with filaments of sizes about 20-300 pc, velocities of ˜200-500 km s-1, densities in the range 10-1-10 cm-3, embedded in a hot, low-density gas with a density smaller than 10-2 cm-3 and a velocity of ˜2000 km s-1. Otherwise, a `superbubble-like' wind develops, with very poor or no cold filamentary structures. Also, the numerical evolution of the SN ejecta has allowed us to obtain the abundance distribution over the first ˜3 kpc extension of the wind and we have found that the SN explosions change significantly the metallicity only of the hot, low-density wind component for which we obtained abundances ˜5-10 Z⊙ in fair consistency with the observations. Moreover, we have found that the SN-driven wind transports to outside the disc large amounts of energy, momentum and gas, but the more massive high-density component reaches only intermediate altitudes smaller than 1.5 kpc. Therefore, no significant amounts of gas mass are lost to the intergalactic medium and the mass evolution of the galaxy is not much affected by the starburst events occurring in the nuclear region.
3D conductive coupling for efficient generation of prominent Fano resonances in metamaterials
Liu, Zhiguang; Liu, Zhe; Li, Jiafang; Li, Wuxia; Li, Junjie; Gu, Changzhi; Li, Zhi-Yuan
2016-01-01
We demonstrate a 3D conductive coupling mechanism for the efficient generation of prominent and robust Fano resonances in 3D metamaterials (MMs) formed by integrating vertical U-shape split-ring resonators (SRRs) or vertical rectangular plates along a planar metallic hole array with extraordinary optical transmission (EOT). In such a configuration, intensified vertical E-field is induced along the metallic holes and naturally excites the electric resonances of the vertical structures, which form non-radiative “dark” modes. These 3D conductive “dark” modes strongly interfere with the “bright” resonance mode of the EOT structure, generating significant Fano resonances with both prominent destructive and constructive interferences. The demonstrated 3D conductive coupling mechanism is highly universal in that both 3D MMs with vertical SRRs and vertical plates exhibit the same prominent Fano resonances despite their dramatic structural difference, which is conceptually different from conventional capacitive and inductive coupling mechanisms that degraded drastically upon small structural deviations. PMID:27296109
3D conductive coupling for efficient generation of prominent Fano resonances in metamaterials
NASA Astrophysics Data System (ADS)
Liu, Zhiguang; Liu, Zhe; Li, Jiafang; Li, Wuxia; Li, Junjie; Gu, Changzhi; Li, Zhi-Yuan
2016-06-01
We demonstrate a 3D conductive coupling mechanism for the efficient generation of prominent and robust Fano resonances in 3D metamaterials (MMs) formed by integrating vertical U-shape split-ring resonators (SRRs) or vertical rectangular plates along a planar metallic hole array with extraordinary optical transmission (EOT). In such a configuration, intensified vertical E-field is induced along the metallic holes and naturally excites the electric resonances of the vertical structures, which form non-radiative “dark” modes. These 3D conductive “dark” modes strongly interfere with the “bright” resonance mode of the EOT structure, generating significant Fano resonances with both prominent destructive and constructive interferences. The demonstrated 3D conductive coupling mechanism is highly universal in that both 3D MMs with vertical SRRs and vertical plates exhibit the same prominent Fano resonances despite their dramatic structural difference, which is conceptually different from conventional capacitive and inductive coupling mechanisms that degraded drastically upon small structural deviations.
NASA Astrophysics Data System (ADS)
Sonoi, T.; Samadi, R.; Belkacem, K.; Ludwig, H.-G.; Caffau, E.; Mosser, B.
2015-11-01
Context. The CoRoT and Kepler space-borne missions have provided us with a wealth of high-quality observational data that allows for seismic inferences of stellar interiors. This requires the computation of precise and accurate theoretical frequencies, but imperfect modeling of the uppermost stellar layers introduces systematic errors. To overcome this problem, an empirical correction has been introduced by Kjeldsen et al. (2008, ApJ, 683, L175) and is now commonly used for seismic inferences. Nevertheless, we still lack a physical justification allowing for the quantification of the surface-effect corrections. Aims: Our aim is to constrain the surface-effect corrections across the Hertzsprung-Russell (HR) diagram using a set of 3D hydrodynamical simulations. Methods: We used a grid of these simulations computed with the CO5BOLD code to model the outer layers of solar-like stars. Upper layers of the corresponding 1D standard models were then replaced by the layers obtained from the horizontally averaged 3D models. The frequency differences between these patched models and the 1D standard models were then calculated using the adiabatic approximation and allowed us to constrain the Kjeldsen et al. power law, as well as a Lorentzian formulation. Results: We find that the surface effects on modal frequencies depend significantly on both the effective temperature and the surface gravity. We further provide the variation in the parameters related to the surface-effect corrections using their power law as well as a Lorentzian formulation. Scaling relations between these parameters and the elevation (related to the Mach number) is also provided. The Lorentzian formulation is shown to be more robust for the whole frequency spectrum, while the power law is not suitable for the frequency shifts in the frequency range above νmax. Finally, we show that, owing to turbulent pressure, the elevation of the uppermost layers modifies the location of the hydrogen ionization zone and
NASA Astrophysics Data System (ADS)
Angerer, Andreas; Astner, Thomas; Wirtitsch, Daniel; Sumiya, Hitoshi; Onoda, Shinobu; Isoya, Junichi; Putz, Stefan; Majer, Johannes
2016-07-01
We design and implement 3D-lumped element microwave cavities that spatially focus magnetic fields to a small mode volume. They allow coherent and uniform coupling to electron spins hosted by nitrogen vacancy centers in diamond. We achieve large homogeneous single spin coupling rates, with an enhancement of more than one order of magnitude compared to standard 3D cavities with a fundamental resonance at 3 GHz. Finite element simulations confirm that the magnetic field distribution is homogeneous throughout the entire sample volume, with a root mean square deviation of 1.54%. With a sample containing 1017 nitrogen vacancy electron spins, we achieve a collective coupling strength of Ω = 12 MHz, a cooperativity factor C = 27, and clearly enter the strong coupling regime. This allows to interface a macroscopic spin ensemble with microwave circuits, and the homogeneous Rabi frequency paves the way to manipulate the full ensemble population in a coherent way.
The Transient 3-D Transport Coupled Code TORT-TD/ATTICA3D for High-Fidelity Pebble-Bed HTGR Analyses
NASA Astrophysics Data System (ADS)
Seubert, Armin; Sureda, Antonio; Lapins, Janis; Bader, Johannes; Laurien, Eckart
2012-01-01
This article describes the 3D discrete ordinates-based coupled code system TORT-TD/ATTICA3D that aims at steady state and transient analyses of pebble-bed high-temperature gas cooled reactors. In view of increasing computing power, the application of time-dependent neutron transport methods becomes feasible for best estimate evaluations of safety margins. The calculation capabilities of TORT-TD/ATTICA3D are presented along with the coupling approach, with focus on the time-dependent neutron transport features of TORT-TD. Results obtained for the OECD/NEA/NSC PBMR-400 benchmark demonstrate the transient capabilities of TORT-TD/ATTICA3D.
NASA Astrophysics Data System (ADS)
Kempka, Thomas; Nakaten, Benjamin; De Lucia, Marco; Nakaten, Natalie; Otto, Christopher; Pohl, Maik; Tillner, Elena; Kühn, Michael
2016-04-01
Utilization of the geological subsurface for production and storage of hydrocarbons, chemical energy and heat as well as for waste disposal requires the quantification and mitigation of environmental impacts as well as the improvement of georesources utilization in terms of efficiency and sustainability. The development of tools for coupled process simulations is essential to tackle these challenges, since reliable assessments are only feasible by integrative numerical computations. Coupled processes at reservoir to regional scale determine the behaviour of reservoirs, faults and caprocks, generally demanding for complex 3D geological models to be considered besides available monitoring and experimenting data in coupled numerical simulations. We have been developing a flexible numerical simulation framework that provides efficient workflows for integrating the required data and software packages to carry out coupled process simulations considering, e.g., multiphase fluid flow, geomechanics, geochemistry and heat. Simulation results are stored in structured data formats to allow for an integrated 3D visualization and result interpretation as well as data archiving and its provision to collaborators. The main benefits in using the flexible simulation framework are the integration of data geological and grid data from any third party software package as well as data export to generic 3D visualization tools and archiving formats. The coupling of the required process simulators in time and space is feasible, while different spatial dimensions in the coupled simulations can be integrated, e.g., 0D batch with 3D dynamic simulations. User interaction is established via high-level programming languages, while computational efficiency is achieved by using low-level programming languages. We present three case studies on the assessment of geological subsurface utilization based on different process coupling approaches and numerical simulations.
Flux limiters in the coupling of radiation and hydrodynamic models
NASA Astrophysics Data System (ADS)
Seaid, M.; Klar, A.; Dubroca, B.
2004-07-01
Two numerical approximations to radiative heat transfer problem based on asymptotic and entropy approaches are proposed for hydrodynamics radiation coupling. We compare the radiative fluxes between the two approaches and we show that the coupling based on the entropy approach is flux limited, while the other approach does not preserve this condition. Relaxation schemes are considered for the hydrodynamic part, and an iterative procedure is used for radiation. The new splitting algorithm avoids the use of Riemann solvers and Newton iterations. Numerical examples are carried out on two and three dimensional problems.
Effects of the LBV Primary's Mass-loss Rate on the 3D Hydrodynamics of eta Carinae's Colliding Winds
NASA Technical Reports Server (NTRS)
Madura, Thomas I.; Gull, Theodore R.; Cocoran, M.; Okazaki, A.; Owocki, S.; Russell, C.; Hamaguchi, K.; Clementel, N; Groh, J.; Hillier, D. J.
2013-01-01
At the heart of eta Carinae's spectacular "Homunculus" nebula lies an extremely luminous (L(sub Total) greater than approximately 5 × 10(exp 6) solar luminosity) colliding wind binary with a highly eccentric (e approximately 0.9), 5.54-year orbit (Figure 1). The primary of the system, a Luminous Blue Variable (LBV), is our closest (D approximately 2.3 kpc) and best example of a pre-hypernova or pre-gamma ray burst environment. The remarkably consistent and periodic RXTE X-ray light curve surprisingly showed a major change during the system's last periastron in 2009, with the X-ray minimum being approximately 50% shorter than the minima of the previous two cycles1. Between 1998 and 2011, the strengths of various broad stellar wind emission lines (e.g. Halpha, Fe II) in line-of-sight (l.o.s.) also decreased by factors of 1.5 - 3 relative to the continuum2. The current interpretation for these changes is that they are due to a gradual factor of 2 - 4 drop in the primary's mass-loss rate over the last approximately 15 years1, 2. However, while a secular change is seen for a direct view of the central source, little to no change is seen in profiles at high stellar latitudes or reflected off of the dense, circumbinary material known as the "Weigelt blobs"2, 3. Moreover, model spectra generated with CMFGEN predict that a factor of 2 - 4 drop in the primary's mass-loss rate should lead to huge changes in the observed spectrum, which thus far have not been seen. Here we present results from large- (plus or minus 1620 AU) and small- (plus or minus 162 AU) domain, full 3D smoothed particle hydrodynamics (SPH) simulations of eta Car's massive binary colliding winds for three different primary-star mass-loss rates (2.4, 4.8, and 8.5 × 10(exp -4) solar mass/yr). The goal is to investigate how the mass-loss rate affects the 3D geometry and dynamics of eta Car's optically-thick wind and spatially-extended wind-wind collision (WWC) regions, both of which are known sources of
Hydrodynamics of Normal Atomic Gases with Spin-orbit Coupling
Hou, Yan-Hua; Yu, Zhenhua
2015-01-01
Successful realization of spin-orbit coupling in atomic gases by the NIST scheme opens the prospect of studying the effects of spin-orbit coupling on many-body physics in an unprecedentedly controllable way. Here we derive the linearized hydrodynamic equations for the normal atomic gases of the spin-orbit coupling by the NIST scheme with zero detuning. We show that the hydrodynamics of the system crucially depends on the momentum susceptibilities which can be modified by the spin-orbit coupling. We reveal the effects of the spin-orbit coupling on the sound velocities and the dipole mode frequency of the gases by applying our formalism to the ideal Fermi gas. We also discuss the generalization of our results to other situations. PMID:26483090
Quantification of avian embryonic cardiac outflow hemodynamics through 3D-0D coupling
NASA Astrophysics Data System (ADS)
Lindsey, Stephanie; Vignon-Clementel, Irene; Butcher, Jonathan
2014-11-01
Outflow malformations account for over 20% of CHDs in the US. While the etiology of these malformations is poorly understood, most can be traced back to perturbations in the patterning of the pharyngeal arch arteries (PAAs), the precursors to the great vessels. Here, we examine the effects of normal and aberrant PAA flow, through the use of two computational models. A 0D electric analog model allows for rapid computation and global tuning of the embryo's vasculature relative to the arches. A second 3D-0D model replaces the electric analog representation of the arches with a 3D reconstruction, thereby leading to more extensive pressure and flow characterization. We obtain 3D arch artery reconstructions from nano-CT stacks and couple them to 0D outlets. In contrast to standard boundary conditions, such coupling maintains the physiologically desired cranial-caudal flow split in control embryos and predicts how this will change with vessel occlusion. We use flow inputs from Doppler velocity tracings to compute 0D and 3D-0D pulsatile hemodynamic simulations in HH18 (day 3), HH24 (day 4), and HH26 (day 5) geometries. We then calculate flow distributions and wall shear stress maps for control embryos. From here, we modify HH18 geometries to simulate varying levels of PAA occlusion. Pulsatile simulations are run in each geometry and results compared to that of controls. Results serve as a basis for examining flow-mediated growth and adaptation in cardiac outflow morphogenesis.
Laplace-domain waveform modeling and inversion for the 3D acoustic-elastic coupled media
NASA Astrophysics Data System (ADS)
Shin, Jungkyun; Shin, Changsoo; Calandra, Henri
2016-06-01
Laplace-domain waveform inversion reconstructs long-wavelength subsurface models by using the zero-frequency component of damped seismic signals. Despite the computational advantages of Laplace-domain waveform inversion over conventional frequency-domain waveform inversion, an acoustic assumption and an iterative matrix solver have been used to invert 3D marine datasets to mitigate the intensive computing cost. In this study, we develop a Laplace-domain waveform modeling and inversion algorithm for 3D acoustic-elastic coupled media by using a parallel sparse direct solver library (MUltifrontal Massively Parallel Solver, MUMPS). We precisely simulate a real marine environment by coupling the 3D acoustic and elastic wave equations with the proper boundary condition at the fluid-solid interface. In addition, we can extract the elastic properties of the Earth below the sea bottom from the recorded acoustic pressure datasets. As a matrix solver, the parallel sparse direct solver is used to factorize the non-symmetric impedance matrix in a distributed memory architecture and rapidly solve the wave field for a number of shots by using the lower and upper matrix factors. Using both synthetic datasets and real datasets obtained by a 3D wide azimuth survey, the long-wavelength component of the P-wave and S-wave velocity models is reconstructed and the proposed modeling and inversion algorithm are verified. A cluster of 80 CPU cores is used for this study.
Ui, Atsushi; Miyaji, Takamasa
2004-10-15
The best-estimate coupled three-dimensional (3-D) core and thermal-hydraulic code system TRAC-BF1/COS3D has been developed. COS3D, based on a modified one-group neutronic model, is a 3-D core simulator used for licensing analyses and core management of commercial boiling water reactor (BWR) plants in Japan. TRAC-BF1 is a plant simulator based on a two-fluid model. TRAC-BF1/COS3D is a coupled system of both codes, which are connected using a parallel computing tool. This code system was applied to the OECD/NRC BWR Turbine Trip Benchmark. Since the two-group cross-section tables are provided by the benchmark team, COS3D was modified to apply to this specification. Three best-estimate scenarios and four hypothetical scenarios were calculated using this code system. In the best-estimate scenario, the predicted core power with TRAC-BF1/COS3D is slightly underestimated compared with the measured data. The reason seems to be a slight difference in the core boundary conditions, that is, pressure changes and the core inlet flow distribution, because the peak in this analysis is sensitive to them. However, the results of this benchmark analysis show that TRAC-BF1/COS3D gives good precision for the prediction of the actual BWR transient behavior on the whole. Furthermore, the results with the modified one-group model and the two-group model were compared to verify the application of the modified one-group model to this benchmark. This comparison shows that the results of the modified one-group model are appropriate and sufficiently precise.
RF study and 3-D simulations of a side-coupling thermionic RF-gun
NASA Astrophysics Data System (ADS)
Rimjaem, S.; Kusoljariyakul, K.; Thongbai, C.
2014-02-01
A thermionic RF-gun for generating ultra-short electron bunches was optimized, developed and used as a source at a linac-based THz radiation research laboratory of the Plasma and Beam Physics Research Facility, Chiang Mai University, Thailand. The RF-gun is a π/2-mode standing wave structure, which consists of two S-band accelerating cells and a side-coupling cavity. The 2856 MHz RF wave is supplied from an S-band klystron to the gun through the waveguide input-port at the cylindrical wall of the second cell. A fraction of the RF power is coupled from the second cell to the first one via a side-coupling cavity. Both the waveguide input-port and the side-coupling cavity lead to an asymmetric geometry of the gun. RF properties and electromagnetic field distributions inside the RF-gun were studied and numerically simulated by using computer codes SUPERFISH 7.19 and CST Microwave Studio 2012©. RF characterizations and tunings of the RF-gun were performed to ensure the reliability of the gun operation. The results from 3D simulations and measurements are compared and discussed in this paper. The influence of asymmetric field distributions inside the RF-gun on the electron beam properties was investigated via 3D beam dynamics simulations. A change in the coupling-plane of the side-coupling cavity is suggested to improve the gun performance.
NASA Astrophysics Data System (ADS)
Wu, Congjun; Li, Yi; Zhou, Xiangfa
2013-03-01
We study unconventional condensations of two-component bosons in a harmonic trap subject to the 3D σ --> . p --> -type spin-orbit (SO) coupling. The topology of condensate wavefunctions manifests in the quaternionic representation. The spatial distributions of the S3 quaternionic phase exhibit 3D skyrmion configurations, while those of the S2 spin orientation possess non-zero Hopf invariants. As increasing SO coupling strength, spin textures evolve from concentric distributions to lattice structures at weak interactions. Strong interactions change condensates into spin-polarized plane-wave states, or, superpositions of two plane-waves exhibiting helical spin spirals. This project is supported by AFOSR FA9550-11-1- 0067(YIP) and NSF-DMR-1105945.
NASA Technical Reports Server (NTRS)
Madura, T. I.; Gull, T. R.; Okazaki, A. T.; Russell, C. M. P.; Owocki, S. P.; Groh, J. H.; Corcoran, M. F.; Hamaguchi, K.; Teodoro, M.
2013-01-01
Recent work suggests that the mass-loss rate of the primary star Eta-A in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present result from large- (+/- 1545 au) and small- (+/- 155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of Eta Car's colliding winds for three Eta-A mass-loss rates ( (dot-M(sub Eta-A) = 2.4, 4.8 and 8.5 × 10(exp -4) M(solar)/ yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot-M Eta-A greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star Eta-B switches from the adiabatic to the radiative-cooling regime during periastron passage (Phi approx.= 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot-M Eta-A and is caused by the encroachment of the wind of Eta-A into the acceleration zone of Eta-B's wind, plus radiative inhibition of Eta-B's wind by Eta-A. The SPH simulations together with 1D radiative transfer models of Eta-A's spectra reveal that a factor of 2 or more drop in dot-M EtaA should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot- M Eta-A was likely by a factor of approx. < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot-M Eta-A may be responsible, but changes in the wind/stellar parameter of Eta-B, while less likely, cannot yet be fully ruled out. We suggest observations during Eta-Car's next periastron in 2014 to further
NASA Astrophysics Data System (ADS)
Madura, T. I.; Gull, T. R.; Okazaki, A. T.; Russell, C. M. P.; Owocki, S. P.; Groh, J. H.; Corcoran, M. F.; Hamaguchi, K.; Teodoro, M.
2013-12-01
Recent work suggests that the mass-loss rate of the primary star ηA in the massive colliding wind binary η Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (±1545 au) and small- (±155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of η Car's colliding winds for three ηA mass-loss rates (dot{M}_{η A}= 2.4, 4.8 and 8.5 × 10-4 M⊙ yr-1), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot{M}_{η A} greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star ηB switches from the adiabatic to the radiative-cooling regime during periastron passage (φ ≈ 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot{M}_{η A} and is caused by the encroachment of the wind of ηA into the acceleration zone of ηB's wind, plus radiative inhibition of ηB's wind by ηA. The SPH simulations together with 1D radiative transfer models of ηA's spectra reveal that a factor of 2 or more drop in dot{M}_{η A} should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot{M}_{η A} was likely by a factor of ≲2 and occurred after 2004. We speculate that most of the recent observed changes in η Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot{M}_{η A} may be responsible, but changes in the wind/stellar parameters of ηB, while less likely, cannot yet be fully ruled out. We suggest observations during η Car's next periastron in 2014 to further test for decreases in dot{M}_{η A}. If
Nonlinear Generalized Hydrodynamic Wave Equations in Strongly Coupled Dusty Plasmas
Veeresha, B. M.; Sen, A.; Kaw, P. K.
2008-09-07
A set of nonlinear equations for the study of low frequency waves in a strongly coupled dusty plasma medium is derived using the phenomenological generalized hydrodynamic (GH) model and is used to study the modulational stability of dust acoustic waves to parallel perturbations. Dust compressibility contributions arising from strong Coulomb coupling effects are found to introduce significant modifications in the threshold and range of the instability domain.
3D neutronic codes coupled with thermal-hydraulic system codes for PWR, and BWR and VVER reactors
Langenbuch, S.; Velkov, K.; Lizorkin, M.
1997-07-01
This paper describes the objectives of code development for coupling 3D neutronics codes with thermal-hydraulic system codes. The present status of coupling ATHLET with three 3D neutronics codes for VVER- and LWR-reactors is presented. After describing the basic features of the 3D neutronic codes BIPR-8 from Kurchatov-Institute, DYN3D from Research Center Rossendorf and QUABOX/CUBBOX from GRS, first applications of coupled codes for different transient and accident scenarios are presented. The need of further investigations is discussed.
3D Hydrodynamic Simulations with Yguazú-A Code to Model a Jet in a Galaxy Cluster
NASA Astrophysics Data System (ADS)
Haro-Corzo, S. A. R.; Velazquez, P.; Diaz, A.
2009-05-01
We present preliminary results for a galaxy's jet expanding into an intra-cluster medium (ICM). We attempt to model the jet-gas interaction and the evolution of a extragalactic collimated jet placed at center of computational grid, which it is modeled as a cylinder ejecting gas in the z-axis direction with fixed velocity. It has precession motion around z-axis (period of 10^5 sec.) and orbital motion in XY-plane (period of 500 yr.). This jet is embedded in the ICM, which is modeled as surrounding wind in the XZ plane. We carried out 3D hydrodynamical simulations using Yguazú-A code. This simulation do not include radiative losses. In order to compare the numerical results with observations, we generated synthetic X-ray emission images. X-ray observations with high-resolution of rich cluster of galaxies show diffuse emission with filamentary structure (sometimes called as cooling flow or X-ray filament). Radio observations show a jet-like emission of the central region of the cluster. Joining these observations, in this work we explore the possibility that the jet-ambient gas interaction leads to a filamentary morphology in the X-ray domain. We have found that simulation considering orbital motion offers the possibility to explain the diffuse emission observed in the X-ray domain. The circular orbital motion, additional to precession motion, contribute to disperse the shocked gas and the X-ray appearance of the 3D simulation reproduce some important details of Abel 1795 X-ray emission (Rodriguez-Martinez et al. 2006, A&A, 448, 15): A bright bow-shock at north (spot), where interact directly the jet and the ICM and which is observed in the X-ray image. Meanwhile, in the south side there is no bow-shock X-ray emission, but the wake appears as a X-ray source. This wake is part of the diffuse shocked ambient gas region.
Modelling of river plume dynamics in Öre estuary (Baltic Sea) with Telemac-3D hydrodynamic model
NASA Astrophysics Data System (ADS)
Sokolov, Alexander
2016-04-01
The main property of river plumes is their buoyancy, fresh water discharged by rivers is less dense than the receiving, saline waters. To study the processes of plume formation in case of river discharge into a brackish estuary where salinity is low (3.5 - 5 psu) a three dimensional hydrodynamic model was applied to the Öre estuary in the Baltic Sea. This estuary is a small fjord-like bay in the north part of the Baltic Sea. Size of the bay is about 8 by 8 km with maximum depth of 35 metres. River Öre has a small average freshwater discharge of 35 m3/s. But in spring during snowmelt the discharge can be many times higher. For example, in April 2015 the discharge increased from 8 m3/s to 160 m3/s in 18 days. To study river plume dynamics a finite element based three dimensional baroclinic model TELEMAC - 3D is used. The TELEMAC modelling suite is developed by the National Laboratory of Hydraulics and Environment (LNHE) of Electricité de France (EDF). Modelling domain was approximated by an unstructured mesh with element size varies from 50 to 500 m. In vertical direction a sigma-coordinate with 20 layers was used. Open sea boundary conditions were obtained from the Baltic Sea model HIROMB-BOOS using COPERNICUS marine environment monitoring service. Comparison of modelling results with observations obtained by BONUS COCOA project's field campaign in Öre estuary in 2015 shows that the model plausible simulate river plume dynamics. Modelling of age of freshwater is also discussed. This work resulted from the BONUS COCOA project was supported by BONUS (Art 185), funded jointly by the EU and the Swedish Research Council Formas.
User Guide for the R5EXEC Coupling Interface in the RELAP5-3D Code
Forsmann, J. Hope; Weaver, Walter L.
2015-04-01
This report describes the R5EXEC coupling interface in the RELAP5-3D computer code from the users perspective. The information in the report is intended for users who want to couple RELAP5-3D to other thermal-hydraulic, neutron kinetics, or control system simulation codes.
Isotropization and Hydrodynamization in Weakly Coupled Heavy-Ion Collisions.
Kurkela, Aleksi; Zhu, Yan
2015-10-30
We numerically solve the (2+1)-dimensional effective kinetic theory of weak coupling QCD under longitudinal expansion, relevant for early stages of heavy-ion collisions. We find agreement with viscous hydrodynamics and classical Yang-Mills simulations in the regimes where they are applicable. By choosing initial conditions that are motivated by a color-glass-condensate framework, we find that for Q_{s}=2 GeV and α_{s}=0.3 the system is approximately described by viscous hydrodynamics well before τ≲1.0 fm/c. PMID:26565462
3D modelling of coupled mass and heat transfer of a convection-oven roasting process.
Feyissa, Aberham Hailu; Gernaey, Krist V; Adler-Nissen, Jens
2013-04-01
A 3D mathematical model of coupled heat and mass transfer describing oven roasting of meat has been developed from first principles. The proposed mechanism for the mass transfer of water is modified and based on a critical literature review of the effect of heat on meat. The model equations are based on a conservation of mass and energy, coupled through Darcy's equations of porous media - the water flow is mainly pressure-driven. The developed model together with theoretical and experimental assessments were used to explain the heat and water transport and the effect of the change in microstructure (permeability, water binding capacity and elastic modulus) that occur during the meat roasting process. The developed coupled partial differential equations were solved by using COMSOL Multiphysics®3.5 and state variables are predicted as functions of both position and time. The proposed mechanism was partially validated by experiments in a convection oven where temperatures were measured online. PMID:23305831
Radiation Coupling with the FUN3D Unstructured-Grid CFD Code
NASA Technical Reports Server (NTRS)
Wood, William A.
2012-01-01
The HARA radiation code is fully-coupled to the FUN3D unstructured-grid CFD code for the purpose of simulating high-energy hypersonic flows. The radiation energy source terms and surface heat transfer, under the tangent slab approximation, are included within the fluid dynamic ow solver. The Fire II flight test, at the Mach-31 1643-second trajectory point, is used as a demonstration case. Comparisons are made with an existing structured-grid capability, the LAURA/HARA coupling. The radiative surface heat transfer rates from the present approach match the benchmark values within 6%. Although radiation coupling is the focus of the present work, convective surface heat transfer rates are also reported, and are seen to vary depending upon the choice of mesh connectivity and FUN3D ux reconstruction algorithm. On a tetrahedral-element mesh the convective heating matches the benchmark at the stagnation point, but under-predicts by 15% on the Fire II shoulder. Conversely, on a mixed-element mesh the convective heating over-predicts at the stagnation point by 20%, but matches the benchmark away from the stagnation region.
Automated design of coupled RF cavities using 2-D and 3-D codes
Smith, Peter; Christiansen, D. W.; Greninger, P. T.; Spalek, G.
2001-01-01
Coupled RF cavities in the Accelerator Production of Tritium Project have been designed using a procedure in which a 2-D code (CCT) searches for a design that meets frequency and coupling requirements, while a 3-D code (HFSS) is used to obtain empirical factors used by CCT to characterize the coupling slot between cavities. Using assumed values of the empirical factors, CCT runs the Superfish code iteratively to solve for a trial cavity design that has a specified frequency and coupling. The frequency shifts and the coupling constant k of the slot are modeled in CCT using a perturbation theory, the results of which are adjusted using the empirical factors. Given a trial design, HFSS is run using periodic boundary conditions to obtain a mode spectrum. The mode spectrum is processed using the DISPER code to obtain values of the coupling and the frequencies with slots. These results are used to calculate a new set of empirical factors, which are fed back into CCT for another design iteration. Cold models have been fabricated and tested to validate the codes, and results will be presented.
Coupling frontal photopolymerization and surface instabilities for a novel 3D patterning technology
NASA Astrophysics Data System (ADS)
Vitale, Alessandra; Hennessy, Matthew; Matar, Omar; Douglas, Jack; Cabral, João
2015-03-01
Patterning of soft matter provides an exceptional route for the generation of micro/nanostructured and functional surfaces. We describe a new 3D fabrication process based on coupling frontal photopolymerization (FPP) with precisely controlled, yet spontaneous, interfacial wrinkling. FPP is a complex spatio-temporal process that can lead to well-defined propagating fronts of network formation, both stable and unstable. We investigate this process focusing on the interfacial monomer-to-polymer conversion profile and its wave propagation. A simple coarse-grained model is found to describe remarkably well the planar frontal logarithmic kinetics, capturing the effects of UV light exposure time (or dose) and temperature, as well as the front position. In defined conditions, surface instabilities are introduced and interfere with wave planarity, resulting in the formation of ``minimal'' surfaces with complex 3D geometries. Building on this understanding on the propagation of wavefronts of network formation during photopolymerization, we demonstrate the design and fabrication of 3D patterned polymer materials with tunable shapes with optical and surface functionality.
Dynamic regimes of hydrodynamically coupled self-propelling particles
NASA Astrophysics Data System (ADS)
Llopis, I.; Pagonabarraga, I.
2006-09-01
We analyze the collective dynamics of self-propelling particles (spps) which move at small Reynolds numbers including the hydrodynamic coupling to the suspending solvent through numerical simulations. The velocity distribution functions show marked deviations from Gaussian behavior at short times, and the mean-square displacement at long times shows a transition from diffusive to ballistic motion for appropriate driving mechanism at low concentrations. We discuss the structures the spps form at long times and how they correlate to their dynamic behavior.
Room temperature strong light-matter coupling in 3D THz meta-atoms (Conference Presentation)
NASA Astrophysics Data System (ADS)
Paulillo, Bruno; Manceau, Jean-Michel; Li, Lianhe; Linfield, Edmund; Colombelli, Raffaele
2016-04-01
We demonstrate strong light-matter coupling at room temperature in the terahertz (THz) spectral region using 3D meta-atoms with extremely sub-wavelength volumes. Using an air-bridge fabrication scheme, we have implemented sub-wavelength 3D THz micro-resonators that rely on suspended loop antennas connected to semiconductor-filled patch cavities. We have experimentally shown that they possess the functionalities of lumped LC resonators: their frequency response can be adjusted by independently tuning the inductance associated the antenna element or the capacitance provided by the metal-semiconductor-metal cavity. Moreover, the radiation coupling and efficiency can be engineered acting on the design of the loop antenna, similarly to conventional RF antennas. Here we take advantage of this rich playground in the context of cavity electrodynamics/intersubband polaritonics. In the strong light-matter coupling regime, a cavity and a two-level system exchange energy coherently at a characteristic rate called the vacuum Rabi frequency ΩR which is dominant with respect to all other loss mechanisms involved. The signature, in the frequency domain, is the appearance of a splitting between the bare cavity and material system resonances: the new states are called upper and a lower polariton branches. So far, most experimental demonstrations of strong light-matter interaction between an intersubband transition and a deeply sub-wavelength mode in the THz or mid-infrared ranges rely on wavelength-scale or larger resonators such as photonic crystals, diffractive gratings, dielectric micro-cavities or patch cavities. Lately, planar metamaterials have been used to enhance the light-matter interaction and strongly reduce the interaction volume by engineering the electric and magnetic resonances of the individual subwavelength constituents. In this contribution we provide evidence of strong coupling between a THz intersubband transition and an extremely sub-wavelength mode (≈λ/10
NASA Astrophysics Data System (ADS)
Gorby, M.; Schwadron, N.; Torok, T.; Downs, C.; Lionello, R.; Linker, J.; Titov, V. S.; Mikic, Z.; Riley, P.; Desai, M. I.; Dayeh, M. A.
2014-12-01
Recent work on the coupling between the Energetic Particle Radiation Environment Module (EPREM, a 3D energetic particle model) and Magnetohydrodynamics Around a Sphere (MAS, an MHD code developed at Predictive Science, Inc.) has demonstrated the efficacy of compression regions around fast coronal mass ejections (CMEs) for particle acceleration low in the corona (˜ 3 - 6 solar radii). These couplings show rapid particle acceleration over a broad longitudinal extent (˜ 80 degrees) resulting from the pile-up of magnetic flux in the compression regions and their subsequent expansion. The challenge for forming large SEP events in such compression-acceleration scenarios is to have enhanced scattering within the acceleration region while also allowing for efficient escape of accelerated particles downstream (away from the Sun) from the compression region. We present here the most recent simulation results including energetic particle and CME plasma profiles, the subsequent flux and dosages at 1AU, and an analysis of the compressional regions as efficient accelerators.
Linearly perturbed MHD equilibria and 3D eddy current coupling via the control surface method
NASA Astrophysics Data System (ADS)
Portone, A.; Villone, F.; Liu, Y.; Albanese, R.; Rubinacci, G.
2008-08-01
In this paper, a coupling strategy based on the control surface concept is used to self-consistently couple linear MHD solvers to 3D codes for the eddy current computation of eddy currents in the metallic structures surrounding the plasma. The coupling is performed by assuming that the plasma inertia (and, with it, all Alfven wave-like phenomena) can be neglected on the time scale of interest, which is dictated by the relevant electromagnetic time of the metallic structures. As is shown, plasma coupling with the metallic structures results in perturbations to the inductance matrix operator. In particular, by adopting the Fourier decomposition in poloidal and toroidal modes, it turns out that each toroidal mode can be associated with a matrix (additively) perturbing the inductance matrix that commonly describes the magnetic coupling of currents in vacuum. In this way, the treatment of resistive wall modes instabilities of various toroidal mode numbers and their possible cross-talk through the currents induced in the metallic structures can be easily studied.
Seo, Sujin; Zhou, Xiangfei; Liu, Gang Logan
2016-07-01
Plasmonic substrates have fixed sensitivity once the geometry of the structure is defined. In order to improve the sensitivity, significant research effort has been focused on designing new plasmonic structures, which involves high fabrication costs; however, a method is reported for improving sensitivity not by redesigning the structure but by simply assembling plasmonic nanoparticles (NPs) near the evanescent field of the underlying 3D plasmonic nanostructure. Here, a nanoscale Lycurgus cup array (nanoLCA) is employed as a base colorimetric plasmonic substrate and an assembly template. Compared to the nanoLCA, the NP assembled nanoLCA (NP-nanoLCA) exhibits much higher sensitivity for both bulk refractive index sensing and biotin-streptavidin binding detection. The limit of detection of the NP-nanoLCA is at least ten times smaller when detecting biotin-streptavidin conjugation. The numerical calculations confirm the importance of the additive plasmon coupling between the NPs and the nanoLCA for a denser and stronger electric field in the same 3D volumetric space. Tunable sensitivity is accomplished by controlling the number of NPs in each nanocup, or the number density of the hot spots. This simple yet scalable and cost-effective method of using additive heterogeneous plasmon coupling effects will benefit various chemical, medical, and environmental plasmon-based sensors. PMID:27206214
A 3D coupled hydro-mechanical granular model for the prediction of hot tearing formation
NASA Astrophysics Data System (ADS)
Sistaninia, M.; Phillion, A. B.; Drezet, J.-M.; Rappaz, M.
2012-07-01
A new 3D coupled hydro-mechanical granular model that simulates hot tearing formation in metallic alloys is presented. The hydro-mechanical model consists of four separate 3D modules. (I) The Solidification Module (SM) is used for generating the initial solid-liquid geometry. Based on a Voronoi tessellation of randomly distributed nucleation centers, this module computes solidification within each polyhedron using a finite element based solute diffusion calculation for each element within the tessellation. (II) The Fluid Flow Module (FFM) calculates the solidification shrinkage and deformation-induced pressure drop within the intergranular liquid. (III) The Semi-solid Deformation Module (SDM) is used to simulate deformation of the granular structure via a combined finite element / discrete element method. In this module, deformation of the solid grains is modeled using an elasto-viscoplastic constitutive law. (IV) The Failure Module (FM) is used to simulate crack initiation and propagation with the fracture criterion estimated from the overpressure required to overcome the capillary forces at the liquid-gas interface. The FFM, SDM, and FM are coupled processes since solid deformation, intergranular flow, and crack initiation are deeply linked together. The granular model predictions have been validated against bulk data measured experimentally and calculated with averaging techniques.
McConnell, Amber C.; Fishman, Randy Scott; Miller, Joel S.
2012-01-01
Mean field expressions based on the simple Heisenberg model were derived to correlate the inter- and intralayer exchange coupling to the critical temperatures, Tc, for several TCNE (tetracyanoethylene) based magnets with extended 2- and 3-D structure types. These expressions were used to estimate the exchange coupling, J, for 2-D ferrimagnetic [MII(TCNE)(NCMe)2]+ (M = Mn, Fe), 3-D antiferromagnetic MnII(TCNE)[C4(CN)8]1/2, and 3-D ferrimagnetic MnII(TCNE)3/2(I3)1/2. The sign and magnitude of the exchange coupling are in accord with previously reported magnetic data.
Studies of coupled cavity LINAC (CCL) accelerating structures with 3-D codes
Spalek, G.
2000-08-01
The cw CCL being designed for the Accelerator Production of Tritium (APT) project accelerates protons from 96 MeV to 211 MeV. It consists of 99 segments each containing up to seven accelerating cavities. Segments are coupled by intersegment coupling cavities and grouped into supermodules. The design method needs to address not only basic cavity sizing for a given coupling and pi/2 mode frequency, but also the effects of high power densities on the cavity frequency, mechanical stresses, and the structure's stop band during operation. On the APT project, 3-D RF (Ansoft Corp.'s HFSS) and coupled RF/structural (Ansys Inc.'s ANSYS) codes are being used. to develop tools to address the above issues and guide cooling channel design. The code's predictions are being checked against available low power Aluminum models. Stop band behavior under power will be checked once the tools are extended to CCDTL structures that have been tested at high power. A summary of calculations made to date and agreement with measured results will be presented.
NASA Astrophysics Data System (ADS)
Sirjacobs, D.; Grégoire, M.; Delhez, E.; Nihoul, J.
2003-04-01
Within the context of the EU INCO-COPERNICUS program "Desertification in the Aral Sea Region: A study of the Natural and Anthropogenic Impacts" (Contract IAC2-CT-2000-10023), a large-scale 3D hydrodynamic model was adapted to address specifically the macroscale processes affecting the Aral Sea water circulation and ventilation. The particular goal of this research is to simulate the effect of lasting negative water balance on the 3D seasonal circulation, temperature, salinity and water-mixing fields of the Aral Sea. The original Aral Sea seasonal hydrodynamism is simulated with the average seasonal forcings corresponding to the period from 1956 to 1960. This first investigation concerns a period of relative stability of the water balance, before the beginning of the drying process. The consequences of the drying process on the hydrodynamic of the Sea will be studied by comparing this first results with the simulation representing the average situation for the years 1981 to 1985, a very low river flow period. For both simulation periods, the forcing considered are the seasonal fluctuations of wind fields, precipitation, evaporation, river discharge and salinity, cloud cover, air temperature and humidity. The meteorological forcings were adapted to the common optimum one-month temporal resolution of the available data sets. Monthly mean kinetic energy flux and surface tensions were calculated from daily ECMWF wind data. Monthly in situ precipitation, surface air temperature and humidity fields were interpolated from data obtained from the Russian Hydrological and Meteorological Institute. Monthly water discharge and average salinity of the river water were considered for both Amu Darya and Syr Darya river over each simulation periods. The water mass conservation routines allowed the simulation of a changing coastline by taking into account local drying and flooding events of particular grid points. Preliminary barotropic runs were realised (for the 1951
NASA Astrophysics Data System (ADS)
Baraka, S. M.; Ben-Jaffel, L. B.
2014-12-01
We use particle-in-cell PIC 3D Electromagnetic, relativistic global code to address large-scale problems in magnetosphere electrodynamics. Terrestrial bow shock is simulated as an example. 3D Magnetohydrodynamics model ,MHD GUMICS in CCMC project, have been used in parallel with PIC under same scaled Solar wind (SW) and IMF conditions. We report new results from the coupling between the two models. Further investigations are required for confirmations of these results. In both codes the Earth's bow shock position is found at ~14.8 RE along the Sun-Earth line, and ~29 RE on the dusk side which is consistent with past in situ observation. Both simulations reproduce the theoretical jump conditions at the shock. However, PIC code density and temperature distributions are inflated and slightly shifted sunward when compared to MHD results. Reflected ions upstream of the bow shock may cause this sunward shift for density and temperature. Distribution of reflected ions and electrons are shown in the foreshock region, within the transition of the shock and in the downstream. The current version of PIC code can be run under modest computing facilities and resources. Additionally, existing MHD simulations should be useful to calibrate scaled properties of plasma resulting from PIC simulations for comparison with observations. Similarities and drawbacks of the results obtained by the two models are listed. The ultimate goal of using these different models in a complimentary manner rather than competitive is to better understand the macrostructure of the magnetosphere
Finite element modeling of a 3D coupled foot-boot model.
Qiu, Tian-Xia; Teo, Ee-Chon; Yan, Ya-Bo; Lei, Wei
2011-12-01
Increasingly, musculoskeletal models of the human body are used as powerful tools to study biological structures. The lower limb, and in particular the foot, is of interest because it is the primary physical interaction between the body and the environment during locomotion. The goal of this paper is to adopt the finite element (FE) modeling and analysis approaches to create a state-of-the-art 3D coupled foot-boot model for future studies on biomechanical investigation of stress injury mechanism, foot wear design and parachute landing fall simulation. In the modeling process, the foot-ankle model with lower leg was developed based on Computed Tomography (CT) images using ScanIP, Surfacer and ANSYS. Then, the boot was represented by assembling the FE models of upper, insole, midsole and outsole built based on the FE model of the foot-ankle, and finally the coupled foot-boot model was generated by putting together the models of the lower limb and boot. In this study, the FE model of foot and ankle was validated during balance standing. There was a good agreement in the overall patterns of predicted and measured plantar pressure distribution published in literature. The coupled foot-boot model will be fully validated in the subsequent works under both static and dynamic loading conditions for further studies on injuries investigation in military and sports, foot wear design and characteristics of parachute landing impact in military. PMID:21676642
Automatic segmentation of bladder and prostate using coupled 3D deformable models.
Costa, María Jimena; Delingette, Hervé; Novellas, Sébastien; Ayache, Nicholas
2007-01-01
In this paper, we propose a fully automatic method for the coupled 3D localization and segmentation of lower abdomen structures. We apply it to the joint segmentation of the prostate and bladder in a database of CT scans of the lower abdomen of male patients. A flexible approach on the bladder allows the process to easily adapt to high shape variation and to intensity inhomogeneities that would be hard to characterize (due, for example, to the level of contrast agent that is present). On the other hand, a statistical shape prior is enforced on the prostate. We also propose an adaptive non-overlapping constraint that arbitrates the evolution of both structures based on the availability of strong image data at their common boundary. The method has been tested on a database of 16 volumetric images, and the validation process includes an assessment of inter-expert variability in prostate delineation, with promising results. PMID:18051066
3D position determination in monolithic crystals coupled to SiPMs for PET.
Etxebeste, Ane; Barrio, John; Muñoz, Enrique; Oliver, Josep F; Solaz, Carles; Llosá, Gabriela
2016-05-21
The interest in using continuous monolithic crystals in positron emission tomography (PET) has grown in the last years. Coupled to silicon photomultipliers (SiPMs), the detector can combine high sensitivity and high resolution, the two main factors to be maximized in a positron emission tomograph. In this work, the position determination capability of a detector comprised of a [Formula: see text] mm(3) LYSO crystal coupled to an [Formula: see text]-pixel array of SiPMs is evaluated. The 3D interaction position of γ-rays is estimated using an analytical model of the light distribution including reflections on the facets of the crystal. Monte Carlo simulations have been performed to evaluate different crystal reflectors and geometries. The method has been characterized and applied to different cases. Intrinsic resolution obtained with the position estimation method used in this work, applied to experimental data, achieves sub-millimetre resolution values. Average resolution over the detector surface for 5 mm thick crystal is ∼0.9 mm FWHM and ∼1.2 mm FWHM for 10 mm thick crystal. Depth of interaction resolution is close to 2 mm FWHM in both cases, while the FWTM is ∼5.3 mm for 5 mm thick crystal and ∼9.6 mm for 10 mm thick crystal. PMID:27119737
NASA Astrophysics Data System (ADS)
Guo, Junyan; Zhang, Yong-Wei; Narayanaswamy, Sridhar
2012-02-01
Fabrication of quantum dots (QDs) with high density may be realized by self-assembly via heteroepitaxial growth of thin films. Since the electronic and optoelectronic properties of QDs are sensitive to size, morphology, strain and especially composition, it is of great importance to control their composition profiles and morphology, and engineer the strain in them. Since the growth is a dynamic process, which carries out via surface diffusion driven primarily by strain relaxation and entropy change due to chemical intermixing, a strong coupling between morphological and composition evolutions during this process leads to a rather complex dynamics, which has not been fully understood. In this work, a 3-D finite element model is developed, which is capable of modeling the formation, self-assembly and coarsening of hetero-epitaxial alloy islands by considering the coupling of morphological and compositional evolution. Several interesting experimental observations, such as fast coarsening kinetics; asymmetries in composition profile and island shape; lateral motion of alloy islands have been observed in our simulations. Our model predictions have painted a rather complete picture for the entire dynamic evolution during the growth of nanoscale heteroepitaxial islands.
Accelerating DynEarthSol3D on tightly coupled CPU-GPU heterogeneous processors
NASA Astrophysics Data System (ADS)
Ta, Tuan; Choo, Kyoshin; Tan, Eh; Jang, Byunghyun; Choi, Eunseo
2015-06-01
DynEarthSol3D (Dynamic Earth Solver in Three Dimensions) is a flexible, open-source finite element solver that models the momentum balance and the heat transfer of elasto-visco-plastic material in the Lagrangian form using unstructured meshes. It provides a platform for the study of the long-term deformation of earth's lithosphere and various problems in civil and geotechnical engineering. However, the continuous computation and update of a very large mesh poses an intolerably high computational burden to developers and users in practice. For example, simulating a small input mesh containing around 3000 elements in 20 million time steps would take more than 10 days on a high-end desktop CPU. In this paper, we explore tightly coupled CPU-GPU heterogeneous processors to address the computing concern by leveraging their new features and developing hardware-architecture-aware optimizations. Our proposed key optimization techniques are three-fold: memory access pattern improvement, data transfer elimination and kernel launch overhead minimization. Experimental results show that our proposed implementation on a tightly coupled heterogeneous processor outperforms all other alternatives including traditional discrete GPU, quad-core CPU using OpenMP, and serial implementations by 67%, 50%, and 154% respectively even though the embedded GPU in the heterogeneous processor has significantly less number of cores than high-end discrete GPU.
3D position determination in monolithic crystals coupled to SiPMs for PET
NASA Astrophysics Data System (ADS)
Etxebeste, Ane; Barrio, John; Muñoz, Enrique; Oliver, Josep F.; Solaz, Carles; Llosá, Gabriela
2016-05-01
The interest in using continuous monolithic crystals in positron emission tomography (PET) has grown in the last years. Coupled to silicon photomultipliers (SiPMs), the detector can combine high sensitivity and high resolution, the two main factors to be maximized in a positron emission tomograph. In this work, the position determination capability of a detector comprised of a 12× 12× 10 mm3 LYSO crystal coupled to an 8× 8 -pixel array of SiPMs is evaluated. The 3D interaction position of γ-rays is estimated using an analytical model of the light distribution including reflections on the facets of the crystal. Monte Carlo simulations have been performed to evaluate different crystal reflectors and geometries. The method has been characterized and applied to different cases. Intrinsic resolution obtained with the position estimation method used in this work, applied to experimental data, achieves sub-millimetre resolution values. Average resolution over the detector surface for 5 mm thick crystal is ∼0.9 mm FWHM and ∼1.2 mm FWHM for 10 mm thick crystal. Depth of interaction resolution is close to 2 mm FWHM in both cases, while the FWTM is ∼5.3 mm for 5 mm thick crystal and ∼9.6 mm for 10 mm thick crystal.
Applications of 3D hydrodynamic and particle tracking models in the San Francisco bay-delta estuary
Smith, P.E.; Donovan, J.M.; Wong, H.F.N.
2005-01-01
Three applications of three-dimensional hydrodynamic and particle-tracking models are currently underway by the United States Geological Survey in the San Francisco Bay-Delta Estuary. The first application is to the San Francisco Bay and a portion of the coastal ocean. The second application is to an important, gated control channel called the Delta Cross Channel, located within the northern portion of the Sacramento-San Joaquin River Delta. The third application is to a reach of the San Joaquin River near Stockton, California where a significant dissolved oxygen problem exists due, in part, to conditions associated with the deep-water ship channel for the Port of Stockton, California. This paper briefly discusses the hydrodynamic and particle tracking models being used and the three applications. Copyright ASCE 2005.
Morgan, F. Dale; Sogade, John
2004-12-14
This project was designed as a broad foundational study of spectral induced polarization (SIP) for characterization of contaminated sites. It encompassed laboratory studies of the effects of chemistry on induced polarization, development of 3D forward modeling and inversion codes, and investigations of inductive and capacitive coupling problems. In the laboratory part of the project a physico-chemical model developed in this project was used to invert laboratory IP spectra for the grain size and the effective grain size distribution of the sedimentary rocks as well as the formation factor, porosity, specific surface area, and the apparent fractal dimension. Furthermore, it was established that the IP response changed with the solution chemistry, the concentration of a given solution chemistry, valence of the constituent ions, and ionic radius. In the field part of the project, a 3D complex forward and inverse model was developed. It was used to process data acquired at two frequencies (1/16 Hz and 1/ 4Hz) in a cross-borehole configuration at the A-14 outfall area of the Savannah River Site (SRS) during March 2003 and June 2004. The chosen SRS site was contaminated with Tetrachloroethylene (TCE) and Trichloroethylene (PCE) that were disposed in this area for several decades till the 1980s. The imaginary conductivity produced from the inverted 2003 data correlated very well with the log10 (PCE) concentration derived from point sampling at 1 ft spacing in five ground-truth boreholes drilled after the data acquisition. The equivalent result for the 2004 data revealed that there were significant contaminant movements during the period March 2003 and June 2004, probably related to ground-truth activities and nearby remediation activities. Therefore SIP was successfully used to develop conceptual models of volume distributions of PCE/TCE contamination. In addition, the project developed non-polarizing electrodes that can be deployed in boreholes for years. A total of 28
Coupling Hydrologic and Hydrodynamic Models to Estimate PMF
NASA Astrophysics Data System (ADS)
Felder, G.; Weingartner, R.
2015-12-01
Most sophisticated probable maximum flood (PMF) estimations derive the PMF from the probable maximum precipitation (PMP) by applying deterministic hydrologic models calibrated with observed data. This method is based on the assumption that the hydrological system is stationary, meaning that the system behaviour during the calibration period or the calibration event is presumed to be the same as it is during the PMF. However, as soon as a catchment-specific threshold is reached, the system is no longer stationary. At or beyond this threshold, retention areas, new flow paths, and changing runoff processes can strongly affect downstream peak discharge. These effects can be accounted for by coupling hydrologic and hydrodynamic models, a technique that is particularly promising when the expected peak discharge may considerably exceed the observed maximum discharge. In such cases, the coupling of hydrologic and hydraulic models has the potential to significantly increase the physical plausibility of PMF estimations. This procedure ensures both that the estimated extreme peak discharge does not exceed the physical limit based on riverbed capacity and that the dampening effect of inundation processes on peak discharge is considered. Our study discusses the prospect of considering retention effects on PMF estimations by coupling hydrologic and hydrodynamic models. This method is tested by forcing PREVAH, a semi-distributed deterministic hydrological model, with randomly generated, physically plausible extreme precipitation patterns. The resulting hydrographs are then used to externally force the hydraulic model BASEMENT-ETH (riverbed in 1D, potential inundation areas in 2D). Finally, the PMF estimation results obtained using the coupled modelling approach are compared to the results obtained using ordinary hydrologic modelling.
3-D GRMHD Simulations of Disk-Jet Coupling and Emission
NASA Astrophysics Data System (ADS)
Nishikawa, K.-I.; Mizuno, Y.; Fuerst, S.; Wu, K.; Hardee, P.; Richardson, G.; Koide, S.; Shibata, K.; Kudoh, T.; Fishman, G. J.
2005-11-01
We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity ~ 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units (τS = rS/c where rS ≡ 2GM/c2) which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3rS. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfvén waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star. We have also calculated the free-free emission from a disk/outflow near a rotating black hole using our axisymmetric GRMHD simulation using a covariant radiative transfer formulation. Our calculation shows radiation from a shock, and hence the disk-jet coupling region is observable.
An example of neutronic penalizations in reactivity transient analysis using 3D coupled chain HEMERA
Dubois, F.; Normand, B.; Sargeni, A.
2012-07-01
HEMERA (Highly Evolutionary Methods for Extensive Reactor Analyses), is a fully coupled 3D computational chain developed jointly by IRSN and CEA. It is composed of CRONOS2 (core neutronics, cross sections library from APOLLO2), FLICA4 (core thermal-hydraulics) and the system code CATHARE. Multi-level and multi-dimensional models are developed to account for neutronics, core thermal-hydraulics, fuel thermal analysis and system thermal-hydraulics, dedicated to best-estimate, conservative simulations and sensitivity analysis. In IRSN, the HEMERA chain is widely used to study several types of reactivity accidents and for sensitivity studies. Just as an example of the HEMERA possibilities, we present here two types of neutronic penalizations and their impact on a power transient due to a REA (Rod Ejection Accident): in the first one, we studied a bum-up distribution modification and in the second one, a delayed-neutron fraction modification. Both modifications are applied to the whole core or localized in a few assemblies. Results show that it is possible to use global or local changes but 1) in case of bum-up modification, the total core power can increase when assembly peak power decrease so, care has to be taken if the goal is to maximize a local power peak and 2) for delayed-neutron fraction, a local modification can have the same effect as the one on the whole core, provided that it is large enough. (authors)
3D simulations of gas puff effects on edge density and ICRF coupling in ASDEX Upgrade
NASA Astrophysics Data System (ADS)
Zhang, W.; Bobkov, V.; Lunt, T.; Noterdaeme, J.-M.; Coster, D.; Bilato, R.; Jacquet, P.; Brida, D.; Feng, Y.; Wolfrum, E.; Guimarais, L.; the ASDEX Upgrade Team
2016-03-01
In recent experiments, a local gas puff was found to be an effective way to tailor the scrape-off layer (SOL) density and improve the ion cyclotron range of frequency (ICRF) power coupling in tokamaks. In order to quantitatively reproduce these experiments, to understand the corresponding physics and to optimize the gas valve positions and rates, simulations were carried out with the 3D edge plasma transport code EMC3-EIRENE in ASDEX Upgrade. An inter-ELM phase of an H-mode discharge with a moderate gas puff rate (1.2 × 1022 electrons s-1) is used in our simulations. We simulated cases with gas puff in the lower divertor, the outer mid-plane and the top of the machine while keeping other conditions the same. Compared with the lower divertor gas puff, the outer mid-plane gas puff can increase the local density in front of the antennas most effectively, while a toroidally uniform but significantly smaller enhancement is found for the top gas puff. Good agreement between our simulations and experiments is obtained. With further simulations, the mechanisms of SOL density tailoring via local gas puffing and the strategies of gas puff optimization are discussed in the paper.
Najjar, F M; Solberg, J; White, D
2008-04-17
A verification test suite has been assessed with primary focus on low reynolds number flow of liquid metals. This is representative of the interface between the armature and rail in gun applications. The computational multiphysics framework, ALE3D, is used. The main objective of the current study is to provide guidance and gain confidence in the results obtained with ALE3D. A verification test suite based on 2-D cases is proposed and includes the lid-driven cavity and the Couette flow are investigated. The hydro and thermal fields are assumed to be steady and laminar in nature. Results are compared with analytical solutions and previously published data. Mesh resolution studies are performed along with various models for the equation of state.
Hayes, J C; Norman, M
1999-10-28
This report details an investigation into the efficacy of two approaches to solving the radiation diffusion equation within a radiation hydrodynamic simulation. Because leading-edge scientific computing platforms have evolved from large single-node vector processors to parallel aggregates containing tens to thousands of individual CPU's, the ability of an algorithm to maintain high compute efficiency when distributed over a large array of nodes is critically important. The viability of an algorithm thus hinges upon the tripartite question of numerical accuracy, total time to solution, and parallel efficiency.
Zhuang, Guisheng; Jensen, Thomas G; Kutter, Jörg P
2012-07-01
In this paper, we describe a microfluidic device composed of integrated microoptical elements and a two-layer microchannel structure for highly sensitive light scattering detection of micro/submicrometer-sized particles. In the two-layer microfluidic system, a sample flow stream is first constrained in the out-of-plane direction into a narrow sheet, and then focused in-plane into a small core region, obtaining on-chip three-dimensional (3D) hydrodynamic focusing. All the microoptical elements, including waveguides, microlens, and fiber-to-waveguide couplers, and the in-plane focusing channels are fabricated in one SU-8 layer by standard photolithography. The channels for out-of-plane focusing are made in a polydimethylsiloxane (PDMS) layer by a single cast using a SU-8 master. Numerical and experimental results indicate that the device can realize 3D hydrodynamic focusing reliably over a wide range of Reynolds numbers (0.5 < Re < 20). Polystyrene particles of three sizes (2, 1, and 0.5 μm) were measured in the microfluidic device with integrated optics, demonstrating the feasibility of this approach to detect particles in the low micrometer size range by light scattering detection. PMID:22740459
A Generalized Hydrodynamics Model for Strongly Coupled Plasmas
NASA Astrophysics Data System (ADS)
Diaw, Abdourahmane; Murillo, Michael Sean
2015-11-01
Starting with the equations of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, we obtain the density, momentum and stress tensor-moment equations. The closure proceeds in two steps. The first that guarantees an equilibrium state is given by density functional theory. It ensures self consistency in the equation-of-state properties of the plasma. The second involves modifying the two-body distribution function to include collisions in the relaxation of the stress tensor. The resulting generalized hydrodynamics thus includes all impacts of Coulomb coupling, viscous damping, and the high-frequency response. We compare our results with those of several known models, including generalized hydrodynamic theory and models obtained using the Singwi-Tosi-Land-Sjolander approximation and the quasi-localized charge approximation. We find that the viscoelastic response, including both the high-frequency elastic generalization and viscous wave damping, is important for correctly describing ion-acoustic waves. We illustrate this result by considering three very different systems: ultracold plasmas, dusty plasmas, and dense plasmas. The new model is validated by comparing its results with those obtained from molecular-dynamics simulations of Yukawa plasmas, and the agreement is excellent. This work was supported by the Air Force Office of Scientific Research (Grant No. FA9550-12-1-0344).
NASA Astrophysics Data System (ADS)
Li, Y.; Acharya, K.; Chen, D.; Stone, M.; Yu, Z.; Young, M.; Zhu, J.; Shafer, D. S.; Warwick, J. J.
2009-12-01
Sustained drought in the western United States since 2000 has led to a significant drop (about 35 meters) in the water level of Lake Mead, the largest reservoir by volume in United States. The drought combined with rapid urban development in southern Nevada and emergence of invasive species has threatened the water quality and ecological processes in Lake Mead. A three-dimensional hydrodynamic model, Environmental Fluid Dynamics Code (EFDC), was applied to investigate lake circulation and temperature stratification in parts of Lake Mead (Las Vegas Bay and Boulder Basin) under changing water levels. Besides the inflow from Las Vegas Wash and the Colorado River, the model considered atmospheric changes as well as the boundary conditions restricted by the operation of Hoover Dam. The model was calibrated and verified by using observed data including water level, velocity, and temperature from 2003 and 2005. The model was applied to study the hydrodynamic processes at water level 366.8 m (year 2000) and at water level 338.2 m (year 2008). The high-stage simulation described the pre-drought lake hydrodynamic processes while the low-stage simulation highlighted the drawdown impact on such processes. The results showed that both inflow and wind-driven mixing process played major roles in the thermal stratification and lake circulation in both cases. However, the atmospheric boundary played a more important role than inflow temperature on thermal stratification of Lake Mead during water level decline. Further, the thermal stratification regime and flow circulation pattern in shallow lake regions (e.g.., the Boulder Basin area) were most impacted. The temperature of the lake at the high-stage was more sensitive to inflow temperatures than at low-stage. Furthermore, flow velocities decreased with the decreasing water level due to reduction in wind impacts, particularly in shallow areas of the lake. Such changes in temperature and lake current due to present drought have a
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
3-D GRMHD Simulations of Disk-Jet Coupling and Associated Variabilities and Emission
NASA Technical Reports Server (NTRS)
Nishikawa, K. I.
2005-01-01
We have performed the fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3 r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star. We will report initial results of emission calculated based on simulations with Kerr metric using a new ray-tracing method developed by S. Feurst and K. Wu. Images include free-free emission, absorption, and electron scattering. Fluorescent iron line emission and its variabilities from black holes will also be calculated based on simulation results.
Coupled fvGCM-GCE Modeling System, 3D Cloud-Resolving Model and Cloud Library
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo
2005-01-01
Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional singlecolumn models in simulating various types of clouds and cloud systems from Merent geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloudscale model (termed a super-parameterization or multiscale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameteridon NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. A seed fund is available at NASA Goddard to build a MMF based on the 2D Goddard cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM). A prototype MMF in being developed and production nms will be conducted at the beginning of 2005. In this talk, I will present: (1) A brief review on GCE model and its applications on precipitation processes, (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), (3) A cloud library generated by Goddard MMF, and 3D GCE model, and (4) A brief discussion on the GCE model on developing a global cloud simulator.
Coupled fvGCM-GCE Modeling System, 3D Cloud-Resolving Model and Cloud Library
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo
2005-01-01
Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud- resolving models (CRMs) agree with observations better than traditional single-column models in simulating various types of clouds and cloud systems from different geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a super-parameterization or multi-scale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. A seed fund is available at NASA Goddard to build a MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM). A prototype MMF in being developed and production runs will be conducted at the beginning of 2005. In this talk, I will present: (1) A brief review on GCE model and its applications on precipitation processes, ( 2 ) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), (3) A cloud library generated by Goddard MMF, and 3D GCE model, and (4) A brief discussion on the GCE model on developing a global cloud simulator.
3-D Surface Depression Profiling Using High Frequency Focused Air-Coupled Ultrasonic Pulses
NASA Technical Reports Server (NTRS)
Roth, Don J.; Kautz, Harold E.; Abel, Phillip B.; Whalen, Mike F.; Hendricks, J. Lynne; Bodis, James R.
1999-01-01
Surface topography is an important variable in the performance of many industrial components and is normally measured with diamond-tip profilometry over a small area or using optical scattering methods for larger area measurement. This article shows quantitative surface topography profiles as obtained using only high-frequency focused air-coupled ultrasonic pulses. The profiles were obtained using a profiling system developed by NASA Glenn Research Center and Sonix, Inc (via a formal cooperative agreement). (The air transducers are available as off-the-shelf items from several companies.) The method is simple and reproducible because it relies mainly on knowledge and constancy of the sound velocity through the air. The air transducer is scanned across the surface and sends pulses to the sample surface where they are reflected back from the surface along the same path as the incident wave. Time-of-flight images of the sample surface are acquired and converted to depth/surface profile images using the simple relation (d = V*t/2) between distance (d), time-of-flight (t), and the velocity of sound in air (V). The system has the ability to resolve surface depression variations as small as 25 microns, is useable over a 1.4 mm vertical depth range, and can profile large areas only limited by the scan limits of the particular ultrasonic system. (Best-case depth resolution is 0.25 microns which may be achievable with improved isolation from vibration and air currents.) The method using an optimized configuration is reasonably rapid and has all quantitative analysis facilities on-line including 2-D and 3-D visualization capability, extreme value filtering (for faulty data), and leveling capability. Air-coupled surface profilometry is applicable to plate-like and curved samples. In this article, results are shown for several proof-of-concept samples, plastic samples burned in microgravity on the STS-54 space shuttle mission, and a partially-coated cylindrical ceramic
A generic semi-implicit coupling methodology for use in RELAP5-3D{copyright}
Aumiller, D.L.; Tomlinson, E.T.; Weaver, W.L.
2000-09-01
A generic semi-implicit coupling methodology has been developed and implemented in the RELAP5-3D{copyright} computer program. This methodology allows RELAP5-3D{copyright} to be used with other computer programs to perform integrated analyses of nuclear power reactor systems and related experimental facilities. The coupling methodology potentially allows different programs to be used to model different portions of the system. The programs are chosen based on their capability to model the phenomena that are important in the simulation in the various portions of the system being considered. The methodology was demonstrated using a test case in which the test geometry was divided into two parts each of which was solved as a RELAP5-3D{copyright} simulation. This test problem exercised all of the semi-implicit coupling features which were installed in RELAP5-3D0. The results of this verification test case show that the semi-implicit coupling methodology produces the same answer as the simulation of the test system as a single process.
Fully-coupled analysis of jet mixing problems. Three-dimensional PNS model, SCIP3D
NASA Technical Reports Server (NTRS)
Wolf, D. E.; Sinha, N.; Dash, S. M.
1988-01-01
Numerical procedures formulated for the analysis of 3D jet mixing problems, as incorporated in the computer model, SCIP3D, are described. The overall methodology closely parallels that developed in the earlier 2D axisymmetric jet mixing model, SCIPVIS. SCIP3D integrates the 3D parabolized Navier-Stokes (PNS) jet mixing equations, cast in mapped cartesian or cylindrical coordinates, employing the explicit MacCormack Algorithm. A pressure split variant of this algorithm is employed in subsonic regions with a sublayer approximation utilized for treating the streamwise pressure component. SCIP3D contains both the ks and kW turbulence models, and employs a two component mixture approach to treat jet exhausts of arbitrary composition. Specialized grid procedures are used to adjust the grid growth in accordance with the growth of the jet, including a hybrid cartesian/cylindrical grid procedure for rectangular jets which moves the hybrid coordinate origin towards the flow origin as the jet transitions from a rectangular to circular shape. Numerous calculations are presented for rectangular mixing problems, as well as for a variety of basic unit problems exhibiting overall capabilities of SCIP3D.
Turbulence in strongly coupled dusty plasmas using generalized hydrodynamic description
Tiwari, Sanat Kumar; Dharodi, Vikram Singh; Das, Amita; Patel, Bhavesh G.; Kaw, Predhiman
2015-02-15
The properties of decaying turbulence have been studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of strongly coupled dusty plasma medium in two dimensions. The GHD model treats the strongly coupled dusty plasma system as a visco-elastic medium. The incompressible limit of the GHD model is considered here. The studies carried out here are, however, applicable to a wider class of visco-elastic systems, and are not merely confined to the dusty plasma medium. Our simulations studies show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence, where Weissenberg's number also plays an important role on the onset of turbulence. This feature has been observed in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.
HIFI-C: a robust and fast method for determining NMR couplings from adaptive 3D to 2D projections.
Cornilescu, Gabriel; Bahrami, Arash; Tonelli, Marco; Markley, John L; Eghbalnia, Hamid R
2007-08-01
We describe a novel method for the robust, rapid, and reliable determination of J couplings in multi-dimensional NMR coupling data, including small couplings from larger proteins. The method, "High-resolution Iterative Frequency Identification of Couplings" (HIFI-C) is an extension of the adaptive and intelligent data collection approach introduced earlier in HIFI-NMR. HIFI-C collects one or more optimally tilted two-dimensional (2D) planes of a 3D experiment, identifies peaks, and determines couplings with high resolution and precision. The HIFI-C approach, demonstrated here for the 3D quantitative J method, offers vital features that advance the goal of rapid and robust collection of NMR coupling data. (1) Tilted plane residual dipolar couplings (RDC) data are collected adaptively in order to offer an intelligent trade off between data collection time and accuracy. (2) Data from independent planes can provide a statistical measure of reliability for each measured coupling. (3) Fast data collection enables measurements in cases where sample stability is a limiting factor (for example in the presence of an orienting medium required for residual dipolar coupling measurements). (4) For samples that are stable, or in experiments involving relatively stronger couplings, robust data collection enables more reliable determinations of couplings in shorter time, particularly for larger biomolecules. As a proof of principle, we have applied the HIFI-C approach to the 3D quantitative J experiment to determine N-C' RDC values for three proteins ranging from 56 to 159 residues (including a homodimer with 111 residues in each subunit). A number of factors influence the robustness and speed of data collection. These factors include the size of the protein, the experimental set up, and the coupling being measured, among others. To exhibit a lower bound on robustness and the potential for time saving, the measurement of dipolar couplings for the N-C' vector represents a realistic
Energy Science and Technology Software Center (ESTSC)
2013-06-24
Version 07 TART2012 is a coupled neutron-photon Monte Carlo transport code designed to use three-dimensional (3-D) combinatorial geometry. Neutron and/or photon sources as well as neutron induced photon production can be tracked. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART2012 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared tomore » other similar codes. Use of the entire system can save you a great deal of time and energy. TART2012 extends the general utility of the code to even more areas of application than available in previous releases by concentrating on improving the physics, particularly with regard to improved treatment of neutron fission, resonance self-shielding, molecular binding, and extending input options used by the code. Several utilities are included for creating input files and displaying TART results and data. TART2012 uses the latest ENDF/B-VI, Release 8, data. New for TART2012 is the use of continuous energy neutron cross sections, in addition to its traditional multigroup cross sections. For neutron interaction, the data are derived using ENDF-ENDL2005 and include both continuous energy cross sections and 700 group neutron data derived using a combination of ENDF/B-VI, Release 8, and ENDL data. The 700 group structure extends from 10-5 eV up to 1 GeV. Presently nuclear data are only available up to 20 MeV, so that only 616 of the groups are currently used. For photon interaction, 701 point photon data were derived using the Livermore EPDL97 file. The new 701 point structure extends from 100 eV up to 1 GeV, and is currently used over this entire energy range. TART2012 completely supersedes all older versions of TART, and it is strongly recommended that one use only the most recent version of TART2012 and its data files. Check authors homepage for related information: http
Boley, A. C.; Morris, M. A.; Desch, S. J.
2013-10-20
A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H{sub 2} is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ∼few× 10{sup –8} L{sub ☉}. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.
NASA Astrophysics Data System (ADS)
Kajzer, A.; Pozorski, J.; Szewc, K.
2014-08-01
In the paper we present Large-eddy simulation (LES) results of 3D Taylor- Green vortex obtained by the three different computational approaches: Smoothed Particle Hydrodynamics (SPH), Lattice Boltzmann Method (LBM) and Finite Volume Method (FVM). The Smagorinsky model was chosen as a subgrid-scale closure in LES for all considered methods and a selection of spatial resolutions have been investigated. The SPH and LBM computations have been carried out with the use of the in-house codes executed on GPU and compared, for validation purposes, with the FVM results obtained using the open-source CFD software OpenFOAM. A comparative study in terms of one-point statistics and turbulent energy spectra shows a good agreement of LES results for all methods. An analysis of the GPU code efficiency and implementation difficulties has been made. It is shown that both SPH and LBM may offer a significant advantage over mesh-based CFD methods.
Meineke, G; Hermans, M; Klos, J; Lenenbach, A; Noll, R
2016-03-01
In a miniaturised flow switch fluid flows are controlled by reducing the local viscosity via absorption of laser radiation. Through this, the local flow rates are increased to switch the outlet port of a fluid flow carrying the analyte. The microfluidic chip is fabricated using Selective Laser-Induced Etching (SLE). SLE allows novel 3D-hydrodynamic focusing, realising circular shaped channel cross-sections and adapting interaction volume geometries to the profile of the laser radiation for optimised absorption. The performance of the switch is validated experimentally with a dyed analyte and video image processing. The ability to sort particles like cells is demonstrated at 8 Hz using polystyrene beads having a diameter of 8 μm. PMID:26862603
3D thermal modeling of TRISO fuel coupled with neutronic simulation
Hu, Jianwei; Uddin, Rizwan
2010-01-01
The Very High Temperature Gas Reactor (VHTR) is widely considered as one of the top candidates identified in the Next Generation Nuclear Power-plant (NGNP) Technology Roadmap under the U.S . Depanment of Energy's Generation IV program. TRlSO particle is a common element among different VHTR designs and its performance is critical to the safety and reliability of the whole reactor. A TRISO particle experiences complex thermo-mechanical changes during reactor operation in high temperature and high burnup conditions. TRISO fuel performance analysis requires evaluation of these changes on micro scale. Since most of these changes are temperature dependent, 3D thermal modeling of TRISO fuel is a crucial step of the whole analysis package. In this paper, a 3D numerical thermal model was developed to calculate temperature distribution inside TRISO and pebble under different scenarios. 3D simulation is required because pebbles or TRISOs are always subjected to asymmetric thermal conditions since they are randomly packed together. The numerical model was developed using finite difference method and it was benchmarked against ID analytical results and also results reported from literature. Monte-Carlo models were set up to calculate radial power density profile. Complex convective boundary condition was applied on the pebble outer surface. Three reactors were simulated using this model to calculate temperature distribution under different power levels. Two asymmetric boundary conditions were applied to the pebble to test the 3D capabilities. A gas bubble was hypothesized inside the TRISO kernel and 3D simulation was also carried out under this scenario. Intuition-coherent results were obtained and reported in this paper.
NASA Astrophysics Data System (ADS)
Wu, Hong; Li, Peng; Li, Yulong
2016-02-01
This paper describes the calculation method for unsteady state conditions in the secondary air systems in gas turbines. The 1D-3D-Structure coupled method was applied. A 1D code was used to model the standard components that have typical geometric characteristics. Their flow and heat transfer were described by empirical correlations based on experimental data or CFD calculations. A 3D code was used to model the non-standard components that cannot be described by typical geometric languages, while a finite element analysis was carried out to compute the structural deformation and heat conduction at certain important positions. These codes were coupled through their interfaces. Thus, the changes in heat transfer and structure and their interactions caused by exterior disturbances can be reflected. The results of the coupling method in an unsteady state showed an apparent deviation from the existing data, while the results in the steady state were highly consistent with the existing data. The difference in the results in the unsteady state was caused primarily by structural deformation that cannot be predicted by the 1D method. Thus, in order to obtain the unsteady state performance of a secondary air system more accurately and efficiently, the 1D-3D-Structure coupled method should be used.
NASA Astrophysics Data System (ADS)
Le Fouest, Vincent; Chami, Malik; Verney, Romaric
2015-02-01
The export of riverine suspended particulate matter (SPM) in the coastal ocean has major implications for the biogeochemical cycles. In the Mediterranean Sea (France), the Rhone River inputs of SPM into the Gulf of Lion (GoL) are highly variable in time, which severely impedes the assessment of SPM fluxes. The objectives of this study are (i) to investigate the prediction of the land-to-ocean flux of SPM using the complementarity (i.e., synergy) between a hydrodynamic sediment transport model and satellite observations, and (ii) to analyze the spatial distribution of the SPM export. An original approach that combines the MARS-3D model with satellite ocean color data is proposed. Satellite-derived SPM and light penetration depth are used to initialize MARS-3D and to validate its predictions. A sensitivity analysis is performed to quantify the impact of riverine SPM size composition and settling rate on the horizontal export of SPM. The best agreement between the model and the satellite in terms of SPM spatial distribution and export is obtained for two conditions: (i) when the relative proportion of "heavy and fast" settling particles significantly increases relative to the "light and slow" ones, and (ii) when the settling rate of heavy and light SPM increases by fivefold. The synergy between MARS-3D and the satellite data improved the SPM flux predictions by 48% near the Rhone River mouth. Our results corroborate the importance of implementing satellite observations within initialization procedures of ocean models since data assimilation techniques may fail for river floods showing strong seasonal variability.
NASA Astrophysics Data System (ADS)
Teruzzi, Anna; Dobricic, Srdjan; Solidoro, Cosimo; Cossarini, Gianpiero
2014-01-01
Increasing attention is dedicated to the implementation of suitable marine forecast systems for the estimate of the state of the ocean. Within the framework of the European MyOcean infrastructure, the pre-existing short-term Mediterranean Sea biogeochemistry operational forecast system has been upgraded by assimilating remotely sensed ocean color data in the coupled transport-biogeochemical model OPATM-BFM using a 3-D variational data assimilation (3D-VAR) procedure. In the present work, the 3D-VAR scheme is used to correct the four phytoplankton functional groups included in the OPATM-BFM in the period July 2007 to September 2008. The 3D-VAR scheme decomposes the error covariance matrix using a sequence of different operators that account separately for vertical covariance, horizontal covariance, and covariance among biogeochemical variables. The assimilation solution is found in a reduced dimensional space, and the innovation for the biogeochemical variables is obtained by the sequential application of the covariance operators. Results show a general improvement in the forecast skill, providing a correction of the basin-scale bias of surface chlorophyll concentration and of the local-scale spatial and temporal dynamics of typical bloom events. Further, analysis of the assimilation skill provides insights into the functioning of the model. The computational costs of the assimilation scheme adopted are low compared to other assimilation techniques, and its modular structure facilitates further developments. The 3D-VAR scheme results especially suitable for implementation within a biogeochemistry operational forecast system.
Chenot, Jean-Loup; Bay, Francois
2007-04-07
The different stages of metal forming processes often involve - beyond the mechanical deformations processes - other physical coupled problems, such as heat transfer, electromagnetism or metallurgy. The purpose of this paper is to focus on problems involving electromagnetic couplings. After a brief recall on electromagnetic modeling, we shall then focus on induction heating processes and present some results regarding heat transfer, as well as mechanical couplings. A case showing coupling for metallurgic microstructure evolution will conclude this paper.
NASA Astrophysics Data System (ADS)
Huang, Genghua; Shu, Rong; Hou, Libing; Li, Ming
2014-06-01
Photon counting lidar has an ultra-high sensitivity which can be hundreds even thousands of times higher than the linear detection lidar. It can significantly increase the system's capability of detection rang and imaging density, saving size and power consumings in airborne or space-borne applications. Based on Geiger-mode Si avalanche photodiodes (Si-APD), a prototype photon counting lidar which used 8 APDs coupled with a 1×8-pixel fiber array has been made in June, 2011. The experiments with static objects showed that the photon counting lidar could operate in strong solar background with 0.04 receiving photoelectrons on average. Limited by less counting times in moving platforms, the probability of detection and the 3D imaging density would be lower than that in static platforms. In this paper, a latest fiber array coupled multi-channel photon counting, 3D imaging, airborne lidar system is introduced. The correlation range receiver algorithm of photon counting 3D imaging is improved for airborne signal photon events extraction and noise filter. The 3D imaging experiments in the helicopter shows that the false alarm rate is less than 6×10-7, and the correct rate is better than 99.9% with 4 received photoelectrons and 0.7MHz system noise on average.
Moment coupling in the interaction of atoms and their ions with a 3d-electron shell
Kosarim, A. V.; Smirnov, B. M.; Capitelli, M.; Laricchiuta, A.
2011-09-15
The moment coupling of an interacting ion and an atom with a 3d-electron shell is analyzed for the ground state of identical atoms and ions where resonant charge exchange proceeds with transition of a 4s-electron. The interaction of the ion charge with the atom quadrupole moment is important for this system along with the exchange interactions and spin-orbit interactions inside an isolated atom and an ion. The quadrupole moment for 3d-atoms in the ground states is evaluated. The hierarchy of interactions in a molecular ion is analyzed depending on ion-atom distances and is compared with the standard Hund scheme. The resonant charge exchange proceeds effectively at separations corresponding to an intermediate case between cases 'a' and 'c' of the Hund coupling scheme.
3D dust clouds (Yukawa Balls) in strongly coupled dusty plasmas
Melzer, A.; Passvogel, M.; Miksch, T.; Ikkurthi, V. R.; Schneider, R.; Block, D.; Piel, A.
2010-06-16
Three-dimensional finite systems of charged dust particles confined to concentric spherical shells in a dusty plasma, so-called 'Yukawa balls', have been studied with respect to their static and dynamic properties. Here, we review the charging of particles in a dusty plasma discharge by computer simulations and the respective particle arrangements. The normal mode spectrum of Yukawa balls is measured from the 3D thermal Brownian motion of the dust particles around their equilibrium positions.
Gagner, Renata; Lafitte, Helene; Dormeau, Pascal; Stoudt, Roger H.
2004-07-01
Anticipated Transients Without Scram (ATWS) accident analyses make part of the Safety Analysis Report of the European Pressurized water Reactor (EPR), covering Risk Reduction Category A (Core Melt Prevention) events. This paper deals with three of the most penalizing RRC-A sequences of ATWS caused by mechanical blockage of the control/shutdown rods, regarding their consequences on the Reactor Coolant System (RCS) and core integrity. A new 3D code internal coupling calculation method has been introduced. (authors)
NASA Astrophysics Data System (ADS)
Chen, Wei-Qiu
2015-10-01
Significant progress has been made in mixed boundary-value problems associated with three-dimensional (3D) crack and contact analyses of advanced materials featuring more complexities compared to the conventional isotropic elastic materials. These include material anisotropy and multifield coupling, two typical characteristics of most current multifunctional materials. In this paper we try to present a state-of-the-art description of 3D exact/analytical solutions derived for crack and contact problems of elastic solids with both transverse isotropy and multifield coupling in the latest decade by the potential theory method in the spirit of V. I. Fabrikant, whose ingenious breakthrough brings new vigor and vitality to the old research subject of classical potential theory. We are particularly interested in crack and contact problems with certain nonlinear features. Emphasis is also placed on the coupling between the temperature field (or the like) and other physical fields (e.g., elastic, electric, and magnetic fields). We further highlight the practical significance of 3D contact solutions, in particular in applications related to modern scanning probe microscopes.
Forcing Ferromagnetic Coupling Between Rare-Earth-Metal and 3d Ferromagnetic Films
NASA Astrophysics Data System (ADS)
Sanyal, Biplab; Antoniak, Carolin; Burkert, Till; Krumme, Bernhard; Warland, Anne; Stromberg, Frank; Praetorius, Christian; Fauth, Kai; Wende, Heiko; Eriksson, Olle
2010-04-01
Using density functional calculations, we have studied the magnetic properties of nanocomposites composed of rare-earth-metal elements in contact with 3d transition metals (Fe and Cr). We demonstrate the possibility to obtain huge magnetic moments in such nanocomposites, of order 10μB/rare-earth-metalatom, with a potential to reach the maximum magnetic moment of Fe-Co alloys at the top of the so-called Slater-Pauling curve. A first experimental proof of concept is given by thin-film synthesis of Fe/Gd and Fe/Cr/Gd nanocomposites, in combination with x-ray magnetic circular dichroism.
3D multi-scale analysis of coupled heat and moisture transport and its parallel implementation
NASA Astrophysics Data System (ADS)
Kruis, Jaroslav
2016-06-01
Parallel implementation of two-scale model of coupled heat and moisture transport is described. The coupled heat and moisture transport is based on the Künzel model. Motivation for the two-scale analysis comes from the requirement to describe distribution of the relative humidity and temperature in historical masonry structures.
Coupled aeroelastic oscillations of a turbine blade row in 3D transonic flow
NASA Astrophysics Data System (ADS)
Gnesin, Vitaly; Kolodyazhnaya, Lyubov; Rzadkowski, Romuald
2001-10-01
This paper presents the mutual time - marching method to predict the aeroelastic stability of an oscillating blade row in 3D transonic flow. The ideal gas flow through a blade row is governed by the time dependent Euler equations in conservative form which are integrated by using the explicit monotonous second order accurate Godunov-Kolgan finite volume scheme and moving hybrid H-O grid. The structure analysis uses the modal approach and 3D finite element dynamic model of blade. The blade movement is assumed as a linear combination of the first modes of blade natural oscillations with the modal coefficients depending on time. To demonstrate the capability and correctness of the method, two experimentally investigated test cases have been selected, in which the blades had performed tuned harmonic bending or torsional vibrations (The 1st and 4th standard configurations of the “Workshop on Aeroelasticity in Turbomachines” by Bolcs and Fransson, 1986). The calculated results of aeroelastic behaviour of the blade row (4th standard configuration), are presented over a wide frequency range under different start regimes of interblade phase angle.
Progress and Challenges in Coupled Hydrodynamic-Ecological Estuarine Modeling
Numerical modeling has emerged over the last several decades as a widely accepted tool for investigations in environmental sciences. In estuarine research, hydrodynamic and ecological models have moved along parallel tracks with regard to complexity, refinement, computational po...
NASA Astrophysics Data System (ADS)
Kim, Ho Jun; Lee, Hae June
2016-06-01
The wide applicability of capacitively coupled plasma (CCP) deposition has increased the interest in developing comprehensive numerical models, but CCP imposes a tremendous computational cost when conducting a transient analysis in a three-dimensional (3D) model which reflects the real geometry of reactors. In particular, the detailed flow features of reactive gases induced by 3D geometric effects need to be considered for the precise calculation of radical distribution of reactive species. Thus, an alternative inclusive method for the numerical simulation of CCP deposition is proposed to simulate a two-dimensional (2D) CCP model based on the 3D gas flow results by simulating flow, temperature, and species fields in a 3D space at first without calculating the plasma chemistry. A numerical study of a cylindrical showerhead-electrode CCP reactor was conducted for particular cases of SiH4/NH3/N2/He gas mixture to deposit a hydrogenated silicon nitride (SiN x H y ) film. The proposed methodology produces numerical results for a 300 mm wafer deposition reactor which agree very well with the deposition rate profile measured experimentally along the wafer radius.
Coupling between Microstrip Lines Embedded in Polyimide Layers for 3D-MMICs on Si
NASA Technical Reports Server (NTRS)
Ponchak, George E.; Tentzeris, Emmanouil M.; Papapolymerou, John
2001-01-01
Three-dimensional circuits built upon multiple layers of polyimide are required for constructing SilSiGe monolithic microwavdmillimeter-wave integrated circuits on CMOS (low resistivity) Si wafers. However, the closely spaced transmission lines are susceptible to high levels of coupling, which degrades circuit performance. In this paper, Finite Difference Time Domain (FDTD) analysis and measured characteristics of novel shielding structures that significantly reduce coupling between embedded microstrip lines are presented. A discussion of the electric and magnetic field distributions for the coupled microstrip lines is presented to provide a physical rationale for the presented results.
NASA Astrophysics Data System (ADS)
Bieler, Andre; Altwegg, Kathrin; Balsiger, Hans; Berthelier, Jean-Jacques; Calmonte, Ursina; Combi, Michael; De Keyser, Johan; Fiethe, Björn; Fougere, Nicolas; Fuselier, Stephen; Gasc, Sébastien; Gombosi, Tamas; Hansen, Kenneth; Hässig, Myrtha; Huang, Zhenguang; Jäckel, Annette; Jia, Xianzhe; Le Roy, Lena; Mall, Urs A.; Rème, Henri; Rubin, Martin; Tenishev, Valeriy; Tóth, Gábor; Tzou, Chia-Yu; Wurz, Peter
2015-11-01
67P/Churyumov-Gerasimenko (67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P's neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 1026 molecules s-1.
A Coupled Finite-Volume Model for 2-D Surface and 3-D Subsurface Flows
Technology Transfer Automated Retrieval System (TEKTRAN)
Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed; therefore, hydrologic modeling tools should consider these interactions to provide reliable predictions, especially during rainfall-runoff processes. This paper presents a fully implicit coupled...
NASA Astrophysics Data System (ADS)
Capet, Arthur; Meysman, Filip J. R.; Akoumianaki, Ioanna; Soetaert, Karline; Grégoire, Marilaure
2016-05-01
Three-dimensional (3D) ecosystem models of shelf environments should properly account for the biogeochemical cycling within the sea floor. However, a full and explicit representation of sediment biogeochemistry into 3D ocean models is computationally demanding. Here, we describe a simplified approach to include benthic processes in 3D ocean models, which includes a parameterization of the different pathways for organic matter mineralization and allows for organic matter remobilization by bottom currents and waves. This efficient approach enables decadal simulations that resolve the inertial contribution of the sea floor to the biogeochemical cycling in shelf environments. The model was implemented to analyze the benthic-pelagic coupling in the northwestern shelf of the Black Sea. Three distinct biogeochemical provinces were identified on the basis of fluxes and rates associated with benthic-pelagic coupling. Our model simulations suitably capture the seasonal variability of in situ flux data as well as their regional variation, which stresses the importance of incorporating temporally varying sediment biogeochemistry and resuspension/redeposition cycles in shelf ecosystem models.
NASA Astrophysics Data System (ADS)
Goryachev, Maxim; Tobar, Michael E.
2015-02-01
A new electromagnetic cavity structure, a lattice of 3D cavities consisting of an array of posts and gaps is presented. The individual cavity elements are based on the cylindrical re-entrant (or Klystron) cavity. We show that these cavities can also be thought of as 3D split-ring resonators, which is confirmed by applying symmetry transformations, each of which is an electromagnetic resonator with spatially separated magnetic and electric field. The characteristics of the cavity is used to mimic phonon behaviour of a one-dimensional (1D) chain of atoms. It is demonstrated how magnetic field coupling can lead to phonon-like dispersion curves with acoustical and optical branches. The system is able to reproduce a number of effects typical to 1D lattices exhibiting acoustic vibration, such as band gaps, phonon trapping, and effects of impurities. In addition, quasicrystal emulations predict the results expected from this class of ordered structures. The system is easily scalable to simulate two-dimensional and 3D lattices and shows a new way to engineer arrays of coupled microwave resonators with a variety of possible applications to hybrid quantum systems proposed.
Coupling Between Microstrip Lines Embedded in Polyimide Layers for 3D-MMICs on Si
NASA Technical Reports Server (NTRS)
Ponchak, George E.; Tentzeris, Emmanouil M.; Papapolymerou, John
2001-01-01
Three-dimensional circuits built upon multiple layers of polyimide are required for constructing Si/SiGe monolithic microwave/millimeter-wave integrated circuits on CMOS (low resistivity) Si wafers. However, the closely spaced transmission lines are susceptible to high levels of coupling, which degrades circuit performance. In this paper, Finite Difference Time Domain (FDTD) analysis and measured characteristics of novel shielding structures that significantly reduce coupling between embedded microstrip lines are presented.
Wang, Xu; Ding, Jie; Guo, Wan-Qian; Ren, Nan-Qi
2010-12-01
Investigating how a bioreactor functions is a necessary precursor for successful reactor design and operation. Traditional methods used to investigate flow-field cannot meet this challenge accurately and economically. Hydrodynamics model can solve this problem, but to understand a bioreactor in sufficient depth, it is often insufficient. In this paper, a coupled hydrodynamics-reaction kinetics model was formulated from computational fluid dynamics (CFD) code to simulate a gas-liquid-solid three-phase biotreatment system for the first time. The hydrodynamics model is used to formulate prediction of the flow field and the reaction kinetics model then portrays the reaction conversion process. The coupled model is verified and used to simulate the behavior of an expanded granular sludge bed (EGSB) reactor for biohydrogen production. The flow patterns were visualized and analyzed. The coupled model also demonstrates a qualitative relationship between hydrodynamics and biohydrogen production. The advantages and limitations of applying this coupled model are discussed. PMID:20727741
Georgy, Jacques; Noureldin, Aboelmagd
2011-01-01
Satellite navigation systems such as the global positioning system (GPS) are currently the most common technique used for land vehicle positioning. However, in GPS-denied environments, there is an interruption in the positioning information. Low-cost micro-electro mechanical system (MEMS)-based inertial sensors can be integrated with GPS and enhance the performance in denied GPS environments. The traditional technique for this integration problem is Kalman filtering (KF). Due to the inherent errors of low-cost MEMS inertial sensors and their large stochastic drifts, KF, with its linearized models, has limited capabilities in providing accurate positioning. Particle filtering (PF) was recently suggested as a nonlinear filtering technique to accommodate for arbitrary inertial sensor characteristics, motion dynamics and noise distributions. An enhanced version of PF called the Mixture PF is utilized in this study to perform tightly coupled integration of a three dimensional (3D) reduced inertial sensors system (RISS) with GPS. In this work, the RISS consists of one single-axis gyroscope and a two-axis accelerometer used together with the vehicle’s odometer to obtain 3D navigation states. These sensors are then integrated with GPS in a tightly coupled scheme. In loosely-coupled integration, at least four satellites are needed to provide acceptable GPS position and velocity updates for the integration filter. The advantage of the tightly-coupled integration is that it can provide GPS measurement update(s) even when the number of visible satellites is three or lower, thereby improving the operation of the navigation system in environments with partial blockages by providing continuous aiding to the inertial sensors even during limited GPS satellite availability. To effectively exploit the capabilities of PF, advanced modeling for the stochastic drift of the vertically aligned gyroscope is used. In order to benefit from measurement updates for such drift, which are
Georgy, Jacques; Noureldin, Aboelmagd
2011-01-01
Satellite navigation systems such as the global positioning system (GPS) are currently the most common technique used for land vehicle positioning. However, in GPS-denied environments, there is an interruption in the positioning information. Low-cost micro-electro mechanical system (MEMS)-based inertial sensors can be integrated with GPS and enhance the performance in denied GPS environments. The traditional technique for this integration problem is Kalman filtering (KF). Due to the inherent errors of low-cost MEMS inertial sensors and their large stochastic drifts, KF, with its linearized models, has limited capabilities in providing accurate positioning. Particle filtering (PF) was recently suggested as a nonlinear filtering technique to accommodate for arbitrary inertial sensor characteristics, motion dynamics and noise distributions. An enhanced version of PF called the Mixture PF is utilized in this study to perform tightly coupled integration of a three dimensional (3D) reduced inertial sensors system (RISS) with GPS. In this work, the RISS consists of one single-axis gyroscope and a two-axis accelerometer used together with the vehicle's odometer to obtain 3D navigation states. These sensors are then integrated with GPS in a tightly coupled scheme. In loosely-coupled integration, at least four satellites are needed to provide acceptable GPS position and velocity updates for the integration filter. The advantage of the tightly-coupled integration is that it can provide GPS measurement update(s) even when the number of visible satellites is three or lower, thereby improving the operation of the navigation system in environments with partial blockages by providing continuous aiding to the inertial sensors even during limited GPS satellite availability. To effectively exploit the capabilities of PF, advanced modeling for the stochastic drift of the vertically aligned gyroscope is used. In order to benefit from measurement updates for such drift, which are
Correlated 3D Nanoscale Mapping and Simulation of Coupled Plasmonic Nanoparticles
2015-01-01
Electron tomography in combination with electron energy-loss spectroscopy (EELS) experiments and simulations was used to unravel the interplay between structure and plasmonic properties of a silver nanocuboid dimer. The precise 3D geometry of the particles fabricated by means of electron beam lithography was reconstructed through electron tomography, and the full three-dimensional information was used as an input for simulations of energy-loss spectra and plasmon resonance maps. Excellent agreement between experiment and theory was found throughout, bringing the comparison between EELS imaging and simulations to a quantitative and correlative level. In addition, interface mode patterns, normally masked by the projection nature of a transmission microscopy investigation, could be unambiguously identified through tomographic reconstruction. This work overcomes the need for geometrical assumptions or symmetry restrictions of the sample in simulations and paves the way for detailed investigations of realistic and complex plasmonic nanostructures. PMID:26495933
Coupling 2-D cylindrical and 3-D x-y-z transport computations
Abu-Shumays, I.K.; Yehnert, C.E.; Pitcairn, T.N.
1998-06-30
This paper describes a new two-dimensional (2-D) cylindrical geometry to three-dimensional (3-D) rectangular x-y-z splice option for multi-dimensional discrete ordinates solutions to the neutron (photon) transport equation. Of particular interest are the simple transformations developed and applied in order to carry out the required spatial and angular interpolations. The spatial interpolations are linear and equivalent to those applied elsewhere. The angular interpolations are based on a high order spherical harmonics representation of the angular flux. Advantages of the current angular interpolations over previous work are discussed. An application to an intricate streaming problem is provided to demonstrate the advantages of the new method for efficient and accurate prediction of particle behavior in complex geometries.
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.
Plasma Biomarker Discovery Using 3D Protein Profiling Coupled with Label-Free Quantitation
Beer, Lynn A.; Tang, Hsin-Yao; Barnhart, Kurt T.; Speicher, David W.
2011-01-01
In-depth quantitative profiling of human plasma samples for biomarker discovery remains quite challenging. One promising alternative to chemical derivatization with stable isotope labels for quantitative comparisons is direct, label-free, quantitative comparison of raw LC–MS data. But, in order to achieve high-sensitivity detection of low-abundance proteins, plasma proteins must be extensively pre-fractionated, and results from LC–MS runs of all fractions must be integrated efficiently in order to avoid misidentification of variations in fractionation from sample to sample as “apparent” biomarkers. This protocol describes a powerful 3D protein profiling method for comprehensive analysis of human serum or plasma proteomes, which combines abundant protein depletion and high-sensitivity GeLC–MS/MS with label-free quantitation of candidate biomarkers. PMID:21468938
NASA Astrophysics Data System (ADS)
Ceccuzzi, Silvio; Maggiora, Riccardo; Milanesio, Daniele; Mirizzi, Francesco; Panaccione, Luigi
2011-12-01
The present work compares and experimentally validates the results coming out from the following three Lower Hybrid (LH) coupling codes: Brambilla code (M. Brambilla), GRILL3D-U (Mikhail Irzak, A. F. Ioffe Physico-Technical Institute, Russia) and TOPLHA (Politecnico di Torino, Italy). The conventional grill antenna, operating in FTU in different scenarios, is used as benchmark. The validation with experimental data is carried out with respect to the average reflection coefficients at the input of a row of the grill, considering two different phasings between adjacent waveguides: -90 ° and -75 °. A comparison between calculated power spectra is also presented. Good agreement can be observed for all the simulated plasma profiles and waveguide phasings between experimental data and codes, in particular for the most recent numerical tools, namely GRILL3D-U and TOPLHA.
Teruzzi, Anna; Dobricic, Srdjan; Solidoro, Cosimo; Cossarini, Gianpiero
2014-01-01
[1] Increasing attention is dedicated to the implementation of suitable marine forecast systems for the estimate of the state of the ocean. Within the framework of the European MyOcean infrastructure, the pre-existing short-term Mediterranean Sea biogeochemistry operational forecast system has been upgraded by assimilating remotely sensed ocean color data in the coupled transport-biogeochemical model OPATM-BFM using a 3-D variational data assimilation (3D-VAR) procedure. In the present work, the 3D-VAR scheme is used to correct the four phytoplankton functional groups included in the OPATM-BFM in the period July 2007 to September 2008. The 3D-VAR scheme decomposes the error covariance matrix using a sequence of different operators that account separately for vertical covariance, horizontal covariance, and covariance among biogeochemical variables. The assimilation solution is found in a reduced dimensional space, and the innovation for the biogeochemical variables is obtained by the sequential application of the covariance operators. Results show a general improvement in the forecast skill, providing a correction of the basin-scale bias of surface chlorophyll concentration and of the local-scale spatial and temporal dynamics of typical bloom events. Further, analysis of the assimilation skill provides insights into the functioning of the model. The computational costs of the assimilation scheme adopted are low compared to other assimilation techniques, and its modular structure facilitates further developments. The 3D-VAR scheme results especially suitable for implementation within a biogeochemistry operational forecast system. PMID:26213670
Simulating HFIR Core Thermal Hydraulics Using 3D-2D Model Coupling
Travis, Adam R; Freels, James D; Ekici, Kivanc
2013-01-01
A model utilizing interdimensional variable coupling is presented for simulating the thermal hydraulic interactions of the High Flux Isotope Reactor (HFIR) core at Oak Ridge National Laboratory (ORNL). The model s domain consists of a single, explicitly represented three-dimensional fuel plate and a simplified two-dimensional coolant channel slice. In simplifying the coolant channel, and thus the number of mesh points in which the Navier-Stokes equations must be solved, the computational cost and solution time are both greatly reduced. In order for the reduced-dimension coolant channel to interact with the explicitly represented fuel plate, however, interdimensional variable coupling must be enacted along all shared boundaries. The primary focus of this paper is in detailing the collection, storage, passage, and application of variables across this interdimensional interface. Comparisons are made showing the general speed-up associated with this simplified coupled model.
Goda, Ibrahim; Ganghoffer, Jean-François
2015-11-01
The purpose of this paper is to develop a homogeneous, orthotropic couple-stress continuum model as a substitute of the 3D periodic heterogeneous cellular solid model of vertebral trabecular bone. Vertebral trabecular bone is modeled as a porous material with an idealized periodic structure made of 3D open cubic cells, which is effectively orthotropic. The chosen architecture is based on studies of samples taken from the central part of vertebral bodies. The effective properties are obtained based on the response of the representative volume element under prescribed boundary conditions. Mixed boundary conditions comprising both traction and displacement boundary conditions are applied on the structure boundaries. In this contribution, the effective mechanical constants of the effective couple-stress continuum are deduced by an equivalent strain energy method. The characteristic lengths for bending and torsion are identified from the resulting homogenized orthotropic moduli. We conduct this study computationally using a finite element approach. Vertebral trabecular bone is modeled either as a cellular solid or as a two-phase material consisting of bone tissue (stiff phase) forming a trabecular network, and a surrounding soft tissue referring to the bone marrow present in the pores. Both the bone tissue forming the network and the pores are assumed to be homogeneous linear elastic, and isotropic media. The scale effects on the predicted couple stress moduli of these networks are investigated by varying the size of the bone specimens over which the boundary conditions are applied. The analysis using mixed boundary conditions gives results that are independent of unit cell size when computing the first couple stress tensor, while it is dependent on the cell size as to the second couple stress tensor moduli. This study provides overall guidance on how the size of the trabecular specimen influence couple stresses elastic moduli of cellular materials, with focus on bones
Scattering and coupling effects of electromagnetic waves in 3D networks of spheres
NASA Astrophysics Data System (ADS)
Defos Du Rau, M.; Pessan, F.; Ruffie, G.; Vignéras-Lefebvre, V.; Parneix, J. P.
1998-01-01
In this paper, the problem of electromagnetic scattering from a 3D system of spheres is considered and an iterative solution that accounts for multiple scattering is proposed. The Mie formalism used for a single sphere is extended to account for multiple scattered fields between several particles. The translational addition theorems for spherical wave functions are used to express the electromagnetic field scattered by a sphere S_i in terms of an incident field for a sphere S_k in a spherical coordinates system attached to the sphere S_k. In this work, the numerical convergence of the method is discussed and associated computational times are given. Numerical computations including Radar Cross Section (RCS) and radiation patterns for various 3D configurations are presented. Some of them are compared with free-space measurements made in the 8 to 100 GHz frequency band using vectorial network analyzers. 11.55.-m S-matrix theory; analytic structure of amplitudes Cet article étudie la diffusion des ondes électromagnétiques par des réseaux tridimensionnels de sphères et propose une méthode itérative pour prendre en compte les effets de multidiffusion. Le formalisme de Mie utilisé dans le cas d'une sphère est étendu pour calculer les champs "multidiffusés" entre plusieurs particules. Les théorèmes d'addition et de translation des fonctions d'onde sphériques sont utilisés pour exprimer le champ diffusé par une sphère S_i comme étant incident sur une sphère S_k, dans un système de coordonnées sphériques lié au centre de S_k. La convergence numérique de la méthode est discutée et des temps de calcul sont donnés. Des résultats numériques tels que des Surfaces Équivalentes Radar (SER) et des diagrammes de rayonnement pour différentes configurations tridimensionnelles sont montrés. Certains d'entre eux sont comparés à des mesures en espace libre faites à l'aide d'analyseurs de réseaux vectoriels dans la bande de fréquence 8{-}100 GHz.
Using High Frequency Focused Water-Coupled Ultrasound for 3-D Surface Depression Profiling
NASA Technical Reports Server (NTRS)
Roth, Don J.; Whalen, Mike F.; Hendricks, J. Lynne; Bodis, James R.
1999-01-01
Surface topography is an important variable in the performance of many industrial components and is normally measured with diamond-tip profilometry over a small area or using optical scattering methods for larger area measurement. A prior study was performed demonstrating that focused air-coupled ultrasound at 1 MHz was capable of profiling surfaces with 25 micron depth resolution and 400 micron lateral resolution over a 1.4 mm depth range. In this article, the question of whether higher-frequency focused water-coupled ultrasound can improve on these specifications is addressed. 10 and 25 MHz focused ultrasonic transducers were employed in the water-coupled mode. Time-of-flight images of the sample surface were acquired and converted to depth / surface profile images using the simple relation (d = V*t/2) between distance (d), time-of-flight (t), and the velocity of sound in water (V). Results are compared for the two frequencies used and with those from the 1 MHz air-coupled configuration.
Unsteady Analysis of Inlet-Compressor Acoustic Interactions Using Coupled 3-D and 1-D CFD Codes
NASA Technical Reports Server (NTRS)
Suresh, A.; Cole, G. L.
2000-01-01
It is well known that the dynamic response of a mixed compression supersonic inlet is very sensitive to the boundary condition imposed at the subsonic exit (engine face) of the inlet. In previous work, a 3-D computational fluid dynamics (CFD) inlet code (NPARC) was coupled at the engine face to a 3-D turbomachinery code (ADPAC) simulating an isolated rotor and the coupled simulation used to study the unsteady response of the inlet. The main problem with this approach is that the high fidelity turbomachinery simulation becomes prohibitively expensive as more stages are included in the simulation. In this paper, an alternative approach is explored, wherein the inlet code is coupled to a lesser fidelity 1-D transient compressor code (DYNTECC) which simulates the whole compressor. The specific application chosen for this evaluation is the collapsing bump experiment performed at the University of Cincinnati, wherein reflections of a large-amplitude acoustic pulse from a compressor were measured. The metrics for comparison are the pulse strength (time integral of the pulse amplitude) and wave form (shape). When the compressor is modeled by stage characteristics the computed strength is about ten percent greater than that for the experiment, but the wave shapes are in poor agreement. An alternate approach that uses a fixed rise in duct total pressure and temperature (so-called 'lossy' duct) to simulate a compressor gives good pulse shapes but the strength is about 30 percent low.
NASA Astrophysics Data System (ADS)
Wang, Bo; Xiao, Jiaying; Jiang, Huabei
2014-08-01
Objective. Neurovascular coupling in epilepsy is poorly understood; its study requires simultaneous monitoring of hemodynamic changes and neural activity in the brain. Approach. Here for the first time we present a combined real-time 3D photoacoustic tomography (PAT) and electrophysiology/electroencephalography (EEG) system for the study of neurovascular coupling in epilepsy, whose ability was demonstrated with a pentylenetetrazol (PTZ) induced generalized seizure model in rats. Two groups of experiments were carried out with different wavelengths to detect the changes of oxy-hemoglobin (HbO2) and deoxy-hemoglobin (HbR) signals in the rat brain. We extracted the average PAT signals of the superior sagittal sinus (SSS), and compared them with the EEG signal. Main results. Results showed that the seizure process can be divided into three stages. A ‘dip’ lasting for 1-2 min in the first stage and the following hyperfusion in the second stage were observed. The HbO2 signal and the HbR signal were generally negatively correlated. The change of blood flow was also estimated. All the acquired results here were in accordance with other published results. Significance. Compared to other existing functional neuroimaging tools, the method proposed here enables reliable tracking of hemodynamic signal with both high spatial and high temporal resolution in 3D, so it is more suitable for neurovascular coupling study of epilepsy.
NASA Astrophysics Data System (ADS)
Taha, Uday; Shabeeb, Ahmed; dragonetti, giovanna; Lamaddalena, Nicola; Coppola, Antonio
2016-04-01
This work analyzed the variability of sprinkler irrigation application over a bare soil, both in terms of water application efficiency and uniformity, by integrating and comparing the information on the irrigation depth data (ID), as measured by catch cans, soil water storage in the upper root zone, as measured by TDR probes, and a 3D simulations of water flow in soils. Three irrigation tests were performed at three different pressures (2, 3 and 4 bar). A lateral water redistribution was observed and simulated after each irrigation event by comparing spatial distributions of site-specific water application efficiency (AEs), as well as ratios of site-specific actual water storage increase (SWEs) and irrigation depth (IDs) to the water content before irrigation. Because of soil water redistribution processes, distribution uniformity based on soil storages was systematically higher than the catch can uniformity. The obvious consequence of lateral water redistribution processes was that the soil smoothing action on non-uniformity observed at the surface increased both with depth and over time. At a given depth the uniformity of soil water storages always attained the same value, whatever the pressure considered and the catch can-based uniformity coefficient. It was concluded that, for the case of random distribution of ID, the uniformity of water storages is driven by the soil behavior rather than by the irrigation system.
NASA Astrophysics Data System (ADS)
Bauknecht, Andreas; Steinert, Torsten; Spengler, Carsten; Suck, Gerrit
2013-07-01
Thermoelectric (TE) modules with annular geometry are very attractive for waste heat recovery within the automotive world, especially when integrated as stacks into tubular heat exchangers. The required temperature difference is built up between the coolant, which flows inside an inner tube, and the exhaust gas, which flows around an outer tube. The flow pattern of the exhaust gas can be axial or circumferential, which can lead to higher heat transfer coefficients on the outer surface of the tube. However, this multidimensional construction in combination with a complex flow pattern can lead to a nonuniform heat flux. Additionally, the system experiences a nonuniform temperature distribution which consequently leads to complex conditions regarding the electrical potential. The relevant effects are investigated using a three-dimensional (3-D) numerical model implemented in the computational fluid dynamics (CFD) simulation environment Star-CCM+. The model supports temperature-dependent characteristics of the materials, contact resistances, and parasitic effects in the TE module. Furthermore, it involves techniques to quickly find the exact maximum power point of the TE module with the given boundary conditions. Using the validated model the influence of the nonuniform temperature distribution is investigated with emphasis on the electrical output and TE efficiency.
Collective dynamics in a binary mixture of hydrodynamically coupled microrotors.
Yeo, Kyongmin; Lushi, Enkeleida; Vlahovska, Petia M
2015-05-01
We study, numerically, the collective dynamics of self-rotating nonaligning particles by considering a monolayer of spheres driven by constant clockwise or counterclockwise torques. We show that hydrodynamic interactions alter the emergence of large-scale dynamical patterns compared to those observed in dry systems. In dilute suspensions, the flow stirred by the rotors induces clustering of opposite-spin rotors, while at higher densities same-spin rotors phase separate. Above a critical rotor density, dynamic hexagonal crystals form. Our findings underscore the importance of inclusion of the many-body, long-range hydrodynamic interactions in predicting the phase behavior of active particles. PMID:26001020
A 3-D Coupled CFD-DSMC Solution Method With Application to the Mars Sample Return Orbiter
NASA Technical Reports Server (NTRS)
Glass, Christopher E.; Gnoffo, Peter A.
2000-01-01
A method to obtain coupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC), 3-D flow field solutions for highly blunt bodies at low incidence is presented and applied to one concept of the Mars Sample Return Orbiter vehicle as a demonstration of the technique. CFD is used to solve the high-density blunt forebody flow defining an inflow boundary condition for a DSMC solution of the afterbody wake flow. By combining the two techniques in flow regions where most applicable, the entire mixed flow field is modeled in an appropriate manner.
A 3D Model to Compute Lightning and HIRF Coupling Effects on Avionic Equipment of an Aircraft
NASA Astrophysics Data System (ADS)
Perrin, E.; Tristant, F.; Guiffaut, C.; Terrade, F.; Reineix, A.
2012-05-01
This paper describes the 3D FDTD model of an aircraft developed to compute the lightning and HIRF (High Intentity Radiated Fields) coupling effects on avionic equipment and all the wire harness associated. This virtual prototype aims at assisting the aircraft manufacturer during the lightning and HIRF certification processes. The model presented here permits to cover a frequency range from lightning spectrum to the low frequency HIRF domain, i.e. 0 to 100 MHz. Moreover, the entire aircraft, including the frame, the skin, the wire harness and the equipment are taken into account in only one model. Results obtained are compared to measurements on a real aircraft.
NASA Astrophysics Data System (ADS)
Phanikumar, M. S.; McGuire, Jennifer T.
2004-06-01
This paper reports the development and application of a three-dimensional multi-component reactive transport model (BGTK) to simulate a wide range of biogeochemical processes in subsurface environments. The model can handle both equilibrium and kinetically controlled reactions and is based on the well tested modular models RT3D [Clement et al., 1998] and PHREEQC-2 [Parkhurst and Appelo, 1999]. Here we describe the details of the new coupled model and demonstrate its capabilities using test problems involving microbial transport in a laboratory column and redox zonation in a contaminated aquifer.
NASA Astrophysics Data System (ADS)
Uhrig, Matthias P.; Kim, Jin-Yeon; Jacobs, Laurence J.
2016-02-01
This research presents a 3D numerical finite element (FE) model which, previously developed, precisely simulates non-contact, air-coupled measurements of nonlinear Rayleigh wave propagation. The commercial FE-solver ABAQUS is used to perform the simulations. First, frequency dependent pressure wave attenuation is investigated numerically to reconstruct the sound pressure distribution along the active surface of the non-contact receiver. Second, constitutive law and excitation source properties are optimized to match nonlinear ultrasonic experimental data. Finally, the FE-model data are fit with analytical solutions showing a good agreement and thus, indicating the significance of the study performed.
NASA Astrophysics Data System (ADS)
Kessouri, Fayçal; Ulses, Caroline; Estournel, Claude; Marsaleix, Patrick
2015-04-01
The Mediterranean Sea presents a wide variety of trophic regimes since the large and intense spring bloom of the North-Western Mediterranean Sea (NWMS) that follows winter convection to the extreme oligotrophic regions of the South-eastern basin. The Mediterranean Sea displays a strong time variability revealing its high sensitivity to climate and anthropic pressures. In this context, it is crucial to develop tools allowing to understand the evolution of the Mediterranean hydrology and marine ecosystem as a response to external forcing. Numerical coupled hydrodynamic and biogeochemical modelling carefully calibrated in the different regions of the basin is the only tool that can answer this question. However, this important step of calibration is particularly difficult because of the lack of coherent sets of data describing the seasonal evolution of the main parameters characterizing the physical and biogeochemical environment in the different sub-basins. The chlorophyll satellite data from 4km MODIS products, a multiple in situ data from MerMEX MOOSE and DEWEX cruises and Bio-Argo floats from NAOS project are believed to be an opportunity to strongly improve the realism of ecosystem models. The model is a 3D coupled simulation using NemoMed12 for hydrodynamics and ECO 3MS for biogeochemistry and covers the whole Mediterranean Sea and runs at 1/12°. The relevant variables mentioned are phytoplankton, organic and inorganic matters faced to water masses dynamics, over ten years since summer 2003. After a short validation, we will expose two topics: First, through this coupling we quantify the nutrients fluxes across the Mediterranean straits over the years. For example, we found an annual net average around 150 Giga moles NO3 per year at Gibraltar, where we expect low annual fluctuations. In contrast, the Strait of Sicily shows greater annual variability going from 70 to 92 Giga moles NO3 per year. All the fluxes are resumed in a detailed diagram of the transport
Intra-chain superexchange couplings in quasi-1D 3d transition-metal magnetic compounds.
Xiang, Hongping; Tang, Yingying; Zhang, Suyun; He, Zhangzhen
2016-07-13
The electronic structure and magnetic properties of the quasi-1D transition-metal borates PbMBO4 (M = Ti, V, Cr, Mn, Fe, Co) have been investigated by density functional theory, including electronic correlation. The results evidence PbCrBO4 and PbFeBO4 as antiferromagnetic (AFM) semiconductors (intra-chain AFM and inter-chain FM) and PbMnBO4 as a ferromagnetic (FM) semiconductor (both intra- and inter-chain FM) in accordance with experimental observations. For non-synthesized PbTiBO4, PbVBO4, and PbCoBO4, the ground-state magnetic structures are paramagnetic, FM, and paramagnetic, respectively. In this series of compounds, there are two kinds of superexchange couplings dominating their magnetic properties, i.e. the direction M-M delocalization superexchange and indirect M-O-M correlation superexchange. For PbMBO4 with M (3+) d (n) , n ⩽ 3 (M = V and Cr), the main intra-chain spin coupling is the M-M t 2g-t 2g direct delocalization superexchange, while for PbMBO4 with M (3+) d (n) , n > 3 (M = Mn and Fe), the main intra-chain spin coupling is the near 90° M-O-M e g-p-e g indirect correlation superexchange. PMID:27213502
Intra-chain superexchange couplings in quasi-1D 3d transition-metal magnetic compounds
NASA Astrophysics Data System (ADS)
Xiang, Hongping; Tang, Yingying; Zhang, Suyun; He, Zhangzhen
2016-07-01
The electronic structure and magnetic properties of the quasi-1D transition-metal borates PbMBO4 (M = Ti, V, Cr, Mn, Fe, Co) have been investigated by density functional theory, including electronic correlation. The results evidence PbCrBO4 and PbFeBO4 as antiferromagnetic (AFM) semiconductors (intra-chain AFM and inter-chain FM) and PbMnBO4 as a ferromagnetic (FM) semiconductor (both intra- and inter-chain FM) in accordance with experimental observations. For non-synthesized PbTiBO4, PbVBO4, and PbCoBO4, the ground-state magnetic structures are paramagnetic, FM, and paramagnetic, respectively. In this series of compounds, there are two kinds of superexchange couplings dominating their magnetic properties, i.e. the direction M–M delocalization superexchange and indirect M–O–M correlation superexchange. For PbMBO4 with M 3+ d n , n ⩽ 3 (M = V and Cr), the main intra-chain spin coupling is the M–M t 2g–t 2g direct delocalization superexchange, while for PbMBO4 with M 3+ d n , n > 3 (M = Mn and Fe), the main intra-chain spin coupling is the near 90° M–O–M e g–p–e g indirect correlation superexchange.
3D-PIC simulation of an inductively coupled ion source
NASA Astrophysics Data System (ADS)
Henrich, Robert; Muehlich, Nina Sarah; Becker, Michael; Heiliger, Christian
2015-09-01
Inductively coupled ion sources are applied to a wide range of plasma applications, especially surface modifications. The knowledge of the behavior and precise information of the plasma parameters are of main importance. These values are tedious to measure without influencing the discharge. By applying our fully three-dimensional PlasmaPIC tool we are able to reach these plasma parameters with a spatial and temporal resolution which is quite hard to achieve experimentally. PlasmaPIC is used for modeling discharges in arbitrary geometries without limitations to any symmetry. By this means we are able to demonstrate that the plasma density has an irrotational character. Furthermore, we will show the dependence of the plasma parameters of different working conditions. We will show that for gridded inductively coupled ion sources the neutral gas pressure inside the discharge chamber depends on the extraction of ions. This effect is considered in PlasmaPIC by a self-consistent coupling of the neutral gas simulation and the plasma simulation whereas the neutral gas distribution is calculated using the direct simulation Monte Carlo method (DSMC). This work has been supported by the ``Bundesministerium fuer Wirtschaft und Energie.'' Grant 50RS1507.
NASA Astrophysics Data System (ADS)
Nguyen, Minh-Hang; Chu, Thi-Xuan; Nguyen, Long; Nguyen, Hai-Binh; Lee, Chun-Wei; Tseng, Fan-Gang; Chen, Te-Chang; Lee, Ming-Chang
2016-07-01
Fabrication of three-dimensional (3D) SU-8 (an epoxy-based negative photoresist from MicroChem) prisms as low-loss couplers for interconnection between optical components, particularly optical fibers and silicon-on-isolator waveguides (SOI WGs), which have mismatched mode sizes, has been investigated. With an interfacial structure formed by a 3D SU-8 prism partly overlaying an SOI WG end with a portion of buried oxide (BOX) removed under the interface, low-loss coupling is ensured and the transmission efficiency can reach 70%. To fabricate these 3D SU-8 prisms, a simple method with two photolithography steps was used for SU-8 hinges and CYTOP (an amorphous fluoropolymer from AGC Chemicals) prism windows, with mild soft and hard bakes, to define the prism profiles with diluted SU-8 filled in the CYTOP prism windows. A buffered oxide etchant is used to remove BOX parts under the interfaces. Some of the fabricated structures were tested, demonstrating the contribution of overlaying SU-8 prisms to the transmission efficiency of optical interconnections between fibers and SOI WGs.
Modeling partially coupled objects with smooth particle hydrodynamics
Wingate, C.A.
1996-10-01
A very simple phenomenological model is presented to model objects that are partially coupled (i.e. welded or bonded) where usually the coupled interface is weaker than the bulk material. The model works by letting objects fully interact in compression and having the objects only partially interact in tension. A disconnect factor is provided to adjust the tensile interaction to simulate coupling strengths. Three cases of an example impact calculation are shown-no coupling, full coupling and partial coupling.
NASA Astrophysics Data System (ADS)
Hou, Weisheng; Yang, Liang; Deng, Dongcheng; Ye, Jing; Clarke, Keith; Yang, Zhijun; Zhuang, Wenming; Liu, Jianxiong; Huang, Jichun
2016-04-01
Urban underground spaces (UUS), especially those containing natural resources that have not yet been utilized, have been recognized as important for future sustainable development in large cities. One of the key steps in city planning is to estimate the quality of urban underground space resources, since they are major determinants of suitable land use. Yet geological constraints are rarely taken into consideration in urban planning, nor are the uncertainties in the quality of the available assessments. Based on Fuzzy Set theory and the analytic hierarchy process, a 3D stepwise process for the quality assessment of geotechnical properties of natural resources in UUS is presented. The process includes an index system for construction factors; area partitioning; the extraction of geological attributes; the creation of a relative membership grade matrix; the evaluation of subject and destination layers; and indeterminacy analysis. A 3D geological model of the study area was introduced into the process that extracted geological attributes as constraints. This 3D geological model was coupled with borehole data for Foshan City, Guangdong province, South China, and the indeterminacies caused by the cell size and the geological strata constraints were analyzed. The results of the case study show that (1) a relatively correct result can be obtained if the cell size is near to the average sampling distance of the boreholes; (2) the constraints of the 3D geological model have a major role in establishing the UUS quality level and distribution, especially at the boundaries of the geological bodies; and (3) the assessment result is impacted by an interaction between the cell resolution and the geological model used.
The research of 1D / 3D coupling simulation on pump and pipe system
NASA Astrophysics Data System (ADS)
Wu, D. Z.; Liu, Q. L.; Wu, P.; Wang, L. Q.; Paulus, T.; Wang, B. G.; Oesterle, M.
2012-11-01
The research of performances of hydraulic mechanical depends on static complete characteristic curves, which have great difference compared with the actual work condition and have accidents potential. So we need a new way to compute the dynamic system, which is more reasonable. So the method to couple one dimensional simulation and three dimensional CFD analysis based on Flowmaster and Fluent is explored, and the dynamic characteristics and internal flow of the pumping system are analyzed. First, a pipe system model is created in Flowmaster and a pump model is created in Fluent; then VB code and scheme code are used to realize the automated operation for Flowmaster and Fluent; at last, the exchange of data between these two parts is realized by an interface program. In this paper, the interaction between pumps and pipe system are analyzed by coupling one-dimensional and three-dimensional simulations. This study would be helpful to identify the influences of the rapid adjustment process on stability of system and provide guides for design of pump system.
Wu, Kenan; Huai, Ying; Jia, Shuqin; Jin, Yuqi
2011-12-19
Coupled simulation based on intracavity partially coherent light model and 3D CFD model is firstly achieved in this paper. The dynamic equation of partially coherent intracavity field is derived based on partially coherent light theory. A numerical scheme for the coupled simulation as well as a method for computing the intracavity partially coherent field is given. The presented model explains the formation of the sugar scooping phenomenon, and enables studies on the dependence of the spatial mode spectrum on physical parameters of laser cavity and gain medium. Computational results show that as the flow rate of iodine increases, higher order mode components dominate in the partially coherent field. Results obtained by the proposed model are in good agreement with experimental results. PMID:22274214
TART97 a coupled neutron-photon 3-D, combinatorial geometry Monte Carlo transport code
Cullen, D.E.
1997-11-22
TART97 is a coupled neutron-photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy. TART97 is distributed on CD. This CD contains on- line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and its data riles.
NASA Astrophysics Data System (ADS)
Damale, A. V.; Nandurkar, K. N.
2012-10-01
Welding induced distortion is one of the critical defects in the welded structures. Angular distortion is most pronounced which badly affects the welded structures. Non-uniform heating during welding develops this angular distortion. Various methods are available to control/minimize the welding distortions. One of the methods available to control this distortion is restraining, in which clamping pressure is applied on the deforming edges of the plates. In the present study, a 3-dimensional coupled transient thermal analysis is done for simulating the restraining phenomenon of welding. The developed transient thermal heat source was used to simulate the arc welding phenomenon. The element birth and death technique was used for simulating filler material deposition. Thermal model was verified by comparing the macrograph of Finite Element Analysis (FEA) model and the weld and verification of structural model was done by comparing the measured and predicted angular distortions. Experimental and Finite Element (FE) analysis is done for both conventional welding and welding under restraining pressure. Transient thermal and non-linear structural analyses were carried out in order to predict angular distortions. The Finite Element Method analysis and experimental verification is done for manual metal arc welding process.
A Coupled Neutron-Photon 3-D Combinatorial Geometry Monte Carlo Transport Code
Energy Science and Technology Software Center (ESTSC)
1998-06-12
TART97 is a coupled neutron-photon, 3 dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly fast: if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system canmore » save you a great deal of time and energy. TART 97 is distributed on CD. This CD contains on-line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and ist data files.« less
NASA Astrophysics Data System (ADS)
Lee, Min-Gu; Baeg, Seung-Ho; Lee, Ki-Min; Lee, Hae-Seok; Baeg, Moon-Hong; Park, Jong-Ok; Kim, Hong-Ki
2011-06-01
The purpose of this research is to develop a new 3D LIDAR sensor, named KIDAR-B25, for measuring 3D image information with high range accuracy, high speed and compact size. To measure a distance to the target object, we developed a range measurement unit, which is implemented by the direct Time-Of-Flight (TOF) method using TDC chip, a pulsed laser transmitter as an illumination source (pulse width: 10 ns, wavelength: 905 nm, repetition rate: 30kHz, peak power: 20W), and an Si APD receiver, which has high sensitivity and wide bandwidth. Also, we devised a horizontal and vertical scanning mechanism, climbing in a spiral and coupled with the laser optical path. Besides, control electronics such as the motor controller, the signal processing unit, the power distributor and so on, are developed and integrated in a compact assembly. The key point of the 3D LIDAR design proposed in this paper is to use the compact scanning mechanism, which is coupled with optical module horizontally and vertically. This KIDAR-B25 has the same beam propagation axis for emitting pulse laser and receiving reflected one with no optical interference each other. The scanning performance of the KIDAR-B25 has proven with the stable operation up to 20Hz (vertical), 40Hz (horizontal) and the time is about 1.7s to reach the maximum speed. The range of vertical plane can be available up to +/-10 degree FOV (Field Of View) with a 0.25 degree angular resolution. The whole horizontal plane (360 degree) can be also available with 0.125 degree angular resolution. Since the KIDAR-B25 sensor has been planned and developed to be used in mobile robots for navigation, we conducted an outdoor test for evaluating its performance. The experimental results show that the captured 3D imaging data can be usefully applicable to the navigation of the robot for detecting and avoiding the moving objects with real time.
Shu, Yu-Chen; Chern, I-Liang; Chang, Chien C.
2014-10-15
Most elliptic interface solvers become complicated for complex interface problems at those “exceptional points” where there are not enough neighboring interior points for high order interpolation. Such complication increases especially in three dimensions. Usually, the solvers are thus reduced to low order accuracy. In this paper, we classify these exceptional points and propose two recipes to maintain order of accuracy there, aiming at improving the previous coupling interface method [26]. Yet the idea is also applicable to other interface solvers. The main idea is to have at least first order approximations for second order derivatives at those exceptional points. Recipe 1 is to use the finite difference approximation for the second order derivatives at a nearby interior grid point, whenever this is possible. Recipe 2 is to flip domain signatures and introduce a ghost state so that a second-order method can be applied. This ghost state is a smooth extension of the solution at the exceptional point from the other side of the interface. The original state is recovered by a post-processing using nearby states and jump conditions. The choice of recipes is determined by a classification scheme of the exceptional points. The method renders the solution and its gradient uniformly second-order accurate in the entire computed domain. Numerical examples are provided to illustrate the second order accuracy of the presently proposed method in approximating the gradients of the original states for some complex interfaces which we had tested previous in two and three dimensions, and a real molecule ( (1D63)) which is double-helix shape and composed of hundreds of atoms.
NASA Astrophysics Data System (ADS)
Shu, Yu-Chen; Chern, I.-Liang; Chang, Chien C.
2014-10-01
Most elliptic interface solvers become complicated for complex interface problems at those “exceptional points” where there are not enough neighboring interior points for high order interpolation. Such complication increases especially in three dimensions. Usually, the solvers are thus reduced to low order accuracy. In this paper, we classify these exceptional points and propose two recipes to maintain order of accuracy there, aiming at improving the previous coupling interface method [26]. Yet the idea is also applicable to other interface solvers. The main idea is to have at least first order approximations for second order derivatives at those exceptional points. Recipe 1 is to use the finite difference approximation for the second order derivatives at a nearby interior grid point, whenever this is possible. Recipe 2 is to flip domain signatures and introduce a ghost state so that a second-order method can be applied. This ghost state is a smooth extension of the solution at the exceptional point from the other side of the interface. The original state is recovered by a post-processing using nearby states and jump conditions. The choice of recipes is determined by a classification scheme of the exceptional points. The method renders the solution and its gradient uniformly second-order accurate in the entire computed domain. Numerical examples are provided to illustrate the second order accuracy of the presently proposed method in approximating the gradients of the original states for some complex interfaces which we had tested previous in two and three dimensions, and a real molecule (1D63) which is double-helix shape and composed of hundreds of atoms.
NASA Astrophysics Data System (ADS)
Choi, Jiwoong; Yin, Youbing; Hoffman, Eric; Tawhai, Merryn; Lin, Ching-Long
2010-11-01
Correct prediction of regional distribution of inhaled aerosol particles is vital to improve pulmonary medicine. Physiologically consistent regional ventilations of airflow and aerosol particles are simulated with a 3D-1D coupled subject-specific boundary condition (BC). In 3D CT-resolved 7-generation airways, large eddy simulations are performed to capture detailed airflow characteristics and Lagrangian particle simulations are carried to track the particle transport and deposition. Results are compared with two traditional outlet BCs: uniform velocity and uniform pressure. Proposed BC is eligible for physiologically consistent airflow distribution in the lung, while the others are not. The regional ventilation and deposition of particles reflect the regional ventilation of airflow. In this study, two traditional BCs yield up to 98% (334%) over-prediction in lobar particle ventilation (deposition) fraction. Upper to lower particle ventilation ratios of both left and right lungs read ˜0.4 with the proposed BC, while those for the other two BCs vary with the error up to 73%.
NASA Astrophysics Data System (ADS)
Pathak, Himanshu; Singh, Akhilendra; Singh, Indra Vir
2016-06-01
In this work, finite element method (FEM) and element free Galerkin method (EFGM) are coupled for solving 3D crack domains subjected to cyclic thermal load of constant amplitude. Crack growth contours and fatigue life have been obtained for each of the considered numerical examples. Thermo-elastic problems are decoupled into thermal and elastic problems . Firstly, the unknown temperature field is obtained by solving heat conduction equation, then, it is used as the input load in the elastic problem to calculate the displacement and stress fields. The geometrical discontinuity across crack surface is modelled by extrinsically enriched EFGM and the remaining part of the domain is approximated by standard finite element method. At the crack interface, a ramp function based interpolation scheme has been implemented. This coupled approach combines the advantages of both EFGM and FEM. A linear successive crack increment approach is used to model crack growth. The growing crack surface is traced by level set function. Standard Paris law is used for life estimation of the three-dimensional crack models. Different cases of planar and non-planar crack problems have been solved and their results are compared with the results obtained using extended finite element method to check accuracy, efficiency and robustness of the coupled FE-EFG approach implemented in this study.
Tycova, Anna; Prikryl, Jan; Foret, Frantisek
2016-04-01
The use of high quality fused silica capillary nanospray tips is critical for obtaining reliable and reproducible electrospray/MS data; however, reproducible laboratory preparation of such tips is a challenging task. In this work, we report on the design and construction of low-cost grinding device assembled from 3D printed and commercially easily available components. Detailed description and characterization of the grinding device is complemented by freely accessible files in stl and skp format allowing easy laboratory replication of the device. The process of sharpening is aimed at achieving maximal symmetricity, surface smoothness and repeatability of the conus shape. Moreover, the presented grinding device brings possibility to fabricate the nanospray tips of desired dimensions regardless of the commercial availability. On several samples of biological nature (reserpine, rabbit plasma, and the mixture of three aminoacids), performance of fabricated tips is shown on CE coupled to MS analysis. The special interest is paid to the effect of tip sharpness. PMID:26626777
NASA Astrophysics Data System (ADS)
Boissier, J.; Bockelée-Morvan, D.; Rodionov, A. V.; Crifo, J.-F.
2010-02-01
Context. Millimetre line observations of comet C/1995 O1 (Hale-Bopp) close to perihelion, completed using the IRAM Plateau de Bure Interferometer, have detected temporal variations in the CO J(2-1) 230 GHz line shape, and in the position of its maximum emission brightness within the field-of-view, whose heuristic analysis has suggested the presence in the coma of a slowly rotating spiral-shaped enhancement of the CO density. Aims: Here, we reanalyse these data using a physically consistent model of the coma. Methods: We consider a large, rotating, icy nucleus with an arbitrarily aspherical shape and an adhoc rotation mode, and compute its tridimensional, time-dependent (”3D+t”) mixed CO + H2O coma, using a previously developed tridimensional hydrodynamical code (HDC). The line emission of CO is then computed using a molecular excitation and radiation transfer code (ERC). In the present, pioneering phase, the HDC and ERC both contain crude, and not thoroughly mutually consistent approximations. Several alternative CO surface flux distributions are considered, and the resulting CO 230 GHz line spectra and brightness maps are compared with observations. Results: We find that when an uniform surface flux of CO is assumed, the spiral structures created by the nucleus asphericity in the CO coma are too faint to account for the observational data, whereas we confirm earlier conclusions based on a heuristic approach that the assumption of an area of suitable dimensions and localization with increased CO flux leads to results in agreement with a large subset of (but not all) the data. This suggests that the true CO coma production map may be more complex than the presently assumed rather simple-minded one. Refined and mutually consistent HDC and ERC are needed for a more satisfactory interpretation of the present and any similar future data. Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported by INSU/CNRS (France), MPG
NASA Astrophysics Data System (ADS)
Zhou, Nianqing; Westrich, Bernhard; Jiang, Simin; Wang, Yan
2011-01-01
SummaryThe rapid economic development has caused heavy pollution in the Pearl River Delta Region (PRDR) and neighboring regions. In order to investigate the polluted characteristics, a coupled modelling approach has been deployed and the study area in the Pear River Delta is divided into two modelling systems: a river network that is simulated by a 1D hydrodynamic and water quality model, and an estuary that is modelled by a 3D model. The 1D river network model is coupled with the 3D estuary model based on the mass conservation criteria. The hydrodynamic (HD) model is based on 1D hydrostatic shallow water equations, and is solved by the Alternating Direction Implicit Approach. The water quality (WQ) model component is based on advection-diffusion equations. Some physical and bio-chemical processes are included in the WQ model to address the water pollution problems by different pollution substances, in particular, to the chemical oxygen demand (COD). The water quality data of the main rivers from 2004 to 2007 were collected, processed and evaluated. Because no detailed water quality data in the river is available, the COD source pollution intensity in the water quality model is taking the pollutant load from the confluence of the main five outlets of the Pearl River Estuary (PRE) for 1996-2007 as surrogate. The results show that the pollutants from the PRDR have a large impact on the quality of the Hong Kong seawater and much more serious in the wet season with high river discharge. An insight of the flow and transportation mechanism in the river and estuary system has been obtained through the paper, and a reliable basis to control water pollution of the Pearl River Delta is provided in this paper.
NASA Astrophysics Data System (ADS)
Otto, Christopher; Kempka, Thomas; Kapusta, Krzysztof; Stańczyk, Krzysztof
2016-04-01
Underground coal gasification (UCG) has the potential to increase the world-wide coal reserves by utilization of coal deposits not mineable by conventional methods. The UCG process involves combusting coal in situ to produce a high-calorific synthesis gas, which can be applied for electricity generation or chemical feedstock production. Apart from its high economic potentials, UCG may induce site-specific environmental impacts such as fault reactivation, induced seismicity and ground subsidence, potentially inducing groundwater pollution. Changes overburden hydraulic conductivity resulting from thermo-mechanical effects may introduce migration pathways for UCG contaminants. Due to the financial efforts associated with UCG field trials, numerical modeling has been an important methodology to study coupled processes considering UCG performance. Almost all previous UCG studies applied 1D or 2D models for that purpose, that do not allow to predict the performance of a commercial-scale UCG operation. Considering our previous findings, demonstrating that far-field models can be run at a higher computational efficiency by using temperature-independent thermo-mechanical parameters, representative coupled simulations based on complex 3D regional-scale models were employed in the present study. For that purpose, a coupled thermo-mechanical 3D model has been developed to investigate the environmental impacts of UCG based on a regional-scale of the Polish Wieczorek mine located in the Upper Silesian Coal Basin. The model size is 10 km × 10 km × 5 km with ten dipping lithological layers, a double fault and 25 UCG reactors. Six different numerical simulation scenarios were investigated, considering the transpressive stress regime present in that part of the Upper Silesian Coal Basin. Our simulation results demonstrate that the minimum distance between the UCG reactors is about the six-fold of the coal seam thickness to avoid hydraulic communication between the single UCG
NASA Astrophysics Data System (ADS)
Gorby, M.; Schwadron, N.; Lee, M. A.; Booth, A. C.; Spence, H.; Torok, T.; Downs, C.; Lionello, R.; Linker, J.; Titov, V. S.; Mikic, Z.; Riley, P.; Desai, M. I.; Dayeh, M. A.; Kozarev, K. A.
2013-12-01
Recent work on the coupling between 3D energetic particle models (e.g., the Energetic Particle Radiation Environment Model, EPREM) and magnetohydrodynamic (MHD) models of Coronal Mass Ejections (CMEs, e.g., the PSI MAS model) has demonstrated the efficacy of compression regions in front of fast CMEs for particle acceleration from remarkably low in the corona (3-6 solar radii). Typically particle acceleration becomes rapid beyond 3Rs and often in regions where compression regions have not yet formed active shocks. The challenge for forming large SEP events in such compression-acceleration scenarios is to have enhanced scattering within the acceleration region while also allowing for efficient escape of accelerated particles downstream (away from the Sun) from the compression region. Simulations show rapid particle acceleration in the range of 3-8 Rs over a broad longitudinal region (80°) resulting from the pile-up of magnetic flux in the compression and the subsequent expansion of these fields. These results have important implications for multi-instrument observations that will allow Solar Probe Plus and Solar Orbiter to test the developing paradigm of SEP acceleration and transport from coronal compressions. We present here recent coupled simulations for SEP acceleration and transport, including energetic particle and CME plasma profiles. The broadness of the longitudinal profile from such events may be a key observational test of compression acceleration in the low corona.
3D Numerical Simulation on Thermal Flow Coupling Field of Stainless Steel During Twin-Roll Casting
NASA Astrophysics Data System (ADS)
Liu, Lianlian; Liao, Bo; Guo, Jing; Liu, Ligang; Hu, Hongyan; Zhang, Yue; Yang, Qingxiang
2014-01-01
The surface crack and lateral crack of the AISI 304 stainless steel thin strip produced by twin-roll casting were observed. The temperature at the center of outlet during twin-roll-casting process was determined by infrared thermometer. In order to avoid the surface cracks of the casting strip, the thermal flow coupling field of AISI 304 stainless steel during twin-roll-casting process was simulated by a 3D fluid-structure coupling model. According to the simulation result, the effect of the casting speed on thermal flow field was analyzed and the process parameters were optimized. Moreover, by studying heat flux curves, the heat transfer mechanism between molten pool and roll was analyzed. The results show that, with the increase of the casting speed, the temperature of the molten pool increases and the solidification point moves toward the outlet. Meanwhile, the whirlpool above gets larger. Based on the solidification front position, the optimized process parameters are 1500 °C and 0.37 m/s. The heat transfer mechanism between molten pool and roll contains direct contacting heat transfer and air gap heat transfer.
SABAE-HW3D: a Meteor-Hydrological Model Coupling the Land Surface to Groundwater Flow
NASA Astrophysics Data System (ADS)
Loukili, Y.; Woodbury, A. D.
2007-12-01
A new coupled model linking the land surface scheme SABAE-HW with the saturated groundwater flow is introduced. SABAE-HW stands for: Soil Atmosphere Boundary, Accurate Evaluations of Heat and Water. It was presented (Loukili et al., 2007) as a parallel extension to the Canadian LAnd Surface Scheme (CLASS) which permits a free choice of soil column depth and layers. This development is important as it allows for the coupling of the atmosphere, land surface and subsurface zones. The generic domain under study is subdivided into soil columns whose surface areas represent GCM grid squares. The depths of columns extend to the saturated zone where a water table lower boundary condition is prescribed. In each column the unsaturated flow is forced by the corresponding meteorological data, and the moisture (liquid and frozen) and temperature profiles are computed by SABAE-HW using the half hour time step. The horizontal flow in the saturated zone is described by the vertically integrated model incorporating a storage term, and discretized using a finite volume scheme operating on the same GCM quadrilateral mesh. Since the saturated flow model affords larger time steps, bottom drainages from soil columns are summed up and input as cell recharges. Meanwhile, as the water table fluctuates, the individual column's mesh is allowed to deform. This physically based coupling strategy was selected for its superiority in providing stable and consistent results. In fact, the convergence to steady state situations in both the unsaturated and saturated zones is realized in many validation numerical tests. Moreover, our SABAE-HW3D code was successfully benchmarked against the finite element code Seep/W (Geo-Slope International, 2002) for steady and transient groundwater flows through different soil types. Even when handling coarse grids, SABAE-HW3D solutions are free of moisture oscillation and under or overshoot near the water table. The model is used to understand and assess properly
NASA Astrophysics Data System (ADS)
Janssen, G.; Del Val Alonso, L.; Groenendijk, P.; Griffioen, J.
2012-12-01
We developed an on-line coupling between the 1D/quasi-2D nutrient transport model ANIMO and the 3D groundwater transport model code MT3DMS. ANIMO is a detailed, process-oriented model code for the simulation of nitrate leaching to groundwater, N- and P-loads on surface waters and emissions of greenhouse gasses. It is the leading nutrient fate and transport code in the Netherlands where it is used primarily for the evaluation of fertilization related legislation. In addition, the code is applied frequently in international research projects. MT3DMS is probably the most commonly used groundwater solute transport package worldwide. The on-line model coupling ANIMO-MT3DMS combines the state-of-the-art descriptions of the biogeochemical cycles in ANIMO with the advantages of using a 3D approach for the transport through the saturated domain. These advantages include accounting for regional lateral transport, considering groundwater-surface water interactions more explicitly, and the possibility of using MODFLOW to obtain the flow fields. An additional merit of the on-line coupling concept is that it preserves feedbacks between the saturated and unsaturated zone. We tested ANIMO-MT3DMS by simulating nutrient transport for the period 1970-2007 in a Dutch agricultural polder catchment covering an area of 118 km2. The transient groundwater flow field had a temporal resolution of one day and was calculated with MODFLOW-MetaSWAP. The horizontal resolution of the model grid was 100x100m and consisted of 25 layers of varying thickness. To keep computation times manageable, we prepared MT3DMS for parallel computing, which in itself is a relevant development for a large community of groundwater transport modelers. For the parameterization of the soil, we applied a standard classification approach, representing the area by 60 units with unique combinations of soil type, land use and geohydrological setting. For the geochemical parameterization of the deeper subsurface, however, we
Jonkman, J. M.; Sclavounos, P. D.
2006-01-01
Aeroelastic simulation tools are routinely used to design and analyze onshore wind turbines, in order to obtain cost effective machines that achieve favorable performance while maintaining structural integrity. These tools employ sophisticated models of wind-inflow; aerodynamic, gravitational, and inertial loading of the rotor, nacelle, and tower; elastic effects within and between components; and mechanical actuation and electrical responses of the generator and of control and protection systems. For offshore wind turbines, additional models of the hydrodynamic loading in regular and irregular seas, the dynamic coupling between the support platform motions and wind turbine motions, and the dynamic characterization of mooring systems for compliant floating platforms are also important. Hydrodynamic loading includes contributions from hydrostatics, wave radiation, and wave scattering, including free surface memory effects. The integration of all of these models into comprehensive simulation tools, capable of modeling the fully coupled aeroelastic and hydrodynamic responses of floating offshore wind turbines, is presented.
Hydrodynamic coupling of particle inclusions embedded in curved lipid bilayer membranes.
Sigurdsson, Jon Karl; Atzberger, Paul J
2016-08-10
We develop theory and computational methods to investigate particle inclusions embedded within curved lipid bilayer membranes. We consider the case of spherical lipid vesicles where inclusion particles are coupled through (i) intramembrane hydrodynamics, (ii) traction stresses with the external and trapped solvent fluid, and (iii) intermonolayer slip between the two leaflets of the bilayer. We investigate relative to flat membranes how the membrane curvature and topology augment hydrodynamic responses. We show how both the translational and rotational mobility of protein inclusions are effected by the membrane curvature, ratio of intramembrane viscosity to solvent viscosity, and intermonolayer slip. For general investigations of many-particle dynamics, we also discuss how our approaches can be used to treat the collective diffusion and hydrodynamic coupling within spherical bilayers. PMID:27373277
NASA Astrophysics Data System (ADS)
Kataoka, J.; Kishimoto, A.; Nishiyama, T.; Fujita, T.; Takeuchi, K.; Kato, T.; Nakamori, T.; Ohsuka, S.; Nakamura, S.; Hirayanagi, M.; Adachi, S.; Uchiyama, T.; Yamamoto, K.
2013-12-01
The release of radioactive isotopes (mainly 137Cs, 134Cs and 131I) from the crippled Fukushima Daiichi Nuclear Plant remains a serious problem in Japan. To help identify radiation hotspots and ensure effective decontamination operation, we are developing a novel Compton camera weighting only 1 kg and measuring just ∼10 cm2 in size. Despite its compactness, the camera realizes a wide 180° field of vision with a sensitivity about 50 times superior to other cameras being tested in Fukushima. We expect that a hotspot producing a 5 μSv/h dose at a distance of 3 m can be imaged every 10 s, with angular resolution better than 10° (FWHM). The 3D position-sensitive scintillators and thin monolithic MPPC arrays are the key technologies developed here. By measuring the pulse-height ratio of MPPC-arrays coupled at both ends of a Ce:GAGG scintillator block, the depth of interaction (DOI) is obtained for incident gamma rays as well as the usual 2D positions, with accuracy better than 2 mm. By using two identical 10 mm cubic Ce:GAGG scintillators as a scatterer and an absorber, we confirmed that the 3D configuration works well as a high-resolution gamma camera, and also works as spectrometer achieving typical energy resolution of 9.8% (FWHM) for 662 keV gamma rays. We present the current status of the prototype camera (weighting 1.5 kg and measuring 8.5×14×16 cm3 in size) being fabricated by Hamamatsu Photonics K.K. Although the camera still operates in non-DOI mode, angular resolution as high as 14° (FWHM) was achieved with an integration time of 30 s for the assumed hotspot described above.
Coupling poroelasticity and CFD for cerebrospinal fluid hydrodynamics.
Tully, Brett; Ventikos, Yiannis
2009-06-01
This research uses a novel coupling of poroelastic theory and computational fluid dynamics to investigate acute hydrocephalus resulting from stenosis of the cerebral aqueduct. By coupling poroelastic theory with a multidimensional simulation of the cerebral aqueduct we are able to investigate, for the first time, the impact of physically relevant stenosis patterns on ventricular enlargement, accounting for the nonintuitive long time history responses of the ventricular system. Preliminary findings demonstrate clearly the importance that the fluidic-poroelastic coupling plays: ventricular enlargement is significantly smaller with local stenosis patterns and almost all of the observable pressure drop occurs across the stenosis. Short timescale effects [O(heartbeat)] are explored and their contribution to the long timescales interrogated. PMID:19304478
Development of a 3D Soil-Plant-Atmosphere Continuum (SPAC) coupled to a Land Surface Model
NASA Astrophysics Data System (ADS)
Bisht, G.; Riley, W. J.; Lorenzetti, D.; Tang, J.
2015-12-01
Exchange of water between the atmosphere and biosphere via evapotranspiration (ET) influences global hydrological, energy, and biogeochemical cycles. Isotopic analysis has shown that evapotranspiration over the continents is largely dominated by transpiration. Water is taken up from soil by plant roots, transported through the plant's vascular system, and evaporated from the leaves. Yet current Land Surface Models (LSMs) integrated into Earth System Models (ESMs) treat plant roots as passive components. These models distribute the ET sink vertically over the soil column, neglect the vertical pressure distribution along the plant vascular system, and assume that leaves can directly access water from any soil layer within the root zone. Numerous studies have suggested that increased warming due to climate change will lead drought and heat-induced tree mortality. A more mechanistic treatment of water dynamics in the soil-plant-atmosphere continuum (SPAC) is essential for investigating the fate of ecosystems under a warmer climate. In this work, we describe a 3D SPAC model that can be coupled to a LSM. The SPAC model uses the variably saturated Richards equations to simulate water transport. The model uses individual governing equations and constitutive relationships for the various SPAC components (i.e., soil, root, and xylem). Finite volume spatial discretization and backward Euler temporal discretization is used to solve the SPAC model. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is used to numerically integrate the discretized system of equations. Furthermore, PETSc's multi-physics coupling capability (DMComposite) is used to solve the tightly coupled system of equations of the SPAC model. Numerical results are presented for multiple test problems.
Coupling and Hydrodynamic Limit for the Inclusion Process
NASA Astrophysics Data System (ADS)
Opoku, Alex; Redig, Frank
2015-08-01
We show propagation of local equilibrium for the symmetric inclusion process (SIP) after diffusive rescaling of space and time, as well as the local equilibrium property of the non-equilibrium steady state in the boundary driven SIP. The main tool is self-duality and a coupling between SIP particles and independent random walkers.
Fraysse, Marion; Pinazo, Christel; Faure, Vincent Martin; Fuchs, Rosalie; Lazzari, Paolo; Raimbault, Patrick; Pairaud, Ivane
2013-01-01
Terrestrial inputs (natural and anthropogenic) from rivers, the atmosphere and physical processes strongly impact the functioning of coastal pelagic ecosystems. The objective of this study was to develop a tool for the examination of these impacts on the Marseille coastal area, which experiences inputs from the Rhone River and high rates of atmospheric deposition. Therefore, a new 3D coupled physical/biogeochemical model was developed. Two versions of the biogeochemical model were tested, one model considering only the carbon (C) and nitrogen (N) cycles and a second model that also considers the phosphorus (P) cycle. Realistic simulations were performed for a period of 5 years (2007–2011). The model accuracy assessment showed that both versions of the model were able of capturing the seasonal changes and spatial characteristics of the ecosystem. The model also reproduced upwelling events and the intrusion of Rhone River water into the Bay of Marseille well. Those processes appeared to greatly impact this coastal oligotrophic area because they induced strong increases in chlorophyll-a concentrations in the surface layer. The model with the C, N and P cycles better reproduced the chlorophyll-a concentrations at the surface than did the model without the P cycle, especially for the Rhone River water. Nevertheless, the chlorophyll-a concentrations at depth were better represented by the model without the P cycle. Therefore, the complexity of the biogeochemical model introduced errors into the model results, but it also improved model results during specific events. Finally, this study suggested that in coastal oligotrophic areas, improvements in the description and quantification of the hydrodynamics and the terrestrial inputs should be preferred over increasing the complexity of the biogeochemical model. PMID:24324589
Running Faster Together: Huge Speed up of Thermal Ratchets due to Hydrodynamic Coupling
NASA Astrophysics Data System (ADS)
Malgaretti, Paolo; Pagonabarraga, Ignacio; Frenkel, Daan
2012-10-01
We present simulations that reveal a surprisingly large effect of hydrodynamic coupling on the speed of thermal ratchet motors. The model that we use considers particles performing thermal ratchet motion in a hydrodynamic solvent. Using particle-based, mesoscopic simulations that maintain local momentum conservation, we analyze quantitatively how the coupling to the surrounding fluid affects ratchet motion. We find that coupling can increase the mean velocity of the moving particles by almost 2 orders of magnitude, precisely because ratchet motion has both a diffusive and a deterministic component. The resulting coupling also leads to the formation of aggregates at longer times. The correlated motion that we describe increases the efficiency of motor-delivered cargo transport and we speculate that the mechanism that we have uncovered may play a key role in speeding up molecular motor-driven intracellular transport.
Liu, Cai-Ming; Zhang, De-Qing; Zhu, Dao-Ben
2016-04-01
A novel 3D MOF based on dysprosium(iii) oxalate and 1,10-phenanthroline (phen), {[Dy(C2O4)1.5phen]·0.5H2O}n (1), has been hydrothermally synthesized. The Dy(3+) ion acts as a typical Y-shaped node, linking to each other to generate an interesting 3D topology structure. Complex 1 is the first 3D DyMOF displaying both ferromagnetic coupling and field-induced two-step magnetic relaxation. PMID:26961387
The CO2 system in the Mediterranean Sea inferred from a 3D coupled physical-biogeochemical model
NASA Astrophysics Data System (ADS)
Ulses, Caroline; Kessouri, Fayçal; Estournel, Claude; Marsaleix, Patrick; Beuvier, Jonathan; Somot, Samuel; Touratier, Frank; Goyet, Catherine; Coppola, Laurent; Diamond, Emilie; Metzl, Nicolas
2015-04-01
The semi-enclosed Mediterranean Sea characterized by short residence times is considered as a region particularly sensitive to natural and anthropogenic forcing. Due to scarce CO2 measurements in the whole basin, the CO2 system, for instance the air-sea CO2 exchanges and the effects of the increase of atmospheric CO2, are poorly characterized. 3D physical-biogeochemical coupled models are unique tools that can provide integrated view and gain understanding in the temporal and spatial variation of the CO2 system variables (dissolved inorganic carbon, total alkalinity, partial pressure of CO2 and pH). An extended version of the biogeochemical model Eco3m-S (Auger et al., 2014), that describes the cycles of carbon, nitrogen, phosphorus and silica, was forced by a regional circulation model (Beuvier et al., 2012) to investigate the CO2 system in the Mediterranean Sea over a 13-years period (2001-2013). First, the quality of the modelling was evaluated through comparisons with satellite and in situ observations collected in the whole basin over the study period (Touratier and Goyet, 2009; 2011 ; Rivaro et al., 2010 ; Pujo-Pay et al., 2011 ; Alvarez et al, 2014). The model reasonably reproduced the various biological regimes (north-western phytoplanctonic bloom regime, oligotrophic eastern regime, etc.) as well as the recorded spatial distribution and temporal variations of the carbonate system variables. The coupled model was then used to estimate the air-sea pCO2 exchanges and the transport of DIC and TA towards the Atlantic Ocean at the Strait of Gibraltar.
Kohn-Sham Band Structure Benchmark Including Spin-Orbit Coupling for 2D and 3D Solids
NASA Astrophysics Data System (ADS)
Huhn, William; Blum, Volker
2015-03-01
Accurate electronic band structures serve as a primary indicator of the suitability of a material for a given application, e.g., as electronic or catalytic materials. Computed band structures, however, are subject to a host of approximations, some of which are more obvious (e.g., the treatment of the exchange-correlation of self-energy) and others less obvious (e.g., the treatment of core, semicore, or valence electrons, handling of relativistic effects, or the accuracy of the underlying basis set used). We here provide a set of accurate Kohn-Sham band structure benchmarks, using the numeric atom-centered all-electron electronic structure code FHI-aims combined with the ``traditional'' PBE functional and the hybrid HSE functional, to calculate core, valence, and low-lying conduction bands of a set of 2D and 3D materials. Benchmarks are provided with and without effects of spin-orbit coupling, using quasi-degenerate perturbation theory to predict spin-orbit splittings. This work is funded by Fritz-Haber-Institut der Max-Planck-Gesellschaft.
Neiman, Jaclyn A Shepard; Raman, Ritu; Chan, Vincent; Rhoads, Mary G; Raredon, Micha Sam B; Velazquez, Jeremy J; Dyer, Rachel L; Bashir, Rashid; Hammond, Paula T; Griffith, Linda G
2015-04-01
In vitro models that recapitulate the liver's structural and functional complexity could prolong hepatocellular viability and function to improve platforms for drug toxicity studies and understanding liver pathophysiology. Here, stereolithography (SLA) was employed to fabricate hydrogel scaffolds with open channels designed for post-seeding and perfused culture of primary hepatocytes that form 3D structures in a bioreactor. Photopolymerizable polyethylene glycol-based hydrogels were fabricated coupled to chemically activated, commercially available filters (polycarbonate and polyvinylidene fluoride) using a chemistry that permitted cell viability, and was robust enough to withstand perfused culture of up to 1 µL/s for at least 7 days. SLA energy dose, photoinitiator concentrations, and pretreatment conditions were screened to determine conditions that maximized cell viability and hydrogel bonding to the filter. Multiple open channel geometries were readily achieved, and included ellipses and rectangles. Rectangular open channels employed for subsequent studies had final dimensions on the order of 350 µm by 850 µm. Cell seeding densities and flow rates that promoted cell viability were determined. Perfused culture of primary hepatocytes in hydrogel scaffolds in the presence of soluble epidermal growth factor (EGF) prolonged the maintenance of albumin production throughout the 7-day culture relative to 2D controls. This technique of bonding hydrogel scaffolds can be employed to fabricate soft scaffolds for a number of bioreactor configurations and applications. PMID:25384798
Pannala, S; D'Azevedo, E; Zacharia, T
2002-02-26
The goal of the radiation modeling effort was to develop and implement a radiation algorithm that is fast and accurate for the underhood environment. As part of this CRADA, a net-radiation model was chosen to simulate radiative heat transfer in an underhood of a car. The assumptions (diffuse-gray and uniform radiative properties in each element) reduce the problem tremendously and all the view factors for radiation thermal calculations can be calculated once and for all at the beginning of the simulation. The cost for online integration of heat exchanges due to radiation is found to be less than 15% of the baseline CHAD code and thus very manageable. The off-line view factor calculation is constructed to be very modular and has been completely integrated to read CHAD grid files and the output from this code can be read into the latest version of CHAD. Further integration has to be performed to accomplish the same with STAR-CD. The main outcome of this effort is to obtain a highly scalable and portable simulation capability to model view factors for underhood environment (for e.g. a view factor calculation which took 14 hours on a single processor only took 14 minutes on 64 processors). The code has also been validated using a simple test case where analytical solutions are available. This simulation capability gives underhood designers in the automotive companies the ability to account for thermal radiation - which usually is critical in the underhood environment and also turns out to be one of the most computationally expensive components of underhood simulations. This report starts off with the original work plan as elucidated in the proposal in section B. This is followed by Technical work plan to accomplish the goals of the project in section C. In section D, background to the current work is provided with references to the previous efforts this project leverages on. The results are discussed in section 1E. This report ends with conclusions and future scope of work in section F.
Hydrodynamic coupling of two sharp-edged beams vibrating in a viscous fluid
Intartaglia, Carmela; Soria, Leonardo; Porfiri, Maurizio
2014-01-01
In this paper, we study flexural vibrations of two thin beams that are coupled through an otherwise quiescent viscous fluid. While most of the research has focused on isolated beams immersed in placid fluids, inertial and viscous hydrodynamic coupling is ubiquitous across a multitude of engineering and natural systems comprising arrays of flexible structures. In these cases, the distributed hydrodynamic loading experienced by each oscillating structure is not only related to its absolute motion but is also influenced by its relative motion with respect to the neighbouring structures. Here, we focus on linear vibrations of two identical beams for low Knudsen, Keulegan–Carpenter and squeeze numbers. Thus, we describe the fluid flow using unsteady Stokes hydrodynamics and we propose a boundary integral formulation to compute pertinent hydrodynamic functions to study the fluid effect. We validate the proposed theoretical approach through experiments on centimetre-size compliant cantilevers that are subjected to underwater base-excitation. We consider different geometric arrangements, beam interdistances and excitation frequencies to ascertain the model accuracy in terms of the relevant non-dimensional parameters. PMID:24511249
NASA Astrophysics Data System (ADS)
Niasse, Nicolas; Chittenden, Jeremy
2012-10-01
The last few years have seen considerable advances in the application of high performance computing techniques to 3D simulations of wire array Z-pinches. Whilst the intense soft X-ray radiation output is the principle application of wire arrays, the ability to encompass spectrally detailed models of this emission within such 3D calculations was thought to be computationally prohibitive. We have developed a non-LTE atomic and radiation physics model with detailed configuration accounting and n-l splitting which is sufficiently streamlined to run in-line with large scale 3D simulations. In order to handle the volume of data generated by the spectral treatment of the billions of numerical cells, a novel data structure derived from a self-balancing binary search tree was developed, enabling the use of non-LTE DCA calculations within large scale 3D simulations for the first time. A brief description of the model is provided and the application of the simulations to understanding the X-ray generation processes within wire array Z-pinches on the Z generator at Sandia National Laboratory is reported. The contribution of the ion temperature and the motion of the unstable plasma at stagnation to the Doppler widths of the lines is described in detail.
Modeling the coupling of reaction kinetics and hydrodynamics in a collapsing cavity
Mishra, Sudib; Deymier, Pierre; Muralidharan, Krishna; Frantziskonis, G.; Pannala, Sreekanth; Simunovic, Srdjan
2010-01-01
We introduce a model of cavitation based on the multiphase Lattice Boltzmann method (LBM) that allows for coupling between the hydrodynamics of a collapsing cavity and supported solute chemical species. We demonstrate that this model can also be coupled to deterministic or stochastic chemical reactions. In a two-species model of chemical reactions (with a major and a minor specie), the major difference observed between the deterministic and stochastic reactions takes the form of random fluctuations in concentration of the minor species. We demonstrate that advection associated with the hydrodynamics of a collapsing cavity leads to highly inhomogeneous concentration of solutes. In turn these inhomogeneities in concentration may lead to significant increase in concentration-dependent reaction rates and can result in a local enhancement in the production of minor species.
NASA Astrophysics Data System (ADS)
Tang, Xian-Zhu; McDevitt, C. J.; Guo, Zehua; Berk, H. L.
2014-03-01
Inertial confinement fusion requires an imploded target in which a central hot spot is surrounded by a cold and dense pusher. The hot spot/pusher interface can take complicated shape in three dimensions due to hydrodynamic mix. It is also a transition region where the Knudsen and inverse Knudsen layer effect can significantly modify the fusion reactivity in comparison with the commonly used value evaluated with background Maxwellians. Here, we describe a hybrid model that couples the kinetic correction of fusion reactivity to global hydrodynamic implosion simulations. The key ingredient is a non-perturbative treatment of the tail ions in the interface region where the Gamow ion Knudsen number approaches or surpasses order unity. The accuracy of the coupling scheme is controlled by the precise criteria for matching the non-perturbative kinetic model to perturbative solutions in both configuration space and velocity space.
Dynamically Coupled Food-web and Hydrodynamic Modeling with ADH-CASM
NASA Astrophysics Data System (ADS)
Piercy, C.; Swannack, T. M.
2012-12-01
Oysters and freshwater mussels are "ecological engineers," modifying the local water quality by filtering zooplankton and other suspended particulate matter from the water column and flow hydraulics by impinging on the near-bed flow environment. The success of sessile, benthic invertebrates such as oysters depends on environmental factors including but not limited to temperature, salinity, and flow regime. Typically food-web and other types of ecological models use flow and water quality data as direct input without regard to the feedback between the ecosystem and the physical environment. The USACE-ERDC has developed a coupled hydrodynamic-ecological modeling approach that dynamically couples a 2-D hydrodynamic and constituent transport model, Adaptive Hydraulics (ADH), with a bioenergetics food-web model, the Comprehensive Aquatics Systems Model (CASM), which captures the dynamic feedback between aquatic ecological systems and the environment. We present modeling results from restored oyster reefs in the Great Wicomico River on the western shore of the Chesapeake Bay, which quantify ecosystem services such as the influence of the benthic ecosystem on water quality. Preliminary results indicate that while the influence of oyster reefs on bulk flow dynamics is limited due to the localized influence of oyster reefs, large reefs and the associated benthic ecosystem can create measurable changes in the concentrations of nitrogen, phosphorus, and carbon in the areas around reefs. We also present a sensitivity analysis to quantify the relative sensitivity of the coupled ADH-CASM model to both hydrodynamic and ecological parameter choice.
NASA Astrophysics Data System (ADS)
Gutknecht, E.; Dadou, I.; Marchesiello, P.; Cambon, G.; Le Vu, B.; Sudre, J.; Garçon, V.; Machu, E.; Rixen, T.; Kock, A.; Flohr, A.; Paulmier, A.; Lavik, G.
2013-06-01
Eastern boundary upwelling systems (EBUS) are regions of high primary production often associated with oxygen minimum zones (OMZs). They represent key regions for the oceanic nitrogen (N) cycle. By exporting organic matter (OM) and nutrients produced in the coastal region to the open ocean, EBUS can play an important role in sustaining primary production in subtropical gyres. However, losses of fixed inorganic N through denitrification and anammox processes take place in oxygen depleted environments such as EBUS, and can potentially mitigate the role of these regions as a source of N to the open ocean. EBUS can also represent a considerable source of nitrous oxide (N2O) to the atmosphere, affecting the atmospheric budget of N2O. In this paper a 3-D coupled physical/biogeochemical model (ROMS/BioEBUS) is used to investigate the N budget in the Namibian upwelling system. The main processes linked to EBUS and associated OMZs are taken into account. The study focuses on the northern part of the Benguela upwelling system (BUS), especially the Walvis Bay area (between 22° S and 24° S) where the OMZ is well developed. Fluxes of N off the Walvis Bay area are estimated in order to understand and quantify (1) the total N offshore export from the upwelling area, representing a possible N source that sustains primary production in the South Atlantic subtropical gyre; (2) export production and subsequent losses of fixed N via denitrification and anammox under suboxic conditions (O2 < 25 mmol O2 m-3); and (3) the N2O emission to the atmosphere in the upwelling area. In the mixed layer, the total N offshore export is estimated as 8.5 ± 3.9 × 1010 mol N yr-1 at 10° E off the Walvis Bay area, with a mesoscale contribution of 20%. Extrapolated to the whole BUS, the coastal N source for the subtropical gyre corresponds to 0.1 ± 0.04 mol N m-2 yr-1. This N flux represents a major source of N for the gyre compared with other N sources, and contributes 28% of the new primary
NASA Astrophysics Data System (ADS)
Balabanov, Nikolai B.; Peterson, Kirk A.
2006-08-01
Recently developed correlation consistent basis sets for the first row transition metal elements Sc-Zn have been utilized to determine complete basis set (CBS) scalar relativistic electron affinities, ionization potentials, and 4s23dn -2-4s1dn -1 electronic excitation energies with single reference coupled cluster methods [CCSD(T), CCSDT, and CCSDTQ] and multireference configuration interaction with three reference spaces: 3d4s, 3d4s4p, and 3d4s4p3d'. The theoretical values calculated with the highest order coupled cluster techniques at the CBS limit, including extrapolations to full configuration interaction, are well within 1kcal/mol of the corresponding experimental data. For the early transition metal elements (Sc-Mn) the internally contracted multireference averaged coupled pair functional method yielded excellent agreement with experiment; however, the atomic properties for the late transition metals (Mn-Zn) proved to be much more difficult to describe with this level of theory, even with the largest reference function of the present work.
Coupled 2D-3D finite element method for analysis of a skin panel with a discontinuous stiffener
NASA Technical Reports Server (NTRS)
Wang, J. T.; Lotts, C. G.; Davis, D. D., Jr.; Krishnamurthy, T.
1992-01-01
This paper describes a computationally efficient analysis method which was used to predict detailed stress states in a typical composite compression panel with a discontinuous hat stiffener. A global-local approach was used. The global model incorporated both 2D shell and 3D brick elements connected by newly developed transition elements. Most of the panel was modeled with 2D elements, while 3D elements were employed to model the stiffener flange and the adjacent skin. Both linear and geometrically nonlinear analyses were performed on the global model. The effect of geometric nonlinearity induced by the eccentric load path due to the discontinuous hat stiffener was significant. The local model used a fine mesh of 3D brick elements to model the region at the end of the stiffener. Boundary conditions of the local 3D model were obtained by spline interpolation of the nodal displacements from the global analysis. Detailed in-plane and through-the-thickness stresses were calculated in the flange-skin interface near the end of the stiffener.
Garrevoet, Jan; Vekemans, Bart; Tack, Pieter; De Samber, Björn; Schmitz, Sylvia; Brenker, Frank E; Falkenberg, Gerald; Vincze, Laszlo
2014-12-01
A new three-dimensional (3D) micro X-ray fluorescence (μXRF) methodology based on a novel 2D energy dispersive CCD detector has been developed and evaluated at the P06 beamline of the Petra-III storage ring (DESY) in Hamburg, Germany. This method is based on the illumination of the investigated sample cross-section by a horizontally focused beam (vertical sheet beam) while fluorescent X-rays are detected perpendicularly to the sheet beam by a 2D energy dispersive (ED) CCD detector allowing the collection of 2D cross-sectional elemental images of a certain depth within the sample, limited only by signal self-absorption effects. 3D elemental information is obtained by a linear scan of the sample in the horizontal direction across the vertically oriented sheet beam and combining the detected cross-sectional images into a 3D elemental distribution data set. Results of the 3D μXRF analysis of mineral inclusions in natural deep Earth diamonds are presented to illustrate this new methodology. PMID:25346101
Coherent dynamic structure factors of strongly coupled plasmas: A generalized hydrodynamic approach
NASA Astrophysics Data System (ADS)
Luo, Di; Zhao, Bin; Hu, GuangYue; Gong, Tao; Xia, YuQing; Zheng, Jian
2016-05-01
A generalized hydrodynamic fluctuation model is proposed to simplify the calculation of the dynamic structure factor S(ω, k) of non-ideal plasmas using the fluctuation-dissipation theorem. In this model, the kinetic and correlation effects are both included in hydrodynamic coefficients, which are considered as functions of the coupling strength (Γ) and collision parameter (kλei), where λei is the electron-ion mean free path. A particle-particle particle-mesh molecular dynamics simulation code is also developed to simulate the dynamic structure factors, which are used to benchmark the calculation of our model. A good agreement between the two different approaches confirms the reliability of our model.
NASA Astrophysics Data System (ADS)
Xu, H.; Rodgers, A.; Lomov, I.; Petersson, A.; Sjogreen, B.; Vorobiev, O.; Chipman, V.
2011-12-01
We report research being performed to improve underground nuclear explosion (UNE) monitoring by developing capabilities for hydrodynamic modeling of ground motions. This effort involves work along two thrusts: 1) we are coupling hydrodynamic (non-linear shock) and seismic (linear anelastic) wave propagation codes; and 2) we are investigating the effect of source emplacement conditions on ground motions in the near field due to nonlinearity and comparing with the empirical models. For both thrusts we are modeling explosion motions using GEODYN, a fully three-dimensional Eulerian hydrodynamic code developed at LLNL. This code incorporates many important features for modeling shock waves in geologic materials, including non-linear response (e.g. porosity, tensile failure, yielding), topography, gravity, 3D material heterogeneities and adaptive mesh refinement. The calculation accuracy is well validated with the analytical solutions to the Lamb's problem and to the finite dilatational volume source at depth in a linear elastic medium. In order to propagate full waveform solutions from hydrodynamic simulations to distances where seismic measurements are made we are coupling GEODYN to WPP (LLNL's anelastic finite difference code for seismic wave simulation). Complex motions computed by GEODYN for explosions are recorded on a dense grid spanning the ranges where motions become linear (elastic). These wavefield records are processed and embedded into the WPP domain where they are introduced as a boundary driving source and continue to propagate as elastic waves at much lower numerical cost than with nonlinear GEODYN. The coupling scheme is validated by comparing the analytical, direct GEODYN solutions and WPP solutions to the finite dilatational volume source at depth in a linear elastic medium, and also by comparing the direct GEODYN solutions and WPP solutions to a complex 1kt chemical explosion in nonlinear granite at stations beyond the elastic radius. The excellent
Fullwave coupling to a 3D antenna code using Green's function formulation of wave-particle response
NASA Astrophysics Data System (ADS)
Wright, John; Bonoli, P. T.; Bilato, R.; Brambilla, M.; Maggiora, R.; Lancellotti, V.
2006-10-01
Using the fullwave code, TORIC, and the 3D antenna code, TOPICA, we construct a complete linear system for the RF driven plasma. The 3D finite element antenna code, TOPICA, requires an admittance, Y, for the plasma, where B=YE. In this work, TORIC was modified to allow excitation of the (Eη, Eζ) electric field components at the plasma surface, corresponding to a single poloidal and toroidal mode number combination (m,n). This leads to the tensor response: Yn= ( ll Yηη& YηζYζη& Yζζ), where each of the Yn submatrices is Nm in size. It is shown that the admittance matrix is equivalent to a Green's function calculation for the fullwave system and the net work done is less than twice a single fullwave calculation. The admittance calculation is used with loading calculation from TOPICA to construct self consistent plasma and antenna currents.
NASA Astrophysics Data System (ADS)
Lee, H.; Min, D.; Lim, S.; Yang, J.; Kwon, B.; Yoo, H.
2009-12-01
In a conventional marine seismic data analysis, pressure data have been usually interpreted on the basis of acoustic wave equation. The acoustic wave equation, however, only deals with P-wave propagation, and it cannot correctly describe the wave propagation in acoustic-elastic (fluid-solid) coupled media. Recently, in 4C OBC survey (4-component ocean bottom cable), it is possible to acquire both pressure and 3-component displacements (measured at the sea-bottom). Combining pressure and displacement data allows us to interpret subsurface structures more accurately. In order to accurately simulate wave propagation in fluid-solid coupled media, we need an acoustic-elastic coupled modeling algorithm, which deals with displacements in elastic region and pressure in acoustic region. For waveform inversion and reverse-time migration that require a great number of forward modeling, it is essential to develop an efficient scheme that reduces computing time and computer core memory. In this study, we present a 3D time-domain acoustic-elastic coupled modeling algorithm on the basis of the cell-based finite difference method. The cell-based method has proven to properly describe the free-surface boundary, which indicates that it will also properly describe the fluid-solid interface boundaries. In the acoustic-elastic coupled modeling, we first compose cell-based finite differences individually for the 3D acoustic and elastic media, and then combine the differences using the fluid-solid interface boundary conditions. Considering that the 2D acoustic-elastic coupled modeling algorithm gives numerical solutions comparable to analytic solutions, we expect that the 3D acoustic-elastic coupled modeling will correctly describe wave propagation in the fluid-solid coupled media. We apply our algorithm to 3D horizontal two- and three-layer models. Numerical experiments show that the cell-based coupled modeling algorithm properly describes S- and converted waves as well as P-waves. The
Integrated Biogeomorphological Modeling Using Delft3D
NASA Astrophysics Data System (ADS)
Ye, Q.; Jagers, B.
2011-12-01
The skill of numerical morphological models has improved significantly from the early 2D uniform, total load sediment models (with steady state or infrequent wave updates) to recent 3D hydrodynamic models with multiple suspended and bed load sediment fractions and bed stratigraphy (online coupled with waves). Although there remain many open questions within this combined field of hydro- and morphodynamics, we observe an increasing need to include biological processes in the overall dynamics. In riverine and inter-tidal environments, there is often an important influence by riparian vegetation and macrobenthos. Over the past decade more and more researchers have started to extend the simulation environment with wrapper scripts and other quick code hacks to estimate their influence on morphological development in coastal, estuarine and riverine environments. Although one can in this way quickly analyze different approaches, these research tools have generally not been designed with reuse, performance and portability in mind. We have now implemented a reusable, flexible, and efficient two-way link between the Delft3D open source framework for hydrodynamics, waves and morphology, and the water quality and ecology modules. The same link will be used for 1D, 2D and 3D modeling on networks and both structured and unstructured grids. We will describe the concepts of the overall system, and illustrate it with some first results.
Coupled discrete element and smoothed particle hydrodynamics simulations of the die filling process
NASA Astrophysics Data System (ADS)
Breinlinger, Thomas; Kraft, Torsten
2015-08-01
Die filling is an important part of the powder compaction process chain, where defects in the final part can be introduced—or prevented. Simulation of this process is therefore a goal for many part producers and has been studied by some researchers already. In this work, we focus on the influence of the surrounding air on the powder flow. We demonstrate the implementing and coupling of the discrete element method for the granular powder and the smoothed particle hydrodynamics method for the gas flow. Application of the method to the die filling process is demonstrated.
Bérut, A; Imparato, A; Petrosyan, A; Ciliberto, S
2016-02-12
We experimentally study the statistical properties of the energy fluxes between two trapped Brownian particles, interacting through dissipative hydrodynamic coupling, and submitted to an effective temperature difference ΔT, obtained by random forcing the position of one trap. We identify effective heat fluxes between the two particles and show that they satisfy an exchange fluctuation theorem in the stationary state. We also show that after the sudden application of a temperature gradient ΔT, the total hot-cold flux satisfies a transient exchange fluctuation theorem for any integration time, whereas the total cold-hot flux only does it asymptotically for long times. PMID:26919017
Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion
NASA Technical Reports Server (NTRS)
Margolis, S. B.
1999-01-01
For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic
NASA Astrophysics Data System (ADS)
Nguyen, Liliane; Guerin, Antoine; Abellán, Antonio; Carrea, Dario; Derron, Marc-Henri; Jaboyedoff, Michel
2015-04-01
Multiple sources of geological data exist nowadays, most of them are in 2D (e.g. geological maps, geological panoramic sketch), and only a few are in 3D (e.g. block diagram). One of the current challenges in geological mapping consists not only in providing a more consistent 3D data, but also in pursuing a gathering and a harmonisation of the geological information in order to obtain a more consistent interpretations of the 3D geological models. New remote sensing techniques have significantly improved the representation of three-dimensional surfaces during the last decade, especially for steep and inaccessible rockcliffs. Therefore, we present an exploratory study aiming to find a reliable method for carrying out a three-dimensional mapping of geological contacts using a High Resolution Digital Elevation Model (HRDEM) with a 1 meter cell size. To this end, we selected the "Scex Rouge Mountain" as pilot study area. This outcrop, which is located in the Diablerets Massif (Vaud, Swiss Alps), has the particularity to present very distinguishable folded geological boundaries on its Southern face. The Southern slope belongs to the Wildhorn nappe, which is mainly composed of sedimentary rocks. The top-layer is composed of siliceous limestones, the well-visible fold layer is the "Pygurus layer" and consist of sandy limestone. Finally the bottom-layer includes marly schist and clayey limestones. At first, different sources of information has been draped on the HRDEM of the Scex Rouge Mountain, including not only classical geological maps (1:25 000) but also different sources of imagery (e.g. gigapixel panoramas, open access images, etc.). In a second step, several three-dimensional polylines have been drawn following the geological limit on each drapped HRDEM. Then we investigated the accuracy of 2D classical geological maps by comparing these geological limits with the drawn 3D polylines. Furthermore, in order to evaluate the accuracy of the method, a ground truth needs
Hines, J.V.; Landry, S.M.; Varani, G.; Tinoco, I. Jr. Lawrence Berkeley Lab., CA )
1994-06-29
Long range carbon-proton scalar couplings were measured for an RNA hairpin of 12 nucleotides using 3D and [sup 13]C-edited 2D NMR. The large one-bond carbon-proton scalar couplings ([sup 1]J[sub CH]) and small n-bond couplings ([sup 1]J[sub CH]) produce ECOSY type cross-peaks, thus facilitating the determination of the sign and magnitude of the smaller [sup 2]J[sub CH] or [sup 3]J[sub CH]. The UUCGRNA hairpin (5[prime]-rGGACUUCGGUCC-3[prime]), whose structure has been determined by our laboratory, was uniformly [sup 13]C-labeled at 30% isotopic enrichment. The observed [sup 1]J[sub CH] couplings were then correlated to the known structure. The signs of [sup 2]J[sub C4[prime]H5[prime
Fullwave coupling to a 3D antenna code using Green's function formulation of wave-particle response.
NASA Astrophysics Data System (ADS)
Wright, John; Bonoli, Paul; Brambilla, Marco; Lancelloti, Vito; Maggiora, Riccardo; Carter, Mark
2006-04-01
Using the fullwave code, TORIC ,and the 3D antenna code, TOPICA, we construct a complete linear system for the RF driven plasma. The 3D finite element antenna code, TOPICA, requires an admittance, Y, for the plasma, where B=YE. In this work TORIC was modified to allow excitation of the (Eη, Eζ) electric field components at the plasma surface, corresponding to a single poloidal and toroidal mode number combination (m,n). This leads the tensor response: Y=( *20c Yηη & Yηζ Yζη & Yζζ ), where each of the Yn submatrices is Nm in size. It is shown that the admittance matrix is equivalent to a Greens function calculation for the fullwave system and in addition, the net work done in the calculation is on the order of twice a single fullwave calculation. After the admittance calculation is done, the response of a plasma to an antenna driven at a given frequency can be calculated by only running the TOPICA code for a new antenna geometry. In tests of loading, TOPICA has been able reproduce loading of the Alcator D antenna (S12 coefficient accurately.).
3D Dynamics of Magnetic Flux Ropes Across Scales: Solar Eruptions and Sun-Earth Plasma Coupling
NASA Astrophysics Data System (ADS)
Chen, James
2012-10-01
Central to the understanding of the eruptive phenomena on the Sun and their impact on the terrestrial plasma environment is the dynamics of coronal mass ejections (CMEs)---a 3D magnetic flux rope configuration---and the evolution of their magnetic fields. I will discuss the basic physics of CME eruption and solar flare energy release in the context of the analytic erupting flux rope model of CMEs. In this ideal MHD model, a CME is treated as a 3D flux rope with its two stationary footpoints anchored in the Sun. The model structure is non-axisymmetric and embedded in a model corona/solar wind. The initial flux rope is driven out of equilibrium by ``injection'' of poloidal flux and propagates under the Lorentz hoop force from the Sun to 1 AU, across a wide range of spatial and temporal scales. Comparisons of the model results and recent STEREO observations show that the solutions that best fit the observed CME position-time data (to within 1-2% of data) also correctly replicate the temporal profiles of associated flare X-ray emissions (GOES data) and the in situ magnetic field and plasma data of the CME ejecta at 1 AU where such data are available (e.g., ACE and STEREO/IMPAXCT/PLASTIC data), providing a unified basis of understanding CME dynamics and flare energetics.
NASA Technical Reports Server (NTRS)
Ponchak, G. E.; Bushyager, N.; Papapolymerou, J.; Tentzeris, E. M.; Laskar, J.
2002-01-01
Three-dimensional circuits built upon multiple layers of polyimide are required for constructing Si/SiGe monolithic microwave/mm-wave integrated circuits on CMOS (low resistivity) Si wafers. It is expected that these circuits will replace the ones fabricated on GaAs and reduce the overall system cost. However, the closely spaced transmission lines that are required for a high-density circuit environment are susceptible to high levels of cross-coupling, which degrades the overall circuit performance. In this paper, theoretical and experimental results on coupling and ways to reduce it are presented for two types of transmission lines: a) the microstrip line and b) the Finite Ground Coplanar (FGC) line. For microstrip lines it is shown that a fence of metalized via-holes can significantly reduce coupling, especially in the case when both lines are on the same polyimide layer or when the shielding structure extends through several polyimide layers. For closely spaced microstrip lines, coupling is lower for a metal filled trench shield than a via-hole fence. Coupling amongst microstrip lines is dependent on the ratio of line separation to polyimide thickness and is primarily due to magnetic fields. For FGC lines it is shown that they have in general low coupling that can be reduced significantly when there is even a small gap between the ground planes of each line. FGC lines have approximately 8 dB lower coupling than coupled coplanar waveguides (CPW). In addition, forward and backward characteristics of the FGC lines do not resemble those of other transmission lines such as microstrip. Therefore, the coupling mechanism of the FGC lines is different compared to thin film microstrip lines.
Coupling of Smoothed Particle Hydrodynamics with Finite Volume method for free-surface flows
NASA Astrophysics Data System (ADS)
Marrone, S.; Di Mascio, A.; Le Touzé, D.
2016-04-01
A new algorithm for the solution of free surface flows with large front deformation and fragmentation is presented. The algorithm is obtained by coupling a classical Finite Volume (FV) approach, that discretizes the Navier-Stokes equations on a block structured Eulerian grid, with an approach based on the Smoothed Particle Hydrodynamics (SPH) method, implemented in a Lagrangian framework. The coupling procedure is formulated in such a way that each solver is applied in the region where its intrinsic characteristics can be exploited in the most efficient and accurate way: the FV solver is used to resolve the bulk flow and the wall regions, whereas the SPH solver is implemented in the free surface region to capture details of the front evolution. The reported results clearly prove that the combined use of the two solvers is convenient from the point of view of both accuracy and computing time.
NASA Astrophysics Data System (ADS)
Otis, Collin; Ferrero, Pietro; Candler, Graham; Givi, Peyman
2013-11-01
The scalar filtered mass density function (SFMDF) methodology is implemented into the computer code US3D. This is an unstructured Eulerian finite volume hydrodynamic solver and has proven very effective for simulation of compressible turbulent flows. The resulting SFMDF-US3D code is employed for large eddy simulation (LES) on unstructured meshes. Simulations are conducted of subsonic and supersonic flows under non-reacting and reacting conditions. The consistency and the accuracy of the simulated results are assessed along with appraisal of the overall performance of the methodology. The SFMDF-US3D is now capable of simulating high speed flows in complex configurations.
NASA Technical Reports Server (NTRS)
Ponchak, George E.; Dalton, Edan; Tentzeris, Emmanouil M.; Papapolymerou, John; Williams, W. Dan (Technical Monitor)
2001-01-01
Three-dimensional circuits built upon multiple layers of polyimide are required for constructing Si/SiGe monolithic microwave/millimeter-wave integrated circuits on complementary metal oxide semiconductor (CMOS) (low resistivity) Si wafers. Thin film microstrip lines (TFMS) with finite width ground planes embedded in the polyimide are often used. However, the closely spaced TFMS lines are susceptible to high levels of coupling, which degrades circuit performance. In this paper, Finite Difference Time Domain (FDTD) analysis and experimental measurements are used to show that the ground planes must be connected by via holes to reduce coupling in both the forward and backward directions.