NASA Astrophysics Data System (ADS)
Shen, W.; Schulte-Pelkum, V.; Ritzwoller, M. H.
2011-12-01
The joint inversion of surface wave dispersion and receiver functions was proven feasible on a station by station basis more than a decade ago. Joint application to a large number of stations across a broad region such as western US is more challenging, however, because of the different resolutions of the two methods. Improvements in resolution in surface wave studies derived from ambient noise and array-based methods applied to earthquake data now allow surface wave dispersion and receiver functions to be inverted simultaneously across much of the Earthscope/USArray Transportable Array (TA), and we have developed a Monte-Carlo procedure for this purpose. As a proof of concept we applied this procedure to a region containing 186 TA stations in the intermountain west, including a variety of tectonic settings such as the Colorado Plateau, the Basin and Range, the Rocky Mountains, and the Great Plains. This work has now been expanded to encompass all TA stations in the western US. Our approach includes three main components. (1) We enlarge the Earthscope Automated Receiver Survey (EARS) receiver function database by adding more events within a quality control procedure. A back-azimuth-independent receiver function and its associated uncertainties are constructed using a harmonic stripping algorithm. (2) Rayleigh wave dispersion curves are generated from the eikonal tomography applied to ambient noise cross-correlation data and Helmoholtz tomography applied to teleseismic surface wave data to yield dispersion maps from 8 sec to 80 sec period. (3) We apply a Metropolis Monte Carlo algorithm to invert for the average velocity structure beneath each station. Simple kriging is applied to interpolate to the discrete results into a continuous 3-D model. This method has now been applied to over 1,000 TA stations in the western US. We show that the receiver functions and surface wave dispersion data can be reconciled beneath more than 80% of the stations using a smooth
NASA Astrophysics Data System (ADS)
Wilde-Piorko, Monika; Polkowski, Marcin; Grad, Marek
2015-04-01
Geological and seismic structure under area of Poland is well studied by over one hundred thousand boreholes, over thirty deep seismic refraction and wide angle reflection profiles and by vertical seismic profiling, magnetic, gravity, magnetotelluric and thermal methods. Compilation of these studies allowed to create a high-resolution 3D P-wave velocity model down to 60 km depth in the area of Poland (Polkowski et al. 2014). Model also provides details about the geometry of main layers of sediments (Tertiary and Quaternary, Cretaceous, Jurassic, Triassic, Permian, old Paleozoic), consolidated/crystalline crust (upper, middle and lower) and uppermost mantle. This model gives an unique opportunity for calculation synthetic receiver function and compering it with observed receiver function calculated for permanent and temporary seismic stations. Modified ray-tracing method (Langston, 1977) can be used directly to calculate the response of the structure with dipping interfaces to the incoming plane wave with fixed slowness and back-azimuth. So, 3D P-wave velocity model has been interpolated to 2.5D P-wave velocity model beneath each seismic station and back-azimuthal sections of components of receiver function have been calculated. Vp/Vs ratio is assumed to be 1.8, 1.67, 1.73, 1.77 and 1.8 in the sediments, upper/middle/lower consolidated/crystalline crust and uppermost mantle, respectively. Densities were calculated with combined formulas of Berteussen (1977) and Gardner et al. (1974). Additionally, to test a visibility of the lithosphere-asthenosphere boundary phases at receiver function sections models have been extended to 250 km depth based on P4-mantle model (Wilde-Piórko et al., 2010). National Science Centre Poland provided financial support for this work by NCN grant DEC-2011/02/A/ST10/00284 and by NCN grant UMO-2011/01/B/ST10/06653.
Frozen Gaussian approximation for 3-D seismic wave propagation
NASA Astrophysics Data System (ADS)
Chai, Lihui; Tong, Ping; Yang, Xu
2017-01-01
We present a systematic introduction on applying frozen Gaussian approximation (FGA) to compute synthetic seismograms in 3-D earth models. In this method, seismic wavefield is decomposed into frozen (fixed-width) Gaussian functions, which propagate along ray paths. Rather than the coherent state solution to the wave equation, this method is rigorously derived by asymptotic expansion on phase plane, with analysis of its accuracy determined by the ratio of short wavelength over large domain size. Similar to other ray-based beam methods (e.g. Gaussian beam methods), one can use relatively small number of Gaussians to get accurate approximations of high-frequency wavefield. The algorithm is embarrassingly parallel, which can drastically speed up the computation with a multicore-processor computer station. We illustrate the accuracy and efficiency of the method by comparing it to the spectral element method for a 3-D seismic wave propagation in homogeneous media, where one has the analytical solution as a benchmark. As another proof of methodology, simulations of high-frequency seismic wave propagation in heterogeneous media are performed for 3-D waveguide model and smoothed Marmousi model, respectively. The second contribution of this paper is that, we incorporate the Snell's law into the FGA formulation, and asymptotically derive reflection, transmission and free surface conditions for FGA to compute high-frequency seismic wave propagation in high contrast media. We numerically test these conditions by computing traveltime kernels of different phases in the 3-D crust-over-mantle model.
NASA Astrophysics Data System (ADS)
Harvey, R. W. (Bob); Petrov, Yu. V.; Jaeger, E. F.; Berry, L. A.; Bonoli, P. T.; Bader, A.
2015-11-01
A time-dependent simulation of C-Mod pulsed ICRF power is made calculating minority hydrogen ion distribution functions with the CQL3D-Hybrid-FOW finite-orbit-width Fokker-Planck code. ICRF fields are calculated with the AORSA full wave code, and RF diffusion coefficients are obtained from these fields using the DC Lorentz gyro-orbit code. Prior results with a zero-banana-width simulation using the CQL3D/AORSA/DC time-cycles showed a pronounced enhancement of the H distribution in the perpendicular velocity direction compared to results obtained from Stix's quasilinear theory, in general agreement with experiment. The present study compares the new FOW results, including relevant gyro-radius effects, to determine the importance of these effects on the the NPA synthetic diagnostic time-dependence. The new NPA results give increased agreement with experiment, particularly in the ramp-down time after the ICRF pulse. Funded, through subcontract with Massachusetts Institute of Technology, by USDOE sponsored SciDAC Center for Simulation of Wave-Plasma Interactions.
Instability and Wave Propagation in Structured 3D Composites
NASA Astrophysics Data System (ADS)
Kaynia, Narges; Fang, Nicholas X.; Boyce, Mary C.
2014-03-01
Many structured composites found in nature possess undulating and wrinkled interfacial layers that regulate mechanical, chemical, acoustic, adhesive, thermal, electrical and optical functions of the material. This research focused on the complex instability and wrinkling pattern arising in 3D structured composites and the effect of the buckling pattern on the overall structural response. The 3D structured composites consisted of stiffer plates supported by soft matrix on both sides. Compression beyond the critical strain led to complex buckling patterns in the initially straight plates. The motivation of our work is to elaborate the formation of a system of prescribed periodic scatterers (metamaterials) due to buckling, and their effect to interfere wave propagation through the metamaterial structures. Such metamaterials made from elastomers enable large reversible deformation and, as a result, significant changes of the wave propagation properties. We developed analytical and finite element models to capture various aspects of the instability mechanism. Mechanical experiments were designed to further explore the modeling results. The ability to actively alter the 3D composite structure can enable on-demand tunability of many different functions, such as active control of wave propagation to create band-gaps and waveguides.
Investigation of surface wave amplitudes in 3-D velocity and 3-D Q models
NASA Astrophysics Data System (ADS)
Ruan, Y.; Zhou, Y.
2010-12-01
It has been long recognized that seismic amplitudes depend on both wave speed structures and anelasticity (Q) structures. However, the effects of lateral heterogeneities in wave speed and Q structures on seismic amplitudes has not been well understood. We investigate the effects of 3-D wave speed and 3-D anelasticity (Q) structures on surface-wave amplitudes based upon wave propagation simulations of twelve globally-distributed earthquakes and 801 stations in Earth models with and without lateral heterogeneities in wave speed and anelasticity using a Spectral Element Method (SEM). Our tomographic-like 3-D Q models are converted from a velocity model S20RTS using a set of reasonable mineralogical parameters, assuming lateral perturbations in both velocity and Q are due to temperature perturbations. Surface-wave amplitude variations of SEM seismograms are measured in the period range of 50--200 s using boxcar taper, cosine taper and Slepian multi-tapers. We calculate ray-theoretical predictions of surface-wave amplitude perturbations due to elastic focusing, attenuation, and anelastic focusing which respectively depend upon the second spatial derivative (''roughness'') of perturbations in phase velocity, 1/Q, and the roughness of perturbations in 1/Q. Both numerical experiments and theoretical calculations show that (1) for short-period (~ 50 s) surface waves, the effects of amplitude attenuation due to 3-D Q structures are comparable with elastic focusing effects due to 3-D wave speed structures; and (2) for long-period (> 100 s) surface waves, the effects of attenuation become much weaker than elastic focusing; and (3) elastic focusing effects are correlated with anelastic focusing at all periods due to the correlation between velocity and Q models; and (4) amplitude perturbations are depend on measurement techniques and therefore cannot be directly compared with ray-theoretical predictions because ray theory does not account for the effects of measurement
3D Ultrasonic Wave Simulations for Structural Health Monitoring
NASA Technical Reports Server (NTRS)
Campbell, Leckey Cara A/; Miler, Corey A.; Hinders, Mark K.
2011-01-01
Structural health monitoring (SHM) for the detection of damage in aerospace materials is an important area of research at NASA. Ultrasonic guided Lamb waves are a promising SHM damage detection technique since the waves can propagate long distances. For complicated flaw geometries experimental signals can be difficult to interpret. High performance computing can now handle full 3-dimensional (3D) simulations of elastic wave propagation in materials. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate ultrasound scattering from flaws in materials. EFIT results have been compared to experimental data and the simulations provide unique insight into details of the wave behavior. This type of insight is useful for developing optimized experimental SHM techniques. 3D EFIT can also be expanded to model wave propagation and scattering in anisotropic composite materials.
Elastic wave modelling in 3D heterogeneous media: 3D grid method
NASA Astrophysics Data System (ADS)
Jianfeng, Zhang; Tielin, Liu
2002-09-01
We present a new numerical technique for elastic wave modelling in 3D heterogeneous media with surface topography, which is called the 3D grid method in this paper. This work is an extension of the 2D grid method that models P-SV wave propagation in 2D heterogeneous media. Similar to the finite-element method in the discretization of a numerical mesh, the proposed scheme is flexible in incorporating surface topography and curved interfaces; moreover it satisfies the free-surface boundary conditions of 3D topography naturally. The algorithm, developed from a parsimonious staggered-grid scheme, solves the problem using integral equilibrium around each node, instead of satisfying elastodynamic differential equations at each node as in the conventional finite-difference method. The computational cost and memory requirements for the proposed scheme are approximately the same as those used by the same order finite-difference method. In this paper, a mixed tetrahedral and parallelepiped grid method is presented; and the numerical dispersion and stability criteria on the tetrahedral grid method and parallelepiped grid method are discussed in detail. The proposed scheme is successfully tested against an analytical solution for the 3D Lamb problem and a solution of the boundary method for the diffraction of a hemispherical crater. Moreover, examples of surface-wave propagation in an elastic half-space with a semi-cylindrical trench on the surface and 3D plane-layered model are presented.
3D Guided Wave Motion Analysis on Laminated Composites
NASA Technical Reports Server (NTRS)
Tian, Zhenhua; Leckey, Cara; Yu, Lingyu
2013-01-01
Ultrasonic guided waves have proved useful for structural health monitoring (SHM) and nondestructive evaluation (NDE) due to their ability to propagate long distances with less energy loss compared to bulk waves and due to their sensitivity to small defects in the structure. Analysis of actively transmitted ultrasonic signals has long been used to detect and assess damage. However, there remain many challenging tasks for guided wave based SHM due to the complexity involved with propagating guided waves, especially in the case of composite materials. The multimodal nature of the ultrasonic guided waves complicates the related damage analysis. This paper presents results from parallel 3D elastodynamic finite integration technique (EFIT) simulations used to acquire 3D wave motion in the subject laminated carbon fiber reinforced polymer composites. The acquired 3D wave motion is then analyzed by frequency-wavenumber analysis to study the wave propagation and interaction in the composite laminate. The frequency-wavenumber analysis enables the study of individual modes and visualization of mode conversion. Delamination damage has been incorporated into the EFIT model to generate "damaged" data. The potential for damage detection in laminated composites is discussed in the end.
3D volumetric radar using 94-GHz millimeter waves
NASA Astrophysics Data System (ADS)
Takács, Barnabás
2006-05-01
This article describes a novel approach to the real-time visualization of 3D imagery obtained from a 3D millimeter wave scanning radar. The MMW radar system employs a spinning antenna to generate a fan-shaped scanning pattern of the entire scene. The beams formed this way provide all weather 3D distance measurements (range/azimuth display) of objects as they appear on the ground. The beam width of the antenna and its side lobes are optimized to produce the best possible resolution even at distances of up to 15 Kms. To create a full 3D data set the fan-pattern is tilted up and down with the help of a controlled stepper motor. For our experiments we collected data at 0.1 degrees increments while using both bi-static as well as a mono-static antennas in our arrangement. The data collected formed a stack of range-azimuth images in the shape of a cone. This information is displayed using our high-end 3D visualization engine capable of displaying high-resolution volumetric models with 30 frames per second. The resulting 3D scenes can then be viewed from any angle and subsequently processed to integrate, fuse or match them against real-life sensor imagery or 3D model data stored in a synthetic database.
Fast wave current drive antenna performance on D3-D
NASA Astrophysics Data System (ADS)
Mayberry, M. J.; Pinsker, R. I.; Petty, C. C.; Chiu, S. C.; Jackson, G. L.; Lippmann, S. I.; Prater, R.; Porkolab, M.
1991-10-01
Fast wave current drive (FWCD) experiments at 60 MHz are being performed on the D3-D tokamak for the first time in high electron temperature, high (beta) target plasmas. A four-element phased-array antenna is used to launch a directional wave spectrum with the peak n(sub parallel) value (approximately = 7) optimized for strong single-pass electron absorption due to electron Landau damping. For this experiment, high power FW injection (2 MW) must be accomplished without voltage breakdown in the transmission lines or antenna, and without significant impurity influx. In addition, there is the technological challenge of impedance matching a four-element antenna while maintaining equal currents and the correct phasing (90 degrees) in each of the straps for a directional spectrum. We describe the performance of the D3-D FWCD antenna during initial FW electron heating and current drive experiments in terms of these requirements.
Protrusive waves guide 3D cell migration along nanofibers
Guetta-Terrier, Charlotte; Monzo, Pascale; Zhu, Jie; Long, Hongyan; Venkatraman, Lakshmi; Zhou, Yue; Wang, PeiPei; Chew, Sing Yian; Mogilner, Alexander
2015-01-01
In vivo, cells migrate on complex three-dimensional (3D) fibrous matrices, which has made investigation of the key molecular and physical mechanisms that drive cell migration difficult. Using reductionist approaches based on 3D electrospun fibers, we report for various cell types that single-cell migration along fibronectin-coated nanofibers is associated with lateral actin-based waves. These cyclical waves have a fin-like shape and propagate up to several hundred micrometers from the cell body, extending the leading edge and promoting highly persistent directional movement. Cells generate these waves through balanced activation of the Rac1/N-WASP/Arp2/3 and Rho/formins pathways. The waves originate from one major adhesion site at leading end of the cell body, which is linked through actomyosin contractility to another site at the back of the cell, allowing force generation, matrix deformation and cell translocation. By combining experimental and modeling data, we demonstrate that cell migration in a fibrous environment requires the formation and propagation of dynamic, actin based fin-like protrusions. PMID:26553933
Protrusive waves guide 3D cell migration along nanofibers.
Guetta-Terrier, Charlotte; Monzo, Pascale; Zhu, Jie; Long, Hongyan; Venkatraman, Lakshmi; Zhou, Yue; Wang, PeiPei; Chew, Sing Yian; Mogilner, Alexander; Ladoux, Benoit; Gauthier, Nils C
2015-11-09
In vivo, cells migrate on complex three-dimensional (3D) fibrous matrices, which has made investigation of the key molecular and physical mechanisms that drive cell migration difficult. Using reductionist approaches based on 3D electrospun fibers, we report for various cell types that single-cell migration along fibronectin-coated nanofibers is associated with lateral actin-based waves. These cyclical waves have a fin-like shape and propagate up to several hundred micrometers from the cell body, extending the leading edge and promoting highly persistent directional movement. Cells generate these waves through balanced activation of the Rac1/N-WASP/Arp2/3 and Rho/formins pathways. The waves originate from one major adhesion site at leading end of the cell body, which is linked through actomyosin contractility to another site at the back of the cell, allowing force generation, matrix deformation and cell translocation. By combining experimental and modeling data, we demonstrate that cell migration in a fibrous environment requires the formation and propagation of dynamic, actin based fin-like protrusions.
Generation Of 3d Periodic Internal Wave Beams:
NASA Astrophysics Data System (ADS)
Chashechkin, Yuli D.; Vasiliev, Alexey Yu.
We study generation of 2D and 3D periodic internal wave beams in continuously strat- ified viscous liquid basing on a complete set of governing equations and exact bound- ary conditions that is no-slip for velocity and attenuation of all disturbances at infinite distance from the source. The linearized governing equations are solved by an integral transform method. A total set of dispersion equation roots contains terms correspond- ing to internal waves and additional roots describing two kinds of periodic boundary layers. The first one is a viscous boundary layer and has an analogue that is a periodic or Stokes' layer in a homogeneous fluid. Its thickness is defined by a kinematic viscos- ity coefficient and a buoyancy frequency. The second one, that is an internal boundary layer, is a specific feature of stratified flows. Its thickness besides the Stokes' scale contains additional factor depending on relative wave frequency and geometry of the problem that is on the local slope of emitting surface and a direction of the waves propagation. We have constructed exact solutions of linear problems describing gen- eration of 2D waves by a strip and 3D by a rectangular with an arbitrary ratio of sides moving along or normally to a sloping plane. We also calculated the wave pattern gen- erated by a part of a vertical cylinder surface with different ratios of intrinsic scales that is of cylinder radius, thickness of the boundary layer and internal viscous scale. All solutions are regularly matched between themselves in limiting cases. The spatial decay of the waves depends on dimension and geometry of the problem. Non-linear generation of internal waves by the Stokes' boundary layer on a periodically rotating horizontal disk or by interacting boundary layers on an arbitrary moving strip is in- vestigated. We found conditions of generation of the main frequency and its second harmonic. In experiments periodic waves beams from different sources are visualised by the
Simulation of 3D Global Wave Propagation Through Geodynamic Models
NASA Astrophysics Data System (ADS)
Schuberth, B.; Piazzoni, A.; Bunge, H.; Igel, H.; Steinle-Neumann, G.
2005-12-01
This project aims at a better understanding of the forward problem of global 3D wave propagation. We use the spectral element program "SPECFEM3D" (Komatitsch and Tromp, 2002a,b) with varying input models of seismic velocities derived from mantle convection simulations (Bunge et al., 2002). The purpose of this approach is to obtain seismic velocity models independently from seismological studies. In this way one can test the effects of varying parameters of the mantle convection models on the seismic wave field. In order to obtain the seismic velocities from the temperature field of the geodynamical simulations we follow a mineral physics approach. Assuming a certain mantle composition (e.g. pyrolite with CMASF composition) we compute the stable phases for each depth (i.e. pressure) and temperature by system Gibbs free energy minimization. Elastic moduli and density are calculated from the equations of state of the stable mineral phases. For this we use a mineral physics database derived from calorimetric experiments (enthalphy and entropy of formation, heat capacity) and EOS parameters.
High Resolution WENO Simulation of 3D Detonation Waves
2012-02-27
pocket behind the detonation front was not observed in their results because the rotating transverse detonation completely consumed the unburned gas. Dou...three-dimensional detonations We add source terms (functions of x, y, z and t) to the PDE system so that the following functions are exact solutions to... detonation rotates counter-clockwise, opposite to that in [48]. It can be seen that, the triple lines and transverse waves collide with the walls, and strong
3D printing of functional biomaterials for tissue engineering.
Zhu, Wei; Ma, Xuanyi; Gou, Maling; Mei, Deqing; Zhang, Kang; Chen, Shaochen
2016-08-01
3D printing is emerging as a powerful tool for tissue engineering by enabling 3D cell culture within complex 3D biomimetic architectures. This review discusses the prevailing 3D printing techniques and their most recent applications in building tissue constructs. The work associated with relatively well-known inkjet and extrusion-based bioprinting is presented with the latest advances in the fields. Emphasis is put on introducing two relatively new light-assisted bioprinting techniques, including digital light processing (DLP)-based bioprinting and laser based two photon polymerization (TPP) bioprinting. 3D bioprinting of vasculature network is particularly discussed for its foremost significance in maintaining tissue viability and promoting functional maturation. Limitations to current bioprinting approaches, as well as future directions of bioprinting functional tissues are also discussed.
3D-Fun: predicting enzyme function from structure.
von Grotthuss, Marcin; Plewczynski, Dariusz; Vriend, Gert; Rychlewski, Leszek
2008-07-01
The 'omics' revolution is causing a flurry of data that all needs to be annotated for it to become useful. Sequences of proteins of unknown function can be annotated with a putative function by comparing them with proteins of known function. This form of annotation is typically performed with BLAST or similar software. Structural genomics is nowadays also bringing us three dimensional structures of proteins with unknown function. We present here software that can be used when sequence comparisons fail to determine the function of a protein with known structure but unknown function. The software, called 3D-Fun, is implemented as a server that runs at several European institutes and is freely available for everybody at all these sites. The 3D-Fun servers accept protein coordinates in the standard PDB format and compare them with all known protein structures by 3D structural superposition using the 3D-Hit software. If structural hits are found with proteins with known function, these are listed together with their function and some vital comparison statistics. This is conceptually very similar in 3D to what BLAST does in 1D. Additionally, the superposition results are displayed using interactive graphics facilities. Currently, the 3D-Fun system only predicts enzyme function but an expanded version with Gene Ontology predictions will be available soon. The server can be accessed at http://3dfun.bioinfo.pl/ or at http://3dfun.cmbi.ru.nl/.
Towards functional assembly of 3D and 2D nanomaterials
NASA Astrophysics Data System (ADS)
Jacobs, Christopher B.; Wang, Kai; Ievlev, Anton V.; Muckley, Eric S.; Ivanov, Ilia N.
2016-09-01
Functional assemblies of materials can be realized by tuning the work function and band gap of nanomaterials by rational material selection and design. Here we demonstrate the structural assembly of 2D and 3D nanomaterials and show that layering a 2D material monolayer on a 3D metal oxide leads to substantial alteration of both the surface potential and optical properties of the 3D material. A 40 nm thick film of polycrystalline NiO was produced by room temperature rf-sputtering, resulting in a 3D nanoparticle assembly. Chemical vapor deposition (CVD) grown 10-30 μm WS2 flakes (2D material) were placed on the NiO surface using a PDMS stamp transfer technique. The 2D/3D WS2/NiO assembly was characterized using confocal micro Raman spectroscopy to evaluate the vibrational properties and using Kelvin probe force microscopy (KPFM) to evaluate the surface potential. Raman maps of the 2D/3D assembly show spatial non-uniformity of the A1g mode ( 418 cm-1) and the disorder-enhanced longitudinal acoustic mode, 2LA(M) ( 350 cm-1), suggesting that the WS2 exists in a strained condition on when transferred onto 3D polycrystalline NiO. KPFM measurements show that single layer WS2 on SiO2 has a surface potential 75 mV lower than that of SiO2, whereas the surface potential of WS2 on NiO is 15 mV higher than NiO, indicating that WS2 could act as electron donor or acceptor depending on the 3D material it is interfaced with. Thus 2D and 3D materials can be organized into functional assemblies with electron flow controlled by the WS2 either as the electron donor or acceptor.
Embedding objects during 3D printing to add new functionalities.
Yuen, Po Ki
2016-07-01
A novel method for integrating and embedding objects to add new functionalities during 3D printing based on fused deposition modeling (FDM) (also known as fused filament fabrication or molten polymer deposition) is presented. Unlike typical 3D printing, FDM-based 3D printing could allow objects to be integrated and embedded during 3D printing and the FDM-based 3D printed devices do not typically require any post-processing and finishing. Thus, various fluidic devices with integrated glass cover slips or polystyrene films with and without an embedded porous membrane, and optical devices with embedded Corning(®) Fibrance™ Light-Diffusing Fiber were 3D printed to demonstrate the versatility of the FDM-based 3D printing and embedding method. Fluid perfusion flow experiments with a blue colored food dye solution were used to visually confirm fluid flow and/or fluid perfusion through the embedded porous membrane in the 3D printed fluidic devices. Similar to typical 3D printed devices, FDM-based 3D printed devices are translucent at best unless post-polishing is performed and optical transparency is highly desirable in any fluidic devices; integrated glass cover slips or polystyrene films would provide a perfect optical transparent window for observation and visualization. In addition, they also provide a compatible flat smooth surface for biological or biomolecular applications. The 3D printed fluidic devices with an embedded porous membrane are applicable to biological or chemical applications such as continuous perfusion cell culture or biocatalytic synthesis but without the need for any post-device assembly and finishing. The 3D printed devices with embedded Corning(®) Fibrance™ Light-Diffusing Fiber would have applications in display, illumination, or optical applications. Furthermore, the FDM-based 3D printing and embedding method could also be utilized to print casting molds with an integrated glass bottom for polydimethylsiloxane (PDMS) device replication
Protein function prediction using local 3D templates.
Laskowski, Roman A; Watson, James D; Thornton, Janet M
2005-08-19
The prediction of a protein's function from its 3D structure is becoming more and more important as the worldwide structural genomics initiatives gather pace and continue to solve 3D structures, many of which are of proteins of unknown function. Here, we present a methodology for predicting function from structure that shows great promise. It is based on 3D templates that are defined as specific 3D conformations of small numbers of residues. We use four types of template, covering enzyme active sites, ligand-binding residues, DNA-binding residues and reverse templates. The latter are templates generated from the target structure itself and scanned against a representative subset of all known protein structures. Together, the templates provide a fairly thorough coverage of the known structures and ensure that if there is a match to a known structure it is unlikely to be missed. A new scoring scheme provides a highly sensitive means of discriminating between true positive and false positive template matches. In all, the methodology provides a powerful new tool for function prediction to complement those already in use.
3D Full-Wave Simulations of Reflectometry
Valeo, E. J.; Kramer, G. J.; Nazikian, R.
2009-11-26
The characterization of fluctuation amplitudes, spatial correlation lengths, and wave vectors through measurement of the correlation properties of reflected microwave diagnostic signals depends on a quantitative knowledge of propagation in toroidal, magnetized plasma. The disparity between the radiation wavelength (mm) and the plasma size makes full wave computations challenging. We extend a two dimensional model which computes propagation in a poloidal plane to include toroidal variation. The model reduces the computational burden compared to that of solving the full-wave equation everywhere-but retains both diffraction and refraction-by merging a description appropriate to the under dense plasma (paraxial) with the required full-wave description near the reflection layer. Initial results for ITER-like profiles demonstrate the utility of the tool as an aid in specifying antenna positioning and setting sensitivity requirements.
Analysis of wave propagation in periodic 3D waveguides
NASA Astrophysics Data System (ADS)
Schaal, Christoph; Bischoff, Stefan; Gaul, Lothar
2013-11-01
Structural Health Monitoring (SHM) is a growing research field in the realm of civil engineering. SHM concepts are implemented using integrated sensors and actuators to evaluate the state of a structure. Within this work, wave-based techniques are addressed. Dispersion effects for propagating waves in waveguides of different materials are analyzed for various different cross-sections. Since analytical theory is limited, a general approach based on the Waveguide Finite Element Method is applied. Numerical results are verified experimentally.
Robust bioengineered 3D functional human intestinal epithelium
Chen, Ying; Lin, Yinan; Davis, Kimberly M.; Wang, Qianrui; Rnjak-Kovacina, Jelena; Li, Chunmei; Isberg, Ralph R.; Kumamoto, Carol A.; Mecsas, Joan; Kaplan, David L.
2015-01-01
Intestinal functions are central to human physiology, health and disease. Options to study these functions with direct relevance to the human condition remain severely limited when using conventional cell cultures, microfluidic systems, organoids, animal surrogates or human studies. To replicate in vitro the tissue architecture and microenvironments of native intestine, we developed a 3D porous protein scaffolding system, containing a geometrically-engineered hollow lumen, with adaptability to both large and small intestines. These intestinal tissues demonstrated representative human responses by permitting continuous accumulation of mucous secretions on the epithelial surface, establishing low oxygen tension in the lumen, and interacting with gut-colonizing bacteria. The newly developed 3D intestine model enabled months-long sustained access to these intestinal functions in vitro, readily integrable with a multitude of different organ mimics and will therefore ensure a reliable ex vivo tissue system for studies in a broad context of human intestinal diseases and treatments. PMID:26374193
Robust bioengineered 3D functional human intestinal epithelium.
Chen, Ying; Lin, Yinan; Davis, Kimberly M; Wang, Qianrui; Rnjak-Kovacina, Jelena; Li, Chunmei; Isberg, Ralph R; Kumamoto, Carol A; Mecsas, Joan; Kaplan, David L
2015-09-16
Intestinal functions are central to human physiology, health and disease. Options to study these functions with direct relevance to the human condition remain severely limited when using conventional cell cultures, microfluidic systems, organoids, animal surrogates or human studies. To replicate in vitro the tissue architecture and microenvironments of native intestine, we developed a 3D porous protein scaffolding system, containing a geometrically-engineered hollow lumen, with adaptability to both large and small intestines. These intestinal tissues demonstrated representative human responses by permitting continuous accumulation of mucous secretions on the epithelial surface, establishing low oxygen tension in the lumen, and interacting with gut-colonizing bacteria. The newly developed 3D intestine model enabled months-long sustained access to these intestinal functions in vitro, readily integrable with a multitude of different organ mimics and will therefore ensure a reliable ex vivo tissue system for studies in a broad context of human intestinal diseases and treatments.
Mach-wave coherence in 3D media with random heterogeneities
NASA Astrophysics Data System (ADS)
Vyas, Jagdish C.; Mai, P. Martin; Galis, Martin; Dunham, Eric M.; Imperatori, Walter
2016-04-01
We investigate Mach-waves coherence for complex super-shear ruptures embedded in 3D random media that lead to seismic scattering. We simulate Mach-wave using kinematic earthquake sources that include fault-regions over which the rupture propagates at super-shear speed. The local slip rate is modeled with the regularized Yoffe function. The medium heterogeneities are characterized by Von Karman correlation function. We consider various realizations of 3D random media from combinations of different values of correlation length (0.5 km, 2 km, 5 km), standard deviation (5%, 10%, 15%) and Hurst exponent (0.2). Simulations in a homogeneous medium serve as a reference case. The ground-motion simulations (maximum resolved frequency of 5 Hz) are conducted by solving the elasto-dynamic equations of motions using a generalized finite-difference method, assuming a vertical strike-slip fault. The seismic wavefield is sampled at numerous locations within the Mach-cone region to study the properties and evolution of the Mach-waves in scattering media. We find that the medium scattering from random heterogeneities significantly diminishes the coherence of Mach-wave in terms of both amplitude and frequencies. We observe that Mach-waves are considerably scattered at distances RJB > 20 km (and beyond) for random media with standard deviation 10%. The scattering efficiency of the medium for small Hurst exponents (H <= 0.2) is mainly controlled by the standard deviation of the velocity heterogeneities, rather than their correlation length, as both theoretical considerations and numerical experiments indicate. Based on our simulations, we propose that local super-shear ruptures may be more common in nature then reported, but are very difficult to detect due to the strong seismic scattering. We suggest that if an earthquake is recorded within 10-15 km fault perpendicular distance and has high PGA, then inversion should be carried out by allowing rupture speed variations from sub
Multiprocess 3D printing for increasing component functionality.
MacDonald, Eric; Wicker, Ryan
2016-09-30
Layer-by-layer deposition of materials to manufacture parts-better known as three-dimensional (3D) printing or additive manufacturing-has been flourishing as a fabrication process in the past several years and now can create complex geometries for use as models, assembly fixtures, and production molds. Increasing interest has focused on the use of this technology for direct manufacturing of production parts; however, it remains generally limited to single-material fabrication, which can limit the end-use functionality of the fabricated structures. The next generation of 3D printing will entail not only the integration of dissimilar materials but the embedding of active components in order to deliver functionality that was not possible previously. Examples could include arbitrarily shaped electronics with integrated microfluidic thermal management and intelligent prostheses custom-fit to the anatomy of a specific patient. We review the state of the art in multiprocess (or hybrid) 3D printing, in which complementary processes, both novel and traditional, are combined to advance the future of manufacturing.
3D dynamic simulation of crack propagation in extracorporeal shock wave lithotripsy
NASA Astrophysics Data System (ADS)
Wijerathne, M. L. L.; Hori, Muneo; Sakaguchi, Hide; Oguni, Kenji
2010-06-01
Some experimental observations of Shock Wave Lithotripsy(SWL), which include 3D dynamic crack propagation, are simulated with the aim of reproducing fragmentation of kidney stones with SWL. Extracorporeal shock wave lithotripsy (ESWL) is the fragmentation of kidney stones by focusing an ultrasonic pressure pulse onto the stones. 3D models with fine discretization are used to accurately capture the high amplitude shear shock waves. For solving the resulting large scale dynamic crack propagation problem, PDS-FEM is used; it provides numerically efficient failure treatments. With a distributed memory parallel code of PDS-FEM, experimentally observed 3D photoelastic images of transient stress waves and crack patterns in cylindrical samples are successfully reproduced. The numerical crack patterns are in good agreement with the experimental ones, quantitatively. The results shows that the high amplitude shear waves induced in solid, by the lithotriptor generated shock wave, play a dominant role in stone fragmentation.
NASA Astrophysics Data System (ADS)
Harvey, R. W.; Petrov, Yu.; Jaeger, E. F.; Berry, L. A.; Bonoli, P. T.; Bader, A.
2015-12-01
A time-dependent simulation of C-Mod pulsed TCRF power is made obtaining minority hydrogen ion distributions with the CQL3D-Hybrid-FOW finite-orbit-width Fokker-Planck code. Cyclotron-resonant TCRF fields are calculated with the AORSA full wave code. The RF diffusion coefficients used in CQL3D are obtained with the DC Lorentz gyro-orbit code for perturbed particle trajectories in the combined equilibrium and TCRF electromagnetic fields. Prior results with a zero-banana-width simulation using the CQL3D/AORSA/DC time-cycles showed a pronounced enhancement of the H distribution in the perpendicular velocity direction compared to results obtained from Stix's quasilinear theory, and this substantially increased the rampup rate of the observed vertically-viewed neutral particle analyzer (NPA) flux, in general agreement with experiment. However, ramp down of the NPA flux after the pulse, remained long compared to the experiment. The present study compares the new FOW results, including relevant gyro-radius effects, to determine the importance of these new effects on the the NPA time-dependence.
Gene3D: modelling protein structure, function and evolution.
Yeats, Corin; Maibaum, Michael; Marsden, Russell; Dibley, Mark; Lee, David; Addou, Sarah; Orengo, Christine A
2006-01-01
The Gene3D release 4 database and web portal (http://cathwww.biochem.ucl.ac.uk:8080/Gene3D) provide a combined structural, functional and evolutionary view of the protein world. It is focussed on providing structural annotation for protein sequences without structural representatives--including the complete proteome sets of over 240 different species. The protein sequences have also been clustered into whole-chain families so as to aid functional prediction. The structural annotation is generated using HMM models based on the CATH domain families; CATH is a repository for manually deduced protein domains. Amongst the changes from the last publication are: the addition of over 100 genomes and the UniProt sequence database, domain data from Pfam, metabolic pathway and functional data from COGs, KEGG and GO, and protein-protein interaction data from MINT and BIND. The website has been rebuilt to allow more sophisticated querying and the data returned is presented in a clearer format with greater functionality. Furthermore, all data can be downloaded in a simple XML format, allowing users to carry out complex investigations at their own computers.
Tailored complex 3D vortex lattice structures by perturbed multiples of three-plane waves.
Xavier, Jolly; Vyas, Sunil; Senthilkumaran, Paramasivam; Joseph, Joby
2012-04-20
As three-plane waves are the minimum number required for the formation of vortex-embedded lattice structures by plane wave interference, we present our experimental investigation on the formation of complex 3D photonic vortex lattice structures by a designed superposition of multiples of phase-engineered three-plane waves. The unfolding of the generated complex photonic lattice structures with higher order helical phase is realized by perturbing the superposition of a relatively phase-encoded, axially equidistant multiple of three noncoplanar plane waves. Through a programmable spatial light modulator assisted single step fabrication approach, the unfolded 3D vortex lattice structures are experimentally realized, well matched to our computer simulations. The formation of higher order intertwined helices embedded in these 3D spiraling vortex lattice structures by the superposition of the multiples of phase-engineered three-plane waves interference is also studied.
Probability Density Function at the 3D Anderson Transition
NASA Astrophysics Data System (ADS)
Rodriguez, Alberto; Vasquez, Louella J.; Roemer, Rudolf
2009-03-01
The probability density function (PDF) for the wavefunction amplitudes is studied at the metal-insulator transition of the 3D Anderson model, for very large systems up to L^3=240^3. The implications of the multifractal nature of the state upon the PDF are presented in detail. A formal expression between the PDF and the singularity spectrum f(α) is given. The PDF can be easily used to carry out a numerical multifractal analysis and it appears as a valid alternative to the more usual approach based on the scaling law of the general inverse participation rations.
A 3D staggered-grid finite difference scheme for poroelastic wave equation
NASA Astrophysics Data System (ADS)
Zhang, Yijie; Gao, Jinghuai
2014-10-01
Three dimensional numerical modeling has been a viable tool for understanding wave propagation in real media. The poroelastic media can better describe the phenomena of hydrocarbon reservoirs than acoustic and elastic media. However, the numerical modeling in 3D poroelastic media demands significantly more computational capacity, including both computational time and memory. In this paper, we present a 3D poroelastic staggered-grid finite difference (SFD) scheme. During the procedure, parallel computing is implemented to reduce the computational time. Parallelization is based on domain decomposition, and communication between processors is performed using message passing interface (MPI). Parallel analysis shows that the parallelized SFD scheme significantly improves the simulation efficiency and 3D decomposition in domain is the most efficient. We also analyze the numerical dispersion and stability condition of the 3D poroelastic SFD method. Numerical results show that the 3D numerical simulation can provide a real description of wave propagation.
3D mapping of elastic modulus using shear wave optical micro-elastography
NASA Astrophysics Data System (ADS)
Zhu, Jiang; Qi, Li; Miao, Yusi; Ma, Teng; Dai, Cuixia; Qu, Yueqiao; He, Youmin; Gao, Yiwei; Zhou, Qifa; Chen, Zhongping
2016-10-01
Elastography provides a powerful tool for histopathological identification and clinical diagnosis based on information from tissue stiffness. Benefiting from high resolution, three-dimensional (3D), and noninvasive optical coherence tomography (OCT), optical micro-elastography has the ability to determine elastic properties with a resolution of ~10 μm in a 3D specimen. The shear wave velocity measurement can be used to quantify the elastic modulus. However, in current methods, shear waves are measured near the surface with an interference of surface waves. In this study, we developed acoustic radiation force (ARF) orthogonal excitation optical coherence elastography (ARFOE-OCE) to visualize shear waves in 3D. This method uses acoustic force perpendicular to the OCT beam to excite shear waves in internal specimens and uses Doppler variance method to visualize shear wave propagation in 3D. The measured propagation of shear waves agrees well with the simulation results obtained from finite element analysis (FEA). Orthogonal acoustic excitation allows this method to measure the shear modulus in a deeper specimen which extends the elasticity measurement range beyond the OCT imaging depth. The results show that the ARFOE-OCE system has the ability to noninvasively determine the 3D elastic map.
3D mapping of elastic modulus using shear wave optical micro-elastography
Zhu, Jiang; Qi, Li; Miao, Yusi; Ma, Teng; Dai, Cuixia; Qu, Yueqiao; He, Youmin; Gao, Yiwei; Zhou, Qifa; Chen, Zhongping
2016-01-01
Elastography provides a powerful tool for histopathological identification and clinical diagnosis based on information from tissue stiffness. Benefiting from high resolution, three-dimensional (3D), and noninvasive optical coherence tomography (OCT), optical micro-elastography has the ability to determine elastic properties with a resolution of ~10 μm in a 3D specimen. The shear wave velocity measurement can be used to quantify the elastic modulus. However, in current methods, shear waves are measured near the surface with an interference of surface waves. In this study, we developed acoustic radiation force (ARF) orthogonal excitation optical coherence elastography (ARFOE-OCE) to visualize shear waves in 3D. This method uses acoustic force perpendicular to the OCT beam to excite shear waves in internal specimens and uses Doppler variance method to visualize shear wave propagation in 3D. The measured propagation of shear waves agrees well with the simulation results obtained from finite element analysis (FEA). Orthogonal acoustic excitation allows this method to measure the shear modulus in a deeper specimen which extends the elasticity measurement range beyond the OCT imaging depth. The results show that the ARFOE-OCE system has the ability to noninvasively determine the 3D elastic map. PMID:27762276
3-D FDTD simulation of shear waves for evaluation of complex modulus imaging.
Orescanin, Marko; Wang, Yue; Insana, Michael
2011-02-01
The Navier equation describing shear wave propagation in 3-D viscoelastic media is solved numerically with a finite differences time domain (FDTD) method. Solutions are formed in terms of transverse scatterer velocity waves and then verified via comparison to measured wave fields in heterogeneous hydrogel phantoms. The numerical algorithm is used as a tool to study the effects on complex shear modulus estimation from wave propagation in heterogeneous viscoelastic media. We used an algebraic Helmholtz inversion (AHI) technique to solve for the complex shear modulus from simulated and experimental velocity data acquired in 2-D and 3-D. Although 3-D velocity estimates are required in general, there are object geometries for which 2-D inversions provide accurate estimations of the material properties. Through simulations and experiments, we explored artifacts generated in elastic and dynamic-viscous shear modulus images related to the shear wavelength and average viscosity.
3D Modeling of Ultrasonic Wave Interaction with Disbonds and Weak Bonds
NASA Technical Reports Server (NTRS)
Leckey, C.; Hinders, M.
2011-01-01
Ultrasonic techniques, such as the use of guided waves, can be ideal for finding damage in the plate and pipe-like structures used in aerospace applications. However, the interaction of waves with real flaw types and geometries can lead to experimental signals that are difficult to interpret. 3-dimensional (3D) elastic wave simulations can be a powerful tool in understanding the complicated wave scattering involved in flaw detection and for optimizing experimental techniques. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate Lamb wave scattering from realistic flaws. This paper discusses simulation results for an aluminum-aluminum diffusion disbond and an aluminum-epoxy disbond and compares results from the disbond case to the common artificial flaw type of a flat-bottom hole. The paper also discusses the potential for extending the 3D EFIT equations to incorporate physics-based weak bond models for simulating wave scattering from weak adhesive bonds.
NASA Astrophysics Data System (ADS)
Leng, Kuangdai; Nissen-Meyer, Tarje; van Driel, Martin
2016-12-01
We present a new, computationally efficient numerical method to simulate global seismic wave propagation in realistic 3-D Earth models. We characterize the azimuthal dependence of 3-D wavefields in terms of Fourier series, such that the 3-D equations of motion reduce to an algebraic system of coupled 2-D meridian equations, which is then solved by a 2-D spectral element method (SEM). Computational efficiency of such a hybrid method stems from lateral smoothness of 3-D Earth models and axial singularity of seismic point sources, which jointly confine the Fourier modes of wavefields to a few lower orders. We show novel benchmarks for global wave solutions in 3-D structures between our method and an independent, fully discretized 3-D SEM with remarkable agreement. Performance comparisons are carried out on three state-of-the-art tomography models, with seismic period ranging from 34 s down to 11 s. It turns out that our method has run up to two orders of magnitude faster than the 3-D SEM, featured by a computational advantage expanding with seismic frequency.
FGG-NUFFT-Based Method for Near-Field 3-D Imaging Using Millimeter Waves
Kan, Yingzhi; Zhu, Yongfeng; Tang, Liang; Fu, Qiang; Pei, Hucheng
2016-01-01
In this paper, to deal with the concealed target detection problem, an accurate and efficient algorithm for near-field millimeter wave three-dimensional (3-D) imaging is proposed that uses a two-dimensional (2-D) plane antenna array. First, a two-dimensional fast Fourier transform (FFT) is performed on the scattered data along the antenna array plane. Then, a phase shift is performed to compensate for the spherical wave effect. Finally, fast Gaussian gridding based nonuniform FFT (FGG-NUFFT) combined with 2-D inverse FFT (IFFT) is performed on the nonuniform 3-D spatial spectrum in the frequency wavenumber domain to achieve 3-D imaging. The conventional method for near-field 3-D imaging uses Stolt interpolation to obtain uniform spatial spectrum samples and performs 3-D IFFT to reconstruct a 3-D image. Compared with the conventional method, our FGG-NUFFT based method is comparable in both efficiency and accuracy in the full sampled case and can obtain more accurate images with less clutter and fewer noisy artifacts in the down-sampled case, which are good properties for practical applications. Both simulation and experimental results demonstrate that the FGG-NUFFT-based near-field 3-D imaging algorithm can have better imaging performance than the conventional method for down-sampled measurements. PMID:27657066
FGG-NUFFT-Based Method for Near-Field 3-D Imaging Using Millimeter Waves.
Kan, Yingzhi; Zhu, Yongfeng; Tang, Liang; Fu, Qiang; Pei, Hucheng
2016-09-19
In this paper, to deal with the concealed target detection problem, an accurate and efficient algorithm for near-field millimeter wave three-dimensional (3-D) imaging is proposed that uses a two-dimensional (2-D) plane antenna array. First, a two-dimensional fast Fourier transform (FFT) is performed on the scattered data along the antenna array plane. Then, a phase shift is performed to compensate for the spherical wave effect. Finally, fast Gaussian gridding based nonuniform FFT (FGG-NUFFT) combined with 2-D inverse FFT (IFFT) is performed on the nonuniform 3-D spatial spectrum in the frequency wavenumber domain to achieve 3-D imaging. The conventional method for near-field 3-D imaging uses Stolt interpolation to obtain uniform spatial spectrum samples and performs 3-D IFFT to reconstruct a 3-D image. Compared with the conventional method, our FGG-NUFFT based method is comparable in both efficiency and accuracy in the full sampled case and can obtain more accurate images with less clutter and fewer noisy artifacts in the down-sampled case, which are good properties for practical applications. Both simulation and experimental results demonstrate that the FGG-NUFFT-based near-field 3-D imaging algorithm can have better imaging performance than the conventional method for down-sampled measurements.
Numerical Investigation of 3D multichannel analysis of surface wave method
NASA Astrophysics Data System (ADS)
Wang, Limin; Xu, Yixian; Luo, Yinhe
2015-08-01
Multichannel analysis of surface wave (MASW) method is an efficient tool to obtain near-surface S-wave velocity, and it has gained popularity in engineering practice. Up to now, most examples of using the MASW technique are focused on 2D models or data from a 1D linear receiver spread. We propose a 3D MASW scheme. A finite-difference (FD) method is used to investigate the method using linear and fan-shaped receiver spreads. Results show that the 3D topography strongly affects propagation of Rayleigh waves. The energy concentration of dispersion image is distorted and bifurcated because of the influence of free-surface topography. These effects are reduced with the 3D MASW method. Lastly we investigate the relation between the array size and the resolution of dispersion measurement.
Extensions of 1d Bgk Electron Solitary Wave Solutions To 3d Magnetized and Unmagnetized Plasmas
NASA Astrophysics Data System (ADS)
Chen, Li-Jen; Parks, George K.
This paper will compare the key results for BGK electron solitary waves in 3D mag- netized and unmagnetized plasmas. For 3D magnetized plasmas with highly magnetic field-aligned electrons, our results predict that the parallel widths of the solitary waves can be smaller than one Debye length, the solitary waves can be large scale features of the magnetosphere, and the parallel width-amplitude relation has a dependence on the perpendicular size. We can thus obtain an estimate on the typical perpendicular size of the observed solitary waves assuming a series of consecutive solitary waves are in the same flux tude with a particular perpendicular span. In 3D unmagnetized plasma systems such as the neutral sheet and magnetic reconnection sites, our theory indi- cates that although mathematical solutions can be constructed as the time-stationary solutions for the nonlinear Vlasov-Poisson equations, there does not exist a param- eter range for the solutions to be physical. We conclude that single-humped solitary potential pulses cannot be self-consistently supported by charged particles in 3D un- magnetized plasmas.
NASA Astrophysics Data System (ADS)
Yang, Si-Tong; Wei, Jiu-Chuan; Cheng, Jiu-Long; Shi, Long-Qing; Wen, Zhi-Jie
2016-12-01
Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling twodimensional wave fields and therefore cannot accurately simulate three-dimensional (3-D) full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity V x, V y, and V z for the same node in 3-D staggered-grid finite difference models by calculating the average value of V y, and V z of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways
Efficient methods to model the scattering of ultrasonic guided waves in 3D
NASA Astrophysics Data System (ADS)
Moreau, L.; Velichko, A.; Wilcox, P. D.
2010-03-01
The propagation of ultrasonic guided waves and their interaction with a defect is of interest to the nondestructive testing community. There is no general solution to the scattering problem and it is still an ongoing research topic. Due to the complexity of guided wave scattering problems, most existing models are related to the 2D case. However, thanks to the increase in computer calculation power, specific 3D problems can also be studied, with the help of numerical or semi-analytical methods. This paper describes two efficient methods aimed at modeling 3D scattering problems. The first method is the use of the Huygens' principle to reduce the size of finite element models. This principle allows the area of interest to be restricted to the very near field of the defect, for both the generation of the incident field and the modal decomposition of the scattered field. The second method consists of separating the 3D problem into two 2D problems for which the solutions are calculated and used to approximate the 3D solution. This can be used at low frequency-thickness products, where Lamb waves have a similar behavior to bulk waves. These two methods are presented briefly and compared on simple scattering cases.
Intersymbol Interference Investigations Using a 3D Time-Dependent Traveling Wave Tube Model
NASA Technical Reports Server (NTRS)
Kory, Carol L.; Andro, Monty
2002-01-01
For the first time, a time-dependent, physics-based computational model has been used to provide a direct description of the effects of the traveling wave tube amplifier (TWTA) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion; gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves. Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry and operating characteristics of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept- amplitude and/or swept-frequency data. First, the TWT model using the three dimensional (3D) electromagnetic code MAFIA is presented. Then, this comprehensive model is used to investigate approximations made in conventional TWT black-box models used in communication system level simulations. To quantitatively demonstrate the effects these approximations have on digital signal performance predictions, including intersymbol interference (ISI), the MAFIA results are compared to the system level analysis tool, Signal Processing Workstation (SPW), using high order modulation schemes including 16 and 64-QAM.
Intersymbol Interference Investigations Using a 3D Time-Dependent Traveling Wave Tube Model
NASA Technical Reports Server (NTRS)
Kory, Carol L.; Andro, Monty; Downey, Alan (Technical Monitor)
2001-01-01
For the first time, a physics based computational model has been used to provide a direct description of the effects of the TWT (Traveling Wave Tube) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion; gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves. Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept amplitude and/or swept frequency data. The fully three-dimensional (3D), time-dependent, TWT interaction model using the electromagnetic code MAFIA is presented. This model is used to investigate assumptions made in TWT black box models used in communication system level simulations. In addition, digital signal performance, including intersymbol interference (ISI), is compared using direct data input into the MAFIA model and using the system level analysis tool, SPW (Signal Processing Worksystem).
NASA Astrophysics Data System (ADS)
Hong, Seongik
In the 3D printing technology, the research for using various materials has been performing. In this research work, 3D printable high viscous materials are suggested as one of the solutions for problems in the traditional 3D printing technology. First, Cu-Ag coreshell was synthesized as a functional material. In terms of the reaction rate, reaction rate limiting step was defined as a fundamental research, and then prepared Cu-Ag coreshell was printed and analyzed. Second, the high viscous Cu paste was prepared and then metal 3D printed structure was fabricated by using new printing method. In the synthesis of Cu-Ag coreshell, different sizes of Cu particle, 2μm and 100nm were used, and when 2μm Cu was applied, the reaction rate was limited by film diffusion control. However, when 100nm Cu was applied, reaction rate was controlled by CuO film and the rate of the reaction, which includes removing CuO film in the solution, is limited by chemical reaction control. The shape of Cu-Ag particle is spherical in the 2μm Cu condition and dendrite shape in the 100nm Cu condition respectively. The conductivity of Cu-Ag coreshell paste increased as increasing content of coreshell particle in the paste and sintering temperature. In order to print high viscous metal paste, the high viscous Cu paste was printed by using screw extruder, and the viscosity of Cu paste was measured as a fundamental research. As increasing wt.% of Cu in the paste, the viscosity also increased. In addition, the shrinkage factor was reduced by increasing wt.% of Cu in the paste. An optimized printing condition for the high viscous material was obtained, and by using this condition, 3D metal structure was fabricated. The final product was heat treated and polished. Through these processes, a fine quality of metal 3D structure was printed.
Calaminici, Patrizia; Janetzko, Florian; Köster, Andreas M; Mejia-Olvera, Roberto; Zuniga-Gutierrez, Bernardo
2007-01-28
Density functional theory optimized basis sets for gradient corrected functionals for 3d transition metal atoms are presented. Double zeta valence polarization and triple zeta valence polarization basis sets are optimized with the PW86 functional. The performance of the newly optimized basis sets is tested in atomic and molecular calculations. Excitation energies of 3d transition metal atoms, as well as electronic configurations, structural parameters, dissociation energies, and harmonic vibrational frequencies of a large number of molecules containing 3d transition metal elements, are presented. The obtained results are compared with available experimental data as well as with other theoretical data from the literature.
Episcopic 3D Imaging Methods: Tools for Researching Gene Function
Weninger, Wolfgang J; Geyer, Stefan H
2008-01-01
This work aims at describing episcopic 3D imaging methods and at discussing how these methods can contribute to researching the genetic mechanisms driving embryogenesis and tissue remodelling, and the genesis of pathologies. Several episcopic 3D imaging methods exist. The most advanced are capable of generating high-resolution volume data (voxel sizes from 0.5x0.5x1 µm upwards) of small to large embryos of model organisms and tissue samples. Beside anatomy and tissue architecture, gene expression and gene product patterns can be three dimensionally analyzed in their precise anatomical and histological context with the aid of whole mount in situ hybridization or whole mount immunohistochemical staining techniques. Episcopic 3D imaging techniques were and are employed for analyzing the precise morphological phenotype of experimentally malformed, randomly produced, or genetically engineered embryos of biomedical model organisms. It has been shown that episcopic 3D imaging also fits for describing the spatial distribution of genes and gene products during embryogenesis, and that it can be used for analyzing tissue samples of adult model animals and humans. The latter offers the possibility to use episcopic 3D imaging techniques for researching the causality and treatment of pathologies or for staging cancer. Such applications, however, are not yet routine and currently only preliminary results are available. We conclude that, although episcopic 3D imaging is in its very beginnings, it represents an upcoming methodology, which in short terms will become an indispensable tool for researching the genetic regulation of embryo development as well as the genesis of malformations and diseases. PMID:19452045
Importance of a 3D forward modeling tool for surface wave analysis methods
NASA Astrophysics Data System (ADS)
Pageot, Damien; Le Feuvre, Mathieu; Donatienne, Leparoux; Philippe, Côte; Yann, Capdeville
2016-04-01
Since a few years, seismic surface waves analysis methods (SWM) have been widely developed and tested in the context of subsurface characterization and have demonstrated their effectiveness for sounding and monitoring purposes, e.g., high-resolution tomography of the principal geological units of California or real time monitoring of the Piton de la Fournaise volcano. Historically, these methods are mostly developed under the assumption of semi-infinite 1D layered medium without topography. The forward modeling is generally based on Thomson-Haskell matrix based modeling algorithm and the inversion is driven by Monte-Carlo sampling. Given their efficiency, SWM have been transfered to several scale of which civil engineering structures in order to, e.g., determine the so-called V s30 parameter or assess other critical constructional parameters in pavement engineering. However, at this scale, many structures may often exhibit 3D surface variations which drastically limit the efficiency of SWM application. Indeed, even in the case of an homogeneous structure, 3D geometry can bias the dispersion diagram of Rayleigh waves up to obtain discontinuous phase velocity curves which drastically impact the 1D mean velocity model obtained from dispersion inversion. Taking advantages of high-performance computing center accessibility and wave propagation modeling algorithm development, it is now possible to consider the use of a 3D elastic forward modeling algorithm instead of Thomson-Haskell method in the SWM inversion process. We use a parallelized 3D elastic modeling code based on the spectral element method which allows to obtain accurate synthetic data with very low numerical dispersion and a reasonable numerical cost. In this study, we choose dike embankments as an illustrative example. We first show that their longitudinal geometry may have a significant effect on dispersion diagrams of Rayleigh waves. Then, we demonstrate the necessity of 3D elastic modeling as a forward
Pan Xiaoyin; Slamet, Marlina; Sahni, Viraht
2010-04-15
We extend our prior work on the construction of variational wave functions {psi} that are functionals of functions {chi}:{psi}={psi}[{chi}] rather than simply being functions. In this manner, the space of variations is expanded over those of traditional variational wave functions. In this article we perform the constrained search over the functions {chi} chosen such that the functional {psi}[{chi}] satisfies simultaneously the constraints of normalization and the exact expectation value of an arbitrary single- or two-particle Hermitian operator, while also leading to a rigorous upper bound to the energy. As such the wave function functional is accurate not only in the region of space in which the principal contributions to the energy arise but also in the other region of the space represented by the Hermitian operator. To demonstrate the efficacy of these ideas, we apply such a constrained search to the ground state of the negative ion of atomic hydrogen H{sup -}, the helium atom He, and its positive ions Li{sup +} and Be{sup 2+}. The operators W whose expectations are obtained exactly are the sum of the single-particle operators W={Sigma}{sub i}r{sub i}{sup n},n=-2,-1,1,2, W={Sigma}{sub i{delta}}(r{sub i}), W=-(1/2){Sigma}{sub i{nabla}i}{sup 2}, and the two-particle operators W={Sigma}{sub n}u{sup n},n=-2,-1,1,2, where u=|r{sub i}-r{sub j}|. Comparisons with the method of Lagrangian multipliers and of other constructions of wave-function functionals are made. Finally, we present further insights into the construction of wave-function functionals by studying a previously proposed construction of functionals {psi}[{chi}] that lead to the exact expectation of arbitrary Hermitian operators. We discover that analogous to the solutions of the Schroedinger equation, there exist {psi}[{chi}] that are unphysical in that they lead to singular values for the expectations. We also explain the origin of the singularity.
3D genome structure modeling by Lorentzian objective function.
Trieu, Tuan; Cheng, Jianlin
2016-11-29
The 3D structure of the genome plays a vital role in biological processes such as gene interaction, gene regulation, DNA replication and genome methylation. Advanced chromosomal conformation capture techniques, such as Hi-C and tethered conformation capture, can generate chromosomal contact data that can be used to computationally reconstruct 3D structures of the genome. We developed a novel restraint-based method that is capable of reconstructing 3D genome structures utilizing both intra-and inter-chromosomal contact data. Our method was robust to noise and performed well in comparison with a panel of existing methods on a controlled simulated data set. On a real Hi-C data set of the human genome, our method produced chromosome and genome structures that are consistent with 3D FISH data and known knowledge about the human chromosome and genome, such as, chromosome territories and the cluster of small chromosomes in the nucleus center with the exception of the chromosome 18. The tool and experimental data are available at https://missouri.box.com/v/LorDG.
3D analysis of interaction of Lamb waves with defects in loaded steel plates.
Kazys, R; Mazeika, L; Barauskas, R; Raisutis, R; Cicenas, V; Demcenko, A
2006-12-22
The objective of the research presented here is the investigation of the interaction of guided waves with welds, defects and other non-uniformities in steel plates loaded by liquid. The investigation has been performed using numerical simulation for 2D and 3D cases by the finite differences method, finite element method and measurement of 3D distributions of acoustic fields. Propagation of the S(0) mode in a steel plate and its interaction with non-uniformities was investigated. It was shown that using the measured leaky wave signals in the water loading of the steel plate and by application of signal processing, the 3D ultrasonic field structure inside and outside of the plate can be reconstructed. The presence of leaky wave signals over the defect caused by the mode conversion of Lamb waves has been proved using the numerical modelling and experimental investigations. The developed signal and data processing enables to visualise dynamics of ultrasonic fields over the plate, and also to estimate spatial positions of defects inside the steel plates.
Earthquake source tensor inversion with the gCAP method and 3D Green's functions
NASA Astrophysics Data System (ADS)
Zheng, J.; Ben-Zion, Y.; Zhu, L.; Ross, Z.
2013-12-01
We develop and apply a method to invert earthquake seismograms for source properties using a general tensor representation and 3D Green's functions. The method employs (i) a general representation of earthquake potency/moment tensors with double couple (DC), compensated linear vector dipole (CLVD), and isotropic (ISO) components, and (ii) a corresponding generalized CAP (gCap) scheme where the continuous wave trains are broken into Pnl and surface waves (Zhu & Ben-Zion, 2013). For comparison, we also use the waveform inversion method of Zheng & Chen (2012) and Ammon et al. (1998). Sets of 3D Green's functions are calculated on a grid of 1 km3 using the 3-D community velocity model CVM-4 (Kohler et al. 2003). A bootstrap technique is adopted to establish robustness of the inversion results using the gCap method (Ross & Ben-Zion, 2013). Synthetic tests with 1-D and 3-D waveform calculations show that the source tensor inversion procedure is reasonably reliable and robust. As initial application, the method is used to investigate source properties of the March 11, 2013, Mw=4.7 earthquake on the San Jacinto fault using recordings of ~45 stations up to ~0.2Hz. Both the best fitting and most probable solutions include ISO component of ~1% and CLVD component of ~0%. The obtained ISO component, while small, is found to be a non-negligible positive value that can have significant implications for the physics of the failure process. Work on using higher frequency data for this and other earthquakes is in progress.
Gravitational Wave Signals from 2D and 3D Core Collapse Supernova Explosions
NASA Astrophysics Data System (ADS)
Yakunin, Konstantin; Mezzacappa, Anthony; Marronetti, Pedro; Bruenn, Stephen; Hix, W. Raphael; Lentz, Eric J.; Messer, O. E. Bronson; Harris, J. Austin; Endeve, Eirik; Blondin, John
2016-03-01
We study two- and three-dimensional (2D and 3D) core-collapse supernovae (CCSN) using our first-principles CCSN simulations performed with the neutrino hydrodynamics code CHIMERA. The following physics is included: Newtonian hydrodynamics with a nuclear equation of state capable of describing matter in both NSE and non-NSE, MGFLD neutrino transport with realistic neutrino interactions, an effective GR gravitational potential, and a nuclear reaction network. Both our 2D and 3D models achieve explosion, which in turn enables us to determine their complete gravitational wave signals. In this talk, we present them, and we analyze the similarities and differences between the 2D and 3D signals.
A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling
NASA Astrophysics Data System (ADS)
Ren, Zhengyong; Kalscheuer, Thomas; Greenhalgh, Stewart; Maurer, Hansruedi
2013-08-01
We have developed a novel goal-oriented adaptive mesh refinement approach for finite-element methods to model plane wave electromagnetic (EM) fields in 3-D earth models based on the electric field differential equation. To handle complicated models of arbitrary conductivity, magnetic permeability and dielectric permittivity involving curved boundaries and surface topography, we employ an unstructured grid approach. The electric field is approximated by linear curl-conforming shape functions which guarantee the divergence-free condition of the electric field within each tetrahedron and continuity of the tangential component of the electric field across the interior boundaries. Based on the non-zero residuals of the approximated electric field and the yet to be satisfied boundary conditions of continuity of both the normal component of the total current density and the tangential component of the magnetic field strength across the interior interfaces, three a-posterior error estimators are proposed as a means to drive the goal-oriented adaptive refinement procedure. The first a-posterior error estimator relies on a combination of the residual of the electric field, the discontinuity of the normal component of the total current density and the discontinuity of the tangential component of the magnetic field strength across the interior faces shared by tetrahedra. The second a-posterior error estimator is expressed in terms of the discontinuity of the normal component of the total current density (conduction plus displacement current). The discontinuity of the tangential component of the magnetic field forms the third a-posterior error estimator. Analytical solutions for magnetotelluric (MT) and radiomagnetotelluric (RMT) fields impinging on a homogeneous half-space model are used to test the performances of the newly developed goal-oriented algorithms using the above three a-posterior error estimators. A trapezoidal topographical model, using normally incident EM waves
An endoscope designed with 3D measurement functions
NASA Astrophysics Data System (ADS)
Zhu, Qi-hai; Li, Zheng-lin; Wang, Li-qiang; Zhao, Chang-Ming; Xu, Peng
2016-10-01
The endoscopic system is widely used in medical and industrial areas, but how to realize the high-precision three-dimensional measurement in the limited space scale still faces many challenges. A method based on the four-step phase-shifting structured light illumination is proposed in this paper for endoscopic 3D measurements. Structured light of which the adjacent phase shift is 90 degrees is generated by the different parts of the time-sharing lighting stripe grating of the optical fiber bundle; CMOS camera is used to collect four structured light images with the phase shift. Finally, the method of four-step phase-shifting is used to demodulate 3D information from the images, and a relative measurement accuracy of 95% within the range of 15-200mm can be obtained. The endoscope with a field of view of 90 degrees, a image resolution of 1280 * 800 and 3D depth calculation time of 0.2 seconds has advantages of simple structure, large field of view, high accuracy and good real-time measurement.
The Vajont disaster: a 3D numerical simulation for the slide and the waves
NASA Astrophysics Data System (ADS)
Rubino, Angelo; Androsov, Alexey; Vacondio, Renato; Zanchettin, Davide; Voltzinger, Naum
2016-04-01
A very high resolution O(5 m), 3D hydrostatic nonlinear numerical model was used to simulate the dynamics of both the slide and the surface waves produced during the Vajont disaster (north Italy, 1963), one of the major landslide-induced tsunamis ever documented. Different simulated wave phenomena like, e.g., maximum run-up on the opposite shore, maximum height, and water velocity were analyzed and compared with data available in literature, including the results of a fully 3D simulation obtained with a Smoothed Particle Hydrodynamic code. The difference between measured and simulated after-slide bathymetries was calculated and used in an attempt to quantify the relative magnitude and extension of rigid and fluid motion components during the event.
Full-wave Moment Tensor and Tomographic Inversions Based on 3D Strain Green Tensor
2010-01-31
G. Jahnke, Wave propagation in 3D spherical sections: effects of subduction zones , Phys. Earth Planet. Inter., 132, 219-234, 2002. Komastitsch, D...is at scales smaller than the Fresnel zone . For example, a 1-Hz P/Pn wave recorded by a receiver ~1000 km from the source has a Fresnel zone width...approach, Eos Trans. AGU, 89(53), Fall Meet. Suppl., abstract T11E-06 Invited, 2008b. Sigloch, K., N. McQuarrie, G. Nolet, Two-stage subduction
3D Modeling of Antenna Driven Slow Waves Excited by Antennas Near the Plasma Edge
NASA Astrophysics Data System (ADS)
Smithe, David; Jenkins, Thomas
2016-10-01
Prior work with the 3D finite-difference time-domain (FDTD) plasma and sheath model used to model ICRF antennas in fusion plasmas has highlighted the possibility of slow wave excitation at the very low end of the SOL density range, and thus the prudent need for a slow-time evolution model to treat SOL density modifications due to the RF itself. At higher frequency, the DIII-D helicon antenna has much easier access to a parasitic slow wave excitation, and in this case the Faraday screen provides the dominant means of controlling the content of the launched mode, with antenna end-effects remaining a concern. In both cases, the danger is the same, with the slow-wave propagating into a lower-hybrid resonance layer a short distance ( cm) away from the antenna, which would parasitically absorb power, transferring energy to the SOL edge plasma, primarily through electron-neutral collisions. We will present 3D modeling of antennas at both ICRF and helicon frequencies. We've added a slow-time evolution capability for the SOL plasma density to include ponderomotive force driven rarefaction from the strong fields in the vicinity of the antenna, and show initial application to NSTX antenna geometry and plasma configurations. The model is based on a Scalar Ponderomotive Potential method, using self-consistently computed local field amplitudes from the 3D simulation.
NASA Astrophysics Data System (ADS)
Ebrahimi, Farzad; Barati, Mohammad Reza
2016-09-01
This article examines the application of nonlocal strain gradient elasticity theory to wave dispersion behavior of a size-dependent functionally graded (FG) nanobeam in thermal environment. The theory contains two scale parameters corresponding to both nonlocal and strain gradient effects. A quasi-3D sinusoidal beam theory considering shear and normal deformations is employed to present the formulation. Mori-Tanaka micromechanical model is used to describe functionally graded material properties. Hamilton's principle is employed to obtain the governing equations of nanobeam accounting for thickness stretching effect. These equations are solved analytically to find the wave frequencies and phase velocities of the FG nanobeam. It is indicated that wave dispersion behavior of FG nanobeams is significantly affected by temperature rise, nonlocality, length scale parameter and material composition.
NASA Astrophysics Data System (ADS)
Lieske, Mike; Schlurmann, Torsten
2016-04-01
INTRODUCTION & MOTIVATION The design of structures in coastal and offshore areas and their maintenance are key components of coastal protection. Usually, assessments of processes and loads on coastal structures are derived from experiments with flow and wave parameters in separate physical models. However, Peregrin (1976) already points out that processes in natural shallow coastal waters flow and sea state processes do not occur separately, but influence each other nonlinearly. Kemp & Simons (1982) perform 2D laboratory tests and study the interactions between a turbulent flow and following waves. They highlight the significance of wave-induced changes in the current properties, especially in the mean flow profiles, and draw attention to turbulent fluctuations and bottom shear stresses. Kemp & Simons (1983) also study these processes and features with opposing waves. Studies on the wave-current interaction in three-dimensional space for a certain wave height, wave period and water depth were conducted by MacIver et al. (2006). The research focus is set on the investigation of long-crested waves on obliquely opposing and following currents in the new 3D wave-current basin. METHODOLOGY In a first step the flow analysis without waves is carried out and includes measurements of flow profiles in the sweet spot of the basin at predefined measurement positions. Five measuring points in the water column have been delineated in different water depths in order to obtain vertical flow profiles. For the characterization of the undisturbed flow properties in the basin, an uniformly distributed flow was generated in the wave basin. In the second step wave analysis without current, the unidirectional wave propagation and wave height were investigated for long-crested waves in intermediate wave conditions. In the sweet spot of the wave basin waves with three different wave directions, three wave periods and uniform wave steepness were examined. For evaluation, we applied a common
Wave Propagation from Complex 3D Sources using the Representation Theorem
2009-09-30
wavenumber integration. The equations for the Green’s functions for surface waves are given by Bache et al. (1982). The Green’s functions for the...Green’s functions for body waves are generated by a procedure similar to that described by Bache and Harkrider (1976) using a saddle point...931 951. Bache , T. C. and D. G Harkrider (1976). The Body Waves Due to a General Seismic source in a Layered Earth Model, Bull. Seism. Soc. Am. 66
3D resolution tests of two-plane wave approach using synthetic seismograms
NASA Astrophysics Data System (ADS)
Ceylan, S.; Larmat, C. S.; Sandvol, E. A.
2012-12-01
Two-plane wave tomography (TPWT) is becoming a standard approach to obtain fundamental mode Rayleigh wave phase velocities for a variety of tectonic settings. A recent study by Ceylan et al. (2012) has applied this method to eastern Tibet, using data from INDEPTH-IV and Namche-Barwa seismic experiments. The TPWT assumes that distortion of wavefronts at each station can be expressed as the sum of two plane waves. However, there is currently no robust or complete resolution test for TPWT, to address its limitations such as wavefront healing. In this study, we test the capabilities of TPWT and resolution of INDEPTH-IV seismic experiment, by performing 3D resolution tests using synthetic seismograms. Utilizing SPECFEM3D software, we compute synthetic data sets resolving periods down to ~30 s. We implement a checkerboard upper mantle (for depths between 50 and 650 km) with variable cell sizes, superimposed to PREM as the background model. We then calculate fundamental mode surface wave phase velocities using TPWT for periods between 33-143 seconds, using synthetic seismograms computed from our three dimensional hypothetical model. Assuming a constant Poisson's ratio, we use partial derivatives from Saito (1988) to invert for shear wave velocities. We show that the combination of TPWT and Saito (1988) methods is capable of retrieving anomalies down to depths of ~200 km for Rayleigh waves. Below these depths, we observe evidence of both lateral and vertical smearing. We also find that the traditional method for estimating the resolution of TPWT consistently overestimates phase velocity resolutions. Love waves exhibit adequate resolution down to depths of ~100 km. At depths greater than 100 km, smearing is more evident in SH wave results than those of SV waves. Increased smearing of SH waves is most probably due to propagation characteristics and shallower sensitivity of Love waves. Our results imply that TPWT can be applied to Love waves, making future investigations of
Ultra wide band millimeter wave holographic ``3-D`` imaging of concealed targets on mannequins
Collins, H.D.; Hall, T.E.; Gribble, R.P.
1994-08-01
Ultra wide band (chirp frequency) millimeter wave ``3-D`` holography is a unique technique for imaging concealed targets on human subjects with extremely high lateral and depth resolution. Recent ``3-D`` holographic images of full size mannequins with concealed weapons illustrate the efficacy of this technique for airport security. A chirp frequency (24 GHz to 40 GHz) holographic system was used to construct extremely high resolution images (optical quality) using polyrod antenna in a bi-static configuration using an x-y scanner. Millimeter wave chirp frequency holography can be simply described as a multi-frequency detection and imaging technique where the target`s reflected signals are decomposed into discrete frequency holograms and reconstructed into a single composite ``3-D`` image. The implementation of this technology for security at airports, government installations, etc., will require real-time (video rate) data acquisition and computer image reconstruction of large volumetric data sets. This implies rapid scanning techniques or large, complex ``2-D`` arrays and high speed computing for successful commercialization of this technology.
A global 3-D MHD model of the solar wind with Alfven waves
NASA Technical Reports Server (NTRS)
Usmanov, A. V.
1995-01-01
A fully three-dimensional solar wind model that incorporates momentum and heat addition from Alfven waves is developed. The proposed model upgrades the previous one by considering self-consistently the total system consisting of Alfven waves propagating outward from the Sun and the mean polytropic solar wind flow. The simulation region extends from the coronal base (1 R(sub s) out to beyond 1 AU. The fully 3-D MHD equations written in spherical coordinates are solved in the frame of reference corotating with the Sun. At the inner boundary, the photospheric magnetic field observations are taken as boundary condition and wave energy influx is prescribed to be proportional to the magnetic field strength. The results of the model application for several time intervals are presented.
3D P-Wave Velocity Structure of the Deep Galicia Rifted Margin
NASA Astrophysics Data System (ADS)
Bayrakci, Gaye; Minshull, Timothy; Davy, Richard; Sawyer, Dale; Klaeschen, Dirk; Papenberg, Cord; Reston, Timothy; Shillington, Donna; Ranero, Cesar
2015-04-01
The combined wide-angle reflection-refraction and multi-channel seismic (MCS) experiment, Galicia 3D, was carried out in 2013 at the Galicia rifted margin in the northeast Atlantic Ocean, west of Spain. The main geological features within the 64 by 20 km (1280 km²) 3D box investigated by the survey are the peridotite ridge (PR), the fault bounded, rotated basement blocks and the S reflector, which has been interpreted to be a low angle detachment fault. 44 short period four-component ocean bottom seismometers and 28 ocean bottom hydrophones were deployed in the 3D box. 3D MCS profiles sampling the whole box were acquired with two airgun arrays of 3300 cu.in. fired alternately every 37.5 m. We present the results from 3D first-arrival time tomography that constrains the P-wave velocity in the 3D box, for the entire depth sampled by reflection data. Results are validated by synthetic tests and by the comparison with Galicia 3D MCS lines. The main outcomes are as follows: 1- The 3.5 km/s iso-velocity contour mimics the top of the acoustic basement observed on MCS profiles. Block bounding faults are imaged as velocity contrasts and basement blocks exhibit 3D topographic variations. 2- On the southern profiles, the top of the PR rises up to 5.5 km depth whereas, 20 km northward, its basement expression (at 6.5 km depth) nearly disappears. 3- The 6.5 km/s iso-velocity contour matches the topography of the S reflector where the latter is visible on MCS profiles. Within a depth interval of 0.6 km (in average), velocities beneath the S reflector increase from 6.5 km/s to 7 km/s, which would correspond to a decrease in the degree of serpentinization from ~45 % to ~30 % if these velocity variations are caused solely by variations in hydration. At the intersections between the block bounding normal faults and the S reflector, this decrease happens over a larger depth interval (> 1 km), suggesting that faults act as conduit for the water flow in the upper mantle.
3D P-Wave Velocity Structure of the Deep Galicia Rifted Margin
NASA Astrophysics Data System (ADS)
Bayrakci, G.; Minshull, T. A.; Davy, R. G.; Sawyer, D. S.; Klaeschen, D.; Papenberg, C. A.; Reston, T. J.; Shillington, D. J.; Ranero, C. R.
2014-12-01
The combined wide-angle reflection-refraction and multi-channel seismic (MCS) experiment, Galicia 3D, was carried out in 2013 at the Galicia rifted margin in the northeast Atlantic Ocean, west of Spain. The main geological features within the 64 by 20 km (1280 km²) 3D box investigated by the survey are the peridotite ridge (PR), the fault bounded, rotated basement blocks and the S reflector, which has been interpreted to be a low angle detachment fault. 44 short period four-component ocean bottom seismometers and 28 ocean bottom hydrophones were deployed in the 3D box. 3D MCS profiles sampling the whole box were acquired with two airgun arrays of 3300 cu.in. fired alternately every 37.5 m. We present the results from 3D first-arrival time tomography that constrains the P-wave velocity in the 3D box, for the entire depth sampled by reflection data. Results are validated by synthetic tests and by the comparison with Galicia 3D MCS lines. The main outcomes are as follows: 1- The 3.5 km/s iso-velocity contour mimics the top of the acoustic basement observed on MCS profiles. Block bounding faults are imaged as velocity contrasts and basement blocks exhibit 3D topographic variations. 2- On the southern profiles, the top of the PR rises up to 5.5 km depth whereas, 20 km northward, its basement expression (at 6.5 km depth) nearly disappears. 3- The 6.5 km/s iso-velocity contour matches the topography of the S reflector where the latter is visible on MCS profiles. Within a depth interval of 0.6 km (in average), velocities beneath the S reflector increase from 6.5 km/s to 7 km/s, which would correspond to a decrease in the degree of serpentinization from ~45 % to ~30 % if these velocity variations are caused solely by variations in hydration. At the intersections between the block bounding normal faults and the S reflector, this decrease happens over a larger depth interval (> 1 km), suggesting that faults act as conduit for the water flow in the upper mantle.
NASA Astrophysics Data System (ADS)
Montiel, F.; Squire, V. A.
2013-12-01
A new ocean wave/sea-ice interaction model is proposed that simulates how a directional wave spectrum evolves as it travels through a realistic marginal ice zone (MIZ), where wave/ice dynamics are entirely governed by coherent conservative wave scattering effects. Field experiments conducted by Wadhams et al. (1986) in the Greenland Sea generated important data on wave attenuation in the MIZ and, particularly, on whether the wave spectrum spreads directionally or collimates with distance from the ice edge. The data suggest that angular isotropy, arising from multiple scattering by ice floes, occurs close to the edge and thenceforth dominates wave propagation throughout the MIZ. Although several attempts have been made to replicate this finding theoretically, including by the use of numerical models, none have confronted this problem in a 3D MIZ with fully randomised floe distribution properties. We construct such a model by subdividing the discontinuous ice cover into adjacent infinite slabs of finite width parallel to the ice edge. Each slab contains an arbitrary (but finite) number of circular ice floes with randomly distributed properties. Ice floes are modeled as thin elastic plates with uniform thickness and finite draught. We consider a directional wave spectrum with harmonic time dependence incident on the MIZ from the open ocean, defined as a continuous superposition of plane waves traveling at different angles. The scattering problem within each slab is then solved using Graf's interaction theory for an arbitrary incident directional plane wave spectrum. Using an appropriate integral representation of the Hankel function of the first kind (see Cincotti et al., 1993), we map the outgoing circular wave field from each floe on the slab boundaries into a directional spectrum of plane waves, which characterizes the slab reflected and transmitted fields. Discretizing the angular spectrum, we can obtain a scattering matrix for each slab. Standard recursive
NASA Astrophysics Data System (ADS)
Xie, J.; Yang, Y.; Ni, S.; Zhao, K.
2015-12-01
In the past decade, ambient noise tomography (ANT) has become an estimated method to construct the earth's interior structures thanks to its advantage in extracting surface waves from cross-correlations of ambient noise without using earthquake data. However, most of previous ambient noise tomography studies concentrate on short and intermediate periods (<50sec) due to the dominant energy of the microseism at these periods. Studies of long period surface waves from cross-correlation of ambient noise are limited. In this study, we verify the accuracy of the long period (50-250sec) surface wave (Rayleigh wave) from ambient noise by comparing both dispersion curves and seismic waveforms from ambient noise with those from earthquake records quantitatively. After that, we calculate vertical-vertical cross-correlation functions among more than1800 USArray Transportable Array stations and extract high quality interstation phase velocity dispersion curves from them at 10-200 sec periods. Then, we adopt a finite frequency ambient noise tomography method based on Born approximation to obtain high resolution phase velocity maps using the obtained dispersion measurements at 10-150 sec periods. Afterward, we extract local dispersion curves from these dispersion maps and invert them for 1D shear wave velocity profiles at individual grids using a Bayesian Monte Carlo method. Finally, a 3D shear velocity model is constructed by assembling all the 1D Vs profiles. Our 3D model is overall similar to other models constructed using earthquake surface waves and body waves. In summary, we demonstrate that the long period surface waves can be extracted from ambient noise, and the long period dispersion measurements from ambient noise are as accurate as those from earthquake data and can be used to construct 3D lithospheric structure from surface down to lithosphere/asthenosphere depths.
NASA Astrophysics Data System (ADS)
Scheffold, Frank
2014-08-01
To characterize the structural and dynamic properties of soft materials and small particles, information on the relevant mesoscopic length scales is required. Such information is often obtained from traditional static and dynamic light scattering (SLS/DLS) experiments in the single scattering regime. In many dense systems, however, these powerful techniques frequently fail due to strong multiple scattering of light. Here I will discuss some experimental innovations that have emerged over the last decade. New methods such as 3D static and dynamic light scattering (3D LS) as well as diffusing wave spectroscopy (DWS) can cover a much extended range of experimental parameters ranging from dilute polymer solutions, colloidal suspensions to extremely opaque viscoelastic emulsions.
Total-Field Technique for 3-D Modeling of Short Period Teleseismic Waves
NASA Astrophysics Data System (ADS)
Monteiller, V.; Beller, S.; Operto, S.; Nissen-Meyer, T.; Tago Pacheco, J.; Virieux, J.
2014-12-01
The massive development of dense seismic arrays and the rapid increase in computing capacity allow today to consider application of full waveform inversion of teleseismic data for high-resolution lithospheric imaging. We present an hybrid numerical method that allows for the modellingof short period teleseismic waves in 3D lithospheric target with both the discontinuous Galerkin finite elements method and finite difference method, opening the possibility to perform waveform inversion of seismograms recorded by dense regional broadband arrays. However, despite the supercomputer ability, the forward-problem remains expensive at global scale for teleseismic configuration especially when 3D numerical methods are considered. In order to perform the forward problem in a reasonable amount of time, we reduce the computational domain in which full waveform modelling is performed. We define a 3D regional domain located below the seismological network that is embedded in a homogeneous background or axisymmetric model, in which the seismic wavefield can be computed efficiently. The background wavefield is used to compute the full wavefield in the 3D regional domain using the so-called total-field/scattered-field technique. This method relies on the decomposition of the wavefield into a background and a scattered wavefields. The computational domain is subdivided into three sub-domains: an outer domain formed by the perfectly-matched absorbing layers, an intermediate domain in which only the outgoing wavefield scattered by the lithospheric heterogeneities is computed, and the inner domain formed by the lithospheric target in which the full wavefield is computed. In this study, we shall present simulations in realistic lithospheric target when the axisymetric background wavefield is computed with the AxiSEM softwave and the 3D simulation in lithospheric target model is performed with the discontinuous Galerkin or finite difference method.
3D frequency-domain finite-difference modeling of acoustic wave propagation
NASA Astrophysics Data System (ADS)
Operto, S.; Virieux, J.
2006-12-01
We present a 3D frequency-domain finite-difference method for acoustic wave propagation modeling. This method is developed as a tool to perform 3D frequency-domain full-waveform inversion of wide-angle seismic data. For wide-angle data, frequency-domain full-waveform inversion can be applied only to few discrete frequencies to develop reliable velocity model. Frequency-domain finite-difference (FD) modeling of wave propagation requires resolution of a huge sparse system of linear equations. If this system can be solved with a direct method, solutions for multiple sources can be computed efficiently once the underlying matrix has been factorized. The drawback of the direct method is the memory requirement resulting from the fill-in of the matrix during factorization. We assess in this study whether representative problems can be addressed in 3D geometry with such approach. We start from the velocity-stress formulation of the 3D acoustic wave equation. The spatial derivatives are discretized with second-order accurate staggered-grid stencil on different coordinate systems such that the axis span over as many directions as possible. Once the discrete equations were developed on each coordinate system, the particle velocity fields are eliminated from the first-order hyperbolic system (following the so-called parsimonious staggered-grid method) leading to second-order elliptic wave equations in pressure. The second-order wave equations discretized on each coordinate system are combined linearly to mitigate the numerical anisotropy. Secondly, grid dispersion is minimized by replacing the mass term at the collocation point by its weighted averaging over all the grid points of the stencil. Use of second-order accurate staggered- grid stencil allows to reduce the bandwidth of the matrix to be factorized. The final stencil incorporates 27 points. Absorbing conditions are PML. The system is solved using the parallel direct solver MUMPS developed for distributed
3D Numerical Simulation on the Sloshing Waves Excited by the Seismic Shacking
NASA Astrophysics Data System (ADS)
Zhang, Lin; Wu, Tso-Ren
2016-04-01
In the event of 2015 Nepal earthquake, a video clip broadcasted worldwide showed a violent water spilling in a hotel swimming pool. This sloshing phenomenon indicates a potential water loss in the sensitive facilities, e.g. the spent fuel pools in nuclear power plant, has to be taken into account carefully under the consideration of seismic-induced ground acceleration. In the previous studies, the simulation of sloshing mainly focused on the pressure force on the structure by using a simplified Spring-Mass Method developed in the field of solid mechanics. However, restricted by the assumptions of plane water surface and limited wave height, significant error will be made in evaluating the amount of water loss in the tank. In this paper, the computational fluid dynamical model, Splash3D, was adopted for studying the sloshing problem accurately. Splash3D solved 3D Navier-Stokes Equation directly with Large-Eddy Simulation (LES) turbulent closure. The Volume-of-fluid (VOF) method with piecewise linear interface calculation (PLIC) was used to track the complex breaking water surface. The time series acceleration of a design seismic was loaded to excite the water. With few restrictions from the assumptions, the accuracy of the simulation results were improved dramatically. A series model validations were conducted by compared to a 2D theoretical solution, and a 3D experimental data. Good comparisons can be seen. After the validation, we performed the simulation for considering a sloshing case in a rectangular water tank with a dimension of 12 m long, 8 m wide, 8 m deep, which contained water with 7 m in depth. The seismic movement was imported by considering time-series acceleration in three dimensions, which were about 0.5 g to 1.2 g in the horizontal directions, and 0.3 g to 1 g in the vertical direction. We focused the discussions on the kinematics of the water surface, wave breaking, velocity field, pressure field, water force on the side walls, and, most
3D crustal structure of the Alpine belt and foreland basins as imaged by ambient-noise surface wave
NASA Astrophysics Data System (ADS)
Molinari, Irene; Morelli, Andrea; Cardi, Riccardo; Boschi, Lapo; Poli, Piero; Kissling, Edi
2016-04-01
We derive a 3-D crustal structure (S wave velocity) underneath northern Italy and the wider Alpine region, from an extensive data set of measurements of Rayleigh-wave phase- and group-velocities from ambient noise correlation among all seismographic stations available to date in the region, via a constrained tomographic inversion made to honor detailed active source reflection/refraction profiles and other geological information. We first derive a regional-scale surface wave tomography from ambient-noise-based phase- and group- surface wave velocity observations (Verbeke et al., 2012). Our regional 3D model (Molinari et al., 2015) shows the low velocity area beneath the Po Plain and the Molasse basin; the contrast between the low-velocity crust of the Adriatic domain and the high-velocity crust of the Tyrrhenian domain is clearly seen, as well as an almost uniform crystalline crust beneath the Alpine belt. However, higher frequency data can be exploited to achieve higher resolution images of the Po Plain and Alpine foreland 3D crustal structure. We collected and analyze one year of noise records (2011) of ~100 North Italy seismic broadband stations, we derive the Green functions between each couple of stations and we measure the phase- and group-Rayleigh wave velocity. We conduct a suite of linear least squares inversion of both phase- and group-velocity data, resulting in 2-D maps of Rayleigh-wave phase and group velocity at periods between 3 and 40s with a resolution of 0.1x0.1 degrees. The maps are then inverted to get the 3D structure with unprecedented details. We present here our results, we compare them with other studies, and we discuss geological/geodynamical implications. We believe that such a model stands for the most up-to-date seismological information on the crustal structure of the Alpine belt and foreland basins, and it can represent a reliable reference for further, more detailed, studies to come, based on the high seismograph station density
Wang, Tao; Green, Ryan; Nair, Rajesh Ramakrishnan; Howell, Mark; Mohapatra, Subhra; Guldiken, Rasim; Mohapatra, Shyam Sundar
2015-01-01
Detection and quantification of cell viability and growth in two-dimensional (2D) and three-dimensional (3D) cell cultures commonly involve harvesting of cells and therefore requires a parallel set-up of several replicates for time-lapse or dose–response studies. Thus, developing a non-invasive and touch-free detection of cell growth in longitudinal studies of 3D tumor spheroid cultures or of stem cell regeneration remains a major unmet need. Since surface acoustic waves (SAWs) permit mass loading-based biosensing and have been touted due to their many advantages including low cost, small size and ease of assembly, we examined the potential of SAW-biosensing to detect and quantify cell growth. Herein, we demonstrate that a shear horizontal-surface acoustic waves (SH-SAW) device comprising two pairs of resonators consisting of interdigital transducers and reflecting fingers can be used to quantify mass loading by the cells in suspension as well as within a 3D cell culture platform. A 3D COMSOL model was built to simulate the mass loading response of increasing concentrations of cells in suspension in the polydimethylsiloxane (PDMS) well in order to predict the characteristics and optimize the design of the SH-SAW biosensor. The simulated relative frequency shift from the two oscillatory circuit systems (one of which functions as control) were found to be concordant to experimental data generated with RAW264.7 macrophage and A549 cancer cells. In addition, results showed that SAW measurements per se did not affect viability of cells. Further, SH-SAW biosensing was applied to A549 cells cultured on a 3D electrospun nanofiber scaffold that generate tumor spheroids (tumoroids) and the results showed the device's ability to detect changes in tumor spheroid growth over the course of eight days. Taken together, these results demonstrate the use of SH-SAW device for detection and quantification of cell growth changes over time in 2D suspension cultures and in 3D cell
Wang, Tao; Green, Ryan; Nair, Rajesh Ramakrishnan; Howell, Mark; Mohapatra, Subhra; Guldiken, Rasim; Mohapatra, Shyam Sundar
2015-12-19
Detection and quantification of cell viability and growth in two-dimensional (2D) and three-dimensional (3D) cell cultures commonly involve harvesting of cells and therefore requires a parallel set-up of several replicates for time-lapse or dose-response studies. Thus, developing a non-invasive and touch-free detection of cell growth in longitudinal studies of 3D tumor spheroid cultures or of stem cell regeneration remains a major unmet need. Since surface acoustic waves (SAWs) permit mass loading-based biosensing and have been touted due to their many advantages including low cost, small size and ease of assembly, we examined the potential of SAW-biosensing to detect and quantify cell growth. Herein, we demonstrate that a shear horizontal-surface acoustic waves (SH-SAW) device comprising two pairs of resonators consisting of interdigital transducers and reflecting fingers can be used to quantify mass loading by the cells in suspension as well as within a 3D cell culture platform. A 3D COMSOL model was built to simulate the mass loading response of increasing concentrations of cells in suspension in the polydimethylsiloxane (PDMS) well in order to predict the characteristics and optimize the design of the SH-SAW biosensor. The simulated relative frequency shift from the two oscillatory circuit systems (one of which functions as control) were found to be concordant to experimental data generated with RAW264.7 macrophage and A549 cancer cells. In addition, results showed that SAW measurements per se did not affect viability of cells. Further, SH-SAW biosensing was applied to A549 cells cultured on a 3D electrospun nanofiber scaffold that generate tumor spheroids (tumoroids) and the results showed the device's ability to detect changes in tumor spheroid growth over the course of eight days. Taken together, these results demonstrate the use of SH-SAW device for detection and quantification of cell growth changes over time in 2D suspension cultures and in 3D cell
3D Simulation of an Audible Ultrasonic Electrolarynx Using Difference Waves
Mills, Patrick; Zara, Jason
2014-01-01
A total laryngectomy removes the vocal folds which are fundamental in forming voiced sounds that make speech possible. Although implanted prosthetics are commonly used in developed countries, simple handheld vibrating electrolarynxes are still common worldwide. These devices are easy to use but suffer from many drawbacks including dedication of a hand, mechanical sounding voice, and sound leakage. To address some of these drawbacks, we introduce a novel electrolarynx that uses vibro-acoustic interference of dual ultrasonic waves to generate an audible fundamental frequency. A 3D simulation of the principles of the device is presented in this paper. PMID:25401965
A 3D unstructured non-hydrostatic ocean model for internal waves
NASA Astrophysics Data System (ADS)
Ai, Congfang; Ding, Weiye
2016-10-01
A 3D non-hydrostatic model is developed to compute internal waves. A novel grid arrangement is incorporated in the model. This not only ensures the homogenous Dirichlet boundary condition for the non-hydrostatic pressure can be precisely and easily imposed but also renders the model relatively simple in its discretized form. The Perot scheme is employed to discretize horizontal advection terms in the horizontal momentum equations, which is based on staggered grids and has the conservative property. Based on previous water wave models, the main works of the present paper are to (1) utilize a semi-implicit, fractional step algorithm to solve the Navier-Stokes equations (NSE); (2) develop a second-order flux-limiter method satisfying the max-min property; (3) incorporate a density equation, which is solved by a high-resolution finite volume method ensuring mass conservation and max-min property based on a vertical boundary-fitted coordinate system; and (4) validate the developed model by using four tests including two internal seiche waves, lock-exchange flow, and internal solitary wave breaking. Comparisons of numerical results with analytical solutions or experimental data or other model results show reasonably good agreement, demonstrating the model's capability to resolve internal waves relating to complex non-hydrostatic phenomena.
Source Process of the 1923 Kanto Earthquake Using New Fault Geometry and 3-D Green's Functions
NASA Astrophysics Data System (ADS)
Kobayashi, R.; Koketsu, K.
2005-12-01
The September 1, 1923, Kanto earthquake caused severe damage and more than 100,000 fatalities in the Tokyo metropolitan area. This earthquake is an interplate event along the Sagami trough where the Philippine Sea plate is subducting beneath a continental plate. We have investigated the source process of this earthquake using the geodetic, teleseismic, and strong motion data (Kobayashi and Koketsu, 2005). The resultant slip distributions show that two asperities (areas of large slips) are located around the base of the Izu peninsula and the Uraga channel. In 2002 and 2003, four seismic surveys were carried out to determine crustal structures and fault locations in the Kanto region (Sato et al., 2005). The seismic reflections from the surface of the Philippine Sea slab suggested that the slab surface should be shallower than the previous models (e.g., Ishida, 1992; Matsu'ura et al., 1980). The fault model of Kobayashi and Koketsu (2005) was also based on Matsu'ura et al. (1980). In this study, we adopt new fault geometry consistent with the result of the reflection surveys and perform another source process inversion. The new slip distribution showed that the western asperity moved from the Uraga channel to the tip of the Miura peninsula, while the western asperity did not move considerably. Green's functions that Kobayashi and Koketsu (2005) used were calculated in a halfspace for geodetic data or in a 1-D model for strong motions. However, the real structure in the Kanto region is three-dimensionally complex as suggested by the geographical setting and seismic surveys. In fact, Kobayashi and Koketsu (2005) showed that the long coda of the observed seismogram at Hongo, Tokyo, was not reproduced in the synthetic one. The forward modeling with a 3-D structure (Sato et al., 1999) suggested that surface waves excited along the boundary between the Kanto mountains and Kanto basin can explain the large coda. Thus we calculate 3-D Green's functions for the strong motion
Vdovin V.L.
2005-08-15
In this report we describe theory and 3D full wave code description for the wave excitation, propagation and absorption in 3-dimensional (3D) stellarator equilibrium high beta plasma in ion cyclotron frequency range (ICRF). This theory forms a basis for a 3D code creation, urgently needed for the ICRF heating scenarios development for the operated LHD, constructed W7-X, NCSX and projected CSX3 stellarators, as well for re evaluation of ICRF scenarios in operated tokamaks and in the ITER . The theory solves the 3D Maxwell-Vlasov antenna-plasma-conducting shell boundary value problem in the non-orthogonal flux coordinates ({Psi}, {theta}, {var_phi}), {Psi} being magnetic flux function, {theta} and {var_phi} being the poloidal and toroidal angles, respectively. All basic physics, like wave refraction, reflection and diffraction are self consistently included, along with the fundamental ion and ion minority cyclotron resonances, two ion hybrid resonance, electron Landau and TTMP absorption. Antenna reactive impedance and loading resistance are also calculated and urgently needed for an antenna -generator matching. This is accomplished in a real confining magnetic field being varying in a plasma major radius direction, in toroidal and poloidal directions, through making use of the hot dense plasma wave induced currents with account to the finite Larmor radius effects. We expand the solution in Fourier series over the toroidal ({var_phi}) and poloidal ({theta}) angles and solve resulting ordinary differential equations in a radial like {Psi}-coordinate by finite difference method. The constructed discretization scheme is divergent-free one, thus retaining the basic properties of original equations. The Fourier expansion over the angle coordinates has given to us the possibility to correctly construct the ''parallel'' wave number k{sub //}, and thereby to correctly describe the ICRF waves absorption by a hot plasma. The toroidal harmonics are tightly coupled with each
Investigation of 3D surface acoustic waves in granular media with 3-color digital holography
NASA Astrophysics Data System (ADS)
Leclercq, Mathieu; Picart, Pascal; Penelet, Guillaume; Tournat, Vincent
2017-01-01
This paper reports the implementation of digital color holography to investigate elastic waves propagating along a layer of a granular medium. The holographic set-up provides simultaneous recording and measurement of the 3D dynamic displacement at the surface. Full-field measurements of the acoustic amplitude and phase at different excitation frequencies are obtained. It is shown that the experimental data can be used to obtain the dispersion curve of the modes propagating in this granular medium layer. The experimental dispersion curve and that obtained from a finite element modeling of the problem are found to be in good agreement. In addition, full-field images of the interaction of an acoustic wave guided in the granular layer with a buried object are also shown.
Wave optics theory and 3-D deconvolution for the light field microscope.
Broxton, Michael; Grosenick, Logan; Yang, Samuel; Cohen, Noy; Andalman, Aaron; Deisseroth, Karl; Levoy, Marc
2013-10-21
Light field microscopy is a new technique for high-speed volumetric imaging of weakly scattering or fluorescent specimens. It employs an array of microlenses to trade off spatial resolution against angular resolution, thereby allowing a 4-D light field to be captured using a single photographic exposure without the need for scanning. The recorded light field can then be used to computationally reconstruct a full volume. In this paper, we present an optical model for light field microscopy based on wave optics, instead of previously reported ray optics models. We also present a 3-D deconvolution method for light field microscopy that is able to reconstruct volumes at higher spatial resolution, and with better optical sectioning, than previously reported. To accomplish this, we take advantage of the dense spatio-angular sampling provided by a microlens array at axial positions away from the native object plane. This dense sampling permits us to decode aliasing present in the light field to reconstruct high-frequency information. We formulate our method as an inverse problem for reconstructing the 3-D volume, which we solve using a GPU-accelerated iterative algorithm. Theoretical limits on the depth-dependent lateral resolution of the reconstructed volumes are derived. We show that these limits are in good agreement with experimental results on a standard USAF 1951 resolution target. Finally, we present 3-D reconstructions of pollen grains that demonstrate the improvements in fidelity made possible by our method.
Wave optics theory and 3-D deconvolution for the light field microscope
Broxton, Michael; Grosenick, Logan; Yang, Samuel; Cohen, Noy; Andalman, Aaron; Deisseroth, Karl; Levoy, Marc
2013-01-01
Light field microscopy is a new technique for high-speed volumetric imaging of weakly scattering or fluorescent specimens. It employs an array of microlenses to trade off spatial resolution against angular resolution, thereby allowing a 4-D light field to be captured using a single photographic exposure without the need for scanning. The recorded light field can then be used to computationally reconstruct a full volume. In this paper, we present an optical model for light field microscopy based on wave optics, instead of previously reported ray optics models. We also present a 3-D deconvolution method for light field microscopy that is able to reconstruct volumes at higher spatial resolution, and with better optical sectioning, than previously reported. To accomplish this, we take advantage of the dense spatio-angular sampling provided by a microlens array at axial positions away from the native object plane. This dense sampling permits us to decode aliasing present in the light field to reconstruct high-frequency information. We formulate our method as an inverse problem for reconstructing the 3-D volume, which we solve using a GPU-accelerated iterative algorithm. Theoretical limits on the depth-dependent lateral resolution of the reconstructed volumes are derived. We show that these limits are in good agreement with experimental results on a standard USAF 1951 resolution target. Finally, we present 3-D reconstructions of pollen grains that demonstrate the improvements in fidelity made possible by our method. PMID:24150383
Effect of background rotation on the evolution of 3D internal gravity wave beams
NASA Astrophysics Data System (ADS)
Fan, Boyu; Akylas, T. R.
2016-11-01
The effect of background rotation on the 3D propagation of internal gravity wave beams (IGWB) is studied, assuming that variations in the along-beam and transverse directions are of long length scale relative to the beam width. The present study generalizes the asymptotic model of KA (Kataoka & Akylas 2015) who considered the analogous problem in the absence of rotation. It is shown that the role of mean vertical vorticity in the earlier analysis is now taken by the flow mean potential vorticity (MPV). Specifically, 3D variations enable resonant transfer of energy to the flow MPV, resulting in strong nonlinear coupling between a 3D IGWB and its induced mean flow. This coupling mechanism is governed by a system of two nonlinear equations of the same form as those derived in KA. Accordingly, the induced mean flow features a purely inviscid modulational component, as well as a viscous one akin to acoustic streaming; the latter grows linearly with time for a quasi-steady IGWB. On the other hand, owing to background rotation, the induced mean flow in the vicinity of the IGWB is no longer purely horizontal and develops an asymmetric behavior. Supported by NSF.
3D Plenoptic PIV Measurements of a Shock Wave Boundary Layer Interaction
NASA Astrophysics Data System (ADS)
Thurow, Brian; Bolton, Johnathan; Arora, Nishul; Alvi, Farrukh
2016-11-01
Plenoptic particle image velocimetry (PIV) is a relatively new technique that uses the computational refocusing capability of a single plenoptic camera and volume illumination with a double-pulsed light source to measure the instantaneous 3D/3C velocity field of a flow field seeded with particles. In this work, plenoptic PIV is used to perform volumetric velocity field measurements of a shock-wave turbulent boundary layer interaction (SBLI). Experiments were performed in a Mach 2.0 flow with the SBLI produced by an unswept fin at 15°angle of attack. The measurement volume was 38 x 25 x 32 mm3 and illuminated with a 400 mJ/pulse Nd:YAG laser with 1.7 microsecond inter-pulse time. Conventional planar PIV measurements along two planes within the volume are used for comparison. 3D visualizations of the fin generated shock and subsequent SBLI are presented. The growth of the shock foot and separation region with increasing distance from the fin tip is observed and agrees with observations made using planar PIV. Instantaneous images depict 3D fluctuations in the position of the shock foot from one image to the next. The authors acknowledge the support of the Air Force Office of Scientific Research.
3D P-wave Velocity Structure Beneath the Eastern Canadian Shield and Northern Appalachian Region
NASA Astrophysics Data System (ADS)
Villemaire, M.; Darbyshire, F. A.; Bastow, I. D.
2010-12-01
Previous seismic studies of the upper mantle of the Canadian Shield have indicated some low-velocity anomalies within the cratonic lithosphere in the Abitibi-Grenville region. The lack of seismograph station coverage to the east and south-east of the studied area prevented definition of the 3D geometry of these anomalies. Adding new stations from the province of Quebec and from the northeastern United States allows us to carry out new studies of the P-wave velocity structure of the upper mantle, in order to better understand the complexity of the region and the interaction of the lithosphere with possible thermal anomalies in the underlying mantle. We analysed teleseismic P wave arrivals from almost 200 earthquakes, recorded at 45 stations deployed across the provinces of Quebec and Ontario and across the northeastern US. The relative arrival times of teleseismic P waves across the array were measured using the cross-correlation method of VanDecar & Crosson (1990). The travel time data were then inverted to estimate the 3D P-wave velocity structure beneath the region, using the least-squares tomographic inversion code of VanDecar (1991). The model shows some interesting features. We see a diffuse low-velocity structure beneath New-England that extends to at least 500 km depth, and that may be related to the Appalachian Mountain belt. There is also a linear low-velocity structure, flanked by higher velocities, perpendicular to the Grenville Front, and along the Ottawa Valley. We interpret this feature as a mantle signature of the Great Meteor Hotspot track. We have looked for systematic differences between the mantle underlying the Archean Superior craton and the Proterozoic Grenville Province but did not find a significant difference in the upper mantle. We investigate the role of thermal and compositional effects to interpret the velocity models and to relate the patterns of the anomalies to past and present tectonic structures.
Poroelastic Wave Propagation With a 3D Velocity-Stress-Pressure Finite-Difference Algorithm
NASA Astrophysics Data System (ADS)
Aldridge, D. F.; Symons, N. P.; Bartel, L. C.
2004-12-01
Seismic wave propagation within a three-dimensional, heterogeneous, isotropic poroelastic medium is numerically simulated with an explicit, time-domain, finite-difference algorithm. A system of thirteen, coupled, first-order, partial differential equations is solved for the particle velocity vector components, the stress tensor components, and the pressure associated with solid and fluid constituents of the two-phase continuum. These thirteen dependent variables are stored on staggered temporal and spatial grids, analogous to the scheme utilized for solution of the conventional velocity-stress system of isotropic elastodynamics. Centered finite-difference operators possess 2nd-order accuracy in time and 4th-order accuracy in space. Seismological utility is enhanced by an optional stress-free boundary condition applied on a horizontal plane representing the earth's surface. Absorbing boundary conditions are imposed on the flanks of the 3D spatial grid via a simple wavefield amplitude taper approach. A massively parallel computational implementation, utilizing the spatial domain decomposition strategy, allows investigation of large-scale earth models and/or broadband wave propagation within reasonable execution times. Initial algorithm testing indicates that a point force density and/or moment density source activated within a poroelastic medium generates diverging fast and slow P waves (and possibly an S-wave)in accord with Biot theory. Solid and fluid particle velocities are in-phase for the fast P-wave, whereas they are out-of-phase for the slow P-wave. Conversions between all wave types occur during reflection and transmission at interfaces. Thus, although the slow P-wave is regarded as difficult to detect experimentally, its presence is strongly manifest within the complex of waves generated at a lithologic or fluid boundary. Very fine spatial and temporal gridding are required for high-fidelity representation of the slow P-wave, without inducing excessive
3-D Inverse Teleseismic Scattered Wave Imaging using the Kirchhoff Approximation
NASA Astrophysics Data System (ADS)
Liu, K.; Levander, A.
2012-04-01
We have developed a 3-D teleseismic imaging technique for scattered elastic wavefields using the Kirchhoff approximation. Kirchhoff migration/inversion have been well developed in exploration seismology within the inverse scattering framework (e.g. Miller et al., 1987; Beylkin and Burridge, 1990) to image subsurface structure that generates secondary wavefields caused by localized heterogeneities. Application of this method in global seismology has been largely limited to 2-D images made with 1-D reference models due to high computational cost and the lack of adequately dense receiver arrays (Bostock, 2002, Poppeliers and Pavlis, 2003; Frederiksen and Revenaugh, 2004; Cao et al., 2010). The deployment of the USArray Transportable and Flexible arrays in the United States and dense array recordings in other countries motivate developing teleseismic scattered wavefield imaging with the Kirchhoff approximation for 3-D velocity models for both scalar and vector wavefields to improve upper mantle imaging. Following Bostock's development of the 2-D problem (2002), we derive the 3-D P-to-S scattering inversion formula by phrasing the inverse problem in terms of the generalized Radon transform (GRT) and singular functions of discontinuity surfaces. In the forward scattering modeling, we extend the method to utilize a 3-D migration velocity model by calculating 3-D finite-difference traveltimes, backprojected from the receivers using an eikonal solver. To demonstrate the relative accuracy of the inversion, we examine several synthetic cases with a variety of discontinuity surfaces (sinuous, dipping, dome- and crater-shaped discontinuity interfaces, point scatterers, etc.). The Kirchhoff GRT imaging can successfully recover the shapes of these structures very well. We compare our Kirchhoff approximation imaging with the Born-approximate results, as well as the common-conversion point (CCP) stacked receiver function imaging for the various synthetic cases, and show a field
NASA Astrophysics Data System (ADS)
Yoshizawa, K.
2014-10-01
A new radially anisotropic shear wave speed model for the Australasian region is constructed from multi-mode phase dispersion of Love and Rayleigh waves. An automated waveform fitting technique based on a global optimization with the Neighbourhood Algorithm allows the exploitation of large numbers of three-component broad-band seismograms to extract path-specific dispersion curves covering the entire continent. A 3-D shear wave model is constructed including radial anisotropy from a set of multi-mode phase speed maps for both Love and Rayleigh waves. These maps are derived from an iterative inversion scheme incorporating the effects of ray-path bending due to lateral heterogeneity, as well as the finite frequency of the surface waves for each mode. The new S wave speed model exhibits major tectonic features of this region that are in good agreement with earlier shear wave models derived primarily from Rayleigh waves. The lateral variations of depth and thickness of the lithosphere-asthenosphere transition (LAT) are estimated from the isotropic (Voigt average) S wave speed model and its vertical gradient, which reveals correlations between the lateral variations of the LAT and radial anisotropy. The thickness of the LAT is very large beneath the Archean cratons in western Australia, whereas that in south Australia is thinner. The radial anisotropy model shows faster SH wave speed than SV beneath eastern Australia and the Coral Sea at the lithospheric depth. The faster SH anomaly in the lithosphere is also seen in the suture zone between the three cratonic blocks of Australia. One of the most conspicuous features of fast SH anisotropy is found in the asthenosphere beneath the central Australia, suggesting anisotropy induced by shear flow in the asthenosphere beneath the fast drifting Australian continent.
NASA Astrophysics Data System (ADS)
Plotnikov, Illya; Vourlidas, Angelos; Tylka, Allan J.; Pinto, Rui; Rouillard, Alexis; Tirole, Margot
2016-07-01
Identifying the physical mechanisms that produce the most energetic particles is a long-standing observational and theoretical challenge in astrophysics. Strong pressure waves have been proposed as efficient accelerators both in the solar and astrophysical contexts via various mechanisms such as diffusive-shock/shock-drift acceleration and betatron effects. In diffusive-shock acceleration, the efficacy of the process relies on shock waves being super-critical or moving several times faster than the characteristic speed of the medium they propagate through (a high Alfven Mach number) and on the orientation of the magnetic field upstream of the shock front. High-cadence, multipoint imaging using the NASA STEREO, SOHO and SDO spacecrafts now permits the 3-D reconstruction of pressure waves formed during the eruption of coronal mass ejections. Using these unprecedented capabilities, some recent studies have provided new insights on the timing and longitudinal extent of solar energetic particles, including the first derivations of the time-dependent 3-dimensional distribution of the expansion speed and Mach numbers of coronal shock waves. We will review these recent developments by focusing on particle events that occurred between 2011 and 2015. These new techniques also provide the opportunity to investigate the enigmatic long-duration gamma ray events.
Wave Propagation from Complex 3D Sources Using the Representation Theorem
2008-09-30
functions for surface waves are given by Bache et al. (1982). The Green’s functions for the complete seismograms are computed using a ring load source...procedure similar to that described by Bache and Harkrider (1976), using a saddle point approximation to calculate a far-field plane wave for a given takeoff...space, Part II, Bull. Seism. Soc. Am. 73: 931-951. Bache , T. C. and D. G. Harkrider (1976). The body waves due to a general seismic source in a layered
NASA Astrophysics Data System (ADS)
Chen, Kai-Xun; Chen, Po-Fei; Liang, Wen-Tzong; Chen, Li-Wei; Gung, YuanCheng
2015-04-01
The Yilan Plain (YP) in NE Taiwan locates on the western YP of the Okinawa Trough and displays high geothermal gradients with abundant hot springs, likely resulting from magmatism associated with the back-arc spreading as attested by the offshore volcanic island (Kueishantao). YP features NS distinctive characteristics that the South YP exhibits thin top sedimentary layer, high on-land seismicity and significant SE movements, relative those of the northern counterpart. A dense network (~2.5 km station interval) of 89 Texan instruments was deployed in Aug. 2014, covering most of the YP and its vicinity. The ray path coverage density of each 0.015 degree cells are greater than 150 km that could provide the robustness assessment of tomographic results. We analyze ambient noise signals to invert a high-resolution 3D S-wave model for shallow velocity structures in and around YP. The aim is to investigate the velocity anomalies corresponding to geothermal resources and the NS geological distinctions aforementioned. We apply the Welch's method to generate empirical Rayleigh wave Green's functions between two stations records of continuous vertical components. The group velocities of thus derived functions are then obtained by the multiple-filter analysis technique measured at the frequency range between 0.25 and 1 Hz. Finally, we implement a wavelet-based multi-scale parameterization technique to construct 3D model of S-wave velocity. Our first month results exhibit low velocity in the plain, corresponding existing sediments, those of whole YP show low velocity offshore YP and those of high-resolution south YP reveal stark velocity contrast across the Sanshin fault. Key words: ambient seismic noises, Welch's method, S-wave, Yilan Plain
Wang, Jingyu; Chen, Fengling; Liu, Longwei; Qi, Chunxiao; Wang, Bingjie; Yan, Xiaojun; Huang, Chenyu; Hou, Wei; Zhang, Michael Q; Chen, Yang; Du, Yanan
2016-06-01
Accompanied by decreased hepatic functions, epithelial-mesenchymal transition (EMT) was observed in two dimensional (2D) cultured hepatocytes with elongated morphology, loss of polarity and weakened cell-cell interaction, while upgrading to 3D culture has been considered as significant improvement of its 2D counterpart for hepatocyte maintenance. Here we hypothesize that 3D culture enhances hepatic functions through regulating the EMT status. Biomaterial-engineered EMT was achieved by culturing HepaRG as 3D spheroids (SP-3D) or 3D stretched cells (ST-3D) in non-adherent and adherent micro-scaffold respectively. In SP-3D, constrained EMT of HepaRG, a hepatic stem cell line, as represented by increased epithelial markers and decreased mesenchymal markers, was echoed by improved hepatic functions. To investigate the relationship between EMT status and hepatic functions, time-series RNA-Seq and gene network analysis were used for comparing different cell culture models, which identified histone deacetylases (HDACs) as key mediating factors. Protein analysis confirmed that high HDAC activity was correlated with high expression of Cadherin-1 (CDH1) and hepatic function genes, which were decreased upon HDAC inhibitor treatment in SP-3D, suggesting HDACs may play positive role in regulating EMT and hepatic functions. To illustrate the application of 3D micro-scaffold culture in drug safety evaluation, hepatotoxicity and metabolism assays of two hepatotoxins (i.e. N-acetyl-p-aminophenol and Doxorubicin) were performed and SP-3D showed more biomimetic toxicity response, indicating regulation of EMT as a vital consideration in designing 3D hepatocyte culture configuration.
Prediction of Tsunami Waves and Runup Generated by 3d Granular Landslides
NASA Astrophysics Data System (ADS)
Mohammed, F.; Fritz, H. M.
2008-12-01
attenuation. The measured wave speed of the leading wave corresponds well to the theoretical approximation of the solitary wave speed while the trailing waves are considerably slower in nature. Attenuation functions of the leading wave crest amplitude, the wave length and the time period were obtained to study the wave behavior in the near field and far field regions. The experimental data also provided the energy conversion rate between the landslide source and the generated tsunami waves. The slide characteristics measured in the experiment provides the landslide source for numerically modeling these landslide tsunamis. The measured wave data serves the validation and advancement of 3-dimensional numerical landslide tsunami and prediction models.
3D Discontinuous Galerkin elastic seismic wave modeling based upon a grid injection method
NASA Astrophysics Data System (ADS)
Monteiller, V.
2015-12-01
Full waveform inversion (FWI) is a seismic imaging method that estimates thesub-surface physical properties with a spatial resolution of the order of thewavelength. FWI is generally recast as the iterative optimization of anobjective function that measures the distance between modeled and recordeddata. In the framework of local descent methods, FWI requires to perform atleast two seismic modelings per source and per FWI iteration.Due to the resulting computational burden, applications of elastic FWI have been usuallyrestricted to 2D geometries. Despite the continuous growth of high-performancecomputing facilities, application of 3D elastic FWI to real-scale problemsremain computationally too expensive. To perform elastic seismic modeling with a reasonable amount of time, weconsider a reduced computational domain embedded in a larger background modelin which seismic sources are located. Our aim is to compute repeatedly thefull wavefield in the targeted domain after model alteration, once theincident wavefield has been computed once for all in the background model. Toachieve this goal, we use a grid injection method referred to as the Total-Field/Scattered-Field (TF/SF) technique in theelectromagnetic community. We implemented the Total-Field/Scattered-Field approach in theDiscontinuous Galerkin Finite Element method (DG-FEM) that is used to performmodeling in the local domain. We show how to interface the DG-FEM with any modeling engine (analytical solution, finite difference or finite elements methods) that is suitable for the background simulation. One advantage of the Total-Field/Scattered-Field approach is related to thefact that the scattered wavefield instead of the full wavefield enter thePMLs, hence making more efficient the absorption of the outgoing waves at theouter edges of the computational domain. The domain reduction in which theDG-FEM is applied allows us to use modest computational resources opening theway for high-resolution imaging by full
Richter, Benjamin; Hahn, Vincent; Bertels, Sarah; Claus, Tanja K; Wegener, Martin; Delaittre, Guillaume; Barner-Kowollik, Christopher; Bastmeyer, Martin
2017-02-01
The combination of three different photoresists into a single direct laser written 3D microscaffold permits functionalization with two bioactive full-length proteins. The cell-instructive microscaffolds consist of a passivating framework equipped with light activatable constituents featuring distinct protein-binding properties. This allows directed cell attachment of epithelial or fibroblast cells in 3D.
3D P and S Wave Velocity Structure and Tremor Locations in the Parkfield Region
NASA Astrophysics Data System (ADS)
Zeng, X.; Thurber, C. H.; Shelly, D. R.; Bennington, N. L.; Cochran, E. S.; Harrington, R. M.
2014-12-01
We have assembled a new dataset to refine the 3D seismic velocity model in the Parkfield region. The S arrivals from 184 earthquakes recorded by the Parkfield Experiment to Record MIcroseismicity and Tremor array (PERMIT) during 2010-2011 were picked by a new S wave picker, which is based on machine learning. 74 blasts have been assigned to four quarries, whose locations were identified with Google Earth. About 1000 P and S wave arrivals from these blasts at permanent seismic network were also incorporated. Low frequency earthquakes (LFEs) occurring within non-volcanic tremor (NVT) are valuable for improving the precision of NVT location and the seismic velocity model at greater depths. Based on previous work (Shelley and Hardebeck, 2010), waveforms of hundreds of LFEs in same family were stacked to improve signal qualify. In a previous study (McClement et al., 2013), stacked traces of more than 30 LFE families at the Parkfileld Array Seismic Observatory (PASO) have been picked. We expanded our work to include LFEs recorded by the PERMIT array. The time-frequency Phase Weight Stacking (tf-PWS) method was introduced to improve the stack quality, as direct stacking does not produce clear S-wave arrivals on the PERMIT stations. This technique uses the coherence of the instantaneous phase among the stacked signals to enhance the signal-to-noise ratio (SNR) of the stack. We found that it is extremely effective for picking LFE arrivals (Thurber et al., 2014). More than 500 P and about 1000 S arrivals from 58 LFE families were picked at the PERMIT and PASO arrays. Since the depths of LFEs are much deeper than earthquakes, we are able to extend model resolution to lower crustal depths. Both P and S wave velocity structure have been obtained with the tomoDD method. The result suggests that there is a low velocity zone (LVZ) in the lower crust and the location of the LVZ is consistent with the high conductivity zone beneath the southern segment of the Rinconada fault that
NASA Astrophysics Data System (ADS)
Bayrakci, Gaye; Minshull, Timothy A.; Davy, Richard G.; Karplus, Marianne S.; Kaeschen, Dirk; Papenberg, Cord; Krabbenhoeft, Anne; Sawyer, Dale; Reston, Timothy J.; Shillington, Donna J.; Ranero, César R.
2014-05-01
Galicia 3D, a reflection-refraction and long offset seismic experiment was carried out from May through September 2013, at the Galicia rifted margin (in the northeast Atlantic Ocean, west of Spain) as a collaboration between US, UK, German and Spanish groups. The 3D multichannel seismic acquisition conducted by R/V Marcus Langseth covered a 64 km by 20 km (1280 km2) zone where the main geological features are the Peridotite Ridge (PR), composed of serpentinized peridotite and thought be upper mantle exhumed to the seafloor during rifting, and the S reflector which has been interpreted to be a low angle detachment fault overlain by fault bounded, rotated, continental crustal blocks. In the 3D box, two airgun arrays of 3300 cu.in. were fired alternately (in flip-flop configuration) every 37.5 m. All shots are recorded by 44 short period four component ocean bottom seismometers (OBS) and 26 ocean bottom hydrophones (OBH) deployed and recovered by R/V Poseidon, as well as four 6 km hydrophone streamers with 12.5 m channel spacing towed by R/V Marcus Langseth. We present the preliminary results of the first arrival time tomography study which is carried out with a subset of the wide-angle dataset, in order to generate a 3D P-wave velocity volume for the entire depth sampled by the reflection data. After the relocation of OBSs and OBHs, an automatic first-arrival time picking approach is applied to a subset of the dataset, which comprises more than 5.5 million source-receiver pairs. Then, the first-arrival times are checked visually, in 3-dimensions. The a priori model used for the first-arrival time tomography is built up using information from previous seismic surveys carried out at the Galicia margin (e.g. ISE, 1997). The FAST algorithm of Zelt and Barton (1998) is used for the first-arrival time inversion. The 3D P-wave velocity volume can be used in interpreting the reflection dataset, as a starting point for migration, to quantify the thinning of the crustal layers
NASA Astrophysics Data System (ADS)
Tago, J.; Cruz-Atienza, V. M.; Etienne, V.; Virieux, J.; Benjemaa, M.; Sanchez-Sesma, F. J.
2010-12-01
Simulating any realistic seismic scenario requires incorporating physical basis into the model. Considering both the dynamics of the rupture process and the anelastic attenuation of seismic waves is essential to this purpose and, therefore, we choose to extend the hp-adaptive Discontinuous Galerkin finite-element method to integrate these physical aspects. The 3D elastodynamic equations in an unstructured tetrahedral mesh are solved with a second-order time marching approach in a high-performance computing environment. The first extension incorporates the viscoelastic rheology so that the intrinsic attenuation of the medium is considered in terms of frequency dependent quality factors (Q). On the other hand, the extension related to dynamic rupture is integrated through explicit boundary conditions over the crack surface. For this visco-elastodynamic formulation, we introduce an original discrete scheme that preserves the optimal code performance of the elastodynamic equations. A set of relaxation mechanisms describes the behavior of a generalized Maxwell body. We approximate almost constant Q in a wide frequency range by selecting both suitable relaxation frequencies and anelastic coefficients characterizing these mechanisms. In order to do so, we solve an optimization problem which is critical to minimize the amount of relaxation mechanisms. Two strategies are explored: 1) a least squares method and 2) a genetic algorithm (GA). We found that the improvement provided by the heuristic GA method is negligible. Both optimization strategies yield Q values within the 5% of the target constant Q mechanism. Anelastic functions (i.e. memory variables) are introduced to efficiently evaluate the time convolution terms involved in the constitutive equations and thus to minimize the computational cost. The incorporation of anelastic functions implies new terms with ordinary differential equations in the mathematical formulation. We solve these equations using the same order
A pair-conformation-dependent scoring function for evaluating 3D RNA-protein complex structures
Li, Haotian; Huang, Yangyu
2017-01-01
Computational prediction of RNA-protein complex 3D structures includes two basic steps: one is sampling possible structures and another is scoring the sampled structures to pick out the correct one. At present, constructing accurate scoring functions is still not well solved and the performances of the scoring functions usually depend on used benchmarks. Here we propose a pair-conformation-dependent scoring function, 3dRPC-Score, for 3D RNA-protein complex structure prediction by considering the nucleotide-residue pairs having the same energy if their conformations are similar, instead of the distance-only dependence of the most existing scoring functions. Benchmarking shows that 3dRPC-Score has a consistent performance in three test sets. PMID:28358834
A pair-conformation-dependent scoring function for evaluating 3D RNA-protein complex structures.
Li, Haotian; Huang, Yangyu; Xiao, Yi
2017-01-01
Computational prediction of RNA-protein complex 3D structures includes two basic steps: one is sampling possible structures and another is scoring the sampled structures to pick out the correct one. At present, constructing accurate scoring functions is still not well solved and the performances of the scoring functions usually depend on used benchmarks. Here we propose a pair-conformation-dependent scoring function, 3dRPC-Score, for 3D RNA-protein complex structure prediction by considering the nucleotide-residue pairs having the same energy if their conformations are similar, instead of the distance-only dependence of the most existing scoring functions. Benchmarking shows that 3dRPC-Score has a consistent performance in three test sets.
Optimal-tradeoff circular harmonic function filters for 3D target recognition
NASA Astrophysics Data System (ADS)
Vijaya Kumar, Bhagavatula V. K.; Xie, Chunyan; Mahalanobis, Abhijit
2003-09-01
3D target recognition is of significant interest because representing the object in 3D space couuld essentially provide a solution to pose variation and self-occlusion problems that are big challenges in 2D pattern recognition. Correlation filers have been used in a variety of 2D pattern matching applications and many correlation filter designs have been developed to handle problems such as rotations. Correlation filters also offer other benefits such as shift-invariance, graceful degradation and closed-form solutions. The 3D extension of correlation filter is a natural extension to handle 3D pattern recognition problem. In this paper, we propose a 3D correlation filter design method based on cylindrical circular harmonic function (CCHF) and use LADAR imagery to illustrate the good performance of CCHF filters.
Three-dimensional potential flows from functions of a 3D complex variable
NASA Technical Reports Server (NTRS)
Kelly, Patrick; Panton, Ronald L.; Martin, E. D.
1990-01-01
Potential, or ideal, flow velocities can be found from the gradient of an harmonic function. An ordinary complex valued analytic function can be written as the sum of two real valued functions, both of which are harmonic. Thus, 2D complex valued functions serve as a source of functions that describe two-dimensional potential flows. However, this use of complex variables has been limited to two-dimensions. Recently, a new system of three-dimensional complex variables has been developed at the NASA Ames Research Center. As a step toward application of this theory to the analysis of 3D potential flow, several functions of a three-dimensional complex variable have been investigated. The results for two such functions, the 3D exponential and 3D logarithm, are presented in this paper. Potential flows found from these functions are investigated. Important characteristics of these flows fields are noted.
Capurso, Daniel; Bengtsson, Henrik; Segal, Mark R.
2016-01-01
The spatial organization of the genome influences cellular function, notably gene regulation. Recent studies have assessed the three-dimensional (3D) co-localization of functional annotations (e.g. centromeres, long terminal repeats) using 3D genome reconstructions from Hi-C (genome-wide chromosome conformation capture) data; however, corresponding assessments for continuous functional genomic data (e.g. chromatin immunoprecipitation-sequencing (ChIP-seq) peak height) are lacking. Here, we demonstrate that applying bump hunting via the patient rule induction method (PRIM) to ChIP-seq data superposed on a Saccharomyces cerevisiae 3D genome reconstruction can discover ‘functional 3D hotspots’, regions in 3-space for which the mean ChIP-seq peak height is significantly elevated. For the transcription factor Swi6, the top hotspot by P-value contains MSB2 and ERG11 – known Swi6 target genes on different chromosomes. We verify this finding in a number of ways. First, this top hotspot is relatively stable under PRIM across parameter settings. Second, this hotspot is among the top hotspots by mean outcome identified by an alternative algorithm, k-Nearest Neighbor (k-NN) regression. Third, the distance between MSB2 and ERG11 is smaller than expected (by resampling) in two other 3D reconstructions generated via different normalization and reconstruction algorithms. This analytic approach can discover functional 3D hotspots and potentially reveal novel regulatory interactions. PMID:26869583
Anomalously low amplitude of S waves produced by the 3D structures in the lower mantle
NASA Astrophysics Data System (ADS)
To, Akiko; Capdeville, Yann; Romanowicz, Barbara
2016-07-01
Direct S and Sdiff phases with anomalously low amplitudes are recorded for the earthquakes in Papua New Guinea by seismographs in northern America. According to the prediction by a standard 1D model, the amplitudes are the lowest at stations in southern California, at a distance and azimuth of around 95° and 55°, respectively, from the earthquake. The amplitude anomaly is more prominent at frequencies higher than 0.03 Hz. We checked and ruled out the possibility of the anomalies appearing because of the errors in the focal mechanism used in the reference synthetic waveform calculations. The observed anomaly distribution changes drastically with a relatively small shift in the location of the earthquake. The observations indicate that the amplitude reduction is likely due to the 3D shear velocity (Vs) structure, which deflects the wave energy away from the original ray paths. Moreover, some previous studies suggested that some of the S and Sdiff phases in our dataset are followed by a prominent postcursor and show a large travel time delay, which was explained by placing a large ultra-low velocity zone (ULVZ) located on the core-mantle boundary southwest of Hawaii. In this study, we evaluated the extent of amplitude anomalies that can be explained by the lower mantle structures in the existing models, including the previously proposed ULVZ. In addition, we modified and tested some models and searched for the possible causes of low amplitudes. Full 3D synthetic waveforms were calculated and compared with the observations. Our results show that while the existing models explain the trends of the observed amplitude anomalies, the size of such anomalies remain under-predicted especially at large distances. Adding a low velocity zone, which is spatially larger and has less Vs reduction than ULVZ, on the southwest side of ULVZ, contributes to explain the low amplitudes observed at distances larger than 100° from the earthquake. The newly proposed low velocity zone
3dRNAscore: a distance and torsion angle dependent evaluation function of 3D RNA structures
Wang, Jian; Zhao, Yunjie; Zhu, Chunyan; Xiao, Yi
2015-01-01
Model evaluation is a necessary step for better prediction and design of 3D RNA structures. For proteins, this has been widely studied and the knowledge-based statistical potential has been proved to be one of effective ways to solve this problem. Currently, a few knowledge-based statistical potentials have also been proposed to evaluate predicted models of RNA tertiary structures. The benchmark tests showed that they can identify the native structures effectively but further improvements are needed to identify near-native structures and those with non-canonical base pairs. Here, we present a novel knowledge-based potential, 3dRNAscore, which combines distance-dependent and dihedral-dependent energies. The benchmarks on different testing datasets all show that 3dRNAscore are more efficient than existing evaluation methods in recognizing native state from a pool of near-native states of RNAs as well as in ranking near-native states of RNA models. PMID:25712091
Pt nanoparticles functionalized 3D SnO2 nanoflowers for gas sensor application
NASA Astrophysics Data System (ADS)
Liu, Yinglin; Huang, Jing; Yang, Jiedi; Wang, Shurong
2017-04-01
3D SnO2 nanoflowers (NFs) assembled by rod-like nanostructures were synthesized by a facile hydrothermal method only using simple and inexpensive SnCl4·5H2O and NaOH as the starting materials, without using any surfactants or templates. The as-synthesized 3D SnO2 NFs were further functionalized by Pt nanoparticles (NPs) by a simple ammonia precipitate method, and the derived Pt NP-functionalized 3D SnO2 NFs were further investigated for gas sensor application using ethanol as a probe gas. Obtained results showed that the Pt NP-functionalized 3D SnO2 NF sensor exhibited much higher response in comparison with pure SnO2 sensor, altogether with short response/recovery times and good reproducibility. The enhanced gas sensing performances could be attributed to spill-over effect of Pt NPs for promoting gas sensing reactions, the synergic electronic interaction between Pt NPs and SnO2 support, the high surface-to-volume ratio and good electron mobility of the 1D SnO2 nanorod units, and unique 3D hierarchical flower-like nanostructures. It is also expected that the as-prepared 3D SnO2 NFs and Pt NP-functionalized product can be used in other fields such as optoelectronic devices, Li-ion battery and dye sensitized solar cells.
Arai, Kenichi; Yoshida, Toshiko; Okabe, Motonori; Goto, Mitsuaki; Mir, Tanveer Ahmad; Soko, Chika; Tsukamoto, Yoshinari; Akaike, Toshihiro; Nikaido, Toshio; Zhou, Kaixuan; Nakamura, Makoto
2016-09-19
The development of new three-dimensional (3D) cell culture system that maintains the physiologically relevant signals of hepatocytes is essential in drug discovery and tissue engineering research. Conventional two-dimensional (2D) culture yields cell growth, proliferation, and differentiation. However, gene expression and signaling profiles can be different from in vivo environment. Here, we report the fabrication of a 3D culture system using an artificial scaffold and our custom-made inkjet 3D bioprinter as a new strategy for studying liver-specific functions of hepatocytes. We built a 3D culture platform for hepatocytes-attachment and formation of cell monolayer by interacting the galactose chain of galactosylated alginate gel (GA-gel) with asialoglycoprotein receptor (ASGPR) of hepatocytes. The 3D geometrical arrangement of cells was controlled by using 3D bioprinter, and cell polarity was controlled with the galactosylated hydrogels. The fabricated GA-gel was able to successfully promote adhesion of hepatocytes. To observe liver-specific functions and to mimic hepatic cord, an additional parallel layer of hepatocytes was generated using two gel sheets. These results indicated that GA-gel biomimetic matrices can be used as a 3D culture system that could be effective for the engineering of liver tissues. © 2016 Wiley Periodicals, Inc. J Biomed Mater Res Part A, 2016.
Seismic waves in 3-D: from mantle asymmetries to reliable seismic hazard assessment
NASA Astrophysics Data System (ADS)
Panza, Giuliano F.; Romanelli, Fabio
2014-10-01
A global cross-section of the Earth parallel to the tectonic equator (TE) path, the great circle representing the equator of net lithosphere rotation, shows a difference in shear wave velocities between the western and eastern flanks of the three major oceanic rift basins. The low-velocity layer in the upper asthenosphere, at a depth range of 120 to 200 km, is assumed to represent the decoupling between the lithosphere and the underlying mantle. Along the TE-perturbed (TE-pert) path, a ubiquitous LVZ, about 1,000-km-wide and 100-km-thick, occurs in the asthenosphere. The existence of the TE-pert is a necessary prerequisite for the existence of a continuous global flow within the Earth. Ground-shaking scenarios were constructed using a scenario-based method for seismic hazard analysis (NDSHA), using realistic and duly validated synthetic time series, and generating a data bank of several thousands of seismograms that account for source, propagation, and site effects. Accordingly, with basic self-organized criticality concepts, NDSHA permits the integration of available information provided by the most updated seismological, geological, geophysical, and geotechnical databases for the site of interest, as well as advanced physical modeling techniques, to provide a reliable and robust background for the development of a design basis for cultural heritage and civil infrastructures. Estimates of seismic hazard obtained using the NDSHA and standard probabilistic approaches are compared for the Italian territory, and a case-study is discussed. In order to enable a reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered, resulting in a new, very efficient, analytical procedure for computing the broadband seismic wave-field in a 3-D anelastic Earth model.
3D Anisotropic structure of the south-central Mongolia from Rayleigh and Love wave tomography
NASA Astrophysics Data System (ADS)
Yu, D.; Wu, Q.; Montagner, J. P.
2014-12-01
A better understanding of the geodynamics of the crust and mantle below Baikal-Mongolia is required to identify the role of mantle processes versus that of far-field tectonic effects from India-Asia collision. Anisotropy tomography can provide new perspective to the continental growth mechanism. In order to study the 3D anisotropic structure of the upper mantle in the south-central Mongolia, we collected the vertical and transverse components of seismograms recorded at 69 broadband seismic stations. We have measured inter-station phase velocities of 7181 Rayleigh waves and 901 Love waves using the frequency-time analysis of wavelet transformation method for the fundamental mode at period range 10~80s. The lateral phase velocity variations are computed by using a regionalization method. These phase velocities have been inverted to obtain the first anisotropic model including Sv velocities, azimuthal and radial anisotropy. The Middle Gobi is associated with low velocity. Based on the distribution of the Cenozoic basalts in the Middle Gobi, it refers that the low velocity anomaly is related to the Cenozoic volcanism. In the northern domain, near to Baikal zone, the azimuthal anisotropy is normal to the Baikal rift and consistent with the fast direction of previous SKS splitting measurements. In the South Gobi, north to Main Mongolian Lineament, the azimuthal anisotropy is NEE-SWW in the crust and NW-SE in the mantle. It indicates that the crust and mantle are decoupled. We propose that the crustal deformation is related to the far-field effects of India-Asia collision and that the mantle flow is correlated with the Baikal rift activity. Further study in process will provide more evidence and insight to better understand the geodynamics in this region.
NASA Astrophysics Data System (ADS)
Saini, A.; Christenson, C. W.; Khattab, T. A.; Wang, R.; Twieg, R. J.; Singer, K. D.
2017-01-01
In order to achieve a high capacity 3D optical data storage medium, a nonlinear or threshold writing process is necessary to localize data in the axial dimension. To this end, commercial multilayer discs use thermal ablation of metal films or phase change materials to realize such a threshold process. This paper addresses a threshold writing mechanism relevant to recently reported fluorescence-based data storage in dye-doped co-extruded multilayer films. To gain understanding of the essential physics, single layer spun coat films were used so that the data is easily accessible by analytical techniques. Data were written by attenuating the fluorescence using nanosecond-range exposure times from a 488 nm continuous wave laser overlapping with the single photon absorption spectrum. The threshold writing process was studied over a range of exposure times and intensities, and with different fluorescent dyes. It was found that all of the dyes have a common temperature threshold where fluorescence begins to attenuate, and the physical nature of the thermal process was investigated.
Nonhydrostatic granular flow over 3-D terrain: New Boussinesq-type gravity waves?
NASA Astrophysics Data System (ADS)
Castro-Orgaz, Oscar; Hutter, Kolumban; Giraldez, Juan V.; Hager, Willi H.
2015-01-01
granular mass flow is a basic step in the prediction and control of natural or man-made disasters related to avalanches on the Earth. Savage and Hutter (1989) pioneered the mathematical modeling of these geophysical flows introducing Saint-Venant-type mass and momentum depth-averaged hydrostatic equations using the continuum mechanics approach. However, Denlinger and Iverson (2004) found that vertical accelerations in granular mass flows are of the same order as the gravity acceleration, requiring the consideration of nonhydrostatic modeling of granular mass flows. Although free surface water flow simulations based on nonhydrostatic depth-averaged models are commonly used since the works of Boussinesq (1872, 1877), they have not yet been applied to the modeling of debris flow. Can granular mass flow be described by Boussinesq-type gravity waves? This is a fundamental question to which an answer is required, given the potential to expand the successful Boussinesq-type water theory to granular flow over 3-D terrain. This issue is explored in this work by generalizing the basic Boussinesq-type theory used in civil and coastal engineering for more than a century to an arbitrary granular mass flow using the continuum mechanics approach. Using simple test cases, it is demonstrated that the above question can be answered in the affirmative way, thereby opening a new framework for the physical and mathematical modeling of granular mass flow in geophysics, whereby the effect of vertical motion is mathematically included without the need of ad hoc assumptions.
Characterization of an SRF gun: a 3D full wave simulation
Wang, E.; Ben-Zvi, I.; Wang, J.
2011-03-28
We characterized a BNL 1.3GHz half-cell SRF gun is tested for GaAs photocathode. The gun already was simulated several years ago via two-dimensional (2D) numerical codes (i.e., Superfish and Parmela) with and without the beam. In this paper, we discuss our investigation of its characteristics using a three dimensional (3D) full-wave code (CST STUDIO SUITE{trademark}).The input/pickup couplers are sited symmetrically on the same side of the gun at an angle of 180{sup o}. In particular, the inner conductor of the pickup coupler is considerably shorter than that of the input coupler. We evaluated the cross-talk between the beam (trajectory) and the signal on the input coupler compared our findings with published results based on analytical models. The CST STUDIO SUITE{trademark} also was used to predict the field within the cavity; particularly, a combination of transient/eigenmode solvers was employed to accurately construct the RF field for the particles, which also includes the effects of the couplers. Finally, we explored the beam's dynamics with a particle in cell (PIC) simulation, validated the results and compare them with 2D code result.
Wearable 3-D Photoacoustic Tomography for Functional Brain Imaging in Behaving Rats
Tang, Jianbo; Coleman, Jason E.; Dai, Xianjin; Jiang, Huabei
2016-01-01
Understanding the relationship between brain function and behavior remains a major challenge in neuroscience. Photoacoustic tomography (PAT) is an emerging technique that allows for noninvasive in vivo brain imaging at micrometer-millisecond spatiotemporal resolution. In this article, a novel, miniaturized 3D wearable PAT (3D-wPAT) technique is described for brain imaging in behaving rats. 3D-wPAT has three layers of fully functional acoustic transducer arrays. Phantom imaging experiments revealed that the in-plane X-Y spatial resolutions were ~200 μm for each acoustic detection layer. The functional imaging capacity of 3D-wPAT was demonstrated by mapping the cerebral oxygen saturation via multi-wavelength irradiation in behaving hyperoxic rats. In addition, we demonstrated that 3D-wPAT could be used for monitoring sensory stimulus-evoked responses in behaving rats by measuring hemodynamic responses in the primary visual cortex during visual stimulation. Together, these results show the potential of 3D-wPAT for brain study in behaving rodents. PMID:27146026
3-D, bluff body drag estimation using a Green's function/Gram-Charlier series approach.
Barone, Matthew Franklin; De Chant, Lawrence Justin
2004-05-01
In this study, we describe the extension of the 2-d preliminary design bluff body drag estimation tool developed by De Chant to apply for 3-d flows. As with the 2-d method, the 3-d extension uses a combined approximate Green's function/Gram-Charlier series approach to retain the body geometry information. Whereas, the 2-d methodology relied solely upon the use of small disturbance theory for the inviscid flow field associated with the body of interest to estimate the near-field initial conditions, e.g. velocity defect, the 3-d methodology uses both analytical (where available) and numerical inviscid solutions. The defect solution is then used as an initial condition in an approximate 3-d Green's function solution. Finally, the Green's function solution is matched to the 3-d analog of the classical 2-d Gram-Charlier series and then integrated to yield the net form drag on the bluff body. Preliminary results indicate that drag estimates computed are of accuracy equivalent to the 2-d method for flows with large separation, i.e. less than 20% relative error. As was the lower dimensional method, the 3-d concept is intended to be a supplement to turbulent Navier-Stokes and experimental solution for estimating drag coefficients over blunt bodies.
3-D, bluff body drag estimation using a Green's function/Gram-Charlier series approach.
Barone, Matthew Franklin; De Chant, Lawrence Justin
2005-01-01
In this study, we describe the extension of the 2-d preliminary design bluff body drag estimation tool developed by De Chant1 to apply for 3-d flows. As with the 2-d method, the 3-d extension uses a combined approximate Green's function/Gram-Charlier series approach to retain the body geometry information. Whereas, the 2-d methodology relied solely upon the use of small disturbance theory for the inviscid flow field associated with the body of interest to estimate the near-field initial conditions, e.g. velocity defect, the 3-d methodology uses both analytical (where available) and numerical inviscid solutions. The defect solution is then used as an initial condition in an approximate 3-d Green's function solution. Finally, the Green's function solution is matched to the 3-d analog of the classical 2-d Gram-Charlier series and then integrated to yield the net form drag on the bluff body. Preliminary results indicate that drag estimates computed are of accuracy equivalent to the 2-d method for flows with large separation, i.e. less than 20% relative error. As was the lower dimensional method, the 3-d concept is intended to be a supplement to turbulent Navier-Stokes and experimental solution for estimating drag coefficients over blunt bodies.
Shuang, Bo; Wang, Wenxiao; Shen, Hao; Tauzin, Lawrence J.; Flatebo, Charlotte; Chen, Jianbo; Moringo, Nicholas A.; Bishop, Logan D. C.; Kelly, Kevin F.; Landes, Christy F.
2016-01-01
Super-resolution microscopy with phase masks is a promising technique for 3D imaging and tracking. Due to the complexity of the resultant point spread functions, generalized recovery algorithms are still missing. We introduce a 3D super-resolution recovery algorithm that works for a variety of phase masks generating 3D point spread functions. A fast deconvolution process generates initial guesses, which are further refined by least squares fitting. Overfitting is suppressed using a machine learning determined threshold. Preliminary results on experimental data show that our algorithm can be used to super-localize 3D adsorption events within a porous polymer film and is useful for evaluating potential phase masks. Finally, we demonstrate that parallel computation on graphics processing units can reduce the processing time required for 3D recovery. Simulations reveal that, through desktop parallelization, the ultimate limit of real-time processing is possible. Our program is the first open source recovery program for generalized 3D recovery using rotating point spread functions. PMID:27488312
Wearable 3-D Photoacoustic Tomography for Functional Brain Imaging in Behaving Rats.
Tang, Jianbo; Coleman, Jason E; Dai, Xianjin; Jiang, Huabei
2016-05-05
Understanding the relationship between brain function and behavior remains a major challenge in neuroscience. Photoacoustic tomography (PAT) is an emerging technique that allows for noninvasive in vivo brain imaging at micrometer-millisecond spatiotemporal resolution. In this article, a novel, miniaturized 3D wearable PAT (3D-wPAT) technique is described for brain imaging in behaving rats. 3D-wPAT has three layers of fully functional acoustic transducer arrays. Phantom imaging experiments revealed that the in-plane X-Y spatial resolutions were ~200 μm for each acoustic detection layer. The functional imaging capacity of 3D-wPAT was demonstrated by mapping the cerebral oxygen saturation via multi-wavelength irradiation in behaving hyperoxic rats. In addition, we demonstrated that 3D-wPAT could be used for monitoring sensory stimulus-evoked responses in behaving rats by measuring hemodynamic responses in the primary visual cortex during visual stimulation. Together, these results show the potential of 3D-wPAT for brain study in behaving rodents.
Carmona, Guillaume; Perera, Upamali; Gillett, Cheryl; Naba, Alexandra; Law, Ah-Lai; Sharma, Ved P.; Wang, Jian; Wyckoff, Jeffrey; Balsamo, Michele; Mosis, Fuad; De Piano, Mario; Monypenny, James; Woodman, Natalie; McConnell, Russell E.; Mouneimne, Ghassan; Van Hemelrijck, Mieke; Cao, Yihai; Condeelis, John; Hynes, Richard O.; Gertler, Frank B.; Krause, Matthias
2016-01-01
Cancer invasion is a hallmark of metastasis. The mesenchymal mode of cancer cell invasion is mediated by elongated membrane protrusions driven by the assembly of branched F-actin networks. How deregulation of actin regulators promotes cancer cell invasion is still enigmatic. We report that increased expression and membrane localization of the actin regulator Lamellipodin correlates with reduced metastasis-free survival and poor prognosis in breast cancer patients. In agreement we find that Lamellipodin depletion reduced lung metastasis in an orthotopic mouse breast cancer model. Invasive 3D cancer cell migration as well as invadopodia formation, and matrix degradation were impaired upon Lamellipodin depletion. Mechanistically, we show that Lamellipodin promotes invasive 3D cancer cell migration via both actin-elongating Ena/VASP proteins and the Scar/WAVE complex, which stimulates actin branching. In contrast, Lamellipodin interaction with Scar/WAVE but not Ena/VASP is required for random 2D cell migration. We identify a phosphorylation-dependent mechanism that regulates selective recruitment of these effectors to Lamellipodin: Abl-mediated Lamellipodin phosphorylation promotes its association with both Scar/WAVE and Ena/VASP, while Src-dependent phosphorylation enhances binding to Scar/WAVE but not Ena/VASP. Through these selective, regulated interactions Lamellipodin mediates directional sensing of EGF gradients and invasive 3D migration of breast cancer cells. Our findings imply that increased Lamellipodin levels enhance Ena/VASP and Scar/WAVE activities at the plasma membrane to promote 3D invasion and metastasis. PMID:26996666
Development of 3D functionally graded models by laser-assisted coaxial powder injection
NASA Astrophysics Data System (ADS)
Yakovlev, Artem; Bertrand, Ph.; Smurov, Igor Y.
2004-04-01
Relatively new method of producing 3D objects with Functionally Graded Material (FGM) structure is realized by coaxial powder injection with variable composition into the zone of laser beam action. The desired 3-dimensional material distribution is realized by repetitive deposition process. Theoretical analysis and experimental results show essential role of radiation mode and powder granularity as optimization parameters. Applied laser sources are continuous wave Nd:YAG(HAAS 2006D, 2kW), pulse-periodic Nd:YAG(HAAS HL304P, avg. power 300 W), quazi-cw CO2 (Rofin-Sinar, 300 W). Among applied materials are nanostructured WC/Co, CuSn, Stainless steel 316L, 430L, Co-base alloy, nanostructured FeCu, etc. The originality of obtained results is that different gradient types are produced "in situ" and combined within one sample: smooth, sharp or multilayered gradients. The number of samples is produced and examined with metallographical and SEM analysis. The minimal spatial gradient resolution (transition zone between two different materials) is starting from 10 microns and can be varied in a wide range; the surface roughness depends from powder granularity, best value of Ra is about 5 μm, microhardness of differet zones of samples is varied from 120 to 450 HV. The achieved geometry spatial resolution is 200 μm.
Functional Equivalence Acceptance Testing of FUN3D for Entry Descent and Landing Applications
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Wood, William A.; Kleb, William L.; Alter, Stephen J.; Glass, Christopher E.; Padilla, Jose F.; Hammond, Dana P.; White, Jeffery A.
2013-01-01
The functional equivalence of the unstructured grid code FUN3D to the the structured grid code LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. Examples from an existing suite of regression tests are used to demonstrate the functional equivalence, encompassing various thermochemical models and vehicle configurations. Algorithm modifications required for the node-based unstructured grid code (FUN3D) to reproduce functionality of the cell-centered structured code (LAURA) are also documented. Challenges associated with computation on tetrahedral grids versus computation on structured-grid derived hexahedral systems are discussed.
NASA Astrophysics Data System (ADS)
Cho, Yeunwoo
2014-11-01
The shedding phenomena of 3-D viscous gravity-capillary solitary waves generated by a moving air-forcing on the surface of deep water are investigated. Near the resonance where the forcing speed is close to 23 cm/s, two kinds of shedding modes are possible; Anti-symmetric and symmetric modes. A relevant theoretical model equation is numerically solved for the identification of shedding of solitary waves, and is analytically studied in terms of their linear stability to transverse perturbations. Furthermore, by tracing trajectories of shed solitary waves, the decay rate of a 3-D solitary wave due to viscous dissipation is estimated. This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2014R1A1A1002441).
Nonlinear evolution of 3D-inertial Alfvén wave and turbulent spectra in Auroral region
NASA Astrophysics Data System (ADS)
Rinawa, M. L.; Modi, K. V.; Sharma, R. P.
2014-10-01
In the present paper, we have investigated nonlinear interaction of three dimensional (3D) inertial Alfvén wave and perpendicularly propagating magnetosonic wave for low β-plasma ( β≪ m e / m i ). We have developed the set of dimensionless equations in the presence of ponderomotive nonlinearity due to 3D-inertial Alfvén wave in the dynamics of perpendicularly propagating magnetosonic wave. Stability analysis and numerical simulation has been carried out to study the effect of nonlinear coupling on the formation of localized structures and turbulent spectra, applicable to auroral region. The results reveal that the localized structures become more and more complex as the nonlinear interaction progresses. Further, we have studied the turbulent spectrum which follows spectral index (˜ k -3.57) at smaller scales. Relevance of the obtained results has been shown with the observations received by various spacecrafts like FAST, Hawkeye and Heos 2.
NASA Astrophysics Data System (ADS)
Sladen, A.; Monteiller, V.
2014-12-01
Most large earthquakes are generated in subduction zones. To study the complexity of these events, teleseismic body waves offer many advantages over other types of data: they allow to study both the temporal and spatial evolution of slip during the rupture, they don't depend on the presence of nearby land and they allow to study earthquakes regardless of their location. Since the development of teleseismic finite-fault inversion in the 1980th, teleseismic body waves have been simulated using 1D velocity models to take into account propagation effects at the source. Yet, subduction zones are known to be highly heterogeneous: they are characterized by curved and dipping structures, strong seismic velocity contrasts, strong variations of topography and height of the water column. The main reason for relying on a 1D approximation is the computational cost of 3D simulations. And while forward simulations of teleseismic waves in a 3D Earth are only starting to be tractable on modern computers at the frequency range of interest (0.1Hz or shorter), finite-fault source studies require a large number of these simulations. In this work, we present a new and efficient approach to compute 3D teleseismic body waves, in which the full 3D propagation is only computed in a regional domain using discontinuous Galerkin finite-element method, while the rest of the seismic wave field is propagated in a background axisymmetric Earth. The regional and global wave fields are matched using the so-called Total-Field/Scattered-Field technique. This new simulation approach allows us to study the waveform complexities resulting from 3D propagation and investigate how they could improve the resolution and reduce the non-uniqueness of finite-fault inversions.
SHEAR WAVE SEISMIC STUDY COMPARING 9C3D SV AND SH IMAGES WITH 3C3D C-WAVE IMAGES
John Beecherl; Bob A. Hardage
2004-07-01
The objective of this study was to compare the relative merits of shear-wave (S-wave) seismic data acquired with nine-component (9-C) technology and with three-component (3-C) technology. The original proposal was written as if the investigation would be restricted to a single 9-C seismic survey in southwest Kansas (the Ashland survey), on the basis of the assumption that both 9-C and 3-C S-wave images could be created from that one data set. The Ashland survey was designed as a 9-C seismic program. We found that although the acquisition geometry was adequate for 9-C data analysis, the source-receiver geometry did not allow 3-C data to be extracted on an equitable and competitive basis with 9-C data. To do a fair assessment of the relative value of 9-C and 3-C seismic S-wave data, we expanded the study beyond the Ashland survey and included multicomponent seismic data from surveys done in a variety of basins. These additional data were made available through the Bureau of Economic Geology, our research subcontractor. Bureau scientists have added theoretical analyses to this report that provide valuable insights into several key distinctions between 9-C and 3-C seismic data. These theoretical considerations about distinctions between 3-C and 9-C S-wave data are presented first, followed by a discussion of differences between processing 9-C common-midpoint data and 3-C common-conversion-point data. Examples of 9-C and 3-C data are illustrated and discussed in the last part of the report. The key findings of this study are that each S-wave mode (SH-SH, SV-SV, or PSV) involves a different subsurface illumination pattern and a different reflectivity behavior and that each mode senses a different Earth fabric along its propagation path because of the unique orientation of its particle-displacement vector. As a result of the distinct orientation of each mode's particle-displacement vector, one mode may react to a critical geologic condition in a more optimal way than do
Converted-Wave Processing of a 3D-3C Refection Seismic Survey of Soda Lake Geothermal Field
NASA Astrophysics Data System (ADS)
Louie, J. N.; Kent, T.; Echols, J.
2012-12-01
This 3D-3C seismic survey greatly improves the structural model of the Soda Lake, Nevada geothermal system. The picked top of a mudstone interval above reservoir levels reveals a detailed fault map. The geothermal reservoir is within a complex of nested grabens. Determining a "geothermal indicator" for the deeper reservoir in the seismic signal, and processing of the 3D converted-wave data, have been unsuccessful to date. Due to a high near-surface Vp/Vs ratio the shear-wave energy is under-sampled with 220 ft receiver spacing and 550 ft (168 m) line spacing. The 2D converted-wave data that we can image shows encouraging similarity to the deep structural features in the P-wave sections, but have little resolution of shallow structures. Higher-density receivers and a better shallow shear-wave model are needed in conjunction with this deep reflection study to effectively image the 3D converted waves.
NASA Astrophysics Data System (ADS)
Belashov, V. Yu.; Belashova, E. S.
2016-11-01
On the basis of the model of the three-dimensional (3D) generalized Kadomtsev-Petviashvili equation for magnetic field h = B / B the formation, stability, and dynamics of 3D soliton-like structures, such as the beams of fast magnetosonic (FMS) waves generated in ionospheric and magnetospheric plasma at a low-frequency branch of oscillations when β = 4 πnT/ B 2 ≪ 1 and β > 1, are studied. The study takes into account the highest dispersion correction determined by values of the plasma parameters and the angle θ = ( B, k), which plays a key role in the FMS beam propagation at those angles to the magnetic field that are close to π/2. The stability of multidimensional solutions is studied by an investigation of the Hamiltonian boundness under its deformations on the basis of solving of the corresponding variational problem. The evolution and dynamics of the 3D FMS wave beam are studied by the numerical integration of equations with the use of specially developed methods. The results can be interpreted in terms of the self-focusing phenomenon, as the formation of a stationary beam and the scattering and self-focusing of the solitary beam of FMS waves. These cases were studied with a detailed investigation of all evolutionary stages of the 3D FMS wave beams in the ionospheric and magnetospheric plasma.
A time-space domain stereo finite difference method for 3D scalar wave propagation
NASA Astrophysics Data System (ADS)
Chen, Yushu; Yang, Guangwen; Ma, Xiao; He, Conghui; Song, Guojie
2016-11-01
The time-space domain finite difference methods reduce numerical dispersion effectively by minimizing the error in the joint time-space domain. However, their interpolating coefficients are related with the Courant numbers, leading to significantly extra time costs for loading the coefficients consecutively according to velocity in heterogeneous models. In the present study, we develop a time-space domain stereo finite difference (TSSFD) method for 3D scalar wave equation. The method propagates both the displacements and their gradients simultaneously to keep more information of the wavefields, and minimizes the maximum phase velocity error directly using constant interpolation coefficients for different Courant numbers. We obtain the optimal constant coefficients by combining the truncated Taylor series approximation and the time-space domain optimization, and adjust the coefficients to improve the stability condition. Subsequent investigation shows that the TSSFD can suppress numerical dispersion effectively with high computational efficiency. The maximum phase velocity error of the TSSFD is just 3.09% even with only 2 sampling points per minimum wavelength when the Courant number is 0.4. Numerical experiments show that to generate wavefields with no visible numerical dispersion, the computational efficiency of the TSSFD is 576.9%, 193.5%, 699.0%, and 191.6% of those of the 4th-order and 8th-order Lax-Wendroff correction (LWC) method, the 4th-order staggered grid method (SG), and the 8th-order optimal finite difference method (OFD), respectively. Meanwhile, the TSSFD is compatible to the unsplit convolutional perfectly matched layer (CPML) boundary condition for absorbing artificial boundaries. The efficiency and capability to handle complex velocity models make it an attractive tool in imaging methods such as acoustic reverse time migration (RTM).
3D Cell Printing of Functional Skeletal Muscle Constructs Using Skeletal Muscle-Derived Bioink.
Choi, Yeong-Jin; Kim, Taek Gyoung; Jeong, Jonghyeon; Yi, Hee-Gyeong; Park, Ji Won; Hwang, Woonbong; Cho, Dong-Woo
2016-10-01
Engineered skeletal muscle tissues that mimic the structure and function of native muscle have been considered as an alternative strategy for the treatment of various muscular diseases and injuries. Here, it is demonstrated that 3D cell-printing of decellularized skeletal muscle extracellular matrix (mdECM)-based bioink facilitates the fabrication of functional skeletal muscle constructs. The cellular alignment and the shape of the tissue constructs are controlled by 3D cell-printing technology. mdECM bioink provides the 3D cell-printed muscle constructs with a myogenic environment that supports high viability and contractility as well as myotube formation, differentiation, and maturation. More interestingly, the preservation of agrin is confirmed in the mdECM, and significant increases in the formation of acetylcholine receptor clusters are exhibited in the 3D cell-printed muscle constructs. In conclusion, mdECM bioink and 3D cell-printing technology facilitate the mimicking of both the structural and functional properties of native muscle and hold great promise for producing clinically relevant engineered muscle for the treatment of muscular injuries.
Cornette, Raphaël; Baylac, Michel; Souter, Thibaud; Herrel, Anthony
2013-10-01
Morpho-functional patterns are important drivers of phenotypic diversity given their importance in a fitness-related context. Although modularity of the mandible and skull has been studied extensively in mammals, few studies have explored shape co-variation between these two structures. Despite being developmentally independent, the skull and mandible form a functionally integrated unit. In the present paper we use 3D surface geometric morphometric methods allowing us to explore the form of both skull and mandible in its 3D complexity using the greater white-toothed shrew as a model. This approach allows an accurate 3D description of zones devoid of anatomical landmarks that are functionally important. Two-block partial least-squares approaches were used to describe the co-variation of form between skull and mandible. Moreover, a 3D biomechanical model was used to explore the functional consequences of the observed patterns of co-variation. Our results show the efficiency of the method in investigations of complex morpho-functional patterns. Indeed, the description of shape co-variation between the skull and the mandible highlighted the location and the intensity of their functional relationships through the jaw adductor muscles linking these two structures. Our results also demonstrated that shape co-variation in form between the skull and mandible has direct functional consequences on the recruitment of muscles during biting.
Prediction of enzyme function based on 3D templates of evolutionarily important amino acids
Kristensen, David M; Ward, R Matthew; Lisewski, Andreas Martin; Erdin, Serkan; Chen, Brian Y; Fofanov, Viacheslav Y; Kimmel, Marek; Kavraki, Lydia E; Lichtarge, Olivier
2008-01-01
Background Structural genomics projects such as the Protein Structure Initiative (PSI) yield many new structures, but often these have no known molecular functions. One approach to recover this information is to use 3D templates – structure-function motifs that consist of a few functionally critical amino acids and may suggest functional similarity when geometrically matched to other structures. Since experimentally determined functional sites are not common enough to define 3D templates on a large scale, this work tests a computational strategy to select relevant residues for 3D templates. Results Based on evolutionary information and heuristics, an Evolutionary Trace Annotation (ETA) pipeline built templates for 98 enzymes, half taken from the PSI, and sought matches in a non-redundant structure database. On average each template matched 2.7 distinct proteins, of which 2.0 share the first three Enzyme Commission digits as the template's enzyme of origin. In many cases (61%) a single most likely function could be predicted as the annotation with the most matches, and in these cases such a plurality vote identified the correct function with 87% accuracy. ETA was also found to be complementary to sequence homology-based annotations. When matches are required to both geometrically match the 3D template and to be sequence homologs found by BLAST or PSI-BLAST, the annotation accuracy is greater than either method alone, especially in the region of lower sequence identity where homology-based annotations are least reliable. Conclusion These data suggest that knowledge of evolutionarily important residues improves functional annotation among distant enzyme homologs. Since, unlike other 3D template approaches, the ETA method bypasses the need for experimental knowledge of the catalytic mechanism, it should prove a useful, large scale, and general adjunct to combine with other methods to decipher protein function in the structural proteome. PMID:18190718
Moored Observations of Internal Waves in Luzon Strait: 3-D Structure, Dissipation, and Evolution
2013-09-30
variability, it may be due to waves propagating into Luzon strait from remote sources. Lee Waves and Dissipation on Supercritical Slopes A profiling...variability of the internal wave field in the upper 1000 m of the water column. The phase progression of internal waves as they propagate away from their
Functionalized 3D Architected Materials via Thiol-Michael Addition and Two-Photon Lithography.
Yee, Daryl W; Schulz, Michael D; Grubbs, Robert H; Greer, Julia R
2017-02-20
Fabrication of functionalized 3D architected materials is achieved by a facile method using functionalized acrylates synthesized via thiol-Michael addition, which are then polymerized using two-photon lithography. A wide variety of functional groups can be attached, from Boc-protected amines to fluoroalkanes. Modification of surface wetting properties and conjugation with fluorescent tags are demonstrated to highlight the potential applications of this technique.
NASA Astrophysics Data System (ADS)
Yuan, Huaiyu; Romanowicz, Barbara; Fischer, Karen M.; Abt, David
2011-03-01
Using a combination of long period seismic waveforms and SKS splitting measurements, we have developed a 3-D upper-mantle model (SAWum_NA2) of North America that includes isotropic shear velocity, with a lateral resolution of ˜250 km, as well as radial and azimuthal anisotropy, with a lateral resolution of ˜500 km. Combining these results, we infer several key features of lithosphere and asthenosphere structure. A rapid change from thin (˜70-80 km) lithosphere in the western United States (WUS) to thick lithosphere (˜200 km) in the central, cratonic part of the continent closely follows the Rocky Mountain Front (RMF). Changes with depth of the fast axis direction of azimuthal anisotropy reveal the presence of two layers in the cratonic lithosphere, corresponding to the fast-to-slow discontinuity found in receiver functions. Below the lithosphere, azimuthal anisotropy manifests a maximum, stronger in the WUS than under the craton, and the fast axis of anisotropy aligns with the absolute plate motion, as described in the hotspot reference frame (HS3-NUVEL 1A). In the WUS, this zone is confined between 70 and 150 km, decreasing in strength with depth from the top, from the RMF to the San Andreas Fault system and the Juan de Fuca/Gorda ridges. This result suggests that shear associated with lithosphere-asthenosphere coupling dominates mantle deformation down to this depth in the western part of the continent. The depth extent of the zone of increased azimuthal anisotropy below the cratonic lithosphere is not well resolved in our study, although it is peaked around 270 km, a robust result. Radial anisotropy is such that, predominantly, ξ > 1, where ξ= (Vsh/Vsv)2, under the continent and its borders down to ˜200 km, with stronger ξ in the bordering oceanic regions. Across the continent and below 200 km, alternating zones of weaker and stronger radial anisotropy, with predominantly ξ < 1, correlate with zones of small lateral changes in the fast axis direction of
NASA Astrophysics Data System (ADS)
Velioǧlu, Deniz; Cevdet Yalçıner, Ahmet; Zaytsev, Andrey
2016-04-01
Tsunamis are huge waves with long wave periods and wave lengths that can cause great devastation and loss of life when they strike a coast. The interest in experimental and numerical modeling of tsunami propagation and inundation increased considerably after the 2011 Great East Japan earthquake. In this study, two numerical codes, FLOW 3D and NAMI DANCE, that analyze tsunami propagation and inundation patterns are considered. Flow 3D simulates linear and nonlinear propagating surface waves as well as long waves by solving three-dimensional Navier-Stokes (3D-NS) equations. NAMI DANCE uses finite difference computational method to solve 2D depth-averaged linear and nonlinear forms of shallow water equations (NSWE) in long wave problems, specifically tsunamis. In order to validate these two codes and analyze the differences between 3D-NS and 2D depth-averaged NSWE equations, two benchmark problems are applied. One benchmark problem investigates the runup of long waves over a complex 3D beach. The experimental setup is a 1:400 scale model of Monai Valley located on the west coast of Okushiri Island, Japan. Other benchmark problem is discussed in 2015 National Tsunami Hazard Mitigation Program (NTHMP) Annual meeting in Portland, USA. It is a field dataset, recording the Japan 2011 tsunami in Hilo Harbor, Hawaii. The computed water surface elevation and velocity data are compared with the measured data. The comparisons showed that both codes are in fairly good agreement with each other and benchmark data. The differences between 3D-NS and 2D depth-averaged NSWE equations are highlighted. All results are presented with discussions and comparisons. Acknowledgements: Partial support by Japan-Turkey Joint Research Project by JICA on earthquakes and tsunamis in Marmara Region (JICA SATREPS - MarDiM Project), 603839 ASTARTE Project of EU, UDAP-C-12-14 project of AFAD Turkey, 108Y227, 113M556 and 213M534 projects of TUBITAK Turkey, RAPSODI (CONCERT_Dis-021) of CONCERT
NASA Astrophysics Data System (ADS)
Temmer, M.; Veronig, A. M.; Gopalswamy, N.; Yashiro, S.
We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of ˜ 240 km s-1. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of ˜ 750±50 km s-1, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas (Astrophys. J. 652, 763, 2006), we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.
NASA Astrophysics Data System (ADS)
Temmer, M.; Veronig, A. M.; Gopalswamy, N.; Yashiro, S.
2011-11-01
We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of ˜ 240 km s-1. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of ˜ 750±50 km s-1, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas (Astrophys. J. 652, 763, 2006), we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.
NASA Technical Reports Server (NTRS)
Temmer, M.; Veronig, A. M.; Gopalswamy, N.; Yashiro, S.
2011-01-01
We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of approx 240 km/s. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of approx 750 +/- 50 km/s, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.
NASA Astrophysics Data System (ADS)
Harvey, R. W.; Petrov, Yu. V.; Kinsey, J. E.; Liu, D.; Heidbrink, W. W.; Taylor, G.; Bonoli, P. T.
2014-10-01
Ion distribution function calculations with CQL3D have been substantially advanced through implementation of guiding-center-orbit-based Fokker-Planck Coefficients. The resulting finite-orbit-width (FOW) calculations are carried out with a fast CQL3D-Hybrid-FOW option, and in a slower but neoclassically complete (except no Er yet) CQL3D-FOW option. Good comparison between time-dependent Fast Ion Diagnostic FIDA, NPA, and neutron signals resulting from neutral beaminjection(NBI) and high harmonic fast wave (HHFW) power injected into the NSTX spherical tokamak have been simulated with the CQL3D-Hybrid-FOW, using only the FOW effects on QL diffusion, and particle losses, direct and CX. Comparisons are also made with recent CQL3D-FOW results, as well as between the original FIDA calculation code and a recent fortran version. Supported by USDOE Grants SC0006614, ER54744, and ER44649.
NASA Astrophysics Data System (ADS)
Liew, Y. M.; McLaughlin, R. A.; Chan, B. T.; Aziz, Y. F. Abdul; Chee, K. H.; Ung, N. M.; Tan, L. K.; Lai, K. W.; Ng, S.; Lim, E.
2015-04-01
Cine MRI is a clinical reference standard for the quantitative assessment of cardiac function, but reproducibility is confounded by motion artefacts. We explore the feasibility of a motion corrected 3D left ventricle (LV) quantification method, incorporating multislice image registration into the 3D model reconstruction, to improve reproducibility of 3D LV functional quantification. Multi-breath-hold short-axis and radial long-axis images were acquired from 10 patients and 10 healthy subjects. The proposed framework reduced misalignment between slices to subpixel accuracy (2.88 to 1.21 mm), and improved interstudy reproducibility for 5 important clinical functional measures, i.e. end-diastolic volume, end-systolic volume, ejection fraction, myocardial mass and 3D-sphericity index, as reflected in a reduction in the sample size required to detect statistically significant cardiac changes: a reduction of 21-66%. Our investigation on the optimum registration parameters, including both cardiac time frames and number of long-axis (LA) slices, suggested that a single time frame is adequate for motion correction whereas integrating more LA slices can improve registration and model reconstruction accuracy for improved functional quantification especially on datasets with severe motion artefacts.
Analysis of a 3-D system function measured for magnetic particle imaging.
Rahmer, Jürgen; Weizenecker, Jürgen; Gleich, Bernhard; Borgert, Jörn
2012-06-01
Magnetic particle imaging (MPI) is a new tomographic imaging approach that can quantitatively map magnetic nanoparticle distributions in vivo. It is capable of volumetric real-time imaging at particle concentrations low enough to enable clinical applications. For image reconstruction in 3-D MPI, a system function (SF) is used, which describes the relation between the acquired MPI signal and the spatial origin of the signal. The SF depends on the instrumental configuration, the applied field sequence, and the magnetic particle characteristics. Its properties reflect the quality of the spatial encoding process. This work presents a detailed analysis of a measured SF to give experimental evidence that 3-D MPI encodes information using a set of 3-D spatial patterns or basis functions that is stored in the SF. This resembles filling 3-D k-space in magnetic resonance imaging, but is faster since all information is gathered simultaneously over a broad acquisition bandwidth. A frequency domain analysis shows that the finest structures that can be encoded with the presented SF are as small as 0.6 mm. SF simulations are performed to demonstrate that larger particle cores extend the set of basis functions towards higher resolution and that the experimentally observed spatial patterns require the existence of particles with core sizes of about 30 nm in the calibration sample. A simple formula is presented that qualitatively describes the basis functions to be expected at a certain frequency.
Realistic 3D coherent transfer function inverse filtering of complex fields
Cotte, Yann; Toy, Fatih M.; Arfire, Cristian; Kou, Shan Shan; Boss, Daniel; Bergoënd, Isabelle; Depeursinge, Christian
2011-01-01
We present a novel technique for three-dimensional (3D) image processing of complex fields. It consists in inverting the coherent image formation by filtering the complex spectrum with a realistic 3D coherent transfer function (CTF) of a high-NA digital holographic microscope. By combining scattering theory and signal processing, the method is demonstrated to yield the reconstruction of a scattering object field. Experimental reconstructions in phase and amplitude are presented under non-design imaging conditions. The suggested technique is best suited for an implementation in high-resolution diffraction tomography based on sample or illumination rotation. PMID:21833359
Realistic 3D coherent transfer function inverse filtering of complex fields.
Cotte, Yann; Toy, Fatih M; Arfire, Cristian; Kou, Shan Shan; Boss, Daniel; Bergoënd, Isabelle; Depeursinge, Christian
2011-08-01
We present a novel technique for three-dimensional (3D) image processing of complex fields. It consists in inverting the coherent image formation by filtering the complex spectrum with a realistic 3D coherent transfer function (CTF) of a high-NA digital holographic microscope. By combining scattering theory and signal processing, the method is demonstrated to yield the reconstruction of a scattering object field. Experimental reconstructions in phase and amplitude are presented under non-design imaging conditions. The suggested technique is best suited for an implementation in high-resolution diffraction tomography based on sample or illumination rotation.
NASA Astrophysics Data System (ADS)
Hu, Y.; Ji, Y.; Egbert, G. D.
2015-12-01
The fictitious time domain method (FTD), based on the correspondence principle for wave and diffusion fields, has been developed and used over the past few years primarily for marine electromagnetic (EM) modeling. Here we present results of our efforts to apply the FTD approach to land and airborne TEM problems which can reduce the computer time several orders of magnitude and preserve high accuracy. In contrast to the marine case, where sources are in the conductive sea water, we must model the EM fields in the air; to allow for topography air layers must be explicitly included in the computational domain. Furthermore, because sources for most TEM applications generally must be modeled as finite loops, it is useful to solve directly for the impulse response appropriate to the problem geometry, instead of the point-source Green functions typically used for marine problems. Our approach can be summarized as follows: (1) The EM diffusion equation is transformed to a fictitious wave equation. (2) The FTD wave equation is solved with an explicit finite difference time-stepping scheme, with CPML (Convolutional PML) boundary conditions for the whole computational domain including the air and earth , with FTD domain source corresponding to the actual transmitter geometry. Resistivity of the air layers is kept as low as possible, to compromise between efficiency (longer fictitious time step) and accuracy. We have generally found a host/air resistivity contrast of 10-3 is sufficient. (3)A "Modified" Fourier Transform (MFT) allow us recover system's impulse response from the fictitious time domain to the diffusion (frequency) domain. (4) The result is multiplied by the Fourier transformation （FT） of the real source current avoiding time consuming convolutions in the time domain. (5) The inverse FT is employed to get the final full waveform and full time response of the system in the time domain. In general, this method can be used to efficiently solve most time-domain EM
3-D Modelling of Stretched Solitary Waves along Magnetic Field Lines
NASA Astrophysics Data System (ADS)
Muschietti, L.; Roth, I.; Carlson, C. W.; Berthomier, M.
2001-12-01
A model is presented for a new type of fast solitary waves which is observed by FAST in downward current regions of the auroral zone. The three-dimensional, coherent structures are electrostatic, have a positive potential, and move along the ambient magnetic field lines with speeds on the order of the electron drift. Their potential profile in the parallel direction, which can be directly measured, is flat-top whereby it cannot fit to the Gaussian shape used in previous work. Their potential profile in the perpendicular direction can only be inferred from a measured unipolar electric signal. We develop an extended BGK model which includes a flattened potential and an assumed cylindrical symmetry around a centric magnetic field line. The model envisions concentric shells of trapped electrons slowly drifting azimuthally while bouncing back and forth in the parallel direction. The electron dynamics is analysed in terms of three basic motions that occur on different time scales. These are defined by the cyclotron frequency Ω e, the bounce frequency ω b, and the azimuthal drift frequency ω γ , for which explicit analytical expressions are obtained. Subject to the ordering ω γ <<ωb<< Ωe, we calculate self-consistent distribution functions in terms of approximate constants of motion. Constraints on the parameters characterizing the amplitude and shape of the stretched solitary wave are discussed.
Carmona, G; Perera, U; Gillett, C; Naba, A; Law, A-L; Sharma, V P; Wang, J; Wyckoff, J; Balsamo, M; Mosis, F; De Piano, M; Monypenny, J; Woodman, N; McConnell, R E; Mouneimne, G; Van Hemelrijck, M; Cao, Y; Condeelis, J; Hynes, R O; Gertler, F B; Krause, M
2016-09-29
Cancer invasion is a hallmark of metastasis. The mesenchymal mode of cancer cell invasion is mediated by elongated membrane protrusions driven by the assembly of branched F-actin networks. How deregulation of actin regulators promotes cancer cell invasion is still enigmatic. We report that increased expression and membrane localization of the actin regulator Lamellipodin correlate with reduced metastasis-free survival and poor prognosis in breast cancer patients. In agreement, we find that Lamellipodin depletion reduced lung metastasis in an orthotopic mouse breast cancer model. Invasive 3D cancer cell migration as well as invadopodia formation and matrix degradation was impaired upon Lamellipodin depletion. Mechanistically, we show that Lamellipodin promotes invasive 3D cancer cell migration via both actin-elongating Ena/VASP proteins and the Scar/WAVE complex, which stimulates actin branching. In contrast, Lamellipodin interaction with Scar/WAVE but not with Ena/VASP is required for random 2D cell migration. We identified a phosphorylation-dependent mechanism that regulates selective recruitment of these effectors to Lamellipodin: Abl-mediated Lamellipodin phosphorylation promotes its association with both Scar/WAVE and Ena/VASP, whereas Src-dependent phosphorylation enhances binding to Scar/WAVE but not to Ena/VASP. Through these selective, regulated interactions Lamellipodin mediates directional sensing of epidermal growth factor (EGF) gradients and invasive 3D migration of breast cancer cells. Our findings imply that increased Lamellipodin levels enhance Ena/VASP and Scar/WAVE activities at the plasma membrane to promote 3D invasion and metastasis.
A Novel Human Tissue-Engineered 3-D Functional Vascularized Cardiac Muscle Construct
Valarmathi, Mani T.; Fuseler, John W.; Davis, Jeffrey M.; Price, Robert L.
2017-01-01
Organ tissue engineering, including cardiovascular tissues, has been an area of intense investigation. The major challenge to these approaches has been the inability to vascularize and perfuse the in vitro engineered tissue constructs. Attempts to provide oxygen and nutrients to the cells contained in the biomaterial constructs have had varying degrees of success. The aim of this current study is to develop a three-dimensional (3-D) model of vascularized cardiac tissue to examine the concurrent temporal and spatial regulation of cardiomyogenesis in the context of postnatal de novo vasculogenesis during stem cell cardiac regeneration. In order to achieve the above aim, we have developed an in vitro 3-D functional vascularized cardiac muscle construct using human induced pluripotent stem cell-derived embryonic cardiac myocytes (hiPSC-ECMs) and human mesenchymal stem cells (hMSCs). First, to generate the prevascularized scaffold, human cardiac microvascular endothelial cells (hCMVECs) and hMSCs were co-cultured onto a 3-D collagen cell carrier (CCC) for 7 days under vasculogenic culture conditions. In this milieu, hCMVECs/hMSCs underwent maturation, differentiation, and morphogenesis characteristic of microvessels, and formed extensive plexuses of vascular networks. Next, the hiPSC-ECMs and hMSCs were co-cultured onto this generated prevascularized CCCs for further 7 or 14 days in myogenic culture conditions. Finally, the vascular and cardiac phenotypic inductions were analyzed at the morphological, immunological, biochemical, molecular, and functional levels. Expression and functional analyses of the differentiated cells revealed neo-angiogenesis and neo-cardiomyogenesis. Thus, our unique 3-D co-culture system provided us the apt in vitro functional vascularized 3-D cardiac patch that can be utilized for cellular cardiomyoplasty. PMID:28194397
A Novel Human Tissue-Engineered 3-D Functional Vascularized Cardiac Muscle Construct.
Valarmathi, Mani T; Fuseler, John W; Davis, Jeffrey M; Price, Robert L
2017-01-01
Organ tissue engineering, including cardiovascular tissues, has been an area of intense investigation. The major challenge to these approaches has been the inability to vascularize and perfuse the in vitro engineered tissue constructs. Attempts to provide oxygen and nutrients to the cells contained in the biomaterial constructs have had varying degrees of success. The aim of this current study is to develop a three-dimensional (3-D) model of vascularized cardiac tissue to examine the concurrent temporal and spatial regulation of cardiomyogenesis in the context of postnatal de novo vasculogenesis during stem cell cardiac regeneration. In order to achieve the above aim, we have developed an in vitro 3-D functional vascularized cardiac muscle construct using human induced pluripotent stem cell-derived embryonic cardiac myocytes (hiPSC-ECMs) and human mesenchymal stem cells (hMSCs). First, to generate the prevascularized scaffold, human cardiac microvascular endothelial cells (hCMVECs) and hMSCs were co-cultured onto a 3-D collagen cell carrier (CCC) for 7 days under vasculogenic culture conditions. In this milieu, hCMVECs/hMSCs underwent maturation, differentiation, and morphogenesis characteristic of microvessels, and formed extensive plexuses of vascular networks. Next, the hiPSC-ECMs and hMSCs were co-cultured onto this generated prevascularized CCCs for further 7 or 14 days in myogenic culture conditions. Finally, the vascular and cardiac phenotypic inductions were analyzed at the morphological, immunological, biochemical, molecular, and functional levels. Expression and functional analyses of the differentiated cells revealed neo-angiogenesis and neo-cardiomyogenesis. Thus, our unique 3-D co-culture system provided us the apt in vitro functional vascularized 3-D cardiac patch that can be utilized for cellular cardiomyoplasty.
NASA Astrophysics Data System (ADS)
Yoon, Young Dae; Bellan, Paul M.
2016-10-01
A full three-dimensional computer code was developed in order to simulate a 3D-localized magnetic reconnection. We assume an incompressible two-fluid regime where the ions are stationary, and electron inertia and Hall effects are present. We solve a single dimensionless differential equation for perturbed magnetic fields with arbitrary background fields. The code has successfully reproduced both experimental and analytic solutions to resonance and Gendrin mode whistler waves in a uniform background field. The code was then modified to model 3D-localized magnetic reconnection as a 3D-localized perturbation on a hyperbolic-tangent background field. Three-dimensional properties that are asymmetric in the out-of-plane direction have been observed. These properties pertained to magnetic field lines, electron currents and their convection. Helicity and energy have also been examined, as well as the addition of a guide field.
Wave Phase-Sensitive Transformation of 3d-Straining of Mechanical Fields
NASA Astrophysics Data System (ADS)
Smirnov, I. N.; Speranskiy, A. A.
2015-11-01
It is the area of research of oscillatory processes in elastic mechanical systems. Technical result of innovation is creation of spectral set of multidimensional images which reflect time-correlated three-dimensional vector parameters of metrological, and\\or estimated, and\\or design parameters of oscillations in mechanical systems. Reconstructed images of different dimensionality integrated in various combinations depending on their objective function can be used as homeostatic profile or cybernetic image of oscillatory processes in mechanical systems for an objective estimation of current operational conditions in real time. The innovation can be widely used to enhance the efficiency of monitoring and research of oscillation processes in mechanical systems (objects) in construction, mechanical engineering, acoustics, etc. Concept method of vector vibrometry based on application of vector 3D phase- sensitive vibro-transducers permits unique evaluation of real stressed-strained states of power aggregates and loaded constructions and opens fundamental innovation opportunities: conduct of continuous (on-line regime) reliable monitoring of turboagregates of electrical machines, compressor installations, bases, supports, pipe-lines and other objects subjected to damaging effect of vibrations; control of operational safety of technical systems at all the stages of life cycle including design, test production, tuning, testing, operational use, repairs and resource enlargement; creation of vibro-diagnostic systems of authentic non-destructive control of anisotropic characteristics of materials resistance of power aggregates and loaded constructions under outer effects and operational flaws. The described technology is revolutionary, universal and common for all branches of engineering industry and construction building objects.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2011-12-01
averaging elastic coefficients and three averaging densities are necessary to describe the heterogeneous medium with VTI anisotropy. The resulting system is solved with iterative Krylov methods. The developed method will be incorporated in an inversion scheme for joint seismic-electromagnetic imaging. References. Brown, B.M., M. Jais, I.W. Knowles, 2005, A variational approach to an elastic inverse problem: Inverse Problems, 21, 1953-1973. Commer, M., G. Newman, 2008, New advances in three-dimensional controlled-source electromagnetic inversion: Geophysical Journal International, 172, 513-535. Newman, G. A., M. Commer and J.J. Carazzone, 2010, Imaging CSEM data in the presence of electrical anisotropy: Geophysics 75, 51-61 Petrov, P.V., G. A. Newman (2010), Using 3D Simulation of Elastic Wave Propagation in Laplace Domain for Electromagnetic-Seismic Inverse Modeling, Abstract T21A-2140 presented at 2010 Fall Meeting, AGU, San Francisco, Calif., 13-17 Dec. Shin, C. , W. Ha, 2008, A comparison between the behavior of objective functions for waveform inversion in the frequency and Laplace domains: Geophysics, 73, 119-133. Shin, C. , Y. H. Cha, 2008. Waveform inversion in the Laplace domain: Geophysical Journal International, 173, 922-931.
Creating perfused functional vascular channels using 3D bio-printing technology.
Lee, Vivian K; Kim, Diana Y; Ngo, Haygan; Lee, Young; Seo, Lan; Yoo, Seung-Schik; Vincent, Peter A; Dai, Guohao
2014-09-01
We developed a methodology using 3D bio-printing technology to create a functional in vitro vascular channel with perfused open lumen using only cells and biological matrices. The fabricated vasculature has a tight, confluent endothelium lining, presenting barrier function for both plasma protein and high-molecular weight dextran molecule. The fluidic vascular channel is capable of supporting the viability of tissue up to 5 mm in distance at 5 million cells/mL density under the physiological flow condition. In static-cultured vascular channels, active angiogenic sprouting from the vessel surface was observed whereas physiological flow strongly suppressed this process. Gene expression analysis was reported in this study to show the potential of this vessel model in vascular biology research. The methods have great potential in vascularized tissue fabrication using 3D bio-printing technology as the vascular channel is simultaneously created while cells and matrix are printed around the channel in desired 3D patterns. It can also serve as a unique experimental tool for investigating fundamental mechanisms of vascular remodeling with extracellular matrix and maturation process under 3D flow condition.
Structural and functional imaging of 3D microfluidic mixers using optical coherence tomography.
Xi, Chuanwu; Marks, Daniel L; Parikh, Devang S; Raskin, Lutgarde; Boppart, Stephen A
2004-05-18
To achieve high mixing efficiency in microfluidic devices, complex designs are often required. Microfluidic devices have been evaluated with light and confocal microscopy, but fluid-flow characteristics at different depths are difficult to separate from the en face images produced. By using optical coherence tomography (OCT), an imaging modality capable of imaging 3D microstructures at micrometer-scale resolutions over millimeter-size scales, we obtained 3D dynamic functional and structural data for three representative microfluidic mixers: a Y channel mixer, a 3D serpentine mixer, and a vortex mixer. In the serpentine mixer, OCT image analysis revealed that the mixing efficiency was linearly dependent on the Reynolds number, whereas it appeared to have exponential dependence when imaged with light microscopy. The visual overlap of fluid flows in light-microscopy images leads to an overestimation of the mixing efficiency, an effect that was eliminated with OCT imaging. Doppler OCT measurements determined velocity profiles at various points in the serpentine mixer. Mixing patterns in the vortex mixer were compared with light-microscopy and OCT image analysis. These results demonstrate that OCT can significantly improve the characterization of 3D microfluidic device structure and function.
Measuring a hidden coordinate: Rate-exchange kinetics from 3D correlation functions.
Berg, Mark A; Darvin, Jason R
2016-08-07
Nonexponential kinetics imply the existence of at least one slow variable other than the observable, that is, the system has a "hidden" coordinate. We develop a simple, but general, model that allows multidimensional correlation functions to be calculated for these systems. Homogeneous and heterogeneous mechanisms are both included, and slow exchange of the rates is allowed. This model shows that 2D and 3D correlation functions of the observable measure the distribution and kinetics of the hidden coordinate controlling the rate exchange. Both the mean exchange time and the shape of the exchange relaxation are measurable. However, complications arise because higher correlation functions are sums of multiple "pathways," each of which measures different dynamics. Only one 3D pathway involves exchange dynamics. Care must be used to extract exchange dynamics without contamination from other processes.
A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies
2014-01-01
We show that an Ng bridge function modified version of the three-dimensional reference interaction site model (3D-RISM-NgB) solvation free energy method can accurately predict the hydration free energy (HFE) of a set of 504 organic molecules. To achieve this, a single unique constant parameter was adjusted to the computed HFE of single atom Lennard-Jones solutes. It is shown that 3D-RISM is relatively accurate at predicting the electrostatic component of the HFE without correction but requires a modification of the nonpolar contribution that originates in the formation of the cavity created by the solute in water. We use a free energy functional with the Ng scaling of the direct correlation function [Ng, K. C. J. Chem. Phys.1974, 61, 2680]. This produces a rapid, reliable small molecule HFE calculation for applications in drug design. PMID:24634616
Tian, Jie; Liu, Qianqi; Wang, Xi; Xing, Ping; Yang, Zhuowen; Wu, Changjun
2017-01-01
As breast cancer tissues are stiffer than normal tissues, shear wave elastography (SWE) can locally quantify tissue stiffness and provide histological information. Moreover, tissue stiffness can be observed on three-dimensional (3D) colour-coded elasticity maps. Our objective was to evaluate the diagnostic performances of quantitative features in differentiating breast masses by two-dimensional (2D) and 3D SWE. Two hundred ten consecutive women with 210 breast masses were examined with B-mode ultrasound (US) and SWE. Quantitative features of 3D and 2D SWE were assessed, including elastic modulus standard deviation (ESDE) measured on SWE mode images and ESDU measured on B-mode images, as well as maximum elasticity (Emax). Adding quantitative features to B-mode US improved the diagnostic performance (p < 0.05) and reduced false-positive biopsies (p < 0.0001). The area under the receiver operating characteristic curve (AUC) of 3D SWE was similar to that of 2D SWE for ESDE (p = 0.026) and ESDU (p = 0.159) but inferior to that of 2D SWE for Emax (p = 0.002). Compared with ESDU, ESDE showed a higher AUC on 2D (p = 0.0038) and 3D SWE (p = 0.0057). Our study indicates that quantitative features of 3D and 2D SWE can significantly improve the diagnostic performance of B-mode US, especially 3D SWE ESDE, which shows considerable clinical value. PMID:28106134
NASA Astrophysics Data System (ADS)
Obermann, Anne; Planès, Thomas; Hadziioannou, Céline; Campillo, Michel
2016-10-01
In the context of seismic monitoring, recent studies made successful use of seismic coda waves to locate medium changes on the horizontal plane. Locating the depth of the changes, however, remains a challenge. In this paper, we use 3-D wavefield simulations to address two problems: first, we evaluate the contribution of surface- and body-wave sensitivity to a change at depth. We introduce a thin layer with a perturbed velocity at different depths and measure the apparent relative velocity changes due to this layer at different times in the coda and for different degrees of heterogeneity of the model. We show that the depth sensitivity can be modelled as a linear combination of body- and surface-wave sensitivity. The lapse-time-dependent sensitivity ratio of body waves and surface waves can be used to build 3-D sensitivity kernels for imaging purposes. Second, we compare the lapse-time behaviour in the presence of a perturbation in horizontal and vertical slabs to address, for instance, the origin of the velocity changes detected after large earthquakes.
3D finite element modelling of guided wave scattering at delaminations in composites
NASA Astrophysics Data System (ADS)
Murat, Bibi Intan Suraya; Fromme, Paul
2016-02-01
Carbon fiber laminate composites are increasingly used for aerospace structures as they offer a number of advantages including a good strength to weight ratio. However, impact during the operation and servicing of the aircraft can lead to barely visible and difficult to detect damage. Depending on the severity of the impact, delaminations can occur, reducing the load carrying capacity of the structure. Efficient nondestructive testing of composite panels can be achieved using guided ultrasonic waves propagating along the structure. The guided wave (A0 Lamb wave mode) scattering at delaminations was modeled using full three-dimensional Finite Element (FE) simulations. The influence of the delamination size was systematically investigated from a parameter study. A significant influence of the delamination width on the guided wave scattering was found, especially on the angular dependency of the scattered guided wave amplitude. The sensitivity of guided ultrasonic waves for the detection of delamination damage in composite panels is discussed.
Estimating Fiber Orientation Distribution Functions in 3D-Polarized Light Imaging
Axer, Markus; Strohmer, Sven; Gräßel, David; Bücker, Oliver; Dohmen, Melanie; Reckfort, Julia; Zilles, Karl; Amunts, Katrin
2016-01-01
Research of the human brain connectome requires multiscale approaches derived from independent imaging methods ideally applied to the same object. Hence, comprehensible strategies for data integration across modalities and across scales are essential. We have successfully established a concept to bridge the spatial scales from microscopic fiber orientation measurements based on 3D-Polarized Light Imaging (3D-PLI) to meso- or macroscopic dimensions. By creating orientation distribution functions (pliODFs) from high-resolution vector data via series expansion with spherical harmonics utilizing high performance computing and supercomputing technologies, data fusion with Diffusion Magnetic Resonance Imaging has become feasible, even for a large-scale dataset such as the human brain. Validation of our approach was done effectively by means of two types of datasets that were transferred from fiber orientation maps into pliODFs: simulated 3D-PLI data showing artificial, but clearly defined fiber patterns and real 3D-PLI data derived from sections through the human brain and the brain of a hooded seal. PMID:27147981
3C3D VSP Imaging of Salt Flanks Using Converted Waves in the Gulf of Mexico
NASA Astrophysics Data System (ADS)
Li, Y.; Doherty, F.; Jackson, J.
2005-05-01
Locating salt boundary and imaging updip sediment structures flanking the salt domes are very important tasks for exploration in the Gulf of Mexico since major petroleum reserves are often trapped underneath overhangs of diapiric salt domes. Although the top of salt and less steep structures can be well imaged using current surface seismic methods, the steep sides of a salt dome with irregularly shapes are hard to image with adequate accuracy. Thus, Vertical Seismic Profiling (VSP) surveys with three-component (3C) receivers in wells are usually requested for improving images of subsurface structures. Conventional multi-offset VSP (OVSP) and refraction salt proximity (SP) surveys are widely applied in the Gulf of Mexico to improve images of slat interfaces, sub-salt and salt flank structures using P waves. In this paper, we will focus on using converted waves to image the steep salt-sediment boundary. A VSP dataset, including multi-OVSP and a SP survey, acquired in the Gulf of Mexico was used in this study. We analyzed 3C OVSP data to identify and separate converted waves, such as PS, P-SP, P-SS, generated at a salt boundary. Then both PP wave and converted waves were 3C3D depth migrated to generate images of the steep salt-sediment interface. Both transmitted P-P and P-S converted waves from the SP survey were used to calculate 3D salt exit points which delineate the steep salt face. The VSP results derived from both methods are abundant and a suitable 3D visualization tool is required for visual integration and interpretation. The image volumes and other available geophysical and geological data were integrated using a 3D visualization tool specially designed for VSP solutions. The migrated images using PP and converted waves provides a precise and complete definition of the steep salt face and reservoir sands flanking the salt dome. This study indicates that both reflection and reflection surveys can result in a consistent location of the steep salt flank
Shaping functional nano-objects by 3D confined supramolecular assembly.
Deng, Renhua; Liang, Fuxin; Li, Weikun; Liu, Shanqin; Liang, Ruijing; Cai, Mingle; Yang, Zhenzhong; Zhu, Jintao
2013-12-20
Nano-objects are generated through 3D confined supramolecular assembly, followed by a sequential disintegration by rupturing the hydrogen bonding. The shape of the nano-objects is tunable, ranging from nano-disc, nano-cup, to nano-toroid. The nano-objects are pH-responsive. Functional materials for example inorganic or metal nanoparticles are easily complexed onto the external surface, to extend both composition and microstructure of the nano-objects.
Adaptive multiconfigurational wave functions
Evangelista, Francesco A.
2014-03-28
A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N{sub 2} and the potential energy curves for the first three singlet states of C{sub 2}. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu{sub 2}O{sub 2}{sup 2+} core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.
Engineering of functional, perfusable 3D microvascular networks on a chip.
Kim, Sudong; Lee, Hyunjae; Chung, Minhwan; Jeon, Noo Li
2013-04-21
Generating perfusable 3D microvessels in vitro is an important goal for tissue engineering, as well as for reliable modelling of blood vessel function. To date, in vitro blood vessel models have not been able to accurately reproduce the dynamics and responses of endothelial cells to grow perfusable and functional 3D vascular networks. Here we describe a microfluidic-based platform whereby we model natural cellular programs found during normal development and angiogenesis to form perfusable networks of intact 3D microvessels as well as tumor vasculatures based on the spatially controlled co-culture of endothelial cells with stromal fibroblasts, pericytes or cancer cells. The microvessels possess the characteristic morphological and biochemical markers of in vivo blood vessels, and exhibit strong barrier function and long-term stability. An open, unobstructed microvasculature allows the delivery of nutrients, chemical compounds, biomolecules and cell suspensions, as well as flow-induced mechanical stimuli into the luminal space of the endothelium, and exhibits faithful responses to physiological shear stress as demonstrated by cytoskeleton rearrangement and increased nitric oxide synthesis. This simple and versatile platform provides a wide range of applications in vascular physiology studies as well as in developing vascularized organ-on-a-chip and human disease models for pharmaceutical screening.
NASA Astrophysics Data System (ADS)
Staple, Bevan; Earhart, R. P.; Slaymaker, Philip A.; Drouillard, Thomas F., II; Mahony, Thomas
2005-05-01
3D imaging LADARs have emerged as the key technology for producing high-resolution imagery of targets in 3-dimensions (X and Y spatial, and Z in the range/depth dimension). Ball Aerospace & Technologies Corp. continues to make significant investments in this technology to enable critical NASA, Department of Defense, and national security missions. As a consequence of rapid technology developments, two issues have emerged that need resolution. First, the terminology used to rate LADAR performance (e.g., range resolution) is inconsistently defined, is improperly used, and thus has become misleading. Second, the terminology does not include a metric of the system"s ability to resolve the 3D depth features of targets. These two issues create confusion when translating customer requirements into hardware. This paper presents a candidate framework for addressing these issues. To address the consistency issue, the framework utilizes only those terminologies proposed and tested by leading LADAR research and standards institutions. We also provide suggestions for strengthening these definitions by linking them to the well-known Rayleigh criterion extended into the range dimension. To address the inadequate 3D image quality metrics, the framework introduces the concept of a Range/Depth Modulation Transfer Function (RMTF). The RMTF measures the impact of the spatial frequencies of a 3D target on its measured modulation in range/depth. It is determined using a new, Range-Based, Slanted Knife-Edge test. We present simulated results for two LADAR pulse detection techniques and compare them to a baseline centroid technique. Consistency in terminology plus a 3D image quality metric enable improved system standardization.
Analysis of non linear partially standing waves from 3D velocity measurements
NASA Astrophysics Data System (ADS)
Drevard, D.; Rey, V.; Svendsen, Ib; Fraunie, P.
2003-04-01
Surface gravity waves in the ocean exhibit an energy spectrum distributed in both frequency and direction of propagation. Wave data collection is of great importance in coastal zones for engineering and scientific studies. In particular, partially standing waves measurements near coastal structures and steep or barred beaches may be a requirement, for instance for morphodynamic studies. The aim of the present study is the analysis of partially standing surface waves icluding non-linear effects. According to 1st order Stokes theory, synchronous measurements of horizontal and vertical velocity components allow calculation of rate of standing waves (Drevard et al, 2003). In the present study, it is demonstrated that for deep water conditions, partially standing 2nd order Stokes waves induced velocity field is still represented by the 1st order solution for the velocity potential contrary to the surface elevation which exhibits harmonic components. For intermediate water depth, harmonic components appear not only in the surface elevation but also in the velocity fields, but their weight remains much smaller, because of the vertical decreasing wave induced motion. For irregular waves, the influence of the spectrum width on the non-linear effects in the analysis is discussed. Keywords: Wave measurements ; reflection ; non-linear effects Acknowledgements: This work was initiated during the stay of Prof. Ib Svendsen, as invited Professor, at LSEET in autumn 2002. This study is carried out in the framework of the Scientific French National Programmes PNEC ART7 and PATOM. Their financial supports are acknowledged References: Drevard, D., Meuret, A., Rey, V. Piazzola, J. And Dolle, A.. (2002). "Partially reflected waves measurements using Acoustic Doppler Velocimeter (ADV)", Submitted to ISOPE 03, Honolulu, Hawaii, May 2003.
WaveQ3D: Fast and accurate acoustic transmission loss (TL) eigenrays, in littoral environments
NASA Astrophysics Data System (ADS)
Reilly, Sean M.
This study defines a new 3D Gaussian ray bundling acoustic transmission loss model in geodetic coordinates: latitude, longitude, and altitude. This approach is designed to lower the computation burden of computing accurate environmental effects in sonar training application by eliminating the need to transform the ocean environment into a collection of Nx2D Cartesian radials. This approach also improves model accuracy by incorporating real world 3D effects, like horizontal refraction, into the model. This study starts with derivations for a 3D variant of Gaussian ray bundles in this coordinate system. To verify the accuracy of this approach, acoustic propagation predictions of transmission loss, time of arrival, and propagation direction are compared to analytic solutions and other models. To validate the model's ability to predict real world phenomena, predictions of transmission loss and propagation direction are compared to at-sea measurements, in an environment where strong horizontal refraction effect have been observed. This model has been integrated into U.S. Navy active sonar training system applications, where testing has demonstrated its ability to improve transmission loss calculation speed without sacrificing accuracy.
NASA Astrophysics Data System (ADS)
Xu, Ting; You, Xue-yi
2017-04-01
A 3D sediment transport model based on the modified environmental fluid dynamics code (EFDC) and the nearshore waves simulation model (SWAN) is developed to study the change of suspended sediment concentration and bottom shear stress under the actions of pure current and wave-current. After being validated by the field measured data, the proposed sediment transport model is applied in the Oujiang River Estuary, China. The results show that the ratios of both bottom shear stress and suspended sediment concentration of pure current to those of wave-current show a gradually increase from shallow nearshore water to deep open sea. The results also show that the proportion of wave contributions on bottom shear stress and sediment concentration are above 60%, approximately 20-30% and less than 10% for the water depth of less than 5 m, 5-10 m and more than 20 m, respectively. For the waters among islands, the proportion of wave contribution to bottom shear stress and sediment concentration is reduced to 10-20% for -5 m water depth and this is more obvious for the waves of large amplitude. The bottom stress and suspended sediment concentration between islands are mainly controlled by tidal current, and the effect of wave is not significant.
Investigation of Parametric Excitation of Whistler Waves Using 3D Particle-In-Cell Simulations
NASA Astrophysics Data System (ADS)
Caplinger, James; Sotnikov, Vladimir; Main, Daniel; Rose, David; Paraschiv, Ioana
2016-10-01
Previous theoretical work has shown that a parametric interaction between quasi-electrostatic lower oblique resonance (LOR) and lower frequency (ω < ωLH) ion acoustic or extremely low frequency (ELF) waves can produce electromagnetic whistler waves in a cold magnetized plasma. It was also demonstrated theoretically that this interaction can more efficiently generate electromagnetic whistler waves than by direct excitation by a conventional loop antenna, operating at a single frequency. For the purpose of numerically validating the above result, a series of particle-in-cell simulations were carried out. We first demonstrate the ability to accurately model whistler wave excitation producing the familiar resonant surfaces which comprise the LOR using a modeled loop antenna. Next we demonstrate the ability to generate ion acoustic waves as well as ELF waves, both of which are shown to agree with the expected linear dispersion relations. Finally, we investigate the existence of any nonlinear interaction which indicates the desired parametric excitation and attempt to analyze the efficiency of this method of excitation and radiated power going into the whistler part of the VLF wave spectrum.
Functional classification of protein 3D structures from predicted local interaction sites.
Parasuram, Ramya; Lee, Joslynn S; Yin, Pengcheng; Somarowthu, Srinivas; Ondrechen, Mary Jo
2010-12-01
A new approach to the functional classification of protein 3D structures is described with application to some examples from structural genomics. This approach is based on functional site prediction with THEMATICS and POOL. THEMATICS employs calculated electrostatic potentials of the query structure. POOL is a machine learning method that utilizes THEMATICS features and has been shown to predict accurate, precise, highly localized interaction sites. Extension to the functional classification of structural genomics proteins is now described. Predicted functionally important residues are structurally aligned with those of proteins with previously characterized biochemical functions. A 3D structure match at the predicted local functional site then serves as a more reliable predictor of biochemical function than an overall structure match. Annotation is confirmed for a structural genomics protein with the ribulose phosphate binding barrel (RPBB) fold. A putative glucoamylase from Bacteroides fragilis (PDB ID 3eu8) is shown to be in fact probably not a glucoamylase. Finally a structural genomics protein from Streptomyces coelicolor annotated as an enoyl-CoA hydratase (PDB ID 3g64) is shown to be misannotated. Its predicted active site does not match the well-characterized enoyl-CoA hydratases of similar structure but rather bears closer resemblance to those of a dehalogenase with similar fold.
Chinnakkannu Vijayakumar, Sivaprasad; Venkatakrishnan, Krishnan; Tan, Bo
2017-02-15
We introduce a "non-noble metal" based SERS active nanobiosensor using a self-assembled 3D hybrid nickel nanonetwork. A tunable biomolecule detector fabricated by a bottom-up approach was functionalized using a multiphoton ionization energy mechanism to create a self-assembled 3D hybrid nickel nanonetwork. The nanonetwork was tested for SERS detection of crystal violet (CV) and glutathione (GSH) at two excitation wavelengths, 532 and 785 nm. The results reveal indiscernible peaks with a limit of detection (LOD) of 1 picomolar (pM) concentration. An enhancement factor (EF) of 9.3 × 10(8) was achieved for the chemical molecule CV and 1.8 × 10(9) for the biomolecule GSH, which are the highest reported values so far. The two results, one being the CV molecule proved that nickel nanonetwork is indeed SERS active and the second being the GSH biomolecule detection at both 532 and 785 nm, confirm that the nanonetwork is a biosensor which has potential for both in vivo and in vitro sensing. In addition, the selectivity and versatility of this biosensor is examined with biomolecules such as l-Cysteine, l-Methionine, and sensing GSH in cell culture medium which mimics the complex biological environment. The functionalized self-assembled 3D hybrid nickel nanonetwork exhibits electromagnetic and charge transfer based SERS activation mechanisms.
Berney, Sandra; Bétrancourt, Mireille; Molinari, Gaëlle; Hoyek, Nady
2015-01-01
The emergence of dynamic visualizations of three-dimensional (3D) models in anatomy curricula may be an adequate solution for spatial difficulties encountered with traditional static learning, as they provide direct visualization of change throughout the viewpoints. However, little research has explored the interplay between learning material presentation formats, spatial abilities, and anatomical tasks. First, to understand the cognitive challenges a novice learner would be faced with when first exposed to 3D anatomical content, a six-step cognitive task analysis was developed. Following this, an experimental study was conducted to explore how presentation formats (dynamic vs. static visualizations) support learning of functional anatomy, and affect subsequent anatomical tasks derived from the cognitive task analysis. A second aim was to investigate the interplay between spatial abilities (spatial visualization and spatial relation) and presentation formats when the functional anatomy of a 3D scapula and the associated shoulder flexion movement are learned. Findings showed no main effect of the presentation formats on performances, but revealed the predictive influence of spatial visualization and spatial relation abilities on performance. However, an interesting interaction between presentation formats and spatial relation ability for a specific anatomical task was found. This result highlighted the influence of presentation formats when spatial abilities are involved as well as the differentiated influence of spatial abilities on anatomical tasks.
Bundled-optode implementation for 3D imaging in functional near-infrared spectroscopy
Nguyen, Hoang-Dung; Hong, Keum-Shik
2016-01-01
The paper presents a functional near-infrared spectroscopy (fNIRS)-based bundled-optode method for detection of the changes of oxy-hemoglobin (HbO) and deoxy-hemoglobin (HbR) concentrations. fNIRS with 32 optodes is utilized to measure five healthy male subjects’ brain-hemodynamic responses to arithmetic tasks. Specifically, the coordinates of 256 voxels in the three-dimensional (3D) volume are computed according to the known probe geometry. The mean path length factor in the Beer-Lambert equation is estimated as a function of the emitter-detector distance, which is utilized for computation of the absorption coefficient. The mean values of HbO and HbR obtained from the absorption coefficient are then applied for construction of a 3D fNIRS image. Our results show that the proposed method, as compared with the conventional approach, can detect brain activity with higher spatial resolution. This method can be extended for 3D fNIRS imaging in real-time applications. PMID:27699115
NASA Astrophysics Data System (ADS)
Tavangar, Amirhossein; Tan, Bo; Venkatakrishnan, Krishnan
2013-08-01
This study proposes a novel technique to synthesize functional carbonaceous three-dimensional (3-D) micro/nanocompounds from agricultural by-products using femtosecond laser irradiation. Biowastes of rice husk and wheat straw are value-engineered to carbonaceous structures in a single-step process under ambient conditions. Our results demonstrate that by controlling the laser fluence, structures with a variety of different morphologies from nanostructures to microstructures can be achieved. Also, the results indicate that altering the laser processing parameters influences the chemical composition of the synthesized structures. This sustainable approach presents an important step towards synthesizing 3-D micro/nanofibrous compounds from biowaste materials. These structures, as-synthesized or as nanocomposite fillers, can have practical uses in electronic, sensing, biological, and environmental applications.
Localization of 3D inertial Alfvén wave and generation of turbulence
NASA Astrophysics Data System (ADS)
Sharma, R. P.; Sharma, Prachi; Yadav, N.
2015-06-01
The present paper deals with the nonlinear interaction of Inertial Alfvén wave (IAW) and fast magnetosonic wave in the low beta plasma, where beta is the ratio of thermal pressure to the background magnetic pressure. In this paper, the localization and turbulent spectra of IAW along with the density dips correlated with the fast magnetosonic wave have been investigated. Variation of parallel electric field along and across the field lines has also been studied. Taking ponderomotive nonlinear effect in the dynamics of fast magnetosonic wave, couple of dimensionless equations has been derived. These coupled equations have been simulated numerically using the pseudo-spectral method. The obtained results reveal that the Kolmogorov scaling is followed by a steeper scaling in magnetic power spectrum, which is consistent with the observations by the FAST and Hawkeye spacecraft in auroral region. The relevance of present investigation has been discussed for auroral plasmas.
Testing hypotheses of bat baculum function with 3D models derived from microCT
Herdina, Anna Nele; Kelly, Diane A; Jahelková, Helena; Lina, Peter H C; Horáček, Ivan; Metscher, Brian D
2015-01-01
The baculum (os penis) has been extensively studied as a taxon-specific character in bats and other mammals but its mechanical function is still unclear. There is a wide consensus in the literature that the baculum is probably a sexually selected character. Using a novel approach combining postmortem manipulation and three-dimensional (3D) imaging, we tested two functional hypotheses in the common noctule bat Nyctalus noctula, the common pipistrelle Pipistrellus pipistrellus, and Nathusius’ pipistrelle Pipistrellus nathusii: (i) whether the baculum can protect the distal urethra and urethral opening from compression during erection and copulation; and (ii) whether the baculum and corpora cavernosa form a functional unit to support both the penile shaft and the more distal glans tip. In freshly dead or frozen and thawed bats, we compared flaccid penises with artificially ‘erect’ penises that were inflated with 10% formalin. Penises were stained with alcoholic iodine and imaged with a lab-based high-resolution x-ray microtomography system. Analysis of the 3D images enabled us to compare the changes in relative positions of the baculum, corpora cavernosa, urethra, and corpus spongiosum with one another between flaccid and ‘erect’ penises. Our results support both functional hypotheses, indicating that the baculum probably performs two different roles during erection. Our approach should prove valuable for comparing and testing the functions of different baculum morphologies in bats and other mammals. Moreover, we have validated an essential component of the groundwork necessary to extend this approach with finite element analysis for quantitative 3D biomechanical modeling of penis function. PMID:25655647
Chang, Mei-Chu; Tsai, Ya-Yi; I, Lin
2013-08-15
We experimentally demonstrate the direct observation of defect mediated wave turbulence with fluctuating defects and low amplitude hole filaments, from a 3D self-excited plane dust acoustic wave in a dusty plasma by reducing dissipation. The waveform undulation is found to be the origin for the amplitude and the phase modulations of the local dust density oscillation, the broadening of the sharp peaks in the frequency spectrum, and the fluctuating defects. The corrugated wave crest surface also causes the observed high and low density patches in the transverse (xy) plane. Low oscillation amplitude spots (holes) share the same positions with the defects. Their trajectories in the xyt space appear in the form of chaotic filaments without long term predictability, through uncertain pair generation, propagation, and pair annihilation.
Moored Observations of Internal Waves in Luzon Strait: 3-D Structure, Dissipation, and Evolution
2016-03-01
advancing the performance of operational and climate models, as well as for understanding local problems such as pollutant dispersal and biological...Y.J. Yang, M.-H. Chang , and Q. Li. 2011. From Luzon Strait to Dongsha Plateau: Stages in the life of an internal wave. Oceanography 24(4):64–77...Knowledge of the general problems of internal waves and ocean mixing are important for advancing the performance of operational and climate models, as well
Propagation of 3D nonlinear waves over complex bathymetry using a High-Order Spectral method
NASA Astrophysics Data System (ADS)
Gouin, Maïté; Ducrozet, Guillaume; Ferrant, Pierre
2016-04-01
Scattering of regular and irregular surface gravity waves propagating over a region of arbitrary three-dimensional varying bathymetry is considered here. The three-dimensional High-Order Spectral method (HOS) with an extension to account for a variable bathymetry is used. The efficiency of the model has been proved to be conserved even with this extension. The method is first applied to a bathymetry consisting of an elliptical lens, as used in the Vincent and Briggs (1989) experiment. Incident waves passing across the lens are transformed and a strong convergence region is observed after the elliptical mound. The wave amplification depends on the incident wave. Numerical results for regular and irregular waves are analysed and compared with other methods and experimental data demonstrating the efficiency and practical applicability of the present approach. Then the method is used to model waves propagating over a real bathymetry: the canyons of Scripps/La Jolla in California. The implementation of this complex bathymetry in the model is presented, as well as the first results achieved. They will be compared to the ones obtained with another numerical model.
Conducting a 3D Converted Shear Wave Project to Reduce Exploration Risk at Wister, CA
Matlick, Skip; Walsh, Patrick; Rhodes, Greg; Fercho, Steven
2015-06-30
Ormat sited 2 full-size exploration wells based on 3D seismic interpretation of fractures, prior drilling results, and temperature anomaly. The wells indicated commercial temperatures (>300 F), but almost no permeability, despite one of the wells being drilled within 820 ft of an older exploration well with reported indications of permeability. Following completion of the second well in 2012, Ormat undertook a lengthy program to 1) evaluate the lack of observed permeability, 2) estimate the likelihood of finding permeability with additional drilling, and 3) estimate resource size based on an anticipated extent of permeability.
Display depth analyses with the wave aberration for the auto-stereoscopic 3D display
NASA Astrophysics Data System (ADS)
Gao, Xin; Sang, Xinzhu; Yu, Xunbo; Chen, Duo; Chen, Zhidong; Zhang, Wanlu; Yan, Binbin; Yuan, Jinhui; Wang, Kuiru; Yu, Chongxiu; Dou, Wenhua; Xiao, Liquan
2016-07-01
Because the aberration severely affects the display performances of the auto-stereoscopic 3D display, the diffraction theory is used to analyze the diffraction field distribution and the display depth through aberration analysis. Based on the proposed method, the display depth of central and marginal reconstructed images is discussed. The experimental results agree with the theoretical analyses. Increasing the viewing distance or decreasing the lens aperture can improve the display depth. Different viewing distances and the LCD with two lens-arrays are used to verify the conclusion.
2.5D real waveform and real noise simulation of receiver functions in 3D models
NASA Astrophysics Data System (ADS)
Schiffer, Christian; Jacobsen, Bo; Balling, Niels
2014-05-01
, the comparison can even be made quantitative, and an iterative inverse 3D model updating would be possible. Furthermore the "2.5D" modelling approach of the in reality 3D problem must be investigated in terms of accuracy of the approximation, in particular with focus on highly 3D structures and multiple phases. References: Schiffer et al. 2013, A fossil subduction zone in the East Greenland Caledonides revealed by a Receiver Function analysis, Geophysical Research Abstracts, Vol. 15, EGU2013-6947, 2013, EGU General Assembly 2013 Schlindwein, V. and Jokat, W., 1999, Structure and evolution of the continental crust of northern east Greenland from integrated geophysical studies: Journal of Geophysical Research, v. 104, no. B7, p. 15227-15,245 Voss, M. and Jokat, W., 2007, Continent-ocean transition and voluminous magmatic underplating derived from P-wave velocity modelling of the East Greenland continental margin: Geophysical Journal International, v.170, no. 2, p. 580-604
Matthew Ward, R.; Venner, Eric; Daines, Bryce; Murray, Stephen; Erdin, Serkan; Kristensen, David M.; Lichtarge, Olivier
2009-01-01
Summary:The Evolutionary Trace Annotation (ETA) Server predicts enzymatic activity. ETA starts with a structure of unknown function, such as those from structural genomics, and with no prior knowledge of its mechanism uses the phylogenetic Evolutionary Trace (ET) method to extract key functional residues and propose a function-associated 3D motif, called a 3D template. ETA then searches previously annotated structures for geometric template matches that suggest molecular and thus functional mimicry. In order to maximize the predictive value of these matches, ETA next applies distinctive specificity filters—evolutionary similarity, function plurality and match reciprocity. In large scale controls on enzymes, prediction coverage is 43% but the positive predictive value rises to 92%, thus minimizing false annotations. Users may modify any search parameter, including the template. ETA thus expands the ET suite for protein structure annotation, and can contribute to the annotation efforts of metaservers. Availability:The ETA Server is a web application available at http://mammoth.bcm.tmc.edu/eta/. Contact: lichtarge@bcm.edu PMID:19307237
SAFE-3D analysis of a piezoelectric transducer to excite guided waves in a rail web
NASA Astrophysics Data System (ADS)
Ramatlo, Dineo A.; Long, Craig S.; Loveday, Philip W.; Wilke, Daniel N.
2016-02-01
Our existing Ultrasonic Broken Rail Detection system detects complete breaks and primarily uses a propagating mode with energy concentrated in the head of the rail. Previous experimental studies have demonstrated that a mode with energy concentrated in the head of the rail, is capable of detecting weld reflections at long distances. Exploiting a mode with energy concentrated in the web of the rail would allow us to effectively detect defects in the web of the rail and could also help to distinguish between reflections from welds and cracks. In this paper, we will demonstrate the analysis of a piezoelectric transducer attached to the rail web. The forced response at different frequencies is computed by the Semi-Analytical Finite Element (SAFE) method and compared to a full three-dimensional finite element method using ABAQUS. The SAFE method only requires the rail track cross-section to be meshed using two-dimensional elements. The ABAQUS model in turn requires a full three-dimensional discretisation of the rail track. The SAFE approach can yield poor predictions at cut-on frequencies associated with other modes in the rail. Problematic frequencies are identified and a suitable frequency range identified for transducer design. The forced response results of the two methods were found to be in good agreement with each other. We then use a previously developed SAFE-3D method to analyse a practical transducer over the selected frequency range. The results obtained from the SAFE-3D method are in good agreement with experimental measurements.
Real-time 3D millimeter wave imaging based FMCW using GGD focal plane array as detectors
NASA Astrophysics Data System (ADS)
Levanon, Assaf; Rozban, Daniel; Kopeika, Natan S.; Yitzhaky, Yitzhak; Abramovich, Amir
2014-03-01
Millimeter wave (MMW) imaging systems are required for applications in medicine, communications, homeland security, and space technology. This is because there is no known ionization hazard for biological tissue, and atmospheric attenuation in this range of the spectrum is relatively low. The lack of inexpensive room temperature imaging systems makes it difficult to give a suitable MMW system for many of the above applications. 3D MMW imaging system based on chirp radar was studied previously using a scanning imaging system of a single detector. The system presented here proposes to employ a chirp radar method with a Glow Discharge Detector (GDD) Focal Plane Array (FPA) of plasma based detectors. Each point on the object corresponds to a point in the image and includes the distance information. This will enable 3D MMW imaging. The radar system requires that the millimeter wave detector (GDD) will be able to operate as a heterodyne detector. Since the source of radiation is a frequency modulated continuous wave (FMCW), the detected signal as a result of heterodyne detection gives the object's depth information according to value of difference frequency, in addition to the reflectance of the image. In this work we experimentally demonstrate the feasibility of implementing an imaging system based on radar principles and FPA of GDD devices. This imaging system is shown to be capable of imaging objects from distances of at least 10 meters.
NASA Astrophysics Data System (ADS)
Bisrat, Shishay T.; DeShon, Heather R.; Pesicek, Jeremy; Thurber, Clifford
2014-01-01
A three-dimensional (3-D), high-resolution P wave seismic attenuation model for the New Madrid Seismic Zone (NMSZ) is determined using P wave path attenuation (t*) values of small-magnitude earthquakes (MD < 3.9). Events were recorded at 89 broadband and short-period seismometers of the Cooperative New Madrid Seismic Zone Network and 40 short-period seismometers of the Portable Array for Numerical Data Acquisition experiment. The amplitude spectra of all the earthquakes are simultaneously inverted for source, path (t*), and site parameters. The t* values are inverted for QP using local earthquake tomography methods and a known 3-D P wave velocity model for the region. The four major seismicity arms of the NMSZ exhibit reduced QP (higher attenuation) than the surrounding crust. The highest attenuation anomalies coincide with areas of previously reported high swarm activity attributed to fluid-rich fractures along the southeast extension of the Reelfoot fault. The QP results are consistent with previous attenuation studies in the region, which showed that active fault zones and fractured crust in the NMSZ are highly attenuating.
NASA Astrophysics Data System (ADS)
Ayers, J.; Owens, C. T.; Liu, K. C.; Swenson, E.; Ghoshal, A.; Weiss, V.
2013-01-01
The application of guided waves to interrogate remote areas of structural components has been researched extensively in characterizing damage. However, there exists a sparsity of work in using piezoelectric transducer-generated guided waves as a method of assessing stress intensity factors (SIF). This quantitative information enables accurate estimation of the remaining life of metallic structures exhibiting cracks, such as military and commercial transport vehicles. The proposed full wavefield approach, based on 3D laser vibrometry and piezoelectric transducer-generated guided waves, provides a practical means for estimation of dynamic stress intensity factors (DSIF) through local strain energy mapping via the J-integral. Strain energies and traction vectors can be conveniently estimated from wavefield data recorded using 3D laser vibrometry, through interpolation and subsequent spatial differentiation of the response field. Upon estimation of the Jintegral, it is possible to obtain the corresponding DSIF terms. For this study, the experimental test matrix consists of aluminum plates with manufactured defects representing canonical elliptical crack geometries under uniaxial tension that are excited by surface mounted piezoelectric actuators. The defects' major to minor axes ratios vary from unity to approximately 133. Finite element simulations are compared to experimental results and the relative magnitudes of the J-integrals are examined.
Schilling, Kurt; Janve, Vaibhav; Gao, Yurui; Stepniewska, Iwona; Landman, Bennett A; Anderson, Adam W
2016-01-01
The ability of diffusion MRI (dMRI) fiber tractography to non-invasively map three-dimensional (3D) anatomical networks in the human brain has made it a valuable tool in both clinical and research settings. However, there are many assumptions inherent to any tractography algorithm that can limit the accuracy of the reconstructed fiber tracts. Among them is the assumption that the diffusion-weighted images accurately reflect the underlying fiber orientation distribution (FOD) in the MRI voxel. Consequently, validating dMRI’s ability to assess the underlying fiber orientation in each voxel is critical for its use as a biomedical tool. Here, using post-mortem histology and confocal microscopy, we present a method to perform histological validation of orientation functions in 3D, which has previously been limited to two-dimensional analysis of tissue sections. We demonstrate the ability to extract the 3D FOD from confocal z-stacks, and quantify the agreement between the MRI estimates of orientation information obtained using constrained spherical deconvolution (CSD) and the true geometry of the fibers. We find an orientation error of approximately 6° in voxels containing nearly parallel fibers, and 10-11° in crossing fiber regions, and note that CSD was unable to resolve fibers crossing at angles below 60° in our dataset. This is the first time the 3D white matter orientation distribution is calculated from histology and compared to dMRI. Thus, this technique serves as a gold standard for dMRI validation studies - providing the ability to determine the extent to which the dMRI signal is consistent with the histological FOD, and to establish how well different dMRI models can predict the ground truth FOD. PMID:26804781
NASA Astrophysics Data System (ADS)
Lee, Dong Hun; Lee, Jung Mo; Cho, Hyun-Moo; Kang, Tae-Seob
2016-10-01
A temporary seismic array was in operation between October 2010 and March 2013 in the Gyeongju area of Korea. Teleseismic records of the seismic array appropriate for receiver function analysis were collected, and selected seismograms were split into five groups based on epicenters-the Banda-Molucca, Sumatra, Iran, Aleutian, and Vanuatu groups. 1D velocity structures beneath each seismic station were estimated by inverting the stacked receiver functions for possible groups. The inversion was done by applying a genetic algorithm, whereas surface wave dispersion data were used as constraints to avoid non-uniqueness in the inversion. The composite velocity structure was constructed by averaging the velocity structures weighted by the number of receiver functions used in stacking. The uncertainty analysis for the velocity structures showed that the average of 95% confidence intervals was ± 0.1 km/s. The 3D velocity structure was modeled through interpolation of 1D composite velocity structures. Moho depths were determined in each composite velocity structure based on the AK135-F S-wave velocity model, and the depths were similar to the H-κ analysis results. The deepest Moho depth in the study area was found to be 31.9 km, and the shallowest, was 25.9 km. The Moho discontinuity dips in a southwestward direction beneath the area. A low velocity layer was also detected between 4 and 14 km depth. Adakitic intrusions and/or a high geothermal gradient appear to be the causes of this low velocity layer. The 3D velocity structure can be used to reliably assess seismic hazards in this area.
3D finite element model of the chinchilla ear for characterizing middle ear functions.
Wang, Xuelin; Gan, Rong Z
2016-10-01
Chinchilla is a commonly used animal model for research of sound transmission through the ear. Experimental measurements of the middle ear transfer function in chinchillas have shown that the middle ear cavity greatly affects the tympanic membrane (TM) and stapes footplate (FP) displacements. However, there is no finite element (FE) model of the chinchilla ear available in the literature to characterize the middle ear functions with the anatomical features of the chinchilla ear. This paper reports a recently completed 3D FE model of the chinchilla ear based on X-ray micro-computed tomography images of a chinchilla bulla. The model consisted of the ear canal, TM, middle ear ossicles and suspensory ligaments, and the middle ear cavity. Two boundary conditions of the middle ear cavity wall were simulated in the model as the rigid structure and the partially flexible surface, and the acoustic-mechanical coupled analysis was conducted with these two conditions to characterize the middle ear function. The model results were compared with experimental measurements reported in the literature including the TM and FP displacements and the middle ear input admittance in chinchilla ear. An application of this model was presented to identify the acoustic role of the middle ear septa-a unique feature of chinchilla middle ear cavity. This study provides the first 3D FE model of the chinchilla ear for characterizing the middle ear functions through the acoustic-mechanical coupled FE analysis.
A 3D MPI-Parallel GPU-accelerated framework for simulating ocean wave energy converters
NASA Astrophysics Data System (ADS)
Pathak, Ashish; Raessi, Mehdi
2015-11-01
We present an MPI-parallel GPU-accelerated computational framework for studying the interaction between ocean waves and wave energy converters (WECs). The computational framework captures the viscous effects, nonlinear fluid-structure interaction (FSI), and breaking of waves around the structure, which cannot be captured in many potential flow solvers commonly used for WEC simulations. The full Navier-Stokes equations are solved using the two-step projection method, which is accelerated by porting the pressure Poisson equation to GPUs. The FSI is captured using the numerically stable fictitious domain method. A novel three-phase interface reconstruction algorithm is used to resolve three phases in a VOF-PLIC context. A consistent mass and momentum transport approach enables simulations at high density ratios. The accuracy of the overall framework is demonstrated via an array of test cases. Numerical simulations of the interaction between ocean waves and WECs are presented. Funding from the National Science Foundation CBET-1236462 grant is gratefully acknowledged.
Summary of work on shock wave feature extraction in 3-D datasets
NASA Technical Reports Server (NTRS)
Hesselink, Lambertus (Principal Investigator)
1996-01-01
A method for extracting and visualizing shock waves from three dimensional data-sets is discussed. Issues concerning computation time, robustness to numerical perturbations, and noise introduction are considered and compared with other methods. Finally, results using this method are discussed.
Luquet, David; Marchiano, Régis; Coulouvrat, François
2015-10-28
Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D
NASA Astrophysics Data System (ADS)
Luquet, David; Marchiano, Régis; Coulouvrat, François
2015-10-01
Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D
NASA Astrophysics Data System (ADS)
Lin, C. W.; Wu, T. R.; Chuang, M. H.; Tsai, Y. L.
2015-12-01
The wind in Taiwan Strait is strong and stable which offers an opportunity to build offshore wind farms. However, frequently visited typhoons and strong ocean current require more attentions on the wave force and local scour around the foundation of the turbine piles. In this paper, we introduce an in-house, multi-phase CFD model, Splash3D, for solving the flow field with breaking wave, strong turbulent, and scour phenomena. Splash3D solves Navier-Stokes Equation with Large-Eddy Simulation (LES) for the fluid domain, and uses volume of fluid (VOF) with piecewise linear interface reconstruction (PLIC) method to describe the break free-surface. The waves were generated inside the computational domain by internal wave maker with a mass-source function. This function is designed to adequately simulate the wave condition under observed extreme events based on JONSWAP spectrum and dispersion relationship. Dirichlet velocity boundary condition is assigned at the upper stream boundary to induce the ocean current. At the downstream face, the sponge-layer method combined with pressure Dirichlet boundary condition is specified for dissipating waves and conducting current out of the domain. Numerical pressure gauges are uniformly set on the structure surface to obtain the force distribution on the structure. As for the local scour around the foundation, we developed Discontinuous Bi-viscous Model (DBM) for the development of the scour hole. Model validations were presented as well. The force distribution under observed irregular wave condition was extracted by the irregular-surface force extraction (ISFE) method, which provides a fast and elegant way to integrate the force acting on the surface of irregular structure. From the Simulation results, we found that the total force is mainly induced by the impinging waves, and the force from the ocean current is about 2 order of magnitude smaller than the wave force. We also found the dynamic pressure, wave height, and the
Laminar optical tomography: high-resolution 3D functional imaging of superficial tissues
NASA Astrophysics Data System (ADS)
Hillman, Elizabeth M. C.; Devor, Anna; Dunn, Andrew K.; Boas, David A.
2006-03-01
Laminar Optical Tomography (LOT) is a new medical imaging modality for high-resolution, depth-resolved, functional imaging of superficial tissue such as rodent cortex, skin and the retina. LOT uses visible laser light to image to depths of >2mm (far deeper than microscopy) and is highly sensitive to absorption and fluorescence contrast, enabling spectroscopic functional information such as hemoglobin oxygenation to be imaged with 100-200 micron resolution. LOT has been used to image the hemodynamic response to stimulus in the somatosensory cortex of rats. The resulting three-dimensional (3D) images through the depth of the cortex can be used to delineate the arterial, capillary and venous responses, revealing new information about the intricacies of the oxygenation and blood flow dynamics related to neuronal activation. Additional applications of LOT are being explored, including the integration of 3D Voltage Sensitive Dye fluorescence imaging. LOT imaging uses a system similar to a confocal microscope, quickly scanning a focused beam of light over the surface of the tissue (~8Hz frame rate). Light is detected from both the focus of the scanning beam, and also at increasing distances from the beam's focus. This scattered light has penetrated more deeply into the tissue, and allows features at different depths to be distinguished. An algorithm that includes photon migration modeling of light scattering converts the raw data into 3D images. The motivation for functional optical imaging will be outlined, the basic principles of LOT imaging will be described, and the latest in-vivo results will be presented.
Spherical cavity-expansion forcing function in PRONTO 3D for application to penetration problems
Warren, T.L.; Tabbara, M.R.
1997-05-01
In certain penetration events the primary mode of deformation of the target can be approximated by known analytical expressions. In the context of an analysis code, this approximation eliminates the need for modeling the target as well as the need for a contact algorithm. This technique substantially reduces execution time. In this spirit, a forcing function which is derived from a spherical-cavity expansion analysis has been implemented in PRONTO 3D. This implementation is capable of computing the structural and component responses of a projectile due to three dimensional penetration events. Sample problems demonstrate good agreement with experimental and analytical results.
Low temperature assembly of functional 3D DNA-PNA-protein complexes.
Flory, Justin D; Simmons, Chad R; Lin, Su; Johnson, Trey; Andreoni, Alessio; Zook, James; Ghirlanda, Giovanna; Liu, Yan; Yan, Hao; Fromme, Petra
2014-06-11
Proteins have evolved to carry out nearly all the work required of living organisms within complex inter- and intracellular environments. However, systematically investigating the range of interactions experienced by a protein that influence its function remains challenging. DNA nanostructures are emerging as a convenient method to arrange a broad range of guest molecules. However, flexible methods are needed for arranging proteins in more biologically relevant 3D geometries under mild conditions that preserve protein function. Here we demonstrate how peptide nucleic acid (PNA) can be used to control the assembly of cytochrome c (12.5 kDa, pI 10.5) and azurin (13.9 kDa, pI 5.7) proteins into separate 3D DNA nanocages, in a process that maintains protein function. Toehold-mediated DNA strand displacement is introduced as a method to purify PNA-protein conjugates. The PNA-proteins were assembled within 2 min at room temperature and within 4 min at 11 °C, and hybridize with even greater efficiency than PNA conjugated to a short peptide. Gel electrophoresis and steady state and time-resolved fluorescence spectroscopy were used to investigate the effect of protein surface charge on its interaction with the negatively charged DNA nanocage. These data were used to generate a model of the DNA-PNA-protein complexes that show the negatively charged azurin protein repelled away from the DNA nanocage while the positively charged cytochrome c protein remains within and closely interacts with the DNA nanocage. When conjugated to PNA and incorporated into the DNA nanocage, the cytochrome c secondary structure and catalytic activity were maintained, and its redox potential was reduced modestly by 20 mV possibly due to neutralization of some positive surface charges. This work demonstrates a flexible new approach for using 3D nucleic acid (PNA-DNA) nanostructures to control the assembly of functional proteins, and facilitates further investigation of protein interactions as well
Zhang, Yimeng; Li, Xiong; Samonds, Jason M.
2015-01-01
Bayesian theory has provided a compelling conceptualization for perceptual inference in the brain. Central to Bayesian inference is the notion of statistical priors. To understand the neural mechanisms of Bayesian inference, we need to understand the neural representation of statistical regularities in the natural environment. In this paper, we investigated empirically how statistical regularities in natural 3D scenes are represented in the functional connectivity of disparity-tuned neurons in the primary visual cortex of primates. We applied a Boltzmann machine model to learn from 3D natural scenes, and found that the units in the model exhibited cooperative and competitive interactions, forming a “disparity association field”, analogous to the contour association field. The cooperative and competitive interactions in the disparity association field are consistent with constraints of computational models for stereo matching. In addition, we simulated neurophysiological experiments on the model, and found the results to be consistent with neurophysiological data in terms of the functional connectivity measurements between disparity-tuned neurons in the macaque primary visual cortex. These findings demonstrate that there is a relationship between the functional connectivity observed in the visual cortex and the statistics of natural scenes. They also suggest that the Boltzmann machine can be a viable model for conceptualizing computations in the visual cortex and, as such, can be used to predict neural circuits in the visual cortex from natural scene statistics. PMID:26712581
Fluorescence fluctuation microscopy to reveal 3D architecture and function in the cell nucleus.
Lenser, Thorsten; Weisshart, Klaus; Ulbricht, Tobias; Klement, Karolin; Hemmerich, Peter
2010-01-01
The three-dimensional (3D) architecture of the cell nucleus is determined not only by the presence of subnuclear domains, such as the nuclear envelope, chromosome territories, and nuclear bodies, but also by smaller domains which form in response to specific functions, such as RNA transcription, DNA replication, and DNA repair. Since both stable and dynamic structures contribute to nuclear morphology, it is important to study the biophysical principles of the formation of macromolecular assemblies within the nucleus. For this purpose, a variety of fluorescence fluctuation microscopy techniques can be applied. Here, we summarize our current knowledge on the 3D architecture of the mammalian cell nucleus and describe in detail how the assembly of functional nuclear protein complexes can be analyzed in living cells using fluorescence bleaching techniques, fluorescence correlation spectroscopy, raster image correlation spectroscopy, and mathematical modeling. In conclusion, the application of all these techniques in combination is a powerful tool to assess the full spectrum of nuclear protein dynamics and to understand the biophysical principles underlying nuclear structure and function.
Zhang, Yimeng; Li, Xiong; Samonds, Jason M; Lee, Tai Sing
2016-03-01
Bayesian theory has provided a compelling conceptualization for perceptual inference in the brain. Central to Bayesian inference is the notion of statistical priors. To understand the neural mechanisms of Bayesian inference, we need to understand the neural representation of statistical regularities in the natural environment. In this paper, we investigated empirically how statistical regularities in natural 3D scenes are represented in the functional connectivity of disparity-tuned neurons in the primary visual cortex of primates. We applied a Boltzmann machine model to learn from 3D natural scenes, and found that the units in the model exhibited cooperative and competitive interactions, forming a "disparity association field", analogous to the contour association field. The cooperative and competitive interactions in the disparity association field are consistent with constraints of computational models for stereo matching. In addition, we simulated neurophysiological experiments on the model, and found the results to be consistent with neurophysiological data in terms of the functional connectivity measurements between disparity-tuned neurons in the macaque primary visual cortex. These findings demonstrate that there is a relationship between the functional connectivity observed in the visual cortex and the statistics of natural scenes. They also suggest that the Boltzmann machine can be a viable model for conceptualizing computations in the visual cortex and, as such, can be used to predict neural circuits in the visual cortex from natural scene statistics.
NASA Astrophysics Data System (ADS)
Garth, T.; Rietbrock, A.
2011-12-01
Dispersion of body wave arrivals observed in the fore-arc have been attributed to high frequency guided waves being retained and delayed by a low velocity layer (LVL) in the subducted crust. Lower frequency seismic waves travel at higher velocities in the surrounding mantle. These subduction zone guided waves have the potential to offer unique insights into subducting oceanic crust. Two and three dimensional finite difference (FD) wave propagation models are used to investigate the factors controlling guided wave dispersion and to test which features of the subducted crust can be resolved by guided waves. Other factors that may affect the frequency content of arrivals in the fore-arc such as elevated attenuation are also investigated. Modeling results are compared to observed guided wave dispersion in the Japan, Aleutian and Central American subduction zones. Modeling has shown that trade-offs occur between the velocity contrast and the thickness of the waveguide, with both parameters potentially affecting the frequency content that is delayed. We combine amplitude spectra plots with displacement spectrograms so that the relative amplitudes and relative arrival times of different frequencies can be compared. This allows the specific effects of given parameters to be understood. The effect of elevated attenuation on the frequency content of arrivals in the fore-arc is investigated using a visco-elastic FD wave propagation model (Bohlen 2002). The sensitivity of observed dispersion to variations in the Vp/Vs ratio of the waveguide material is also investigated. Understanding the relative dispersion of P and S waves as well as the relative importance of attenuation in the subduction system may allow us to understand more about the hydrous conditions in subduction zones. Systematic variations in the contrast between the LVL and the surrounding material are investigated. Modeling is designed to test if guided wave dispersion can resolve down dip velocity changes in the
Moored Observations of Internal Waves in Luzon Strait: 3-D Structure, Dissipation, and Evolution
2014-09-30
the performance of operational and climate models, as well as for understanding local problems such as pollutant dispersal and biological productivity...substantially improves both our understanding and predictive ability of linear internal tides and NLIWs in Luzon Strait and the South China Sea...westward into the northeastern South China Sea (SCS). • To better understand generation and propagation of internal waves in a strongly sheared
Full 3D dispersion curve solutions for guided waves in generally anisotropic media
NASA Astrophysics Data System (ADS)
Hernando Quintanilla, F.; Lowe, M. J. S.; Craster, R. V.
2016-02-01
Dispersion curves of guided waves provide valuable information about the physical and elastic properties of waves propagating within a given waveguide structure. Algorithms to accurately compute these curves are an essential tool for engineers working in non-destructive evaluation and for scientists studying wave phenomena. Dispersion curves are typically computed for low or zero attenuation and presented in two or three dimensional plots. The former do not always provide a clear and complete picture of the dispersion loci and the latter are very difficult to obtain when high values of attenuation are involved and arbitrary anisotropy is considered in single or multi-layered systems. As a consequence, drawing correct and reliable conclusions is a challenging task in the modern applications that often utilize multi-layered anisotropic viscoelastic materials. These challenges are overcome here by using a spectral collocation method (SCM) to robustly find dispersion curves in the most complicated cases of high attenuation and arbitrary anisotropy. Solutions are then plotted in three-dimensional frequency-complex wavenumber space, thus gaining much deeper insight into the nature of these problems. The cases studied range from classical examples, which validate this approach, to new ones involving materials up to the most general triclinic class for both flat and cylindrical geometry in multi-layered systems. The apparent crossing of modes within the same symmetry family in viscoelastic media is also explained and clarified by the results. Finally, the consequences of the centre of symmetry, present in every crystal class, on the solutions are discussed.
Well-posedness of linearized motion for 3-D water waves far from equilibrium
Hou, T.Y.; Zhen-huan Teng; Pingwen Zhang
1996-12-31
In this paper, we study the motion of a free surface separating two different layers of fluid in three dimensions. We assume the flow to be inviscid, irrotational, and incompressible. In this case, one can reduce the entire motion by variables on the surface alone. In general, without additional regularizing effects such as surface alone. In general, without additional regularizing effects such as surface tension or viscosity, the flow can be subject to Rayleigh-Taylor or Kelvin-Helmholtz instabilities which will lead to unbounded growth in high frequency wave numbers. In this case, the problem is not well-posed in the Hadamard sense. The problem of water wave with no fluid above is a special case. It is well-known that such motion is well-posed when the free surface is sufficiently close to equilibrium. Beale, Hous and Lowengrub derived a general condition which ensures well-posedness of the linearization about a presumed time-dependent motion in two dimensional case. The linearized equations, when formulated in a proper coordinate system are found to have a qualitative structure surprisingly like that for the simple case of linear waves near equilbrium. Such an analysis is essential in analyzing stability of boundary integral methods for computing free interface problems. 19 refs.
NASA Astrophysics Data System (ADS)
Benhamouche, Mehdi; Bernard, Laurent; Serhir, Mohammed; Pichon, Lionel; Lesselier, Dominique
2013-11-01
This paper proposes a criterion for locating obstacles by time reversal (TR) of electromagnetic (EM) waves based on the analysis of the density of EM energy map in time domain. Contrarily to a monochromatic study of the TR, the wide-band approach requires to determine the instant of the wave focus. This enables us to locate the focal spots that are indicative of the positions. The criterion proposed is compared to the inverse of the minimum entropy criterion as used in the literature [X. Xu, E.L. Miller, C.M. Rappaport, IEEE Trans. Geosci. Remote Sens. 41, 1804 (2003)]. An application for the localization of 3D metal targets is proposed using finite integration technique (FIT) as computational tool at the modeling stage. An experimental validation is presented for canonical three-dimensional configurations with two kinds of metal objects. Contribution to the Topical Issue "Numelec 2012", Edited by Adel Razek.
Exploring the surface reactivity of 3d metal endofullerenes: a density-functional theory study.
Estrada-Salas, Rubén E; Valladares, Ariel A
2009-09-24
Changes in the preferential sites of electrophilic, nucleophilic, and radical attacks on the pristine C60 surface with endohedral doping using 3d transition metal atoms were studied via two useful reactivity indices, namely the Fukui functions and the molecular electrostatic potential. Both of these were calculated at the density functional BPW91 level of theory with the DNP basis set. Our results clearly show changes in the preferential reactivity sites on the fullerene surface when it is doped with Mn, Fe, Co, or Ni atoms, whereas there are no significant changes in the preferential reactivity sites on the C60 surface upon endohedral doping with Cu and Zn atoms. Electron affinities (EA), ionization potentials (IP), and HOMO-LUMO gaps (Eg) were also calculated to complete the study of the endofullerene's surface reactivity. These findings provide insight into endofullerene functionalization, an important issue in their application.
NASA Astrophysics Data System (ADS)
Zhang, Jianzhong; Huang, Yueqin; Song, Lin-Ping; Liu, Qing-Huo
2011-03-01
We propose a new ray tracing technique in a 3-D heterogeneous isotropic media based on bilinear traveltime interpolation and the wave front group marching. In this technique, the media is discretized into a series of rectangular cells. There are two steps to be carried out: one is a forward step where wave front expansion is evolved from sources to whole computational domain and the subsequent one is a backward step where ray paths are calculated for any source-receiver configuration as desired. In the forward step, we derive a closed-form expression to calculate traveltime at an arbitrary point in a cell using a bilinear interpolation of the known traveltimes on the cell's surface. Then the group marching method (GMM), a fast wave front advancing method, is applied to expand the wave front from the source to all girds. In the backward step, ray paths starting from receivers are traced by finding the intersection points of potential ray propagation vectors with the surfaces of relevant cells. In this step, the same TI scheme is used to compute the candidate intersection points on all surfaces of each relevant cell. In this process, the point with the minimum traveltime is selected as a ray point from which the similar step is continued until sources. A number of numerical experiments demonstrate that our 3-D ray tracing technique is able to achieve very accurate computation of traveltimes and ray paths and meanwhile take much less computer time in comparison with the existing popular ones like the finite-difference-based GMM method, which is combined with the maximum gradient ray tracing, and the shortest path method.
3-D upper mantle shear wave speed structure beneath the South Pacific Superswell by a BBOBS array
NASA Astrophysics Data System (ADS)
Isse, T.; Suetsugu, D.; Shiobara, H.; Sugioka, H.; Yoshizawa, K.; Kanazawa, T.; Fukao, Y.
2005-12-01
Previous seismic tomography studies show a broad low velocity anomaly in the lower mantle, so-called superplume, beneath the South Pacific and there are hotspot chains and large scale topographic high at surface of this region. However, the resolution of seismic tomography is poor, especially in the upper mantle, because of limited spatial distribution of seismic stations. To improve the station coverage, we deployed an array of long-term broadband ocean bottom seismometers (BBOBS) in this region. The quality of the vertical component of seismograms recorded by the BBOBS array is comparable with those by island seismic stations. This observation has enabled us to obtain a more precise 3-D shear wave speed structure in the upper mantle of this region by analyzing Rayleigh waves. We employed a two-station method to determine phase velocity of fundamental mode Rayleigh wave recorded by the BBOBS array and island stations in the Pacific Ocean. We obtained 1025 path-average phase velocity dispersion curves including 188 dispersion curves using the BBOBS data in a period range between 40 and 140 seconds. We then inverted them to a 3-D shear wave speed structure down to a depth of 200 km. At shallow depths the eastern part of the French Polynesia region is in general slower than the western part, which indicates an age-dependence of seismic structure of the uppermost mantle. Slow speed anomalies corresponding to the hotspots are apparently superposed on this age-dependence: Slow speed anomalies can be seen from the surface to a depth of 200 km beneath the Society, Pitcairn, and Macdonald hotspots, but they are limited only to the deep part beneath the Samoa hotspot. The slow speed anomalies beneath the Pitcairn and Society hotspots apparently coalesce at a depth of 100 km, where a single anomaly extending upward from below seems to branch into two directions. A resolution analysis indicates that the BBOBS array data has improved the spatial resolution substantially.
Imaging of 3D Ocean Turbulence Microstructure Using Low Frequency Acoustic Waves
NASA Astrophysics Data System (ADS)
Minakov, Alexander; Kolyukhin, Dmitriy; Keers, Henk
2015-04-01
In the past decade the technique of imaging the ocean structure with low-frequency signal (Hz), produced by air-guns and typically employed during conventional multichannel seismic data acquisition, has emerged. The method is based on extracting and stacking the acoustic energy back-scattered by the ocean temperature and salinity micro- and meso-structure (1 - 100 meters). However, a good understanding of the link between the scattered wavefield utilized by the seismic oceanography and physical processes in the ocean is still lacking. We describe theory and the numerical implementation of a 3D time-dependent stochastic model of ocean turbulence. The velocity and temperature are simulated as homogeneous Gaussian isotropic random fields with the Kolmogorov-Obukhov energy spectrum in the inertial subrange. Numerical modeling technique is employed for sampling of realizations of random fields with a given spatial-temporal spectral tensor. The model used is shown to be representative for a wide range of scales. Using this model, we provide a framework to solve the forward and inverse acoustic scattering problem using marine seismic data. Our full-waveform inversion method is based on the ray-Born approximation which is specifically suitable for the modelling of small velocity perturbations in the ocean. This is illustrated by showing a good match between synthetic seismograms computed using ray-Born and synthetic seismograms produced with a more computationally expensive finite-difference method.
NASA Astrophysics Data System (ADS)
Matsuda, Takashi S.; Nakamura, Takuji; Murphy, Damian; Tsutsumi, Masaki; Moffat-Griffin, Tracy; Zhao, Yucheng; Pautet, Pierre-Dominique; Ejiri, Mitsumu K.; Taylor, Michael
2016-07-01
ANGWIN (Antarctic Gravity Wave Imaging/Instrument Network) is an international airglow imager/instrument network in the Antarctic, which commenced observations in 2011. It seeks to reveal characteristics of mesospheric gravity waves, and to study sources, propagation, breaking of the gravity waves over the Antarctic and the effects on general circulation and upper atmosphere. In this study, we compared distributions of horizontal phase velocity of the gravity waves at around 90 km altitude observed in the mesospheric airglow imaging over different locations using our new statistical analysis method of 3-D Fourier transform, developed by Matsuda et al. (2014). Results from the airglow imagers at four stations at Syowa (69S, 40E), Halley (76S, 27W), Davis (69S, 78E) and McMurdo (78S, 156E) out of the ANGWIN imagers have been compared, for the observation period between April 6 and May 21 in 2013. In addition to the horizontal distribution of propagation and phase speed, gravity wave energies have been quantitatively compared, indicating a smaller GW activity in higher latitude stations. We further investigated frequency dependence of gravity wave propagation direction, as well as nightly variation of the gravity wave direction and correlation with the background wind variations. We found that variation of propagation direction is partly due to the effect of background wind in the middle atmosphere, but variation of wave sources could play important role as well. Secondary wave generation is also needed to explain the observed results.
3D functional ultrasound imaging of the cerebral visual system in rodents.
Gesnik, Marc; Blaize, Kevin; Deffieux, Thomas; Gennisson, Jean-Luc; Sahel, José-Alain; Fink, Mathias; Picaud, Serge; Tanter, Mickaël
2017-02-03
3D functional imaging of the whole brain activity during visual task is a challenging task in rodents due to the complex tri-dimensional shape of involved brain regions and the fine spatial and temporal resolutions required to reveal the visual tract. By coupling functional ultrasound (fUS) imaging with a translational motorized stage and an episodic visual stimulation device, we managed to accurately map and to recover the activity of the visual cortices, the Superior Colliculus (SC) and the Lateral Geniculate Nuclei (LGN) in 3D. Cerebral Blood Volume (CBV) responses during visual stimuli were found to be highly correlated with the visual stimulus time profile in visual cortices (r=0.6), SC (r=0.7) and LGN (r=0.7). These responses were found dependent on flickering frequency and contrast, and optimal stimulus parameters for largest CBV increases were obtained. In particular, increasing the flickering frequency higher than 7Hz revealed a decrease of visual cortices response while the SC response was preserved. Finally, cross-correlation between CBV signals exhibited significant delays (d=0.35s +/-0.1s) between blood volume response in SC and visual cortices in response to our visual stimulus. These results emphasize the interest of fUS imaging as a whole brain neuroimaging modality for brain vision studies in rodent models.
NASA Astrophysics Data System (ADS)
Belashov, Vasily
We study the formation, structure, stability and dynamics of the multidimensional soliton-like beam structures forming on the low-frequency branch of oscillation in the ionospheric and magnetospheric plasma for cases when beta=4pinT/B(2) <<1 and beta>1. In first case with the conditions omega
Mapping 3-D functional capillary geometry in rat skeletal muscle in vivo
Milkovich, Stephanie; Goldman, Daniel; Ellis, Christopher G.
2012-01-01
We have developed a novel mapping software package to reconstruct microvascular networks in three dimensions (3-D) from in vivo video images for use in blood flow and O2 transport modeling. An intravital optical imaging system was used to collect video sequences of blood flow in microvessels at different depths in the tissue. Functional images of vessels were produced from the video sequences and were processed using automated edge tracking software to yield location and geometry data for construction of the 3-D network. The same video sequences were analyzed for hemodynamic and O2 saturation data from individual capillaries in the network. Simple user-driven commands allowed the connection of vessel segments at bifurcations, and semiautomated registration enabled the tracking of vessels across multiple focal planes and fields of view. The reconstructed networks can be rotated and manipulated in 3-D to verify vessel connections and continuity. Hemodynamic and O2 saturation measurements made in vivo can be indexed to corresponding vessels and visualized using colorized maps of the vascular geometry. Vessels in each reconstruction are saved as text-based files that can be easily imported into flow or O2 transport models with complete geometry, hemodynamic, and O2 transport conditions. The results of digital morphometric analysis of seven microvascular networks showed mean capillary diameters and overall capillary density consistent with previous findings using histology and corrosion cast techniques. The described mapping software is a valuable tool for the quantification of in vivo microvascular geometry, hemodynamics, and oxygenation, thus providing rich data sets for experiment-based computational models. PMID:22140042
3-D joint inversion of the magnetotelluric phase tensor and vertical magnetic transfer functions
NASA Astrophysics Data System (ADS)
Tietze, Kristina; Ritter, Oliver; Egbert, Gary D.
2015-11-01
With advancing computational resources, 3-D inversion techniques have become feasible in recent years and are now a more widely used tool for magnetotelluric (MT) data interpretation. Galvanic distortion caused by small-scale near-surface inhomogeneities remains an obstacle for 3-D MT inversion which so far has experienced little attention. If not considered properly, the effect on 3-D inversion can be immense and result in erroneous subsurface models and interpretations. To tackle the problem we implemented inversion of the distortion-free phase tensor into the ModEM inversion package. The dimensionless phase tensor components describe only variations of the conductivity structure. When inverting these data, particular care has to be taken of the conductivity structure in the a priori model, which provides the reference frame when transferring the information from phase tensors into absolute conductivity values. Our results obtained with synthetic data show that phase tensor inversion can recover the regional conductivity structure in presence of galvanic distortion if the a priori model provides a reasonable assumption for the regional resistivity average. Joint inversion of phase tensor data and vertical magnetic transfer functions improves recovery of the absolute resistivity structure and is less dependent on the prior model. We also used phase tensor inversion for a data set of more than 250 MT sites from the central San Andreas fault, California, where a number of sites showed significant galvanic distortion. We find the regional structure of the phase tensor inversion results compatible with previously obtained models from impedance inversion. In the vicinity of distorted sites, phase tensor inversion models exhibit more homogeneous/smoother conductivity structures.
3D bioprinting of functional human skin: production and in vivo analysis.
Cubo, Nieves; Garcia, Marta; Del Cañizo, Juan F; Velasco, Diego; Jorcano, Jose L
2016-12-05
Significant progress has been made over the past 25 years in the development of in vitro-engineered substitutes that mimic human skin, either to be used as grafts for the replacement of lost skin, or for the establishment of in vitro human skin models. In this sense, laboratory-grown skin substitutes containing dermal and epidermal components offer a promising approach to skin engineering. In particular, a human plasma-based bilayered skin generated by our group, has been applied successfully to treat burns as well as traumatic and surgical wounds in a large number of patients in Spain. There are some aspects requiring improvements in the production process of this skin; for example, the relatively long time (three weeks) needed to produce the surface required to cover an extensive burn or a large wound, and the necessity to automatize and standardize a process currently performed manually. 3D bioprinting has emerged as a flexible tool in regenerative medicine and it provides a platform to address these challenges. In the present study, we have used this technique to print a human bilayered skin using bioinks containing human plasma as well as primary human fibroblasts and keratinocytes that were obtained from skin biopsies. We were able to generate 100 cm(2), a standard P100 tissue culture plate, of printed skin in less than 35 min (including the 30 min required for fibrin gelation). We have analysed the structure and function of the printed skin using histological and immunohistochemical methods, both in 3D in vitro cultures and after long-term transplantation to immunodeficient mice. In both cases, the generated skin was very similar to human skin and, furthermore, it was indistinguishable from bilayered dermo-epidermal equivalents, handmade in our laboratories. These results demonstrate that 3D bioprinting is a suitable technology to generate bioengineered skin for therapeutical and industrial applications in an automatized manner.
3D Cell Entrapment as a Function of the Weight Percent of Peptide-Amphiphile Hydrogels
Scott, Carolyn M.; Forster, Colleen L.; Kokkoli, Efrosini
2015-01-01
The design of scaffolds which mimic the stiffness, nanofiber structure, and biochemistry of the native extra-cellular matrix (ECM) has been a major objective for the tissue engineering field. Furthermore, mimicking the innate three dimensional (3D) environment of the ECM has been shown to significantly alter cellular response compared to traditional two dimensional (2D) culture. We report the development of a self-assembling, fibronectin-mimetic, peptide-amphiphile nanofiber scaffold for 3D cell culture. To form such a scaffold, 5 mol% of a bioactive PR_g fibronectin-mimetic peptide-amphiphile was mixed with 95 mol% of a diluent peptide-amphiphile (E2) whose purpose was to neutralize electrostatic interactions, increase the gelation kinetics and promote cell survival. Atomic force microscopy verified the fibrilar structure of the gels and the mechanical properties were characterized for various weight percent (wt%) formulations of the 5 mol% PR_g - 95 mol% E2 peptide-amphiphile mixture. The 0.5 wt% formulations had an elastic modulus of 429.0 ± 21.3 Pa while the 1.0 wt% peptide-amphiphile hydrogels had an elastic modulus of 808.6 ± 38.1 Pa. The presence of entrapped cells in the gels decreased the elastic modulus and the decrease was a function of the cell loading. While both formulations supported cell proliferation, the 0.5 wt% gels supported significantly greater NIH3T3/GFP fibroblast cell proliferation throughout the gels than the 1.0 wt% gels. However, compared to the 0.5 wt% formulations, the 1.0 wt% hydrogels promoted greater increase in mRNA expression and production of fibronectin and type IV collagen ECM proteins. This study suggests that this fibronectin-mimetic scaffold holds great promise in the advance of 3D culture applications and cell therapies. PMID:25970351
Understanding the core-halo relation of quantum wave dark matter from 3D simulations.
Schive, Hsi-Yu; Liao, Ming-Hsuan; Woo, Tak-Pong; Wong, Shing-Kwong; Chiueh, Tzihong; Broadhurst, Tom; Hwang, W-Y Pauchy
2014-12-31
We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter, described by the Schrödinger-Poisson (SP) equation with a particle mass ∼10(-22) eV. A distinct gravitationally self-bound solitonic core is found at the center of every halo, with a profile quite different from cores modeled in the warm or self-interacting dark matter scenarios. Furthermore, we show that each solitonic core is surrounded by an extended halo composed of large fluctuating dark matter granules which modulate the halo density on a scale comparable to the diameter of the solitonic core. The scaling symmetry of the SP equation and the uncertainty principle tightly relate the core mass to the halo specific energy, which, in the context of cosmological structure formation, leads to a simple scaling between core mass (Mc) and halo mass (Mh), Mc∝a(-1/2)Mh(1/3), where a is the cosmic scale factor. We verify this scaling relation by (i) examining the internal structure of a statistical sample of virialized halos that form in our 3D cosmological simulations and by (ii) merging multiple solitons to create individual virialized objects. Sufficient simulation resolution is achieved by adaptive mesh refinement and graphic processing units acceleration. From this scaling relation, present dwarf satellite galaxies are predicted to have kiloparsec-sized cores and a minimum mass of ∼10(8)M⊙, capable of solving the small-scale controversies in the cold dark matter model. Moreover, galaxies of 2×10(12)M⊙ at z=8 should have massive solitonic cores of ∼2×10(9)M⊙ within ∼60 pc. Such cores can provide a favorable local environment for funneling the gas that leads to the prompt formation of early stellar spheroids and quasars.
In Vitro Model of the Epidermis: Connecting Protein Function to 3D Structure.
Arnette, Christopher; Koetsier, Jennifer L; Hoover, Paul; Getsios, Spiro; Green, Kathleen J
2016-01-01
Much of our understanding of the biological processes that underlie cellular functions in humans, such as cell-cell communication, intracellular signaling, and transcriptional and posttranscriptional control of gene expression, has been acquired from studying cells in a two-dimensional (2D) tissue culture environment. However, it has become increasingly evident that the 2D environment does not support certain cell functions. The need for more physiologically relevant models prompted the development of three-dimensional (3D) cultures of epithelial, endothelial, and neuronal tissues (Shamir & Ewald, 2014). These models afford investigators with powerful tools to study the contribution of spatial organization, often in the context of relevant extracellular matrix and stromal components, to cellular and tissue homeostasis in normal and disease states.
A novel time dependent gamma evaluation function for dynamic 2D and 3D dose distributions.
Podesta, Mark; Persoon, Lucas C G G; Verhaegen, Frank
2014-10-21
Modern external beam radiotherapy requires detailed verification and quality assurance so that confidence can be placed on both the delivery of a single treatment fraction and on the consistency of delivery throughout the treatment course. To verify dose distributions, a comparison between prediction and measurement must be made. Comparisons between two dose distributions are commonly performed using a Gamma evaluation which is a calculation of two quantities on a pixel by pixel basis; the dose difference, and the distance to agreement. By providing acceptance criteria (e.g. 3%, 3 mm), the function will find the most appropriate match within its two degrees of freedom. For complex dynamic treatments such as IMRT or VMAT it is important to verify the dose delivery in a time dependent manner and so a gamma evaluation that includes a degree of freedom in the time domain via a third parameter, time to agreement, is presented here. A C++ (mex) based gamma function was created that could be run on either CPU and GPU computing platforms that would allow a degree of freedom in the time domain. Simple test cases were created in both 2D and 3D comprising of simple geometrical shapes with well-defined boundaries varying over time. Changes of varying magnitude in either space or time were introduced and repeated gamma analyses were performed varying the criteria. A clinical VMAT case was also included, artificial air bubbles of varying size were introduced to a patient geometry, along with shifts of varying magnitude in treatment time. For all test cases where errors in distance, dose or time were introduced, the time dependent gamma evaluation could accurately highlight the errors.The time dependent gamma function presented here allows time to be included as a degree of freedom in gamma evaluations. The function allows for 2D and 3D data sets which are varying over time to be compared using appropriate criteria without penalising minor offsets of subsequent radiation fields
Skin-Friction Measurements in a 3-D, Supersonic Shock-Wave/Boundary-Layer Interaction
NASA Technical Reports Server (NTRS)
Wideman, J. K.; Brown, J. L.; Miles, J. B.; Ozcan, O.
1994-01-01
The experimental documentation of a three-dimensional shock-wave/boundary-layer interaction in a nominal Mach 3 cylinder, aligned with the free-stream flow, and 20 deg. half-angle conical flare offset 1.27 cm from the cylinder centerline. Surface oil flow, laser light sheet illumination, and schlieren were used to document the flow topology. The data includes surface-pressure and skin-friction measurements. A laser interferometric skin friction data. Included in the skin-friction data are measurements within separated regions and three-dimensional measurements in highly-swept regions. The skin-friction data will be particularly valuable in turbulence modeling and computational fluid dynamics validation.
Pseudo 3-D P wave refraction seismic monitoring of permafrost in steep unstable bedrock
NASA Astrophysics Data System (ADS)
Krautblatter, Michael; Draebing, Daniel
2014-02-01
permafrost in steep rock walls can cause hazardous rock creep and rock slope failure. Spatial and temporal patterns of permafrost degradation that operate at the scale of instability are complex and poorly understood. For the first time, we used P wave seismic refraction tomography (SRT) to monitor the degradation of permafrost in steep rock walls. A 2.5-D survey with five 80 m long parallel transects was installed across an unstable steep NE-SW facing crestline in the Matter Valley, Switzerland. P wave velocity was calibrated in the laboratory for water-saturated low-porosity paragneiss samples between 20°C and -5°C and increases significantly along and perpendicular to the cleavage by 0.55-0.66 km/s (10-13%) and 2.4-2.7 km/s (>100%), respectively, when freezing. Seismic refraction is, thus, technically feasible to detect permafrost in low-porosity rocks that constitute steep rock walls. Ray densities up to 100 and more delimit the boundary between unfrozen and frozen bedrock and facilitate accurate active layer positioning. SRT shows monthly (August and September 2006) and annual active layer dynamics (August 2006 and 2007) and reveals a contiguous permafrost body below the NE face with annual changes of active layer depth from 2 to 10 m. Large ice-filled fractures, lateral onfreezing of glacierets, and a persistent snow cornice cause previously unreported permafrost patterns close to the surface and along the crestline which correspond to active seasonal rock displacements up to several mm/a. SRT provides a geometrically highly resolved subsurface monitoring of active layer dynamics in steep permafrost rocks at the scale of instability.
3D printed fluidics with embedded analytic functionality for automated reaction optimisation
Capel, Andrew J; Wright, Andrew; Harding, Matthew J; Weaver, George W; Li, Yuqi; Harris, Russell A; Edmondson, Steve; Goodridge, Ruth D
2017-01-01
Additive manufacturing or ‘3D printing’ is being developed as a novel manufacturing process for the production of bespoke micro- and milliscale fluidic devices. When coupled with online monitoring and optimisation software, this offers an advanced, customised method for performing automated chemical synthesis. This paper reports the use of two additive manufacturing processes, stereolithography and selective laser melting, to create multifunctional fluidic devices with embedded reaction monitoring capability. The selectively laser melted parts are the first published examples of multifunctional 3D printed metal fluidic devices. These devices allow high temperature and pressure chemistry to be performed in solvent systems destructive to the majority of devices manufactured via stereolithography, polymer jetting and fused deposition modelling processes previously utilised for this application. These devices were integrated with commercially available flow chemistry, chromatographic and spectroscopic analysis equipment, allowing automated online and inline optimisation of the reaction medium. This set-up allowed the optimisation of two reactions, a ketone functional group interconversion and a fused polycyclic heterocycle formation, via spectroscopic and chromatographic analysis. PMID:28228852
Blouin, Christian; Boucher, Yan; Roger, Andrew J.
2003-01-01
Comparative sequence analysis has been used to study specific questions about the structure and function of proteins for many years. Here we propose a knowledge-based framework in which the maximum likelihood rate of evolution is used to quantify the level of constraint on the identity of a site. We demonstrate that site-rate mapping on 3D structures using datasets of rhodopsin-like G-protein receptors and α- and β-tubulins provides an excellent tool for pinpointing the functional features shared between orthologous and paralogous proteins. In addition, functional divergence within protein families can be inferred by examining the differences in the site rates, the differences in the chemical properties of the side chains or amino acid usage between aligned sites. Two novel analytical methods are introduced to characterize rate- independent functional divergence. These are tested using a dataset of two classes of HMG-CoA reductases for which only one class can perform both the forward and reverse reaction. We show that functionally divergent sites occur in a cluster of sites interacting with the catalytic residues and that this information should facilitate the design of experimental strategies to directly test functional properties of residues. PMID:12527789
Blouin, Christian; Boucher, Yan; Roger, Andrew J
2003-01-15
Comparative sequence analysis has been used to study specific questions about the structure and function of proteins for many years. Here we propose a knowledge-based framework in which the maximum likelihood rate of evolution is used to quantify the level of constraint on the identity of a site. We demonstrate that site-rate mapping on 3D structures using datasets of rhodopsin-like G-protein receptors and alpha- and beta-tubulins provides an excellent tool for pinpointing the functional features shared between orthologous and paralogous proteins. In addition, functional divergence within protein families can be inferred by examining the differences in the site rates, the differences in the chemical properties of the side chains or amino acid usage between aligned sites. Two novel analytical methods are introduced to characterize rate- independent functional divergence. These are tested using a dataset of two classes of HMG-CoA reductases for which only one class can perform both the forward and reverse reaction. We show that functionally divergent sites occur in a cluster of sites interacting with the catalytic residues and that this information should facilitate the design of experimental strategies to directly test functional properties of residues.
ProFunc: a server for predicting protein function from 3D structure.
Laskowski, Roman A; Watson, James D; Thornton, Janet M
2005-07-01
ProFunc (http://www.ebi.ac.uk/thornton-srv/databases/ProFunc) is a web server for predicting the likely function of proteins whose 3D structure is known but whose function is not. Users submit the coordinates of their structure to the server in PDB format. ProFunc makes use of both existing and novel methods to analyse the protein's sequence and structure identifying functional motifs or close relationships to functionally characterized proteins. A summary of the analyses provides an at-a-glance view of what each of the different methods has found. More detailed results are available on separate pages. Often where one method has failed to find anything useful another may be more forthcoming. The server is likely to be of most use in structural genomics where a large proportion of the proteins whose structures are solved are of hypothetical proteins of unknown function. However, it may also find use in a comparative analysis of members of large protein families. It provides a convenient compendium of sequence and structural information that often hold vital functional clues to be followed up experimentally.
NASA Astrophysics Data System (ADS)
Greenhalgh, Stewart; Zhou, Bing; Maurer, Hansruedi
2010-05-01
We have developed a modified version of the spectral element method (SEM), called the Gaussian Quadrature Grid (GQG) approach, for frequency domain 3D seismic modelling in arbitrary heterogeneous, anisotropic media. The model may incorporate an arbitrary free-surface topography and irregular subsurface interfaces. Unlike the SEM ,it does not require a powerful mesh generator such as the Delauney Triangular or TetGen. Rather, the GQG approach replaces the element mesh with Gaussian quadrature abscissae to directly sample the physical properties of the model parameters and compute the weighted residual or variational integral. This renders the model discretisation simple and easily matched to the model topography, as well as direct control of the model paramterisation for subsequent inversion. In addition, it offers high accuracy in numerical modelling provided that an appropriate density of the Gaussian quadrature abscissae is employed. The second innovation of the GQG is the incorporation of a new implementation of perfectly matched layers to suppress artificial reflections from the domain margins. We employ PML model parameters (specified complex valued density and elastic moduli) rather than explicitly solving the governing wave equation with a complex co-ordinate system as in conventional approaches. Such an implementation is simple, general, effective and easily extendable to any class of anisotropy and other numerical modelling methods. The accuracy of the GQG approach is controlled by the number of Gaussian quadrature points per minimum wavelength, the so-called sampling density. The optimal sampling density should be the one which enables high definition of geological characteristics and high precision of the variational integral evaluation and spatial differentiation. Our experiments show that satisfactory results can be obtained using sampling densities of 5 points per minimum wavelength. Efficiency of the GQG approach mainly depends on the linear
NASA Astrophysics Data System (ADS)
Moulik, P.; Ekstrom, G.
2012-12-01
We have developed a framework that can be used to investigate anisotropic velocity, density and anelastic heterogeneity in the Earth's mantle using a wide spectrum (0.3-50 mHz) of seismological observables. We start with the extensive dataset of surface-wave phase anomalies, long-period waveforms, and body-wave travel times collected by Kustowski et al. (2008) for the development of the global model S362ANI. The additional data included in our analysis are splitting functions of spheroidal and toroidal modes, which are analogous to phase velocity maps at low frequencies. We include in this set of observations a new dataset containing the splitting functions of 56 spheroidal fundamental modes and overtones, measured by Deuss et al. (2011, 2012) using data from large recent earthquakes. Apart from providing unique constraints on the long-wavelength elastic and density structure in the mantle, the overtone splitting data are especially sensitive to the velocity (and anisotropic) structure in the transition zone and in the deeper mantle. The detection of anisotropy, a marker of flow, in the transition zone has implications for our understanding of mantle convection. Our forward modeling of the splitting functions, like the other types of data, includes the effects of radial anisotropy (Mochizuki, 1986). We show that the upper-mantle shear-wave anisotropy of S362ANI generates a clear contribution to the splitting functions of the modes that are sensitive to the upper-mantle structure. We explore the tradeoffs between fitting the mode splitting functions and the travel-times of body waves that turn in the transition zone or in the lower mantle (e.g. SS), while observing that the waveforms and the surface wave phase-anomalies provide complementary information about the mantle. Our experiments suggest that the splitting data are sufficiently sensitive to the anisotropy in the mantle such that their inclusion may provide a better depth resolution of the anisotropic shear
Capturing atmospheric effects on 3D millimeter wave radar propagation patterns
NASA Astrophysics Data System (ADS)
Cook, Richard D.; Fiorino, Steven T.; Keefer, Kevin J.; Stringer, Jeremy
2016-05-01
Traditional radar propagation modeling is done using a path transmittance with little to no input for weather and atmospheric conditions. As radar advances into the millimeter wave (MMW) regime, atmospheric effects such as attenuation and refraction become more pronounced than at traditional radar wavelengths. The DoD High Energy Laser Joint Technology Offices High Energy Laser End-to-End Operational Simulation (HELEEOS) in combination with the Laser Environmental Effects Definition and Reference (LEEDR) code have shown great promise simulating atmospheric effects on laser propagation. Indeed, the LEEDR radiative transfer code has been validated in the UV through RF. Our research attempts to apply these models to characterize the far field radar pattern in three dimensions as a signal propagates from an antenna towards a point in space. Furthermore, we do so using realistic three dimensional atmospheric profiles. The results from these simulations are compared to those from traditional radar propagation software packages. In summary, a fast running method has been investigated which can be incorporated into computational models to enhance understanding and prediction of MMW propagation through various atmospheric and weather conditions.
Intensity images and statistics from numerical simulation of wave propagation in 3-D random media.
Martin, J M; Flatté, S M
1988-06-01
An extended random medium is modeled by a set of 2-D thin Gaussian phase-changing screens with phase power spectral densities appropriate to the natural medium being modeled. Details of the algorithm and limitations on its application to experimental conditions are discussed, concentrating on power-law spectra describing refractive-index fluctuations of the neutral atmosphere. Inner and outer scale effects on intensity scintillation spectra and intensity variance are also included. Images of single realizations of the intensity field at the observing plane are presented, showing that under weak scattering the small-scale Fresnel length structure of the medium dominates the intensity scattering pattern. As the strength of scattering increases, caustics and interference fringes around focal regions begin to form. Finally, in still stronger scatter, the clustering of bright regions begins to reflect the large-scale structure of the medium. For plane waves incident on the medium, physically reasonable inner scales do not produce the large values of intensity variance observed in the focusing region during laser propagation experiments over kilometer paths in the atmosphere. Values as large as experimental observations have been produced in the simulations, but they require inner scales of the order of 10 cm. Inclusion of an outer scale depresses the low-frequency end of the intensity spectrum and reduces the maximum of the intensity variance. Increasing the steepness of the power law also slightly increases the maximum value of intensity variance.
3D mapping of somatotopic reorganization with small animal functional MRI
Yu, Xin; Wang, Shumin; Chen, Der-Yow; Dodd, Stephen; Goloshevsky, Artem; Koretsky, Alan P.
2009-01-01
There are few in vivo noninvasive methods to study neuroplasticity in animal brains. Functional MRI (fMRI) has been developed for animal brain mapping, but few fMRI studies have analyzed functional alteration due to plasticity in animal models. One major limitation is that fMRI maps are characterized by statistical parametric mapping making the apparent boundary dependent on the statistical threshold used. Here, we developed a method to characterize the location of center-of-mass in fMRI maps that is shown not to be sensitive to statistical threshold. Utilizing centers-of-mass as anchor points to fit the spatial distribution of the BOLD response enabled quantitative group analysis of altered boundaries of functional somatosensory maps. This approach was used to study cortical reorganization in the rat primary somatosensory cortex (S1) after sensory deprivation to the barrel cortex by follicle ablation (F.A.). FMRI demonstrated an enlarged nose S1 representation in the 3D somatotopic functional maps. This result clearly demonstrates that fMRI enables the spatial mapping of functional changes that can characterize multiple regions of S1 cortex and still be sensitive to changes due to plasticity. PMID:19770051
NASA Technical Reports Server (NTRS)
Zhao, W.; Newman, J. C., Jr.; Sutton, M. A.; Wu, X. R.; Shivakumar, K. N.
1995-01-01
Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.
Detection of hidden objects using a real-time 3-D millimeter-wave imaging system
NASA Astrophysics Data System (ADS)
Rozban, Daniel; Aharon, Avihai; Levanon, Assaf; Abramovich, Amir; Yitzhaky, Yitzhak; Kopeika, N. S.
2014-10-01
Millimeter (mm)and sub-mm wavelengths or terahertz (THz) band have several properties that motivate their use in imaging for security applications such as recognition of hidden objects, dangerous materials, aerosols, imaging through walls as in hostage situations, and also in bad weather conditions. There is no known ionization hazard for biological tissue, and atmospheric degradation of THz radiation is relatively low for practical imaging distances. We recently developed a new technology for the detection of THz radiation. This technology is based on very inexpensive plasma neon indicator lamps, also known as Glow Discharge Detector (GDD), that can be used as very sensitive THz radiation detectors. Using them, we designed and constructed a Focal Plane Array (FPA) and obtained recognizable2-dimensional THz images of both dielectric and metallic objects. Using THz wave it is shown here that even concealed weapons made of dielectric material can be detected. An example is an image of a knife concealed inside a leather bag and also under heavy clothing. Three-dimensional imaging using radar methods can enhance those images since it can allow the isolation of the concealed objects from the body and environmental clutter such as nearby furniture or other people. The GDDs enable direct heterodyning between the electric field of the target signal and the reference signal eliminating the requirement for expensive mixers, sources, and Low Noise Amplifiers (LNAs).We expanded the ability of the FPA so that we are able to obtain recognizable 2-dimensional THz images in real time. We show here that the THz detection of objects in three dimensions, using FMCW principles is also applicable in real time. This imaging system is also shown here to be capable of imaging objects from distances allowing standoff detection of suspicious objects and humans from large distances.
A coupled wave-3-D hydrodynamics model of the Taranto Sea (Italy): a multiple-nesting approach
NASA Astrophysics Data System (ADS)
Gaeta, Maria Gabriella; Samaras, Achilleas G.; Federico, Ivan; Archetti, Renata; Maicu, Francesco; Lorenzetti, Giuliano
2016-09-01
The present work describes an operational strategy for the development of a multiscale modeling system, based on a multiple-nesting approach and open-source numerical models. The strategy was applied and validated for the Gulf of Taranto in southern Italy, scaling large-scale oceanographic model results to high-resolution coupled wave-3-D hydrodynamics simulations for the area of Mar Grande in the Taranto Sea. The spatial and temporal high-resolution simulations were performed using the open-source TELEMAC suite, forced by wind data from the COSMO-ME database, boundary wave spectra from the RON buoy at Crotone and results from the Southern Adriatic Northern Ionian coastal Forecasting System (SANIFS) regarding sea levels and current fields. Model validation was carried out using data collected in the Mar Grande basin from a fixed monitoring station and during an oceanographic campaign in October 2014. The overall agreement between measurements and model results in terms of waves, sea levels, surface currents, circulation patterns and vertical velocity profiles is deemed to be satisfactory, and the methodology followed in the process can constitute a useful tool for both research and operational applications in the same field and as support of decisions for management and design of infrastructures.
NASA Astrophysics Data System (ADS)
Freniere, Cole; Pathak, Ashish; Raessi, Mehdi
2016-11-01
Ocean Wave Energy Converters (WECs) are devices that convert energy from ocean waves into electricity. To aid in the design of WECs, an advanced computational framework has been developed which has advantages over conventional methods. The computational framework simulates the performance of WECs in a virtual wave tank by solving the full Navier-Stokes equations in 3D, capturing the fluid-structure interaction, nonlinear and viscous effects. In this work, we present simulations of the performance of pitching cylinder-type WECs and compare against experimental data. WECs are simulated at both model and full scales. The results are used to determine the role of the Keulegan-Carpenter (KC) number. The KC number is representative of viscous drag behavior on a bluff body in an oscillating flow, and is considered an important indicator of the dynamics of a WEC. Studying the effects of the KC number is important for determining the validity of the Froude scaling and the inviscid potential flow theory, which are heavily relied on in the conventional approaches to modeling WECs. Support from the National Science Foundation is gratefully acknowledged.
Pollitz, F.F.
2002-01-01
I present a new algorithm for calculating seismic wave propagation through a three-dimensional heterogeneous medium using the framework of mode coupling theory originally developed to perform very low frequency (f < ???0.01-0.05 Hz) seismic wavefield computation. It is a Greens function approach for multiple scattering within a defined volume and employs a truncated traveling wave basis set using the locked mode approximation. Interactions between incident and scattered wavefields are prescribed by mode coupling theory and account for the coupling among surface waves, body waves, and evanescent waves. The described algorithm is, in principle, applicable to global and regional wave propagation problems, but I focus on higher frequency (typically f ??????0.25 Hz) applications at regional and local distances where the locked mode approximation is best utilized and which involve wavefields strongly shaped by propagation through a highly heterogeneous crust. Synthetic examples are shown for P-SV-wave propagation through a semi-ellipsoidal basin and SH-wave propagation through a fault zone.
Senthivel, Vivek Raj; Sturrock, Marc; Piedrafita, Gabriel; Isalan, Mark
2016-01-01
Nonlinear responses to signals are widespread natural phenomena that affect various cellular processes. Nonlinearity can be a desirable characteristic for engineering living organisms because it can lead to more switch-like responses, similar to those underlying the wiring in electronics. Steeper functions are described as ultrasensitive, and can be applied in synthetic biology by using various techniques including receptor decoys, multiple co-operative binding sites, and sequential positive feedbacks. Here, we explore the inherent non-linearity of a biological signaling system to identify functions that can potentially be exploited using cell genome engineering. For this, we performed genome-wide transcription profiling to identify genes with ultrasensitive response functions to Hepatocyte Growth Factor (HGF). We identified 3,527 genes that react to increasing concentrations of HGF, in Madin-Darby canine kidney (MDCK) cells, grown as cysts in 3D collagen cell culture. By fitting a generic Hill function to the dose-responses of these genes we obtained a measure of the ultrasensitivity of HGF-responsive genes, identifying a subset with higher apparent Hill coefficients (e.g. MMP1, TIMP1, SNORD75, SNORD86 and ERRFI1). The regulatory regions of these genes are potential candidates for future engineering of synthetic mammalian gene circuits requiring nonlinear responses to HGF signalling. PMID:27982133
3D Functional Elements Deep Inside Silicon with Nonlinear Laser Lithography
NASA Astrophysics Data System (ADS)
Tokel, Onur; Turnali, Ahmet; Ergecen, Emre; Pavlov, Ihor; Ilday, Fatih Omer
Functional optical and electrical elements fabricated on silicon (Si) constitute fundamental building blocks of electronics and Si-photonics. However, since the highly successful established lithography are geared towards surface processing, elements embedded inside Si simply do not exist. Here, we present a novel direct-laser writing method for positioning buried functional elements inside Si wafers. This new phenomenon is distinct from previous work, in that the surface of Si is not modified. By exploiting nonlinear interactions of a focused laser, permanent refractive index changes are induced inside Si. The imprinted index contrast is then used to demonstrate a plethora of functional elements and capabilities embedded inside Si. In particular, we demonstrate the first functional optical element inside Si, the first information-storage capability inside Si, creation of high-resolution subsurface holograms, buried multilevel structures, and complex 3D architectures in Si, none of which is currently possible with other methods. This new approach complements available techniques by taking advantage of the real estate under Si, and therefore can pave the way for creating entirely new multilevel devices through electronic-photonic integration.
NASA Astrophysics Data System (ADS)
Pilz, Marco; Parolai, Stefano; Woith, Heiko
2017-01-01