3D Elastic Seismic Wave Propagation Code
Energy Science and Technology Software Center (ESTSC)
1998-09-23
E3D is capable of simulating seismic wave propagation in a 3D heterogeneous earth. Seismic waves are initiated by earthquake, explosive, and/or other sources. These waves propagate through a 3D geologic model, and are simulated as synthetic seismograms or other graphical output.
Fullwave coupling to a 3D antenna code using Green's function formulation of wave-particle response
NASA Astrophysics Data System (ADS)
Wright, John; Bonoli, P. T.; Bilato, R.; Brambilla, M.; Maggiora, R.; Lancellotti, V.
2006-10-01
Using the fullwave code, TORIC, and the 3D antenna code, TOPICA, we construct a complete linear system for the RF driven plasma. The 3D finite element antenna code, TOPICA, requires an admittance, Y, for the plasma, where B=YE. In this work, TORIC was modified to allow excitation of the (Eη, Eζ) electric field components at the plasma surface, corresponding to a single poloidal and toroidal mode number combination (m,n). This leads to the tensor response: Yn= ( ll Yηη& YηζYζη& Yζζ), where each of the Yn submatrices is Nm in size. It is shown that the admittance matrix is equivalent to a Green's function calculation for the fullwave system and the net work done is less than twice a single fullwave calculation. The admittance calculation is used with loading calculation from TOPICA to construct self consistent plasma and antenna currents.
Numerical simulation of 3D breaking waves
NASA Astrophysics Data System (ADS)
Fraunie, Philippe; Golay, Frederic
2015-04-01
Numerical methods dealing with two phase flows basically can be classified in two ways : the "interface tracking" methods when the two phases are resolved separately including boundary conditions fixed at the interface and the "interface capturing" methods when a single flow is considered with variable density. Physical and numerical properties of the two approaches are discussed, based on some numerical experiments performed concerning 3D breaking waves. Acknowledgements : This research was supported by the Modtercom program of Region PACA.
Fullwave coupling to a 3D antenna code using Green's function formulation of wave-particle response.
NASA Astrophysics Data System (ADS)
Wright, John; Bonoli, Paul; Brambilla, Marco; Lancelloti, Vito; Maggiora, Riccardo; Carter, Mark
2006-04-01
Using the fullwave code, TORIC ,and the 3D antenna code, TOPICA, we construct a complete linear system for the RF driven plasma. The 3D finite element antenna code, TOPICA, requires an admittance, Y, for the plasma, where B=YE. In this work TORIC was modified to allow excitation of the (Eη, Eζ) electric field components at the plasma surface, corresponding to a single poloidal and toroidal mode number combination (m,n). This leads the tensor response: Y=( *20c Yηη & Yηζ Yζη & Yζζ ), where each of the Yn submatrices is Nm in size. It is shown that the admittance matrix is equivalent to a Greens function calculation for the fullwave system and in addition, the net work done in the calculation is on the order of twice a single fullwave calculation. After the admittance calculation is done, the response of a plasma to an antenna driven at a given frequency can be calculated by only running the TOPICA code for a new antenna geometry. In tests of loading, TOPICA has been able reproduce loading of the Alcator D antenna (S12 coefficient accurately.).
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.
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.
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.
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.
Wave-CAIPI for Highly Accelerated 3D Imaging
Bilgic, Berkin; Gagoski, Borjan A.; Cauley, Stephen F.; Fan, Audrey P.; Polimeni, Jonathan R.; Grant, P. Ellen; Wald, Lawrence L.; Setsompop, Kawin
2014-01-01
Purpose To introduce the Wave-CAIPI (Controlled Aliasing in Parallel Imaging) acquisition and reconstruction technique for highly accelerated 3D imaging with negligible g-factor and artifact penalties. Methods The Wave-CAIPI 3D acquisition involves playing sinusoidal gy and gz gradients during the readout of each kx encoding line, while modifying the 3D phase encoding strategy to incur inter-slice shifts as in 2D-CAIPI acquisitions. The resulting acquisition spreads the aliasing evenly in all spatial directions, thereby taking full advantage of 3D coil sensitivity distribution. By expressing the voxel spreading effect as a convolution in image space, an efficient reconstruction scheme that does not require data gridding is proposed. Rapid acquisition and high quality image reconstruction with Wave-CAIPI is demonstrated for high-resolution magnitude and phase imaging and Quantitative Susceptibility Mapping (QSM). Results Wave-CAIPI enables full-brain gradient echo (GRE) acquisition at 1 mm isotropic voxel size and R=3×3 acceleration with maximum g-factors of 1.08 at 3T, and 1.05 at 7T. Relative to the other advanced Cartesian encoding strategies 2D-CAIPI and Bunched Phase Encoding, Wave-CAIPI yields up to 2-fold reduction in maximum g-factor for 9-fold acceleration at both field strengths. Conclusion Wave-CAIPI allows highly accelerated 3D acquisitions with low artifact and negligible g-factor penalties, and may facilitate clinical application of high-resolution volumetric imaging. PMID:24986223
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
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.
Subduction zone guided waves: 3D modelling and attenuation effects
NASA Astrophysics Data System (ADS)
Garth, T.; Rietbrock, A.
2013-12-01
Waveform modelling is an important tool for understanding complex seismic structures such as subduction zone waveguides. These structures are often simplified to 2D structures for modelling purposes to reduce computational costs. In the case of subduction zone waveguide affects, 2D models have shown that dispersed arrivals are caused by a low velocity waveguide, inferred to be subducted oceanic crust and/or hydrated outer rise normal faults. However, due to the 2D modelling limitations the inferred seismic properties such as velocity contrast and waveguide thickness are still debated. Here we test these limitations with full 3D waveform modelling. For waveguide effects to be observable the waveform must be accurately modelled to relatively high frequencies (> 2 Hz). This requires a small grid spacing due to the high seismic velocities present in subduction zones. A large area must be modelled as well due to the long propagation distances (400 - 600 km) of waves interacting with subduction zone waveguides. The combination of the large model area and small grid spacing required means that these simulations require a large amount of computational resources, only available at high performance computational centres like the UK National super computer HECTOR (used in this study). To minimize the cost of modelling for such a large area, the width of the model area perpendicular to the subduction trench (the y-direction) is made as small as possible. This reduces the overall volume of the 3D model domain. Therefore the wave field is simulated in a model ';corridor' of the subduction zone velocity structure. This introduces new potential sources of error particularly from grazing wave side reflections in the y-direction. Various dampening methods are explored to reduce these grazing side reflections, including perfectly matched layers (PML) and more traditional exponential dampening layers. Defining a corridor model allows waveguide affects to be modelled up to at least 2
Acetylcholinesterase: From 3D Structure to Function
Dvir, Hay; Silman, Israel; Harel, Michal; Rosenberry, Terrone L.; Sussman, Joel L.
2010-01-01
By rapid hydrolysis of the neurotransmitter, acetylcholine, acetylcholinesterase terminates neurotransmission at cholinergic synapses. Acetylcholinesterase is a very fast enzyme, functioning at a rate approaching that of a diffusion-controlled reaction. The powerful toxicity of organophosphate poisons is attributed primarily to their potent inhibition of acetylcholinesterase. Acetylcholinesterase inhibitors are utilized in the treatment of various neurological disorders, and are the principal drugs approved thus far by the FDA for management of Alzheimer’s disease. Many organophosphates and carbamates serve as potent insecticides, by selectively inhibiting insect acetylcholinesterase. The determination of the crystal structure of Torpedo californica acetylcholinesterase permitted visualization, for the first time, at atomic resolution, of a binding pocket for acetylcholine. It also allowed identification of the active site of acetylcholinesterase, which, unexpectedly, is located at the bottom of a deep gorge lined largely by aromatic residues. The crystal structure of recombinant human acetylcholinesterase in its apo-state is similar in its overall features to that of the Torpedo enzyme; however, the unique crystal packing reveals a novel peptide sequence which blocks access to the active-site gorge. PMID:20138030
Nonlinear dynamics of Airy-vortex 3D wave packets: emission of vortex light waves.
Driben, Rodislav; Meier, Torsten
2014-10-01
The dynamics of 3D Airy-vortex wave packets is studied under the action of strong self-focusing Kerr nonlinearity. Emissions of nonlinear 3D waves out of the main wave packets with the topological charges were demonstrated. Because of the conservation of the total angular momentum, charges of the emitted waves are equal to those carried by the parental light structure. The rapid collapse imposes a severe limitation on the propagation of multidimensional waves in Kerr media. However, the structure of the Airy beam carrier allows the coupling of light from the leading, most intense peak into neighboring peaks and consequently strongly postpones the collapse. The dependence of the critical input amplitude for the appearance of a fast collapse on the beam width is studied for wave packets with zero and nonzero topological charges. Wave packets carrying angular momentum are found to be much more resistant to the rapid collapse. PMID:25360922
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. PMID:27043763
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
Improving 3D Genome Reconstructions Using Orthologous and Functional Constraints
Diament, Alon; Tuller, Tamir
2015-01-01
The study of the 3D architecture of chromosomes has been advancing rapidly in recent years. While a number of methods for 3D reconstruction of genomic models based on Hi-C data were proposed, most of the analyses in the field have been performed on different 3D representation forms (such as graphs). Here, we reproduce most of the previous results on the 3D genomic organization of the eukaryote Saccharomyces cerevisiae using analysis of 3D reconstructions. We show that many of these results can be reproduced in sparse reconstructions, generated from a small fraction of the experimental data (5% of the data), and study the properties of such models. Finally, we propose for the first time a novel approach for improving the accuracy of 3D reconstructions by introducing additional predicted physical interactions to the model, based on orthologous interactions in an evolutionary-related organism and based on predicted functional interactions between genes. We demonstrate that this approach indeed leads to the reconstruction of improved models. PMID:26000633
A Self-Consistent Beam Loaded Travelling Wave Accelerator Model for use in TRACE-3D
NASA Astrophysics Data System (ADS)
Lampel, M. C.
1997-05-01
An optics model of a constant gradient traveling wave (CGTW) accelerator structure has been implemented for TRACE-3D. TRACE-3D is an envelope code including space charge that is used to model bunched beams in magnetic transport systems and radio frequency (rf) accelerators when the effects of beam current might be significant. The new matrix model has been developed to allow incorporation of particle beam loading (current) effects on the accelerator gradient and the accelerator structure's beam focusing properties in a self-consistent manner. The beam loaded electric field for a CGTW accelerator structure is constant for only a particular design current (e.g., 0 current), otherwise it can be written as a function of accelerator attenuation and axial position along the structure. The variation of the electric field through the structure has been taken into account in the new model. CGTW structures differ substantially in focusing properties and beam loading properties from standing wave structures. Examples will be presented using the new TW model, propagating electron beams with different currents through the Stanford Linear Accelerator Center's 3 m structure. The results will be compared to the zero current TW structure model in TRANSPORT and the Tank model (a standing wave structure model) in TRACE-3D. A computer demonstration of the code with the new element will also be presented.
Solving tolerancing and 3D beam shaping problems by multifunctional wave optical design
NASA Astrophysics Data System (ADS)
Buehling, Sven; Wyrowski, Frank
2000-10-01
A strategy for designing optical systems that are optimized for multiple optical functions on the basis of wave optics is presented. Each optical function is composed of an input field, a set of fixed system parameters, and a merit function. A design algorithm is proposed which is applicable for optical systems consisting of an transmission operator followed by an arbitrary linear operator. The goal is to find the transmission operator which is optimal for all optical functions simultaneously. In later design steps, the found transmission operator can be transformed to real optical elements, for instance by using the thin element approximation. It is shown that the algorithm is efficiently applicable by investigating two sample applications for multifunctional wave optical design: the design of tolerant systems and 3D beam shaping.
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
3D RNA and Functional Interactions from Evolutionary Couplings.
Weinreb, Caleb; Riesselman, Adam J; Ingraham, John B; Gross, Torsten; Sander, Chris; Marks, Debora S
2016-05-01
Non-coding RNAs are ubiquitous, but the discovery of new RNA gene sequences far outpaces the research on the structure and functional interactions of these RNA gene sequences. We mine the evolutionary sequence record to derive precise information about the function and structure of RNAs and RNA-protein complexes. As in protein structure prediction, we use maximum entropy global probability models of sequence co-variation to infer evolutionarily constrained nucleotide-nucleotide interactions within RNA molecules and nucleotide-amino acid interactions in RNA-protein complexes. The predicted contacts allow all-atom blinded 3D structure prediction at good accuracy for several known RNA structures and RNA-protein complexes. For unknown structures, we predict contacts in 160 non-coding RNA families. Beyond 3D structure prediction, evolutionary couplings help identify important functional interactions-e.g., at switch points in riboswitches and at a complex nucleation site in HIV. Aided by increasing sequence accumulation, evolutionary coupling analysis can accelerate the discovery of functional interactions and 3D structures involving RNA. PMID:27087444
Simulation of 3D Seismic Wave Propagation with Volcano Topography
NASA Astrophysics Data System (ADS)
Ripperger, J.; Igel, H.; Wassermann, J.
2001-12-01
We investigate the possibilities of using three-dimensional finite difference (FD) methods for numerical simulation of the seismic wave field at active volcanoes. We put special emphasis on the implementation of the boundary conditions for free surface topography. We compare two different approaches to solve the free surface boundary conditions. The algorithms are implemented on parallel hardware and have been tested for correctness and stability. We apply them to smooth artificial topographies and to the real topography of Mount Merapi, Indonesia. We conclude, that grid stretching type methods (e.g. Hestholm & Ruud, 1994) are not well suited for realistic volcano topography as they tend to become unstable for large topographic gradients. The representation of topography through staircase shaped grids (Ohminato & Chouet, 1997) results in stable calculations, while demanding very fine gridding. The simulations show the effects of a three-dimensional surface topography on elastic wave propagation. Ground motion at the surface is severely affected by topography. If neglected, this may jeopardize attempts to determine source location by analyzing particle motion. Numerical studies like this can help to understand wave propagation phenomena observed on field recordings in volcano seismology. Future studies will aim at separating the wave effects of internal scattering, topography and sources (tremors, tectonic events, pyroclastic flows).
Simulating Seismic Wave Propagation in 3-D Structure: A Case Study For Istanbul City
NASA Astrophysics Data System (ADS)
Yelkenci, Seda; Aktar, Mustafa
2013-04-01
Investigation of the wave propagation around the Marmara Sea, in particular for the city of Istanbul is critical because this target area is identified as one of the megacities with the highest seismic risk in the world. This study makes an attempt for creating an integrated 3D seismic/geologic model and precise understanding of 3-D wave propagation in the city of Istanbul. The approach is based on generating synthetic seismograms using realistic velocity structures as well as accurate location, focal mechanism and source parameters of reference earthquakes. The modarate size reference earthquakes occured in the Marmara Sea and were recorded by the National Seismic Network of Turkey as well as the network of Istanbul Early Warning and Rapid Response System. The seismograms are simulated by means of a 3-D finite difference method operated on parallel processing environment. In the content of creating a robust velocity model; 1D velocity models which are derived fom previous crustal studies of Marmara region such as refraction seismic and receiver functions have been conducted firstly for depths greater than 1km. Velocity structure in shallower part of the study region is then derived from recent geophysical and geotechnical surveys. To construct 3-D model from the obtained 1-D model data, a variety of interpolation methods are considered. According to the observations on amplitude and arrival time based on comparison of simulated seismograms, the considered velocity model is refined the way that S delay times are compensated. Another important task of this work is an application of the finite difference method to estimate three-dimensional seismic responses for a specified basin structure including soft sediments with low shear velocities in respect of the surrounded area in the Asian part of Istanbul. The analysis performed both in the time and frequency domain, helps in understanding of the comprehensive wave propagation characteristics and the distribution of
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
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.
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.
Joint inversion of 3D crustal structure with ambient noise and earthquake body wave travel time
NASA Astrophysics Data System (ADS)
Li, Z.; Ni, S.; Chong, J.; Wang, X.
2012-12-01
Surface wave tomography based on the noise correlation function of seismic ambient noise has been widely used in studies of crustal and mantle structure . However, the periods of surface wave dispersions in the ambient noise tomography are typically less than 40 s, which limits its resolution on the lower crust. Travel times of earthquake body waves, such as Sg and SmS, could provide additional constraints to the crustal structure, especially to the lower crust due to the ray paths of SmS traveling through the lower crust twice. Here, we proposed a joint inversion method for 3D crustal structure with ambient noise and earthquake body wave travel time data, with the goal of providing better constraints and resolutions on the whole crust. We constructed the linear equations for joint inversion of crustal S velocity structure with the surface wave dispersion and body wave travel time data, and solved the equations with LSQR algorithm. Different weighting and damping factors, together with smoothing constraints, are adopted for surface wave dispersion and body wave travel time data to fit both dataset simultaneously. Synthetics experiments showed that the joint inversion could resolve the crust structure better than sole tomography of ambient noise or body wave travel time. We conducted the joint inversion around the Yangtze block in the eastern China. Rayleigh wave dispersions are extracted from the seismic ambient noise tomography by Zheng et al (2011) in this area. The body waves (e.g., Sg, SmS, Sn) are coherent to be identified and their travel times are measured with accuracy from high quality waveforms of some recent local earthquakes in this area. In order to minimize the travel time uncertainties, the focal depth and epicenter of these local earthquakes were resolved by depth phases and temporary aftershock observations. The result from joint inversion suggests that the crustal velocity structure, especially the lower crust, was well improved, which not only
A new algorithm for 3D reconstruction from support functions.
Gardner, Richard J; Kiderlen, Markus
2009-03-01
We introduce a new algorithm for reconstructing an unknown shape from a finite number of noisy measurements of its support function. The algorithm, based on a least squares procedure, is very easy to program in standard software such as Matlab, and it works for both 2D and 3D reconstructions (in fact, in principle, in any dimension). Reconstructions may be obtained without any pre- or post-processing steps and with no restriction on the sets of measurement directions except their number, a limitation dictated only by computing time. An algorithm due to Prince and Willsky was implemented earlier for 2D reconstructions, and we compare the performance of their algorithm and ours. But our algorithm is the first that works for 3D reconstructions with the freedom stated in the previous paragraph. Moreover, under mild conditions, theory guarantees that outputs of the new algorithm will converge to the input shape as the number of measurements increases. In addition we offer a linear program version of the new algorithm that is much faster and better, or at least comparable, in performance at low levels of noise and reasonably small numbers of measurements. Another modification of the algorithm, suitable for use in a "focus of attention" scheme, is also described. PMID:19147881
Characterizing the propagation of gravity waves in 3D nonlinear simulations of solar-like stars
NASA Astrophysics Data System (ADS)
Alvan, L.; Strugarek, A.; Brun, A. S.; Mathis, S.; Garcia, R. A.
2015-09-01
Context. The revolution of helio- and asteroseismology provides access to the detailed properties of stellar interiors by studying the star's oscillation modes. Among them, gravity (g) modes are formed by constructive interferences between progressive internal gravity waves (IGWs), propagating in stellar radiative zones. Our new 3D nonlinear simulations of the interior of a solar-like star allows us to study the excitation, propagation, and dissipation of these waves. Aims: The aim of this article is to clarify our understanding of the behavior of IGWs in a 3D radiative zone and to provide a clear overview of their properties. Methods: We use a method of frequency filtering that reveals the path of individual gravity waves of different frequencies in the radiative zone. Results: We are able to identify the region of propagation of different waves in 2D and 3D, to compare them to the linear raytracing theory and to distinguish between propagative and standing waves (g-modes). We also show that the energy carried by waves is distributed in different planes in the sphere, depending on their azimuthal wave number. Conclusions: We are able to isolate individual IGWs from a complex spectrum and to study their propagation in space and time. In particular, we highlight in this paper the necessity of studying the propagation of waves in 3D spherical geometry, since the distribution of their energy is not equipartitioned in the sphere.
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. PMID:21342824
3D surface-wave tomography in the central Baltic Shield
NASA Astrophysics Data System (ADS)
Bruneton, M.; Pedersen, H. A.; Farra, V.; Svekalapko Seismic Tomography Working Group
2003-04-01
The main objective of the SVEKALAPKO deep seismic experiment was to image in details the lithosphere-asthenosphere system of the central Baltic Shield, therefore enhancing our knowledge of the structure and evolution of cratonic lithosphere. During the experiment a regular 2D grid of 46 broad-band stations covered the southern part of Finland. This exceptional stations distribution made it possible to undertake a high precision surface-wave tomography. We developed a technique based on paraxial ray tracing to obtain 2D phase-velocity maps as a function of frequency which can subsequently be inverted for the 3D structure. The major improvement is that we jointly inverted for the velocity model under the array and the shape of incoming wave fronts, therefore reducing artifacts due to structure outside the study region. The data set included carefully selected fundamental mode Rayleigh wave arrival times of 69 teleseismic events, computed using Wiener filtering. An average dispersion curve was obtained imposing the phase-velocity to be quasi constant. It leads to shear-wave velocities for the lithospheric mantle 4% faster than standard Earth model ak135. The inversion of the same data set was also conducted using weaker constraints to obtain the lateral variations of the phase-velocity at each frequency and subsequently of the shear-wave velocity as a function of depth. Three Vs profiles were computed respectively in the Karelian Archean province, in the Proterozoic Svekofennia, and at the suture between the two domains. They showed significant variations, the higher lithospheric velocities were seen in the proterozoic domain, a low velocity zone was necessary only in the suture zone. Our results showed that chemical changes are maintained within the lithosphere over extended periods of time.
3D modeling of ultrasonic wave interaction with disbonds and weak bonds
NASA Astrophysics Data System (ADS)
Leckey, C.; Hinders, M.
2012-05-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.
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)
Tong, Ping; Komatitsch, Dimitri; Tseng, Tai-Lin; Hung, Shu-Huei; Chen, Chin-Wu; Basini, Piero; Liu, Qinya
2014-10-01
We present a three-dimensional (3-D) hybrid method that interfaces the spectral-element method (SEM) with the frequency-wave number (FK) technique to model the propagation of teleseismic plane waves beneath seismic arrays. The accuracy of the resulting 3-D SEM-FK hybrid method is benchmarked against semianalytical FK solutions for 1-D models. The accuracy of 2.5-D modeling based on 2-D SEM-FK hybrid method is also investigated through comparisons to this 3-D hybrid method. Synthetic examples for structural models of the Alaska subduction zone and the central Tibet crust show that this method is capable of accurately capturing interactions between incident plane waves and local heterogeneities. This hybrid method presents an essential tool for the receiver function and scattering imaging community to verify and further improve their techniques. These numerical examples also show the promising future of the 3-D SEM-FK hybrid method in high-resolution regional seismic imaging based on waveform inversions of converted/scattered waves recorded by seismic array.
3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Wright, G. B.; Yuen, D.
2008-12-01
In the past 25 years a wide variety of numerical methods, such as finite-difference, finite-volume , finite- elements, and pseudospectral methods have been employed to study the problem of 3-D mantle convection. All have specialized strengths but also serious weaknesses. The first three methods are generally considered low-order and can involve high algorithmic complexity (as in triangular elements). Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. Here, we introduce a new grid/mesh-free approach using radial basis functions (RBFs). It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and naturally permits local node refinement. It has been shown for shallow-water equations and vortex flows that RBFs outperform other numerical methods in the sense that they obtain a much higher accuracy for the same spatial resolution while being able to take unusually large time steps. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer( planetary surface ) and inner shell ( core-mantle boundary ). The velocity and stress components are expressed in terms of the scalar potential approach (Zebib and Schubert, 1982) and the other remaining variable is the perturbed temperature field. We have studied the problem from the onset of convection to a modest nonlinear regime.
Field structure of collapsing wave packets in 3D strong Langmuir turbulence
NASA Technical Reports Server (NTRS)
Newman, D. L.; Robinson, P. A.; Goldman, M. V.
1989-01-01
A simple model is constructed for the electric fields in the collapsing wave packets found in 3D simulations of driven and damped isotropic strong Langmuir turbulence. This model, based on a spherical-harmonic decomposition of the electrostatic potential, accounts for the distribution of wave-packet shapes observed in the simulations, particularly the predominance of oblate wave packets. In contrast with predictions for undamped and undriven subsonic collapse of scalar fields, oblate vector-field wave packets do not flatten during collapse but, instead, remain approximately self-similar and rigid.
Patterned illumination for analysing neuronal function in 3D
NASA Astrophysics Data System (ADS)
Go, Mary Ann; Stricker, Christian; Redman, Steve; Bachor, Hans; Daria, Vincent R.
2012-06-01
We use patterned 3D multi-spot illumination to perform neuronal multi-site stimulation in rat brain tissue. Using a spatial light modulator, we holograpically project 3D light fields for multi-site two-photon photolysis of caged neurotransmitters to generate synaptic inputs to a neuron. Controlled photostimulation of multiple synapses from various locations in the dendritic tree provides a way to analyze how neurons integrate multiple inputs. Our holographic projection setup is incorporated into a two-photon 3D imaging microscope for visualization and for accurate positioning of specific uncaging sites along the neuron's dendritic tree. We show two-photon images and the neuron's response to holographic photostimulation of synapses along dendrites.
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).
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.
3D WKB solution for fast magnetoacoustic wave behaviour around an X-line
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; Botha, G. J. J.; Régnier, S.; Spoors, D. L.
2016-06-01
Context. We study the propagation of a fast magnetoacoustic wave in a 3D magnetic field created from two magnetic dipoles. The magnetic topology contains an X-line. Aims: We aim to contribute to the overall understanding of MHD wave propagation within inhomogeneous media, specifically around X-lines. Methods: We investigate the linearised, 3D MHD equations under the assumptions of ideal and cold plasma. We utilise the WKB approximation and Charpit's method during our investigation. Results: It is found that the behaviour of the fast magnetoacoustic wave is entirely dictated by the local, inhomogeneous, equilibrium Alfvén speed profile. All parts of the wave experience refraction during propagation, where the magnitude of the refraction effect depends on the location of an individual wave element within the inhomogeneous magnetic field. The X-line, along which the Alfvén speed is identically zero, acts as a focus for the refraction effect. There are two main types of wave behaviour: part of the wave is either trapped by the X-line or escapes the system, and there exists a critical starting region around the X-line that divides these two types of behaviour. For the set-up investigated, it is found that 15.5% of the fast wave energy is trapped by the X-line. Conclusions: We conclude that linear, β = 0 fast magnetoacoustic waves can accumulate along X-lines and thus these will be specific locations of fast wave energy deposition and thus preferential heating. The work here highlights the importance of understanding the magnetic topology of a system. We also demonstrate how the 3D WKB technique described in this paper can be applied to other magnetic configurations.
Isotropic 3D Super-resolution Imaging with a Self-bending Point Spread Function
Jia, Shu; Vaughan, Joshua C.; Zhuang, Xiaowei
2014-01-01
Airy beams maintain their intensity profiles over a large propagation distance without substantial diffraction and exhibit lateral bending during propagation1-5. This unique property has been exploited for micromanipulation of particles6, generation of plasma channels7 and guidance of plasmonic waves8, but has not been explored for high-resolution optical microscopy. Here, we introduce a self-bending point spread function (SB-PSF) based on Airy beams for three-dimensional (3D) super-resolution fluorescence imaging. We designed a side-lobe-free SB-PSF and implemented a two-channel detection scheme to enable unambiguous 3D localization of fluorescent molecules. The lack of diffraction and the propagation-dependent lateral bending make the SB-PSF well suited for precise 3D localization of molecules over a large imaging depth. Using this method, we obtained super-resolution imaging with isotropic 3D localization precision of 10-15 nm over a 3 μm imaging depth from ∼2000 photons per localization. PMID:25383090
Fast Wave Trains Associated with Solar Eruptions: Insights from 3D Thermodynamic MHD Simulations
NASA Astrophysics Data System (ADS)
Downs, C.; Liu, W.; Torok, T.; Linker, J.; Mikic, Z.; Ofman, L.
2015-12-01
EUV imaging observations during the SDO/AIA era have provided new insights into a variety of wave phenomena occurring in the low solar corona. One example is the observation of quasi-periodic, fast-propagating wave trains that are associated with solar eruptions, including flares and CMEs. While there has been considerable progress in understanding such waves from both an observational and theoretical perspective, it remains a challenge to pin down their physical origin. In this work, we detail our results from a case-study 3D thermodynamic MHD simulation of a coronal mass ejection where quasi-periodic wave trains are generated during the simulated eruption. We find a direct correlation between the onset of non-steady reconnection in the flare current sheet and the generation of quasi-periodic wave train signatures when patchy, collimated downflows interact with the flare arcade. Via forward modeling of SDO/AIA observables, we explore how the appearance of the wave trains is affected by line-of-sight integration and the multi-thermal nature of the coronal medium. We also examine how the wave trains themselves are channeled by natural waveguides formed in 3D by the non-uniform background magnetic field. While the physical association of the reconnection dynamics to the generation of quasi-periodic wave trains appears to be a compelling result, unanswered questions posed from recent observations as well as future prospects will be discussed.
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.
Modeling Recent Large Earthquakes Using the 3-D Global Wave Field
NASA Astrophysics Data System (ADS)
Hjörleifsdóttir, V.; Kanamori, H.; Tromp, J.
2003-04-01
We use the spectral-element method (SEM) to accurately compute waveforms at periods of 40 s and longer for three recent large earthquakes using 3D Earth models and finite source models. The M_w~7.6, Jan~26, 2001, Bhuj, India event had a small rupture area and is well modeled at long periods with a point source. We use this event as a calibration event to investigate the effects of 3-D Earth models on the waveforms. The M_w~7.9, Nov~11, 2001, Kunlun, China, event exhibits a large directivity (an asymmetry in the radiation pattern) even at periods longer than 200~s. We used the source time function determined by Kikuchi and Yamanaka (2001) and the overall pattern of slip distribution determined by Lin et al. to guide the wave-form modeling. The large directivity is consistent with a long fault, at least 300 km, and an average rupture speed of 3±0.3~km/s. The directivity at long periods is not sensitive to variations in the rupture speed along strike as long as the average rupture speed is constant. Thus, local variations in rupture speed cannot be ruled out. The rupture speed is a key parameter for estimating the fracture energy of earthquakes. The M_w~8.1, March~25, 1998, event near the Balleny Islands on the Antarctic Plate exhibits large directivity in long period surface waves, similar to the Kunlun event. Many slip models have been obtained from body waves for this earthquake (Kuge et al. (1999), Nettles et al. (1999), Antolik et al. (2000), Henry et al. (2000) and Tsuboi et al. (2000)). We used the slip model from Henry et al. to compute SEM waveforms for this event. The synthetic waveforms show a good fit to the data at periods from 40-200~s, but the amplitude and directivity at longer periods are significantly smaller than observed. Henry et al. suggest that this event comprised two subevents with one triggering the other at a distance of 100 km. To explain the observed directivity however, a significant amount of slip is required between the two subevents
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
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.
NASA Astrophysics Data System (ADS)
Wu, Chensheng; Nelson, William; Davis, Christopher C.
2014-10-01
Plenoptic functions are functions that preserve all the necessary light field information of optical events. Theoretical work has demonstrated that geometric based plenoptic functions can serve equally well in the traditional wave propagation equation known as the "scalar stochastic Helmholtz equation". However, in addressing problems of 3D turbulence simulation, the dominant methods using phase screen models have limitations both in explaining the choice of parameters (on the transverse plane) in real-world measurements, and finding proper correlations between neighboring phase screens (the Markov assumption breaks down). Though possible corrections to phase screen models are still promising, the equivalent geometric approach based on plenoptic functions begins to show some advantages. In fact, in these geometric approaches, a continuous wave problem is reduced to discrete trajectories of rays. This allows for convenience in parallel computing and guarantees conservation of energy. Besides the pairwise independence of simulated rays, the assigned refractive index grids can be directly tested by temperature measurements with tiny thermoprobes combined with other parameters such as humidity level and wind speed. Furthermore, without loss of generality one can break the causal chain in phase screen models by defining regional refractive centers to allow rays that are less affected to propagate through directly. As a result, our work shows that the 3D geometric approach serves as an efficient and accurate method in assessing relevant turbulence problems with inputs of several environmental measurements and reasonable guesses (such as Cn 2 levels). This approach will facilitate analysis and possible corrections in lateral wave propagation problems, such as image de-blurring, prediction of laser propagation over long ranges, and improvement of free space optic communication systems. In this paper, the plenoptic function model and relevant parallel algorithm computing
ATHENA 3D: A finite element code for ultrasonic wave propagation
NASA Astrophysics Data System (ADS)
Rose, C.; Rupin, F.; Fouquet, T.; Chassignole, B.
2014-04-01
The understanding of wave propagation phenomena requires use of robust numerical models. 3D finite element (FE) models are generally prohibitively time consuming. However, advances in computing processor speed and memory allow them to be more and more competitive. In this context, EDF R&D developed the 3D version of the well-validated FE code ATHENA2D. The code is dedicated to the simulation of wave propagation in all kinds of elastic media and in particular, heterogeneous and anisotropic materials like welds. It is based on solving elastodynamic equations in the calculation zone expressed in terms of stress and particle velocities. The particularity of the code relies on the fact that the discretization of the calculation domain uses a Cartesian regular 3D mesh while the defect of complex geometry can be described using a separate (2D) mesh using the fictitious domains method. This allows combining the rapidity of regular meshes computation with the capability of modelling arbitrary shaped defects. Furthermore, the calculation domain is discretized with a quasi-explicit time evolution scheme. Thereby only local linear systems of small size have to be solved. The final step to reduce the computation time relies on the fact that ATHENA3D has been parallelized and adapted to the use of HPC resources. In this paper, the validation of the 3D FE model is discussed. A cross-validation of ATHENA 3D and CIVA is proposed for several inspection configurations. The performances in terms of calculation time are also presented in the cases of both local computer and computation cluster use.
Quantitative analysis of accuracy of seismic wave-propagation codes in 3D random scattering media
NASA Astrophysics Data System (ADS)
Galis, Martin; Imperatori, Walter; Mai, P. Martin
2013-04-01
Several recent verification studies (e.g. Day et al., 2001; Bielak et al., 2010, Chaljub et al., 2010) have demonstrated the importance of assessing the accuracy of available numerical tools at low frequency in presence of large-scale features (basins, topography, etc.). The fast progress in high-performance computing, including efficient optimization of numerical codes on petascale supercomputers, has permitted the simulation of 3D seismic wave propagation at frequencies of engineering interest (up to 10Hz) in highly heterogeneous media (e.g. Hartzell et al, 2010; Imperatori and Mai, 2013). However, high frequency numerical simulations involving random scattering media, characterized by small-scale heterogeneities, are much more challenging for most numerical methods, and their verification may therefore be even more crucial than in the low-frequency case. Our goal is to quantitatively compare the accuracy and the behavior of three different numerical codes for seismic wave propagation in 3D random scattering media at high frequency. We deploy a point source with omega-squared spectrum, and focus on the near-source region, being of great interest in strong motion seismology. We use two codes based on finite-difference method (FD1 and FD2) and one code based on support-operator method (SO). Both FD1 and FD2 are 4-th order staggered-grid finite-difference codes (for FD1 see Olsen et al., 2009; for FD2 see Moczo et al., 2007). The FD1 and FD2 codes are characterized by slightly different medium representations, since FD1 uses point values of material parameters in each FD-cell, while FD2 uses the effective material parameters at each grid-point (Moczo et al., 2002). SO is 2-nd order support-operator method (Ely et al., 2008). We considered models with random velocity perturbations described by van Karman correlation function with different correlation lengths and different standard deviations. Our results show significant variability in both phase and amplitude as
The small data solutions of general 3-D quasilinear wave equations. II
NASA Astrophysics Data System (ADS)
Ding, Bingbing; Witt, Ingo; Yin, Huicheng
2016-07-01
This paper is a continuation of the work in [8], where the authors established the global existence of smooth small data solutions to the general 3-D quasilinear wave equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0 when the weak null condition holds. In the present paper, we show that the smooth small data solutions of equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0 will blow up in finite time when the weak null condition does not hold and a generic nondegenerate condition on the initial data is satisfied, moreover, a precise blowup time is completely determined. Therefore, collecting the main results in this paper and [8], we have given a basically complete study on the blowup or global existence of small data solutions to the 3-D quasilinear wave equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0.
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.
NASA Astrophysics Data System (ADS)
Cheng, Cheng; Bodin, Thomas; Allen, Richard M.
2016-02-01
We present a novel 3D pre-stack Kirchhoff depth migration (PKDM) method for teleseismic receiver functions. The proposed algorithm considers the effects of diffraction, scattering, and travel time alteration caused by 3D volumetric heterogeneities. It is therefore particularly useful for imaging complex 3D structures such as dipping discontinuities, which is hard to accomplish with traditional methods. The scheme is based on the acoustic wave migration principle, where at each time step of the receiver function, the energy is migrated back to the ensemble of potential conversion points in the image, given a smooth 3D reference model. Travel times for P and S waves are computed with an efficient Eikonal solver, the Fast Marching Method. We also consider elastic scattering patterns, where the amplitude of converted S waves depends on the angle between the incident P wave, and the scattered S wave. Synthetic experiments demonstrate the validity of the method for a variety of dipping angle discontinuities. Comparison with the widely used Common Conversion Point (CCP) stacking method reveals that our migration shows considerable improvement. For example, the effect of multiple reflections that usually produce apparent discontinuities is avoided. The proposed approach is practical, computationally efficient, and is therefore a potentially powerful alternative to standard CCP methods for imaging large-scale continental structure under dense networks.
Upscaling small heterogeneities for seismic wave propagation in 3D complex media
NASA Astrophysics Data System (ADS)
Cupillard, P.; Capdeville, Y.
2012-04-01
Seismic waves propagating in the Earth are affected by different sizes of heterogeneities. When modelling these waves (using numerical methods such as the SEM), taking into account heterogeneities that are much smaller than the minimum wavelength is a challenge because meshing small heterogeneities often requires important efforts and leads to high numerical costs. In this work, we present a technique which allows to upscale the small heterogeneities that can lie in an elastic medium. This technique yields a smooth effective medium and effective equations. We describe its implementation in the 3D case and we show relevant examples.
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
Spatial parallelism of a 3D finite difference, velocity-stress elastic wave propagation code
Minkoff, S.E.
1999-12-01
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately, finite difference simulations for 3D elastic wave propagation are expensive. The authors model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MPI library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speedup. Because I/O is handled largely outside of the time-step loop (the most expensive part of the simulation) the authors have opted for straight-forward broadcast and reduce operations to handle I/O. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ghost cells. When this communication is balanced against computation by allocating subdomains of reasonable size, they observe excellent scaled speedup. Allocating subdomains of size 25 x 25 x 25 on each node, they achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Spatial Parallelism of a 3D Finite Difference, Velocity-Stress Elastic Wave Propagation Code
MINKOFF,SUSAN E.
1999-12-09
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately. finite difference simulations for 3D elastic wave propagation are expensive. We model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MP1 library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speed up. Because i/o is handled largely outside of the time-step loop (the most expensive part of the simulation) we have opted for straight-forward broadcast and reduce operations to handle i/o. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ''ghost cells''. When this communication is balanced against computation by allocating subdomains of reasonable size, we observe excellent scaled speed up. Allocating subdomains of size 25 x 25 x 25 on each node, we achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
3D simulation of seismic wave propagation around a tunnel using the spectral element method
NASA Astrophysics Data System (ADS)
Lambrecht, L.; Friederich, W.
2010-05-01
We model seismic wave propagation in the environment of a tunnel for later application to reconnaissance. Elastic wave propagation can be simulated by different numerical techniques such as finite differences and pseudospectral methods. Their disadvantage is the lack of accuracy on free surfaces, numerical dispersion and inflexibility of the mesh. Here we use the software package SPECFEM3D_SESAME in an svn development version, which is based on the spectral element method (SEM) and can handle complex mesh geometries. A weak form of the elastic wave equation leads to a linear system of equations with a diagonal mass matrix, where the free surface boundary of the tunnel can be treated under realistic conditions and can be effectively implemented in parallel. We have designed a 3D external mesh including a tunnel and realistic features such as layers and holes to simulate elastic wave propagation in the zone around the tunnel. The source is acting at the tunnel surface so that we excite Rayleigh waves which propagate to the front face of the tunnel. A conversion takes place and a high amplitude S-wave is radiated in the direction of the tunnel axis. Reflections from perturbations in front of the tunnel can be measured by receivers implemented on the tunnel face. For a shallow tunnel the land surface has high influence on the wave propagation. By implementing additional receivers at this surface we intent to improve the prediction. It shows that the SEM is very capable to handle the complex geometry of the model and especially incorporates the free surfaces of the model.
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
Multiple-mode Lamb wave scattering simulations using 3D elastodynamic finite integration technique.
Leckey, Cara A C; Rogge, Matthew D; Miller, Corey A; Hinders, Mark K
2012-02-01
We have implemented three-dimensional (3D) elastodynamic finite integration technique (EFIT) simulations to model Lamb wave scattering for two flaw-types in an aircraft-grade aluminum plate, a rounded rectangle flat-bottom hole and a disbond of the same shape. The plate thickness and flaws explored in this work include frequency-thickness regions where several Lamb wave modes exist and sometimes overlap in phase and/or group velocity. For the case of the flat-bottom hole the depth was incrementally increased to explore progressive changes in multiple-mode Lamb wave scattering due to the damage. The flat-bottom hole simulation results have been compared to experimental data and are shown to provide key insight for this well-defined experimental case by explaining unexpected results in experimental waveforms. For the rounded rectangle disbond flaw, which would be difficult to implement experimentally, we found that Lamb wave behavior differed significantly from the flat-bottom hole flaw. Most of the literature in this field is restricted to low frequency-thickness regions due to difficulties in interpreting data when multiple modes exist. We found that benchmarked 3D EFIT simulations can yield an understanding of scattering behavior for these higher frequency-thickness regions and in cases that would be difficult to set up experimentally. Additionally, our results show that 2D simulations would not have been sufficient for modeling the complicated scattering that occurred. PMID:21908011
NASA Astrophysics Data System (ADS)
Vidal, A.; San-Blas, A. A.; Quesada-Pereira, F. D.; Pérez-Soler, J.; Gil, J.; Vicente, C.; Gimeno, B.; Boria, V. E.
2015-07-01
A novel technique for the full-wave analysis of 3-D complex waveguide devices is presented. This new formulation, based on the Boundary Integral-Resonant Mode Expansion (BI-RME) method, allows the rigorous full-wave electromagnetic characterization of 3-D arbitrarily shaped metallic structures making use of extremely low CPU resources (both time and memory). The unknown electric current density on the surface of the metallic elements is represented by means of Rao-Wilton-Glisson basis functions, and an algebraic procedure based on a singular value decomposition is applied to transform such functions into the classical solenoidal and nonsolenoidal basis functions needed by the original BI-RME technique. The developed tool also provides an accurate computation of the electromagnetic fields at an arbitrary observation point of the considered device, so it can be used for predicting high-power breakdown phenomena. In order to validate the accuracy and efficiency of this novel approach, several new designs of band-pass waveguides filters are presented. The obtained results (S-parameters and electromagnetic fields) are successfully compared both to experimental data and to numerical simulations provided by a commercial software based on the finite element technique. The results obtained show that the new technique is specially suitable for the efficient full-wave analysis of complex waveguide devices considering an integrated coaxial excitation, where the coaxial probes may be in contact with the metallic insets of the component.
3-D P Wave Velocity Structure of Marmara Region Using Local Earthquake Tomography
NASA Astrophysics Data System (ADS)
Işık, S. E.; Gurbuz, C.
2014-12-01
The 3D P wave velocity model of upper and lower crust of the Marmara Region between 40.200- 41.200N and 26.500- 30.500E is obtained by tomographic inversion (Simulps) of 47034 P wave arrivals of local earthquakes recorded at 90 land stations between October 2009 and December 2012 and 30 OBO stations and 14162 shot arrivals recorded at 35 OBO stations (Seismarmara Survey, 2001). We first obtained a 1D minimum model with Velest code in order to obtain an initial model for 3D inversion with 648 well located earthquakes located within the study area. After several 3D inversion trials we decided to create a more adequate initial model for 3D inversion. Choosing the initial model we estimated the 3D P wave velocity model representing the whole region both for land and sea. The results are tested by making Checkerboard , Restoring Resolution and Characteristic Tests, and the reliable areas of the resulting model is defined in terms of RDE, DWS, SF and Hit count distributions. By taking cross sections from the resulting model we observed the vertical velocity change along profiles crossing both land and sea. All the profiles crossing the basins showed that the high velocities of lower crust make extensions towards the basin area which looks like the force that gives a shape to the basins. These extensions of lower crust towards the basins appeared with an average velocity of 6.3 km/s which might be the result of the deformation due the shearing in the region. It is also interpreted that the development of these high velocities coincide with the development of the basins. Thus, both the basins and the high velocity zones around them might be resulted from the entrance of the NAF into the Marmara Sea and at the same time a shear regime was dominated due to the resistance of the northern Marmara Region (Yılmaz, 2010). The seismicity is observed between 5 km and 15 km after the 3D location of the earthquakes. The locations of the earthquakes improved and the seismogenic zone
Study of nonlinear 3-D evolution of kinetic Alfvén wave and fluctuation spectra
NASA Astrophysics Data System (ADS)
Sharma, Prachi; Yadav, Nitin; Sharma, R. P.
2015-11-01
Waves and instabilities play a very crucial role in astrophysical plasmas e.g. solar wind, Geospace etc. The main objective of current study is to investigate the importance of nonlinear processes associated with kinetic Alfvén waves (KAWs) in order to understand the physical mechanism behind the magnetopause turbulence. Numerical simulation of the coupled equations guiding the dynamics of three dimensionally propagating kinetic Alfvén wave (KAW) and slow magnetosonic wave has been performed for intermediate beta plasma (i.e. me/mi ≪ β < 1, where β is thermal to magnetic pressure ratio) applicable to the magnetopause. A simplified semi-analytical model based on paraxial approach has also been developed. We have examined the field localization and associated power spectrum of 3-D kinetic Alfvén wave for this nonlinear interaction. Governing dynamical equations of KAW and slow magnetosonic wave get coupled when the ponderomotive force arising due to pump KAW is taken into account while studying the slow magnetosonic wave dynamics. The numerical prediction of power law scaling is just consistent with the observation of THEMIS spacecraft in the magnetopause.
NASA Astrophysics Data System (ADS)
Wilder, F. D.; Ergun, R.; Goodrich, K.; Malaspina, D.; Eriksson, S.; Stawarz, J. E.; Sturner, A. P.; Holmes, J.; Burch, J. L.; Torbert, R. B.; Phan, T.; Le Contel, O.; Goldman, M. V.; Newman, D. L.; Lindqvist, P. A.; Khotyaintsev, Y. V.; Strangeway, R. J.; Russell, C. T.; Giles, B. L.; Pollock, C. J.
2015-12-01
The phenomenon of magnetic reconnection, especially at electron scales, is still poorly understood. One process that warrants further investigation is the role of wave phenomenon in mediating magnetic reconnection. Previous observations have shown the presence of electrostatic solitary waves (ESWs) as well as whistler mode waves near the dayside reconnection site. Additionally, recent simulations have suggested that whistler waves might be generated by electron phase space holes associated with ESWs as they propagate along the magnetic separatrix towards the diffusion region. Other observations have shown ESWs with distinct speeds and time scales, suggesting that different instabilities generate the ESWs. NASA's recently launched Magnetospheric Multiscale (MMS) mission presents a unique opportunity to investigate the roles of wave phenomena, such as ESWs and whistlers, in asymmetric reconnection at the dayside magnetopause. We will present 3-D electric and magnetic field data from magnetopause crossings by MMS during its first dayside science phase. Burst mode wave data and electron distributions from all four spacecraft will be analyzed to investigate the origin of these wave phenomena, as well as their impact on the reconnection electric field.
Numerical investigation of wave attenuation by vegetation using a 3D RANS model
NASA Astrophysics Data System (ADS)
Marsooli, Reza; Wu, Weiming
2014-12-01
Vegetation has been recognized as an important natural shoreline protection against storm surges and waves. Understanding of wave-vegetation interaction is essential for assessing the ability of vegetation patches, such as wetlands, to mitigate storm damages. In this study the wave attenuation by vegetation is investigated numerically using a 3-D model which solves the Reynolds-Averaged Navier-Stokes equations (RANS) by means of a finite-volume method based on collocated hexahedron mesh. A mixing length model is used for turbulence closure of the RANS equations. The water surface boundary is tracked using the Volume-of-Fluid (VOF) method with the Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM) to solve the VOF advection equation. The presence of vegetation is taken into account by adding the vegetation drag and inertia forces to the momentum equations. The model is validated by several laboratory experiments of short wave propagation through vegetation over flat and sloping beds. The comparisons show good agreement between the measured data and calculated results, but the swaying motion of flexible vegetation which is neglected in this study can influence the accuracy of the wave height predictions. The model is then applied to one of the validation tests with different vegetation properties, revealing that the wave height attenuation by vegetation depends not only on the wave conditions, but also the vegetation characteristics such as vegetation height and density.
Scattering and coupling effects of electromagnetic waves in 3D networks of spheres
NASA Astrophysics Data System (ADS)
Defos Du Rau, M.; Pessan, F.; Ruffie, G.; Vignéras-Lefebvre, V.; Parneix, J. P.
1998-01-01
In this paper, the problem of electromagnetic scattering from a 3D system of spheres is considered and an iterative solution that accounts for multiple scattering is proposed. The Mie formalism used for a single sphere is extended to account for multiple scattered fields between several particles. The translational addition theorems for spherical wave functions are used to express the electromagnetic field scattered by a sphere S_i in terms of an incident field for a sphere S_k in a spherical coordinates system attached to the sphere S_k. In this work, the numerical convergence of the method is discussed and associated computational times are given. Numerical computations including Radar Cross Section (RCS) and radiation patterns for various 3D configurations are presented. Some of them are compared with free-space measurements made in the 8 to 100 GHz frequency band using vectorial network analyzers. 11.55.-m S-matrix theory; analytic structure of amplitudes Cet article étudie la diffusion des ondes électromagnétiques par des réseaux tridimensionnels de sphères et propose une méthode itérative pour prendre en compte les effets de multidiffusion. Le formalisme de Mie utilisé dans le cas d'une sphère est étendu pour calculer les champs "multidiffusés" entre plusieurs particules. Les théorèmes d'addition et de translation des fonctions d'onde sphériques sont utilisés pour exprimer le champ diffusé par une sphère S_i comme étant incident sur une sphère S_k, dans un système de coordonnées sphériques lié au centre de S_k. La convergence numérique de la méthode est discutée et des temps de calcul sont donnés. Des résultats numériques tels que des Surfaces Équivalentes Radar (SER) et des diagrammes de rayonnement pour différentes configurations tridimensionnelles sont montrés. Certains d'entre eux sont comparés à des mesures en espace libre faites à l'aide d'analyseurs de réseaux vectoriels dans la bande de fréquence 8{-}100 GHz.
Characteristics of Oceanic Waves Caused by Landslides: Insights from 3D-Hydrocode Modeling
NASA Astrophysics Data System (ADS)
Elbeshausen, D.; Wünnemann, K.; Weiss, R.
2008-12-01
In an experimental framework, the generation of tsunami waves can be considered as a two-dimensional or three-dimensional problem by finding the respective geometry of the experimental set up. In nature, of course, it is a fully three-dimensional problem. The generation of tsunami waves caused by landslides (submarine and subaerial) must be approached as a dynamical problem. The understanding of the slide body's dynamics plays a key role in understanding the generated waves. Numerical calculations are a standard tool in tsunami science as the propagation of long waves can be tackled with depth-averaged equations. These classical models have often been used for modeling the propagation and run up of those tsunami waves caused by earthquakes. Tsunamis generated during slide motion are different. They are shorter and have larger amplitudes. It could be demonstrated in respective laboratory experiments and two dimensional numerical studies that in the initial phase the behavior of the waves is very complex, resulting in wave breaking and plunging. We conducted hydrocode simulations to model the generation of tsunami by slide in two dimensions. Results could show the complexity of the initial wave evolution as well as the development of the slide body itself. As an extension to these two-dimensional simulations, we now consider the three-dimensional problem and reveal some differences to the two-dimensional results. For this purpose we are using iSALE-3D, a multi- material, multi-rheology hydrocode capable of studying landslide processes in both two and three dimensions. iSALE-3D has been originally developed to study shock waves and high pressure scenarios like meteorite impacts or explosions and has been successfully validated against theoretical and experimental results as well as other numerical codes. Previous studies revealed differences in the formation and propagation of oceanic waves caused by meteorite impacts at different angles of incidence. Hence, for
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.
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.
Concealed threat detection with the IRAD sub-millimeter wave 3D imaging radar
NASA Astrophysics Data System (ADS)
Robertson, Duncan A.; Cassidy, Scott L.; Jones, Ben; Clark, Anthony
2014-06-01
Sub-millimeter wave 3D imaging radar is a promising technology for the stand-off detection of threats concealed on people. The IRAD 340 GHz 3D imaging radar uses polarization intensity information to identify signatures associated with concealed threats. We report on an extensive trials program which has been carried out involving dozens of individual subjects wearing a variety of different clothing to evaluate the detection of a wide range of threat and benign items. We have developed an automatic algorithm to run on the radar which yields a level of anomaly indication in real time. Statistical analysis of the large volume of recorded data has enabled performance metrics for the radar system to be evaluated.
A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2008-01-01
We apply the method of support operators (SOM) to solve the 3-D, viscoelastic equations of motion for use in earthquake simulations. SOM is a generalized finite-difference method that can utilize meshes of arbitrary structure and incorporate irregular geometry. Our implementation uses a 3-D, logically rectangular, hexahedral mesh. Calculations are second-order in space and time. A correction term is employed for suppression of spurious zero-energy modes (hourglass oscillations). We develop a free surface boundary condition, and an absorbing boundary condition using the method of perfectly matched layers (PML). Numerical tests using a layered material model in a highly deformed mesh show good agreement with the frequency-wavenumber method, for resolutions greater than 10 nodes per wavelength. We also test a vertically incident P wave on a semi-circular canyon, for which results match boundary integral solutions at resolutions greater that 20 nodes per wavelength. We also demonstrate excellent parallel scalability of our code.
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
An efficient flexible-order model for 3D nonlinear water waves
Engsig-Karup, A.P. Bingham, H.B.; Lindberg, O.
2009-04-01
The flexible-order, finite difference based fully nonlinear potential flow model described in [H.B. Bingham, H. Zhang, On the accuracy of finite difference solutions for nonlinear water waves, J. Eng. Math. 58 (2007) 211-228] is extended to three dimensions (3D). In order to obtain an optimal scaling of the solution effort multigrid is employed to precondition a GMRES iterative solution of the discretized Laplace problem. A robust multigrid method based on Gauss-Seidel smoothing is found to require special treatment of the boundary conditions along solid boundaries, and in particular on the sea bottom. A new discretization scheme using one layer of grid points outside the fluid domain is presented and shown to provide convergent solutions over the full physical and discrete parameter space of interest. Linear analysis of the fundamental properties of the scheme with respect to accuracy, robustness and energy conservation are presented together with demonstrations of grid independent iteration count and optimal scaling of the solution effort. Calculations are made for 3D nonlinear wave problems for steep nonlinear waves and a shoaling problem which show good agreement with experimental measurements and other calculations from the literature.
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 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
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-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
Numerical homogenization for seismic wave propagation in 3D geological media
NASA Astrophysics Data System (ADS)
Cupillard, P.; Capdeville, Y.; Botella, A.
2014-12-01
Despite the important increase of the computational power in the last decades, simulating the seismic wave propagation through realistic geological models is still a challenge. By realistic models we here mean 3D media in which a broad variety (in terms of amplitude and extent) of heterogeneities lies, including discontinuities with complex geometry such as faulted and folded horizons, intrusive geological contacts and fault systems. To perform accurate numerical simulations, these discontinuities require complicated meshes which usually contain extremely small elements, yielding large, sometimes prohibitive, computation costs. Fortunately, the recent development of the non-periodic homogenization technique now enables to overcome this problem by computing smooth equivalent models for which a coarse mesh is sufficient to get an accurate wavefield. In this work, we present an efficient implementation of the technique which now allows for the homogenization of large 3D geological models. This implementation relies on a tetrahedral finite-element solution of the elasto-static equation behind the homogenization problem. Because this equation is time-independent, solving it is numerically cheaper than solving the wave equation, but it nevertheless requires some care because of the large size of the stiffness matrix arising from the fine mesh of realistic geological structures. A domain decomposition is therefore adopted. In our strategy, the obtained sub-domains overlap but they are independent so the solution within each of them can be computed either in series or in parallel. In addition, well-balanced loads, efficient search algorithms and multithreading are implemented to speed up the computation. The resulting code enables the homogenization of 3D elastic media in a time that is neglectable with respect to the simulation time of the wave propagation within. This is illustrated through a sub-surface model of the Furfooz karstic region, Belgium.
NASA Astrophysics Data System (ADS)
Woodbury, D.; Kubota, S.; Johnson, I.
2014-10-01
Computer simulations of electromagnetic wave propagation in magnetized plasmas are an important tool for both plasma heating and diagnostics. For active millimeter-wave and microwave diagnostics, accurately modeling the evolution of the beam parameters for launched, reflected or scattered waves in a toroidal plasma requires that calculations be done using the full 3-D geometry. Previously, we reported on the application of GPGPU (General-Purpose computing on Graphics Processing Units) to a 3-D vacuum Maxwell code using the FDTD (Finite-Difference Time-Domain) method. Tests were done for Gaussian beam propagation with a hard source antenna, utilizing the parallel processing capabilities of the NVIDIA K20M. In the current study, we have modified the 3-D code to include a soft source antenna and an induced current density based on the cold plasma approximation. Results from Gaussian beam propagation in an inhomogeneous anisotropic plasma, along with comparisons to ray- and beam-tracing calculations will be presented. Additional enhancements, such as advanced coding techniques for improved speedup, will also be investigated. Supported by U.S. DoE Grant DE-FG02-99-ER54527 and in part by the U.S. DoE, Office of Science, WDTS under the Science Undergraduate Laboratory Internship program.
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
NASA Astrophysics Data System (ADS)
Parisi, Laura; Ferreira, Ana M. G.
2016-04-01
The surface wave full ray theory (FRT) is an efficient tool to calculate synthetic waveforms of surface waves. It combines the concept of local modes with exact ray tracing as a function of frequency, providing a more complete description of surface wave propagation than the widely used great circle approximation (GCA). The purpose of this study is to evaluate the ability of the FRT approach to model teleseismic long-period surface waveforms (T ˜ 45-150 s) in the context of current 3-D Earth models to empirically assess its validity domain and its scope for future studies in seismic tomography. To achieve this goal, we compute vertical and horizontal component fundamental mode synthetic Rayleigh waveforms using the FRT, which are compared with calculations using the highly accurate spectral element method. We use 13 global earth models including 3-D crustal and mantle structure, which are derived by successively varying the strength and lengthscale of heterogeneity in current tomographic models. For completeness, GCA waveforms are also compared with the spectral element method. We find that the FRT accurately predicts the phase and amplitude of long-period Rayleigh waves (T ˜ 45-150 s) for almost all the models considered, with errors in the modelling of the phase (amplitude) of Rayleigh waves being smaller than 5 per cent (10 per cent) in most cases. The largest errors in phase and amplitude are observed for T ˜ 45 s and for the three roughest earth models considered that exhibit shear wave anomalies of up to ˜20 per cent, which is much larger than in current global tomographic models. In addition, we find that overall the GCA does not predict Rayleigh wave amplitudes well, except for the longest wave periods (T ˜ 150 s) and the smoothest models considered. Although the GCA accurately predicts Rayleigh wave phase for current earth models such as S20RTS and S40RTS, FRT's phase errors are smaller, notably for the shortest wave periods considered (T ˜ 45 s and
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
3D rendering of passive millimeter-wave scenes using modified open source software
NASA Astrophysics Data System (ADS)
Murakowski, Maciej; Wilson, John; Murakowski, Janusz; Schneider, Garrett; Schuetz, Christopher; Prather, Dennis
2011-05-01
As millimeter-wave imaging technology becomes more mature, several applications are emerging for which this technology may be useful. However, effectively predicting the nuances of millimeter-wave phenomenology on the usefulness for a given application remains a challenge. To this end, an accurate millimeter-wave scene simulator would have tremendous value in predicting imager requirements for a given application. Herein, we present a passive millimeter-wave scene simulator built on the open-source 3d modeling software Blender. We describe the changes made to the Blender rendering engine to make it suitable for this purpose, including physically accurate reflections at each material interface, volumetric absorption and scattering, and tracking of both s and p polarizations. In addition, we have incorporated a mmW material database and world model that emulates the effects of cold sky profiles for varying weather conditions and frequencies of operation. The images produced by this model have been validated against calibrated experimental imagery captured by a passive scanning millimeter-wave imager for maritime, desert, and standoff detection applications.
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
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. PMID:26994875
NASA Astrophysics Data System (ADS)
Crosta, G.; Imposimato, S.; Roddeman, D.; Frattini, P.
2012-04-01
Fast moving landslides can be originated along slopes in mountainous terrains with natural and artificial lakes, or fjords at the slope foot. This landslides can reach extremely high speed and the impact with the immobile reservoir water can be influenced by the local topography and the landslide mass profile. The impact can generate large impulse waves and landslide tsunami. Initiation, propagation and runup are the three phases that need to be considered. The landslide evolution and the consequent wave can be controlled by the initial mass position (subaerial, partially or completely submerged), the landslide speed, the type of material, the subaerial and subaqueous slope geometry, the landslide depth and length at the impact, and the water depth. Extreme events have been caused by subaerial landslides: the 1963 Vajont rockslide (Italy), the 1958 Lituya Bay event (Alaska), the Tafjord and the Loen multiple events event (Norway), also from volcanic collapses (Hawaii and Canary islands). Various researchers completed a systematic experimental work on 2D and 3D wave generation and propagation (Kamphuis and Bowering, 1970; Huber, 1980; Müller, 1995; Huber and Hager, 1997; Fritz, 2002; Zweifel, 2004; Panizzo et al., 2005; Heller, 2007; Heller and Kinnear, 2010; Sælevik et al., 2009), using both rigid blocks and deformable granular" masses. Model data and results have been used to calibrate and validate numerical modelling tools (Harbitz, 1992; Jiang and LeBlond, 1993; Grilli et al., 2002; Grilli and Watts, 2005; Lynett and Liu, 2005; Tinti et al., 2006; Abadie et al., 2010) generally considering simplified rheologies (e.g. viscous rheologies) for subaerial subaqueous spreading. We use a FEM code (Roddeman, 2011; Crosta et al., 2006, 2009, 2010, 2011) adopting an Eulerian-Lagrangian approach to give accurate results for large deformations. We model both 2D and fully 3D events considering different settings. The material is considered as a fully deformable elasto
Development of a GPU-Accelerated 3-D Full-Wave Code for Reflectometry Simulations
NASA Astrophysics Data System (ADS)
Reuther, K. S.; Kubota, S.; Feibush, E.; Johnson, I.
2013-10-01
1-D and 2-D full-wave codes used as synthetic diagnostics in microwave reflectometry are standard tools for understanding electron density fluctuations in fusion plasmas. The accuracy of the code depends on how well the wave properties along the ignored dimensions can be pre-specified or neglected. In a toroidal magnetic geometry, such assumptions are never strictly correct and ray tracing has shown that beam propagation is inherently a 3-D problem. Previously, we reported on the application of GPGPU's (General-Purpose computing on Graphics Processing Units) to a 2-D FDTD (Finite-Difference Time-Domain) code ported to utilize the parallel processing capabilities of the NVIDIA C870 and C1060. Here, we report on the development of a FDTD code for 3-D problems. Initial tests will use NVIDIA's M2070 GPU and concentrate on the launching and propagation of Gaussian beams in free space. If available, results using a plasma target will also be presented. Performance will be compared with previous generations of GPGPU cards as well as with NVIDIA's newest K20C GPU. Finally, the possibility of utilizing multiple GPGPU cards in a cluster environment or in a single node will also be discussed. Supported by U.S. DoE Grants DE-FG02-99-ER54527 and DE-AC02-09CH11466 and the DoE National Undergraduate Fusion Fellowship.
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.
3D laboratory experiments on a system of low-crested breakwaters under oblique wave attack
NASA Astrophysics Data System (ADS)
Papacharalampous, Georgia; Karantinos, Michalis; Giantsi, Theodora; Moutzouris, Constantinos
2016-04-01
Low-crested breakwaters are being increasingly used for shore protection. Hydrodynamics around coastal structures are complicated and have not been fully understood. A series of large scale (1:40) 3D laboratory experiments were carried out in the Laboratory of Harbour Works, National Technical University of Athens to investigate the wave disturbance around a system of two non-parallel to the shoreline breakwaters. The structures were of the type of low-crested, permeable and attacked by obliquely incident waves. Three different water depths were tested in the basin with a range of various different spectra. The transmission and reflection coefficients were measured in the middle of each breakwater. For this purpose, 1 gauge and 4 gauges (in line) were placed on the landward and seaward side of each breakwater respectively. The effect of diffraction is incorporate at the measured wave heights. The measured coefficients are being compared to their corresponding estimated using existing empirical formulas. Most of those formulas neglect wave obliquity.
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
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.
NASA Astrophysics Data System (ADS)
Jheng, Y.; Hung, S.; Zhou, Y.; Chang, Y.
2012-12-01
Surface wave travel-time tomography has been widely used as a powerful strategy to image shear wave velocity structure of the Earth's crust and upper mantle, providing comparable information other than body wave tomography. Traditionally, lateral variations of dispersive phase velocities are first obtained at multiple frequencies and then used to invert for shear wave velocity with 1-D depth-dependent sensitivity kernels. However, this approach runs short on considering the directional- and depth-dependence of scattering while surface wave propagating through laterally heterogeneous Earth. To refrain from these shortcomings, we here provide a fully 3-D finite-frequency method based on the Born scattering theory formulated with surface wave mode summation, and apply it to regional fundamental Rayleigh wave travel-time tomography in central Tibet. Our data were collected from Project Hi-CLIMB, which deployed an N-S trending linear array of over 100 broadband seismic stations with a large aperture of 800 km and very dense spacing of ~3-8 km across the Lhasa and Qiangtang terranes during 2004-2005. We follow a standard procedure of ambient noise cross correlation to extract empirical Green's functions of fundamental Rayleigh waves at 10-33 s between station pairs. A multi-taper method is employed to measure the phase differences as a function of period between observed and synthetic Rayleigh waves as well as the corresponding sensitivity kernels for the measured phase delays to 3-D shear wave velocity perturbations in a spherically-symmetric model suitable for central Tibet. A wavelet-based, multi-scale parameterization is invoked in the tomographic inversion to deal with the intrinsically multi-scale nature of unevenly distributed data and resolve the structure with data-adaptive spectral and spatial resolutions. The preliminary result shows that to the north of the Banggong-Nujiang suture (BNS), the crustal shear wave velocity beneath the Qiangtang terrane is
System-in-package LTCC platform for 3D RF to millimeter wave
NASA Astrophysics Data System (ADS)
Vähä-Heikkilä, T.; Lahti, M.
2011-04-01
This presentation shows recent trends and results in 3D Low Temperature Co-Fired Ceramics (LTCC) modules in applications from RF to millimeter waves. The system-in-package LTCC platform is a true three dimensional module technology. LTCC is a lightweight multi-layer technology having typically 6-20 ceramic layers and metallizations between. The metallization levels i.e different metal layers can be patterned and connected together with metal vias. Passive devices can also be fabricated on LTCC while active devices and other chips are connected with flip-chip, wire bonding or soldering. In addition to passives directly fabricated to LTCC, several different technologies/ chips can be hybrid integrated to the same module. LTCC platform is also well suited for the realization of antenna arrays for microwave and millimeter wave applications. Potential applications are ranging from short range communications to space and radars. VTT has designed, fabricated and characterized microwave and millimeter wave packages for Radio Frequency (RF) Micro Electro Mechanical Systems (MEMS) as well as active devices. Also, several types of system-in-package modules have been realized containing hybrid integrated CMOS and GaAs MMICs and antenna arrays.
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.
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)
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
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
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
NASA Astrophysics Data System (ADS)
Berssenbrügge, Philipp; Dekiff, Markus; Kemper, Björn; Denz, Cornelia; Dirksen, Dieter
2012-03-01
The increasing importance of optical 3D measurement techniques and the growing number of available methods and systems require a fast and simple method to characterize the measurement accuracy. However, the conventional approach of comparing measured coordinates to known reference coordinates of a test target faces two major challenges: the precise fabrication of the target and - in case of pattern projecting systems - finding the position of the reference points in the obtained point cloud. The modulation transfer function (MTF) on the other hand is an established instrument to describe the resolution characteristics of 2D imaging systems. Here, the MTF concept is applied to two different topometric systems based on fringe and speckle pattern projection to obtain a 3D transfer function. We demonstrate that in the present case fringe projection provides typically 3.5 times the 3D resolution achieved with speckle pattern projection. By combining measurements of the 3D transfer function with 2D MTF measurements the dependency of 2D and 3D resolutions are characterized. We show that the method allows for a simple comparison of the 3D resolution of two 3D sensors using a low cost test target, which is easy to manufacture.
Effect of spatial dispersion on transient acoustic wave propagation in 3D.
Every, A G
2006-12-22
Spatial dispersion is the variation of wave speed with wavelength. It sets in when the acoustic wavelength approaches the natural scale of length of the medium, which could, for example, be the lattice constant of a crystal, the repeat distance in a superlattice, or the grain size in a granular material. In centrosymmetric media, the first onset of dispersion is accommodated by the introduction of fourth order spatial derivatives into the wave equation. These lead to a correction to the phase velocity which is quadratic in the spatial frequency. This paper treats the effect of spatial dispersion on the point force elastodynamic Green's functions of solids. The effects of dispersion are shown to be most pronounced in the vicinity of wave arrivals. These lose their singular form, and are transformed into wave trains known as quasi-arrivals. The step and ramp function wave arrivals are treated, and it is shown that their unfolded quasi-arrival forms can be expressed in terms of integrals involving the Airy function. PMID:16828830
3D microporous base-functionalized covalent organic frameworks for size-selective catalysis.
Fang, Qianrong; Gu, Shuang; Zheng, Jie; Zhuang, Zhongbin; Qiu, Shilun; Yan, Yushan
2014-03-10
The design and synthesis of 3D covalent organic frameworks (COFs) have been considered a challenge, and the demonstrated applications of 3D COFs have so far been limited to gas adsorption. Herein we describe the design and synthesis of two new 3D microporous base-functionalized COFs, termed BF-COF-1 and BF-COF-2, by the use of a tetrahedral alkyl amine, 1,3,5,7-tetraaminoadamantane (TAA), combined with 1,3,5-triformylbenzene (TFB) or triformylphloroglucinol (TFP). As catalysts, both BF-COFs showed remarkable conversion (96% for BF-COF-1 and 98% for BF-COF-2), high size selectivity, and good recyclability in base-catalyzed Knoevenagel condensation reactions. This study suggests that porous functionalized 3D COFs could be a promising new class of shape-selective catalysts. PMID:24604810
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
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
NASA Astrophysics Data System (ADS)
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-08-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.
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-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
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
NASA Astrophysics Data System (ADS)
Polkowski, Marcin; Grad, Marek
2015-04-01
The 3D P-wave seismic velocity model was obtained by combining data from multiple studies during past 50 years. Data sources included refraction seismology, reflection seismology, geological boreholes, vertical seismic profiling, magnetotellurics and gravimetry. Use of many data sources allowed creation of detailed 3D P-wave velocity model that reaches to depth of 60 km and includes 6-layers of sediments and 3-layers of the crust. Purpose of this study is to analyze how 3D model influences local (accuracy of location and source time estimation for local events), regional (identification of wide-angle seismic phases) and global (teleseismic tomography) seismic travel times. Additionally we compare results of forward seismic wave propagation with signals observed on short period and broadband stations. National Science Centre Poland provided financial support for this work by NCN grant DEC-2011/02/A/ST10/00284.
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.
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.
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.
Propagation of Electromagnetic Waves in 3D Opal-based Magnetophotonic Crystals
NASA Astrophysics Data System (ADS)
Pardavi-Horvath, Martha; Makeeva, Galina S.; Golovanov, Oleg A.; Rinkevich, Anatolii B.
2013-03-01
Opals, a class of self-organized 3D nanostructures, are typical representatives of photonic bandgap structures. The voids inside of the opal structure of close packed SiO2 spheres can be infiltrated by a magnetic material, creating magnetically tunable magnetophotonic crystals with interesting and potentially useful properties at GHz and THz frequencies. The propagation of electromagnetic waves at microwave frequencies was investigated numerically in SiO2 opal based magnetic nanostructures, using rigorous mathematical models to solve Maxwell's equations complemented by the Landau-Lifshitz equation with electrodynamic boundary conditions. The numerical approach is based on Galerkin's projection method using the decomposition algorithm on autonomous blocks with Floquet channels. The opal structure consists of SiO2 nanospheres, with inter-sphere voids infiltrated with nanoparticles of Ni-Zn ferrites. Both the opal matrix and the ferrite are assumed to be lossy. A model, taking into account the real structure of the ferrite particles in the opal's voids was developed to simulate the measured FMR lineshape of the ferrite infiltrated opal. The numerical technique shows an excellent agreement when applied to model recent experimental data on similar ferrite opals.
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.
Area and volume coherence for efficient visualization of 3D scalar functions
Max, N. California Univ., Davis, CA ); Hanrahan, P. ); Crawfis, R. )
1990-01-01
We present an algorithm for compositing a combination of density clouds and contour surfaces used to represent a scalar function on a 3-D volume. The volume is divided into convex polyhedra, at whose vertices the function is known, and the polyhedra are sorted in depth before compositing. For data given at scattered 3-D points, we show that this sorting can be done in O(n) time if we chose the tetrahedra in the Delaunay triangulation as the polyhedra. The integrals for cloud opacity and visible cloud intensity along a ray through a convex polyhedron are computed analytically, and this computation is coherent across the polyhedron's area. 33 refs.
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
Extracting chemical information from plane wave calculations by a 3D 'fuzzy atoms' analysis
NASA Astrophysics Data System (ADS)
Bakó, I.; Stirling, A.; Seitsonen, A. P.; Mayer, I.
2013-03-01
Bond order and valence indices have been calculated by the method of the three-dimensional 'fuzzy atoms' analysis, using the numerical molecular orbitals obtained from plane wave DFT calculations, i.e., without introducing any external atom-centered functions. Weight functions of both Hirshfeld and Becke types have been applied. The results are rather close to the similar 'fuzzy atoms' ones obtained by using atom-centered basis sets and agree well with the chemical expectations, stressing the power of the genuine chemical concepts.
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
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.
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)
Tian, Y.; Ritzwoller, M. H.; Shen, W.; Levshin, A. L.; Barmin, M. P.
2014-12-01
The error in the epicentral location of crustal earthquakes across the contiguous US is on the order of 10 km due to the inability of 1D seismic velocity models to capture regional body wave travel time variations. New high resolution 3D models of the crust and uppermost mantle have been constructed recently across the US by inverting surface wave dispersion from ambient noise and earthquakes, receiver functions, and Rayleigh wave H/V ratios using USArray data [e.g., Shen et al., 2013]. These are mostly S-wave models of the lithosphere, however, which are not optimal for predicting regional P-wave travel times. We explore the use of observations of surface waves to improve regional event characterization because the new 3D models are constructed explicitly to model their behavior. In particular, we use measurements of group and phase time delays and the amplitude ratio between different periods of surface waves to estimate the moment tensor, the epicentral location and the earthquake depth. Preliminary estimates of these variables are determined through a simulated annealing algorithm. Afterward, a Bayesian Monte Carlo method is applied to estimate the posterior distribution of all variables in order to assess uncertainties in source characteristics. The reliability and limitations of the location method are tested by systematic relocation of earthquakes across the contiguous US.
Functional Stereology for 3D Particle Size Distributions from 2D Observations: a Practical Approach
NASA Astrophysics Data System (ADS)
Proussevitch, A. A.; Sahagian, D. L.; Jutzeler, M.
2010-12-01
Functional stereology applies known deconvolution techniques to obtain 3D size distributions from 2D cross-section data based on an assumption that both 2D and 3D statistics have known distribution functions with unknown parameters. A new stereological approach enables us to solve this problem by utilizing function minimization to find parameters of the distribution functions. There is no limit to continuous distribution function types that could be used, such as Gaussian, Logistic, Weibull, Gamma, and others. As compared to previously known 2D to 3D conversion methods (e.g. Sahagian and Proussevitch, 1998), functional stereology has much greater practical application to non-spherical particles/objects because it is free of uncontrollable error propagation for all particles shapes. The new practical method of functional stereology has been implemented in Stereonet software adapted for both a) direct logarithmic scales of particle/voids volumes, and b) Phi units of linear dimensions (-log2 of size). Applications of the method include distribution of voids/bubbles in all types of volcanic rocks, pore networks in sedimentary rocks, mineral and grain sizes, volcanic clasts, breccia, and texture features of a wide range of rock formations. Such applications demonstrate utility of this functional stereology approach.
Using CATH-Gene3D to Analyze the Sequence, Structure, and Function of Proteins.
Sillitoe, Ian; Lewis, Tony; Orengo, Christine
2015-01-01
The CATH database is a classification of protein structures found in the Protein Data Bank (PDB). Protein structures are chopped into individual units of structural domains, and these domains are grouped together into superfamilies if there is sufficient evidence that they have diverged from a common ancestor during the process of evolution. A sister resource, Gene3D, extends this information by scanning sequence profiles of these CATH domain superfamilies against many millions of known proteins to identify related sequences. Thus the combined CATH-Gene3D resource provides confident predictions of the likely structural fold, domain organisation, and evolutionary relatives of these proteins. In addition, this resource incorporates annotations from a large number of external databases such as known enzyme active sites, GO molecular functions, physical interactions, and mutations. This unit details how to access and understand the information contained within the CATH-Gene3D Web pages, the downloadable data files, and the remotely accessible Web services. PMID:26087950
3D Hot Test Simulations of a 220 GHz Folded Waveguide Traveling Wave Tube Using a CFDTD PIC Method
NASA Astrophysics Data System (ADS)
Lin, Ming-Chieh; Song, Heather
2015-11-01
Millimeter or sub-THz wave sources centered at 220 GHz is of interest due to the potential for its commercial and military applications including high resolution radar, remote sensing, and high-data-rate communications. It has been demonstrated via 3D cold test finite element method (FEM) simulations that a folded waveguide traveling wave tube (FWTWT) can be designed and optimized at this frequency range with a small signal gain of 18 dB over a comparatively broad (-3 dB) bandwidth of ~ 10%. On the other hand, 3D hot test simulations of a V-band ladder TWT have been successfully demonstrated using a conformal finite-difference time-domain (CFDTD) particle-in-cell (PIC) method for center frequency of 50 GHz. In the present work, the 220 GHz FWTWT designs have been reviewed and studied. 3D Cold test simulations using both the CFDTD and FEM methods have been carried out and compared with each other as basis for 3D hot test CFDTD PIC simulations. The preliminary simulation result shows that the gain-bandwidth features at 220 GHz are achievable while carefully avoiding beam interceptions. Our study shows that the interaction characteristics are very sensitive to the operating beam parameters. Detail simulation results and discussions will be presented.
3D millimeter wave imaging of vertical cracks and its application for the inspection of HDPE pipes
NASA Astrophysics Data System (ADS)
Ghasr, Mohammad Tayeb; Ying, Kuang; Zoughi, Reza
2014-02-01
Robust detection of vertical cracks in high-density polyethylene (HDPE) pipes is a challenging task for the majority of nondestructive testing (NDT) techniques. Vertical cracks are specifically referred to those whose largest planar view is parallel to the signal direction of propagation, leaving very little signal to be scattered for detection. In such pipes this commonly occurs between two pipes sections when thermally or adhesively joined. This work presents the utility and efficacy of three-dimensional (3D) millimeter wave holographical imaging based on synthetic aperture radar (SAR) algorithm for imaging such cracks. Such a 3D millimeter wave image can readily represent the type, size, and location of various flaws within a structure. Two-dimensional (2D) slices of the 3D image, at different orientations, can also be readily produced showing the cross-sectional views of the structure and flaws, further aiding in identifying, and sizing a flaw or vertical crack. Imaging results for planner and curved (pipe section) specimen with machined flaws are presented. These images are produced using a novel field-portable, small, and low-cost wideband millimeter-wave reflectometer capable of rapid 3D image production.
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 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.
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-01-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 5mm 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 were 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. PMID:24965886
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. PMID:24965886
Measuring a hidden coordinate: Rate-exchange kinetics from 3D correlation functions
NASA Astrophysics Data System (ADS)
Berg, Mark A.; Darvin, Jason R.
2016-08-01
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
A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies.
Truchon, Jean-François; Pettitt, B Montgomery; Labute, Paul
2014-03-11
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
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
NASA Astrophysics Data System (ADS)
Gumerov, Nail A.; Karavaev, Alexey V.; Surjalal Sharma, A.; Shao, Xi; Papadopoulos, Konstantinos D.
2011-04-01
Efficient spectral and pseudospectral algorithms for simulation of linear and nonlinear 3D whistler waves in a cold electron plasma are developed. These algorithms are applied to the simulation of whistler waves generated by loop antennas and spheromak-like stationary waves of considerable amplitude. The algorithms are linearly stable and show good stability properties for computations of nonlinear waves over tens of thousands of time steps. Additional speedups by factors of 10-20 (comparing single core CPU and one GPU) are achieved by using graphics processors (GPUs), which enable efficient numerical simulation of the wave propagation on relatively high resolution meshes (tens of millions nodes) in personal computing environment. Comparisons of the numerical results with analytical solutions and experiments show good agreement. The limitations of the codes and the performance of the GPU computing are discussed.
Acrylic-acid-functionalized PolyHIPE scaffolds for use in 3D cell culture.
Hayward, Adam S; Sano, Naoko; Przyborski, Stefan A; Cameron, Neil R
2013-12-01
This study describes the development of a functional porous polymer for use as a scaffold to support 3D hepatocyte culture. A high internal phase emulsion (HIPE) is prepared containing the monomers styrene (STY), divinylbenzene (DVB), and 2-ethylhexyl acrylate (EHA) in the external oil phase and the monomer acrylic acid (Aa) in the internal aqueous phase. Upon thermal polymerization with azobisisobutyronitrile (AIBN), the resulting porous polymer (polyHIPE) is found to have an open-cell morphology and a porosity of 89%, both suitable characteristics for 3D cell scaffold applications. X-ray photo-electron spectroscopy reveals that the polyHIPE surface contained 7.5% carboxylic acid functionality, providing a useful substrate for subsequent surface modifications and bio-conjugations. Initial bio-compatibility assessments with human hepatocytes show that the acid functionality does not have any detrimental effect on cell adhesion. It is therefore believed that this material can be a useful precursor scaffold towards 3D substrates that offer tailored surface functionality for enhanced cell adhesion. PMID:24243821
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.
Internal waves patterns in the wake of a 3D body towed in a two-layer fluid
NASA Astrophysics Data System (ADS)
Lacaze, Laurent; Mercier, Matthieu; Thual, Olivier; Paci, Alexandre
2014-11-01
Stratified flows over obstacles are important features in meteorology and oceanography. The characterization of these flows is crucial in order to propose models of geophysical processes such as mixing and ocean circulation or orographic drag in the atmosphere. For some specific stratification profiles, the energy of internal waves generated by the obstacle can be trapped at a given depth, at the base of the oceanic mixing layer or at the top of the atmospheric boundary layer for instance. This scenario can be modelled by a two-layer stratified fluid for which gravity waves spread at the interface between the two layers. The work presented here focuses on a two-layer flow over a 3D obstacle, or equivalently, an obstacle towed in a fluid at rest. Experiments performed both in the large-scale flume of CNRM-GAME Toulouse (METEO-FRANCE & CNRS) and in a smaller tank apparatus, are presented with a specific attention on the measurement of the 3D wave patterns. A non-hydrostatic linear analysis is used to describe the observed wave patterns. The experiments highlight the strong influence of the Froude number on the generated waves. More specifically, we investigate the nature of the wake angle obtained from the wave pattern, and discuss a transition from Kelvin to Mach angle.
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.
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
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
Functional 3-D cardiac co-culture model using bioactive chitosan nanofiber scaffolds.
Hussain, Ali; Collins, George; Yip, Derek; Cho, Cheul H
2013-02-01
The in vitro generation of a three-dimensional (3-D) myocardial tissue-like construct employing cells, biomaterials, and biomolecules is a promising strategy in cardiac tissue regeneration, drug testing, and tissue engineering applications. Despite significant progress in this field, current cardiac tissue models are not yet able to stably maintain functional characteristics of cardiomyocytes for long-term culture and therapeutic purposes. The objective of this study was to fabricate bioactive 3-D chitosan nanofiber scaffolds using an electrospinning technique and exploring its potential for long-term cardiac function in the 3-D co-culture model. Chitosan is a natural polysaccharide biomaterial that is biocompatible, biodegradable, non-toxic, and cost effective. Electrospun chitosan was utilized to provide structural scaffolding characterized by scale and architectural resemblance to the extracellular matrix (ECM) in vivo. The chitosan fibers were coated with fibronectin via adsorption in order to enhance cellular adhesion to the fibers and migration into the interfibrous milieu. Ventricular cardiomyocytes were harvested from neonatal rats and studied in various culture conditions (i.e., mono- and co-cultures) for their viability and function. Cellular morphology and functionality were examined using immunofluorescent staining for alpha-sarcomeric actin (SM-actin) and gap junction protein, Connexin-43 (Cx43). Scanning electron microscopy (SEM) and light microscopy were used to investigate cellular morphology, spatial organization, and contractions. Calcium indicator was used to monitor calcium ion flux of beating cardiomyocytes. The results demonstrate that the chitosan nanofibers retained their cylindrical morphology in long-term cell cultures and exhibited good cellular attachment and spreading in the presence of adhesion molecule, fibronectin. Cardiomyocyte mono-cultures resulted in loss of cardiomyocyte polarity and islands of non-coherent contractions. However
NASA Astrophysics Data System (ADS)
Uhrig, Matthias P.; Kim, Jin-Yeon; Jacobs, Laurence J.
2016-02-01
This research presents a 3D numerical finite element (FE) model which, previously developed, precisely simulates non-contact, air-coupled measurements of nonlinear Rayleigh wave propagation. The commercial FE-solver ABAQUS is used to perform the simulations. First, frequency dependent pressure wave attenuation is investigated numerically to reconstruct the sound pressure distribution along the active surface of the non-contact receiver. Second, constitutive law and excitation source properties are optimized to match nonlinear ultrasonic experimental data. Finally, the FE-model data are fit with analytical solutions showing a good agreement and thus, indicating the significance of the study performed.
NASA Astrophysics Data System (ADS)
Obermann, Anne; Planès, Thomas; Hadziioannou, Céline; Campillo, Michel
2016-07-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: firstly, 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. Secondly, we compare the lapse-time behavior 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.
A New Global Model for 3-D variations in P Wave Speed in Earth's Mantle
NASA Astrophysics Data System (ADS)
Karason, H.; van der Hilst, R. D.; Li, C.
2003-12-01
In an effort to improve the resolution of mantle structure we have combined complementary data sets of short- and long period (absolute and differential) travel time residuals. Our new model is based on short period P (N\\~7.7x10**6), pP (N\\~2.3x10**5), and PKP (N\\~16x10**4) data from the catalog by Engdahl et al (BSSA, 1998), short-period PKP differential times (N\\~1600) measured by McSweeney & Creager, and long-period differential PP-P times - N\\~20,000 measured by Bolton & Masters and N\\~18,000 by Ritsema - and Pdiff-PKP (N\\~560) measured by Wysession. Inversion tests, spectral analysis, and comparison with geology indicate that the large-scale upper mantle structure is better constrained with the addition of PP-P, whereas the Pdiff and PKP data help constrain deep mantle structure (Karason & Van der Hilst, JGR, 2001). The long period data were measured by cross-correlation. We solved the system of equations using 400 iterations of the iterative algorithm LSQR For the short period (1 Hz) data we use a high frequency approximation and trace rays through a fine grid of constant slowness cells to invert for mantle structure. For low frequency Pdiff and PP data we account for sensitivity to structure away from the optical ray path with 3-D Frechet derivatives (sensitivity kernels) estimated from single forward scattering and projected onto basis functions (constant slowness blocks) used for model parameterization. With such kernels the low frequency data can constrain long wavelength heterogeneity without keeping the short period data from mapping details in densely sampled regions. In addition to finite frequency sensitivity kernels we optimized the localization by using a parameterization that adapts to spatial resolution, with small cells in regions of dense sampling and larger cells in regions where sampling is more sparse (the total number of cells was \\~ 350,000). Finally, we corrected all travel times and surface reflections for lateral variations in
Sabhachandani, P; Motwani, V; Cohen, N; Sarkar, S; Torchilin, V; Konry, T
2016-02-01
Here we describe a robust, microfluidic technique to generate and analyze 3D tumor spheroids, which resembles tumor microenvironment and can be used as a more effective preclinical drug testing and screening model. Monodisperse cell-laden alginate droplets were generated in polydimethylsiloxane (PDMS) microfluidic devices that combine T-junction droplet generation and external gelation for spheroid formation. The proposed approach has the capability to incorporate multiple cell types. For the purposes of our study, we generated spheroids with breast cancer cell lines (MCF-7 drug sensitive and resistant) and co-culture spheroids of MCF-7 together with a fibroblast cell line (HS-5). The device has the capability to house 1000 spheroids on chip for drug screening and other functional analysis. Cellular viability of spheroids in the array part of the device was maintained for two weeks by continuous perfusion of complete media into the device. The functional performance of our 3D tumor models and a dose dependent response of standard chemotherapeutic drug, doxorubicin (Dox) and standard drug combination Dox and paclitaxel (PCT) was analyzed on our chip-based platform. Altogether, our work provides a simple and novel, in vitro platform to generate, image and analyze uniform, 3D monodisperse alginate hydrogel tumors for various omic studies and therapeutic efficiency screening, an important translational step before in vivo studies. PMID:26686985
Layer-by-layer assembly of 3D tissue constructs with functionalized graphene
Shin, Su Ryon; Aghaei-Ghareh-Bolagh, Behnaz; Gao, Xiguang; Nikkhah, Mehdi; Jung, Sung Mi; Dolatshahi-Pirouz, Alireza; Kim, Sang Bok; Kim, Sun Min; Dokmeci, Mehmet R.; Tang, Xiaowu (Shirley); Khademhosseini, Ali
2014-01-01
Carbon-based nanomaterials have been considered as promising candidates to mimic certain structure and function of native extracellular matrix materials for tissue engineering. Significant progress has been made in fabricating carbon nanoparticle-incorporated cell culture substrates, but limited studies have been reported on the development of three-dimensional (3D) tissue constructs using these nanomaterials. Here, we present a novel approach to engineer 3D multi-layered constructs using layer-by-layer (LbL) assembly of cells separated with self-assembled graphene oxide (GO)-based thin films. The GO-based structures are shown to serve as cell adhesive sheets that effectively facilitate the formation of multi-layer cell constructs with interlayer connectivity. By controlling the amount of GO deposited in forming the thin films, the thickness of the multi-layer tissue constructs could be tuned with high cell viability. Specifically, this approach could be useful for creating dense and tightly connected cardiac tissues through the co-culture of cardiomyocytes and other cell types. In this work, we demonstrated the fabrication of stand-alone multi-layer cardiac tissues with strong spontaneous beating behavior and programmable pumping properties. Therefore, this LbL-based cell construct fabrication approach, utilizing GO thin films formed directly on cell surfaces, has great potential in engineering 3D tissue structures with improved organization, electrophysiological function, and mechanical integrity. PMID:25419209
Berney, Sandra; Bétrancourt, Mireille; Molinari, Gaëlle; Hoyek, Nady
2015-01-01
The emergence of dynamic visualizations of three-dimensional (3D) models in anatomy curricula may be an adequate solution for spatial difficulties encountered with traditional static learning, as they provide direct visualization of change throughout the viewpoints. However, little research has explored the interplay between learning material presentation formats, spatial abilities, and anatomical tasks. First, to understand the cognitive challenges a novice learner would be faced with when first exposed to 3D anatomical content, a six-step cognitive task analysis was developed. Following this, an experimental study was conducted to explore how presentation formats (dynamic vs. static visualizations) support learning of functional anatomy, and affect subsequent anatomical tasks derived from the cognitive task analysis. A second aim was to investigate the interplay between spatial abilities (spatial visualization and spatial relation) and presentation formats when the functional anatomy of a 3D scapula and the associated shoulder flexion movement are learned. Findings showed no main effect of the presentation formats on performances, but revealed the predictive influence of spatial visualization and spatial relation abilities on performance. However, an interesting interaction between presentation formats and spatial relation ability for a specific anatomical task was found. This result highlighted the influence of presentation formats when spatial abilities are involved as well as the differentiated influence of spatial abilities on anatomical tasks. PMID:25689057
Layer-by-layer assembly of 3D tissue constructs with functionalized graphene.
Shin, Su Ryon; Aghaei-Ghareh-Bolagh, Behnaz; Gao, Xiguang; Nikkhah, Mehdi; Jung, Sung Mi; Dolatshahi-Pirouz, Alireza; Kim, Sang Bok; Kim, Sun Min; Dokmeci, Mehmet R; Tang, Xiaowu Shirley; Khademhosseini, Ali
2014-10-22
Carbon-based nanomaterials have been considered as promising candidates to mimic certain structure and function of native extracellular matrix materials for tissue engineering. Significant progress has been made in fabricating carbon nanoparticle-incorporated cell culture substrates, but limited studies have been reported on the development of three-dimensional (3D) tissue constructs using these nanomaterials. Here, we present a novel approach to engineer 3D multi-layered constructs using layer-by-layer (LbL) assembly of cells separated with self-assembled graphene oxide (GO)-based thin films. The GO-based structures are shown to serve as cell adhesive sheets that effectively facilitate the formation of multi-layer cell constructs with interlayer connectivity. By controlling the amount of GO deposited in forming the thin films, the thickness of the multi-layer tissue constructs could be tuned with high cell viability. Specifically, this approach could be useful for creating dense and tightly connected cardiac tissues through the co-culture of cardiomyocytes and other cell types. In this work, we demonstrated the fabrication of stand-alone multi-layer cardiac tissues with strong spontaneous beating behavior and programmable pumping properties. Therefore, this LbL-based cell construct fabrication approach, utilizing GO thin films formed directly on cell surfaces, has great potential in engineering 3D tissue structures with improved organization, electrophysiological function, and mechanical integrity. PMID:25419209
RV functional imaging: 3-D echo-derived dynamic geometry and flow field simulations.
Pasipoularides, Ares D; Shu, Ming; Womack, Michael S; Shah, Ashish; Von Ramm, Olaf; Glower, Donald D
2003-01-01
We describe a novel functional imaging approach for quantitative analysis of right ventricular (RV) blood flow patterns in specific experimental animals (or humans) using real-time, three-dimensional (3-D) echocardiography (RT3D). The method is independent of the digital imaging modality used. It comprises three parts. First, a semiautomated segmentation aided by intraluminal contrast medium locates the RV endocardial surface. Second, a geometric scheme for dynamic RV chamber reconstruction applies a time interpolation procedure to the RT3D data to quantify wall geometry and motion at 400 Hz. A volumetric prism method validated the dynamic geometric reconstruction against simultaneous sonomicrometric canine measurements. Finally, the RV endocardial border motion information is used for mesh generation on a computational fluid dynamics solver to simulate development of the early RV diastolic inflow field. Boundary conditions (tessellated endocardial surface nodal velocities) for the solver are directly derived from the endocardial geometry and motion information. The new functional imaging approach may yield important kinematic information on the distribution of instantaneous velocities in the RV diastolic flow field of specific normal or diseased hearts. PMID:12388220
Wan, Alwin Ming-Doug; Inal, Sahika; Williams, Tiffany; Wang, Karin; Leleux, Pierre; Estevez, Luis; Giannelis, Emmanuel P.; Fischbach, Claudia; Malliaras, George G.; Gourdon, Delphine
2015-01-01
We report the fabrication of three dimensional (3D) macroporous scaffolds made from poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) via an ice-templating method. The scaffolds offer tunable pore size and morphology, and are electrochemically active. When a potential is applied to the scaffolds, reversible changes take place in their electrical doping state, which in turn enables precise control over the conformation of adsorbed proteins (e.g., fibronectin). Additionally, the scaffolds support the growth of mouse fibroblasts (3T3-L1) for 7 days, and are able to electrically control cell adhesion and pro-angiogenic capability. These 3D matrix-mimicking platforms offer precise control of protein conformation and major cell functions, over large volumes and long cell culture times. As such, they represent a new tool for biological research with many potential applications in bioelectronics, tissue engineering, and regenerative medicine. PMID:26413300
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.
Tavangar, Amirhossein; Tan, Bo; Venkatakrishnan, Krishnan
2013-01-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. PMID:23924310
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.
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
Radial subsampling for fast cost function computation in intensity-based 3D image registration
NASA Astrophysics Data System (ADS)
Boettger, Thomas; Wolf, Ivo; Meinzer, Hans-Peter; Celi, Juan Carlos
2007-03-01
Image registration is always a trade-off between accuracy and speed. Looking towards clinical scenarios the time for bringing two or more images into registration should be around a few seconds only. We present a new scheme for subsampling 3D-image data to allow for efficient computation of cost functions in intensity-based image registration. Starting from an arbitrary center point voxels are sampled along scan lines which do radially extend from the center point. We analyzed the characteristics of different cost functions computed on the sub-sampled data and compared them to known cost functions with respect to local optima. Results show the cost functions are smooth and give high peaks at the expected optima. Furthermore we investigated capture range of cost functions computed under the new subsampling scheme. Capture range was remarkably better for the new scheme compared to metrics using all voxels or different subsampling schemes and high registration accuracy was achieved as well. The most important result is the improvement in terms of speed making this scheme very interesting for clinical scenarios. We conclude using the new subsampling scheme intensity-based 3D image registration can be performed much faster than using other approaches while maintaining high accuracy. A variety of different extensions of the new approach is conceivable, e.g. non-regular distribution of the scan lines or not to let the scan lines start from a center point only, but from the surface of an organ model for example.
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.
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.
NASA Astrophysics Data System (ADS)
Deng, Y.; Ebert-Uphoff, I.; Chen, J.
2015-12-01
Causal discovery seeks to discover potential cause-effect relationships from observational data. Here we adopt the idea of interpreting large-scale atmospheric dynamical processes, particularly those tied to propagation of large-scale waves, as information flow around the globe, which can then be calculated using causal discovery methods. We apply a well-established causal discovery algorithm - based on constraint-based structure learning of probabilistic graphical models - toward 51 years of 6-hourly, atmospheric isobaric-level geopotential height data to construct the first-ever graphs of 3D information flow in the atmosphere. These graphs are created globally for different seasons and their connection to phase/energy propagation of atmospheric waves are investigated. Specifically, we examine the information flows 1) in the topical region that represent horizontal and vertical propagations of Kelvin and Rossby-gravity waves whose associated momentum transfer are known to play a key role in the Quasi-Biennial Oscillation (QBO), and 2) in the northern extratropics that represent propagations of planetary-scale waves whose heat/momentum fluxes are responsible for vacillations in the polar stratospheric vortex and occurrences of extreme events such as the stratospheric sudden warming. The sensitivity of the constructed graphs of 3D information flow to data resolution and pre-processing methods (e.g., spatial and temporal filtering) will be discussed.
NASA Astrophysics Data System (ADS)
Xu, Hongrui; Luo, Yinhe; Chen, Chao; Xu, Yixian
2016-06-01
Eikonal tomography based on ambient noise data is one of the most effective methods to reveal shallow earth structures. By tracking surface wave phase fronts, constructing travel time surfaces, and computing the gradients of travel time surfaces to generate phase velocity maps, eikonal tomography avoids the ray tracing and matrix construction and inversion in the traditional surface wave tomography methods. In this study, we collect continuous ambient noise data recorded by a dense seismic array in Karamay, Xinjiang to construct a 3D model of shallow structures using eikonal tomography. The seismic array consists of 35 stations with shortest interstation distance close to 1 km. 890 empirical surface wave Green's functions (EGFs) between each station pair are retrieved by cross-correlating one or two months of continuous ambient noise data. From these EGFs, surface wave travel times in the frequency range of 1.8 to 4.0 Hz are measured by a frequency-time analysis technique (FTAN). Then, eikonal tomography is adopted to construct Rayleigh wave phase velocity maps and estimate the phase velocity uncertainties. Finally, we invert the obtained phase velocity dispersion curves for 1D shear velocity profiles and then assemble these 1D profiles to construct a 3D shear velocity model. Major velocity features of our 3D model are correlated well with the known geological features. A shallow east-west velocity discontinuity is observed, which clearly reflects the lithological change between Baogutu formation (C1b) and Xibeikulasi formation (C1x) of lower Carboniferous system. Low shear velocities are observed beneath the location of porphyry copper deposit (V), possibly related to stockwork fracture and hydrothermal brecciation developed during the intrusion of deep magma in forming the deposit.
Schilling, Kurt; Janve, Vaibhav; Gao, Yurui; Stepniewska, Iwona; Landman, Bennett A; Anderson, Adam W
2016-04-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 that 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
Efficient implementation of the rank correlation merit function for 2D/3D registration.
Figl, M; Bloch, C; Gendrin, C; Weber, C; Pawiro, S A; Hummel, J; Markelj, P; Pernus, F; Bergmann, H; Birkfellner, W
2010-10-01
A growing number of clinical applications using 2D/3D registration have been presented recently. Usually, a digitally reconstructed radiograph is compared iteratively to an x-ray image of the known projection geometry until a match is achieved, thus providing six degrees of freedom of rigid motion which can be used for patient setup in image-guided radiation therapy or computer-assisted interventions. Recently, stochastic rank correlation, a merit function based on Spearman's rank correlation coefficient, was presented as a merit function especially suitable for 2D/3D registration. The advantage of this measure is its robustness against variations in image histogram content and its wide convergence range. The considerable computational expense of computing an ordered rank list is avoided here by comparing randomly chosen subsets of the DRR and reference x-ray. In this work, we show that it is possible to omit the sorting step and to compute the rank correlation coefficient of the full image content as fast as conventional merit functions. Our evaluation of a well-calibrated cadaver phantom also confirms that rank correlation-type merit functions give the most accurate results if large differences in the histogram content for the DRR and the x-ray image are present. PMID:20844334
NASA Astrophysics Data System (ADS)
Sidler, Rolf; Carcione, José M.; Holliger, Klaus
2014-02-01
We present a novel approach for the comprehensive, flexible and accurate simulation of poroelastic wave propagation in 3-D cylindrical coordinates. An important application of this method is the realistic modelling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, as of yet largely unresolved, problem in exploration geophysics. To this end, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the centre followed by the mudcake and/or casing, and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction and Fourier expansions in the vertical and azimuthal directions as well as a Runge-Kutta integration scheme for the time evolution. Trigonometric interpolation and a domain decomposition method based on the method of characteristics are used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces as well as to reduce the number of gridpoints in the innermost domain for computational efficiency. We apply this novel modelling approach to the particularly challenging scenario of near-surface borehole environments. To this end, we compare 3-D heterogeneous and corresponding rotationally invariant simulations, assess the sensitivity of Stoneley waves to formation permeability in the presence of a casing and evaluate the effects of an excavation damage zone behind a casing on sonic log recordings. Our results indicate that only first arrival times of fast modes are reasonably well described by rotationally invariant approximations of 3-D heterogenous media. We also find that Stoneley waves are indeed remarkably sensitive to the average permeability behind a perforated PVC casing, and that the presence of an excavation damage zone behind a casing tends to dominate the overall signature of recorded seismograms.
NASA Astrophysics Data System (ADS)
Pereira, J. P.; Duarte, C. A.; Jiao, X.; Guoy, D.
2009-06-01
This paper presents a study of generalized enrichment functions for 3D curved crack fronts. Two coordinate systems used in the definition of singular curved crack front enrichment functions are analyzed. In the first one, a set of Cartesian coordinate systems defined along the crack front is used. In the second case, the geometry of the crack front is approximated by a set of curvilinear coordinate systems. A description of the computation of derivatives of enrichment functions and curvilinear base vectors is presented. The coordinate systems are automatically defined using geometrical information provided by an explicit representation of the crack surface. A detailed procedure to accurately evaluate the surface normal, conormal and tangent vectors along curvilinear crack fronts in explicit crack surface representations is also presented. An accurate and robust definition of orthonormal vectors along crack fronts is crucial for the proper definition of enrichment functions. Numerical experiments illustrate the accuracy and robustness of the proposed approaches.
Skin-friction measurements in a 3-D, supersonic shock-wave/boundary-layer interaction
NASA Astrophysics Data System (ADS)
Wideman, Jeffrey Kenneth
An experimental study has been conducted in a three-dimensional, supersonic shockwave/boundary-layer interaction (3-D SW/BLI) with the intent of providing accurate experimental data for turbulence modeling and computational fluid dynamics (CFD) code validation. The experiment was performed in the High Reynolds Channel 1 (HRCI) wind tunnel at NASA Ames Research Center. The test was conducted at a Mach number of M(sub infinity) = 2.89 and at a Reynolds number of Re = 15 x 106/m. The model consisted of a sting-supported cylinder aligned with the tunnel axis and a 20 deg half-angle conical flare offset 1.27 cm from the cylinder centerline. The generated shock system was verified to be steady by schlieren visualization. The highlight of the study was the acquisition of 3-D skin-friction data by a laser interferometric skin friction (LISF) meter. Surface pressure measurements were obtained in 15 deg intervals around the cylinder and flare. Additional measurements included surface oil flow and laser light sheet illumination which were used to document the flow topology. Skin-friction measurements are proving to be a very challenging test of a CFD code predictive capability. However, at the present time there is a very limited amount of accurate skin-friction data in complex flows such as in 3-D SW/BLI. The LISF technique is advantageous as compared to other skin-friction measurement techniques for application in complex flows like the present since it is non-intrusive and is capable of performing measurements in flows with large shear and pressure gradients where the reliability of other techniques is questionable. Thus, the prevent skin-friction data will prove valuable to turbulence modeling and CFD code validation efforts.
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.
NASA Astrophysics Data System (ADS)
Bennis, A.; Ardhuin, F.; Dumas, F.; Bonneton, P.
2010-12-01
The interaction of waves with three-dimensional current structure is investigated using a two-way coupled modelling system combining MARS3D (Lazure and Dumas 2008) with WAVEWATCH III (Tolman 2008, Ardhuin et al. 2009) , a wave model (NOAA/NCEP, Tolman 2008). After a basic validation in two dimensions, the flow model MARS3D was adapted with three options that solve for the total momentum (Mellor 2003, 2008) or the quasi-Eulerian momentum (Ardhuin et al. 2008b). Adiabatic model results show that, as expected from theory (Ardhuin et al. 2008a), the total momentum fluxes parameterized by Mellor are not self-consistent and can lead to very large errors (Bennis and Ardhuin 2010). We thus use the model option to solve for the quasi-Eulerian momentum, including sources of momentum and turbulent kinetic energy (TKE). The influence of these TKE sources is investigated in the case of the NSTS experiment (Thornton and Guza, 1986). The feedback of the currents on the waves is negligible in this case. The sources of TKE from wave breaking and wave bottom friction are found to have strong influence on the bottom friction, in a way consistent with the parameterizations by Longuet-Higgins (1970) and Mellor (2002). The complete model is then applied to a real case of a large rip current on the South-West coast of France (Bruneau et al., Cont. Shelf Res. 2009). The breaking of waves on the opposed current generates a strong coupling on the rip current that partially controls the strength of the current and it three-dimensional shape.
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.
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.
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. PMID:26712581
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.
Nonlinear Evolution of a 3D Inertial Alfvén Wave and Its Implication in Particle Acceleration
NASA Astrophysics Data System (ADS)
Sharma, Prachi; Yadav, Nitin; Sharma, R. P.
2016-03-01
A simulation based on a pseudo-spectral method has been performed in order to study particle acceleration. A model for the acceleration of charged particles by field localization is developed for the low-β plasma. For this purpose, a fractional diffusion approach has been employed. The nonlinear interaction between a 3D inertial Alfvén wave and a slow magnetosonic wave has been examined, and the dynamical equations of these two waves in the presence of ponderomotive nonlinearity have been solved numerically. The nonlinear evolution of the inertial Alfvén wave in the presence of slow magnetosonic wave undergoes a filamentation instability and results in field intensity localization. The results obtained show the localization and power spectrum of inertial Alfvén wave due to nonlinear coupling. The scaling obtained after the first break point of the magnetic power spectrum has been used to calculate the formation of the thermal tail of energetic particles in the solar corona.
Rigas, Fotis; Sklavounos, Spyros
2005-05-20
Accidental blast wave generation and propagation in the surroundings poses severe threats for people and property. The prediction of overpressure maxima and its change with time at specified distances can lead to useful conclusions in quantitative risk analysis applications. In this paper, the use of a computational fluid dynamics (CFD) code CFX-5.6 on dense explosive detonation events is described. The work deals with the three-dimensional simulation of overpressure wave propagation generated by the detonation of a dense explosive within a small-scale branched tunnel. It also aids at validating the code against published experimental data as well as to study the way that the resulting shock wave propagates in a confined space configuration. Predicted overpressure histories were plotted and compared versus experimental measurements showing a reasonably good agreement. Overpressure maxima and corresponding times were found close to the measured ones confirming that CFDs may constitute a useful tool in explosion hazard assessment procedures. Moreover, it was found that blast wave propagates preserving supersonic speed along the tunnel accompanied by high overpressure levels, and indicating that space confinement favors the formation and maintenance of a shock rather than a weak pressure wave. PMID:15885402
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.
Fatigue crack detection in metallic structures with Lamb waves and 3D laser vibrometry
NASA Astrophysics Data System (ADS)
Staszewski, W. J.; Lee, B. C.; Traynor, R.
2007-03-01
The paper presents the application of ultrasonic guided waves for fatigue crack detection in metallic structures. The study involves a simple fatigue test performed to introduce a crack into an aluminium plate. Lamb waves generated by a low-profile, surface-bonded piezoceramic transducer are sensed using a tri-axis, multi-position scanning laser vibrometer. The results demonstrate the potential of laser vibrometry for simple, rapid and robust detection of fatigue cracks in metallic structures. The method could be used in quality inspection and in-service maintenance of metallic structures in aerospace, civil and mechanical engineering industries.
Non-linear 3D Born Shear Wave Tomography in Southeastern Asia
NASA Astrophysics Data System (ADS)
Cao, A.; Panning, M.; Kim, A.; Romanowicz, B.
2007-12-01
We have developed a 3D radially anisotropic shear velocity model of the upper mantle in southeastern Asia from the inversion of long period seismic multimode waveforms. Our approach is based on normal mode perturbation theory, specifically, on a recent modification of the Born approximation, which we call "N-Born", and which includes a non-linear term that allows the accurate inclusion of accumulated phase shifts which arise when the wavepath traverses a spatially extended region with a smooth velocity anomaly of constant sign. We apply the N-Born approximation in the forward modeling part and calculate linear 3D Born kernels in the inverse part. Our starting model is a 3D radially anisotropic model which we derived from a large dataset of teleseismic multimode long period waveforms in the period range 60 to 400 s, using a finite-frequency 2D approximation (NACT, Li and Romanowicz, 1995). This model covered a larger region of East Asia (longitude 30 to 150 degrees and latitude -10 to 60 degrees), while our N-Born model is restricted to a smaller subregion (longitude 75 to 150 degrees and latitude 0 to 45 degrees) for computational efficiency. In this subregion, our N-Born isotropic and anisotropic models are both parameterized at relatively short wavelengths corresponding to a spherical spline level 6 (~200km). Our N-Born model can fit waveforms as well as the NACT model, with up to ~ 83% variance reduction. While the models agree in general, the N-Born isotropic model shows a stronger fast velocity anomaly beneath the Tibetan plateau in the depth range of 150 km to 250 km, which disappears at greater depth, consistent with other studies. More importantly, the N-Born anisotropic model can recover well the downwelling structure associated with subducted slabs. Beneath the Tibet plateau, radial anisotropy shows VSH>VSV, which is indicative of horizontal rather than vertical flow and may help distinguish between end member models of the tectonics of Tibet.
Effect of fiber diameter on the assembly of functional 3D cardiac patches
NASA Astrophysics Data System (ADS)
Fleischer, Sharon; Miller, Jacob; Hurowitz, Haley; Shapira, Assaf; Dvir, Tal
2015-07-01
The cardiac ECM has a unique 3D structure responsible for tissue morphogenesis and strong contractions. It is divided into three fiber groups with specific roles and distinct dimensions; nanoscale endomysial fibers, perimysial fibers with a diameter of 1 μm, and epimysial fibers, which have a diameter of several micrometers. We report here on our work, where distinct 3D fibrous scaffolds, each of them recapitulating the dimension scales of a single fiber population in the heart matrix, were fabricated. We have assessed the mechanical properties of these scaffolds and the contribution of each fiber population to cardiomyocyte morphogenesis, tissue assembly and function. Our results show that the nanoscale fiber scaffolds were more elastic than the microscale scaffolds, however, cardiomyocytes cultured on microscale fiber scaffolds exhibited enhanced spreading and elongation, both on the single cell and on the engineered tissue levels. In addition, lower fibroblast proliferation rates were observed on these microscale topographies. Based on the collected data we have fabricated composite scaffolds containing micro and nanoscale fibers, promoting superior tissue morphogenesis without compromising tissue contraction. Cardiac tissues, engineered within these composite scaffolds exhibited superior function, including lower excitation threshold and stronger contraction forces than tissue engineered within the single-population fiber scaffolds.
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
Internal wave attractors examined using laboratory experiments and 3D numerical simulations
NASA Astrophysics Data System (ADS)
Brouzet, C.; Sibgatullin, I. N.; Scolan, H.; Ermanyuk, E. V.; Dauxois, T.
2016-04-01
In the present paper, we combine numerical and experimental approaches to study the dynamics of stable and unstable internal wave attractors. The problem is considered in a classic trapezoidal setup filled with a uniformly stratified fluid. Energy is injected into the system at global scale by the small-amplitude motion of a vertical wall. Wave motion in the test tank is measured with the help of conventional synthetic schlieren and PIV techniques. The numerical setup closely reproduces the experimental one in terms of geometry and the operational range of the Reynolds and Schmidt numbers. The spectral element method is used as a numerical tool to simulate the nonlinear dynamics of a viscous salt-stratified fluid. We show that the results of three-dimensional calculations are in excellent qualitative and quantitative agreement with the experimental data, including the spatial and temporal parameters of the secondary waves produced by triadic resonance instability. Further, we explore experimentally and numerically the effect of lateral walls on secondary currents and spanwise distribution of velocity amplitudes in the wave beams. Finally, we test the assumption of a bidimensional flow and estimate the error made in synthetic schlieren measurements due to this assumption.
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.
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 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
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
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
Overview Much of our understanding of the biological processes that underlie cellular functions in humans, such as cell-cell communication, intracellular signaling, and transcriptional and post-transcriptional 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 and 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. PMID:26778564
Miles, A R; Edwards, M J; Greenough, J A
2004-11-08
Perturbations on an interface driven by a strong blast wave grow in time due to a combination of Rayleigh-Taylor, Richtmyer-Meshkov, and decompression effects. In this paper, results from three-dimensional numerical simulations of such a system under drive conditions to be attainable on the National Ignition Facility [E. M. Campbell, Laser Part. Beams, 9(2), 209 (1991)] are presented. Using the multi-physics, adaptive mesh refinement, higher order Godunov Eulerian hydrocode, Raptor [L. H. Howell and J.A. Greenough, J. Comp. Phys. 184, 53 (2003)], the late nonlinear instability evolution, including transition to turbulence, is considered for various multimode perturbation spectra. The 3D post-transition state differs from the 2D result, but the process of transition proceeds similarly in both 2D and 3D. The turbulent mixing transition results in a reduction in the growth rate of the mixing layer relative to its pre-transition value and, in the case of the bubble front, relative to the 2D result. The post-transition spike front velocity is approximately the same in 2D and 3D. Implications for hydrodynamic mixing in core-collapse supernova are discussed.
Weak concentration and wave operator for a 3D coupled nonlinear Schrödinger system
NASA Astrophysics Data System (ADS)
Pastor, Ademir
2015-02-01
Reported in this paper are results concerning the Cauchy problem and the dynamics for a cubic nonlinear Schrödinger system arising in nonlinear optics. A sharp criterion is given concerned with the dichotomy global existence versus finite time blow-up. When a radial solution blows up in finite time, we prove the concentration in the critical Lebesgue space. Sufficient condition for the scattering and the construction of the wave operator in the energy space is also provided.
Generation of field-aligned currents and Alfven waves by 3D magnetic reconnection
Ma, Z.W.; Lee, L.C.; Otto, A.
1995-07-01
The authors have carried out a three-dimensional compressible MHD simulation to study the generation of field-aligned currents (FAC`s) and Alfven waves by magnetic reconnection for locally antiparallel magnetic fields across the current sheet. Reconnection is triggered by a localized resistivity. The results indicate that both FAC`s and Alfven waves are generated by the three-dimensional reconnection process. Two pairs of FAC`s are generated on each side of current sheet. The polarities of the resulting FAC pair in the leading bulge region are opposite to those of a FAC pair in the trailing quasi-steady region. It is further found that a large portion of the FAC`s ({approximately}40%) is located in the closed field line region. They examine the Walen relation between FAC and parallel vorticity and find that Alfven waves are generated and propagate away from the reconnection site. They discuss the relevance of the results to the observed Region 1 FAC`s at noon. 15 refs., 4 figs.
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.
NASA Astrophysics Data System (ADS)
Xu, Feixiang; Zou, Qiushun; Zhou, Quancheng; Wang, Tongbiao; Yu, Tianbao; Liu, Nianhua
2016-05-01
We report that self-imaging effect still can be achieved in photonic quasicrystal waveguides (PtQCWs) just as it does in photonic crystal waveguides. As a possible application of the results, a new kind of compact 3 dB PtQCWs-based power splitters based on this effect for terahertz waves with symmetric interference is presented and analyzed. The finite element method is used to calculate the distributions of stable-state electric field and evaluate transmission efficiency of these structures. The calculated results show that the proposed device provides a new compact model for exporting efficiently THz wave with a broad bandwidth to two channels averagely and can be extended to new designs of PtQCW devices.
Kokeyama, Keiko; Sato, Shuichi; Nishizawa, Atsushi; Kawamura, Seiji; Chen, Yanbei; Sugamoto, Akio
2009-10-23
The displacement- and frequency-noise-free interferometer (DFI) is a multiple laser interferometer array for gravitational-wave detection free from both the displacement noise of optics and laser frequency noise. So far, partial experimental demonstrations of the DFI have been done in 2D table top experiments. In this Letter, we report the complete demonstration of a 3D DFI. The DFI consists of four Mach-Zehnder interferometers with four mirrors and two beam splitters The attained maximum suppression of the displacement noise of both mirrors and beam splitters was 40 dB at about 50 MHz. The nonvanishing DFI response to a gravitational wave was successfully confirmed using multiple electro-optic modulators and computing methods. PMID:19905742
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)
Huang, Qinghua; Li, Zhanhui; Wang, Yanbin
2010-12-01
We presented a parallel 3-D staggered grid pseudospectral time domain (PSTD) method for simulating ground-penetrating radar (GPR) wave propagation. We took the staggered grid method to weaken the global effect in PSTD and developed a modified fast Fourier transform (FFT) spatial derivative operator to eliminate the wraparound effect due to the implicit periodical boundary condition in FFT operator. After the above improvements, we achieved the parallel PSTD computation based on an overlap domain decomposition method without any absorbing condition for each subdomain, which can significantly reduce the required grids in each overlap subdomain comparing with other proposed algorithms. We test our parallel technique for some numerical models and obtained consistent results with the analytical ones and/or those of the nonparallel PSTD method. The above numerical tests showed that our parallel PSTD algorithm is effective in simulating 3-D GPR wave propagation, with merits of saving computation time, as well as more flexibility in dealing with complicated models without losing the accuracy. The application of our parallel PSTD method in applied geophysics and paleoseismology based on GPR data confirmed the efficiency of our algorithm and its potential applications in various subdisciplines of solid earth geophysics. This study would also provide a useful parallel PSTD approach to the simulation of other geophysical problems on distributed memory PC cluster.
NASA Astrophysics Data System (ADS)
Luo, Cong; Friederich, Wolfgang
2016-04-01
Realistic shallow seismic wave propagation simulation is an important tool for studying induced seismicity (e.g., during geothermal energy development). However over a long time, there is a significant problem which constrains computational seismologists from performing a successful simulation conveniently: pre-processing. Conventional pre-processing has often turned out to be inefficient and unrobust because of the miscellaneous operations, considerable complexity and insufficiency of available tools. An integrated web-based platform for shallow seismic wave propagation simulation has been built. It is aiming at providing a user-friendly pre-processing solution, and cloud-based simulation abilities. The main features of the platform for the user include: revised digital elevation model (DEM) retrieving and processing mechanism; generation of multi-layered 3D shallow Earth model geometry (the computational domain) with user specified surface topography based on the DEM; visualization of the geometry before the simulation; a pipeline from geometry to fully customizable hexahedral element mesh generation; customization and running the simulation on our HPC; post-processing and retrieval of the results over cloud. Regarding the computational aspect, currently the widely accepted specfem3D is chosen as the computational package; packages using different types of elements can be integrated as well in the future. According to our trial simulation experiments, this web-based platform has produced accurate waveforms while significantly simplifying and enhancing the pre-processing and improving the simulation success rate.
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
P-wave tomography for 3-D radial and azimuthal anisotropy of Tohoku and Kyushu subduction zones
NASA Astrophysics Data System (ADS)
Wang, Jian; Zhao, Dapeng
2013-06-01
We determined high-resolution P-wave tomography for 3-D radial and azimuthal anisotropy of the Tohoku and Kyushu subduction zones using a large number of high-quality arrival-time data of local earthquakes recorded by the dense seismic network on the Japan Islands. Trench-normal P-wave fast-velocity directions (FVDs) are revealed in the backarc mantle wedge in both Tohoku and Kyushu, which are consistent with the model of slab-driven corner flow. Trench-parallel FVDs with amplitude <4 per cent appear in the forearc mantle wedge under Tohoku and Kyushu, suggesting the existence of B-type olivine fabric there. Trench-parallel FVDs are also visible in the mantle wedge under the volcanic front in Tohoku but not in Kyushu, suggesting that 3-D flow may exist in the mantle wedge under Tohoku and the 3-D flow is affected by the subduction rate of the oceanic plate. Negative radial anisotropy (i.e. vertical velocity being faster than horizontal velocity) is revealed in the low-velocity zones in the mantle wedge under the arc volcanoes in Tohoku and Kyushu as well as in the low-velocity zones below the Philippine Sea slab under Kyushu, which may reflect hot upwelling flows and transitions of olivine fabrics with the presence of water in the upper mantle. Trench-parallel FVDs and positive radial anisotropy (i.e. horizontal velocity being faster than vertical velocity) are revealed in the subducting Pacific slab under Tohoku and the Philippine Sea slab under Kyushu, suggesting that the slabs keep their frozen-in anisotropy formed at the mid-ocean ridge or that the slab anisotropy is induced by the lattice-preferred orientation of the B-type olivine.
NASA Astrophysics Data System (ADS)
Todoriki, M.; Furumura, T.; Maeda, T.
2013-12-01
We have studied the effect of topography and a seawater layer on the propagation of seismic wave propagation towards the realization of a high-resolution 3D FDM simulation of strong ground motions expected from future large subduction zone earthquakes along the Nankai Trough. Although most of the former studies on seismic wave propagation simulation did not consider a seawater layer in their simulation model, some of the recent studies claimed the importance of topography and a seawater layer on the simulation of strong ground motions (e.g., Petukhin et al., 2010; Nakamura, 2012; Maeda et al., 2013). In this study, we examined the effect of these two features on seismic wave propagation by introducing the high-resolution topography with a seawater layer over a wide frequency band. The area of 3D FDM simulation is 1200 km x 1000 km for horizontal directions and 200 km in depth, which covers entirely the area of southwestern Japan centered at 136E and 34.8N. This model was discretized with small grid interval of 0.5 km in horizontal direction and 0.25 km in depth. We used 2400 nodes of the K-computer, which is about 2.9% of its total resources, with a total memory of 1TB. We used a 3D velocity model of Koketsu et al. (2008) and an original source-rupture model from a recent study on the expansion of source-rupture area of the 1707 Hoei earthquake (Furumura et al., 2011). The result of simulation shows that the effect of a seawater layer on ground motion is small in almost all parts of Japan Island with a change of the seismic wave amplitude of less than +-20%. However, around the Northern Kanto area characterized by a belt-shaped anomalous zone, the amplitude of ground motion grows twice as large as that without seawater. This was possibly brought about from amplification of the amplitudes of surface waves generated on the Philippine Sea plate in the Suruga Trough located in the eastern end of the Nankai Trough. It is quite likely that the amplitude of surface wave
A full 3D plane-wave-expansion model for 1-3 piezoelectric composite structures.
Wilm, Mikaël; Ballandras, Sylvain; Laude, Vincent; Pastureaud, Thomas
2002-09-01
The plane-wave-expansion (PWE) approach dedicated to the simulation of periodic devices has been extended to 1-3 connectivity piezoelectric composite structures. The case of simple but actual piezoelectric composite structures is addressed, taking piezoelectricity, acoustic losses, and electrical excitation conditions rigorously into account. The material distribution is represented by using a bidimensional Fourier series and the electromechanical response is simulated using a Bloch-Floquet expansion together with the Fahmy-Adler formulation of the Christoffel problem. Application of the model to 1-3 connectivity piezoelectric composites is reported and compared to previously published analyses of this problem. PMID:12243182
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.
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
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. PMID:25615301
Density functional study of hydrazine N-N bond cleaving on 3d metal surfaces
NASA Astrophysics Data System (ADS)
Fathurrahman, Fadjar; Kasai, Hideaki
2015-11-01
Theoretical calculations based on dispersion-corrected density functional theory (DFT-D2) has been performed to investigate hydrazine adsorption and N-N bond cleaving on closed packed surfaces of 3d metals: Fe(110), Co(0001), Ni(111), Cu(111), and Zn(0001). The activation energies of N-N bond cleaving of hydrazine on each surface are estimated using climbing-image nudged elastic band (CINEB) method. The results showed that the activation energies for this process have increasing trend from Fe(110) to Zn(0001). By examining the electronic structure of the adsorbed hydrazine, it is found that this trend is related to occupation of derived 6σ* orbitals. It is also found that approximate linear relationship between reaction energy and activation energy (the Brønsted-Evans-Polanyi (BEP) relationship) holds for those surfaces.
NASA Astrophysics Data System (ADS)
Salin, B. M.; Salin, M. B.
2015-07-01
Although optical tools for measuring the surface-wave characteristics provide the best spatial and temporal resolutions compared with other methods, they face some difficulties while converting the results of indirect measurements into the absolute levels of waves. In this paper, we propose a combined optical method for measuring the 3D spectral density of the heights and the time realizations of the surface-wave profiles. The method involves, first, synchronous recording of the optical-brightness field on a rough-surface area and the surface-oscillation measurement at one or several points and, second, filtering of the spatial image spectrum, so that the filter parameters are also chosen from the condition of maximum correlation of the reconstructed and measured surface oscillations at one or two points. The second part of this work deals with the results of measuring the multi-dimensional wave spectra on the basis of the proposed method.
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.
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.
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
An assessment of the evidence from ATLAS3D for a variable initial mass function
NASA Astrophysics Data System (ADS)
Clauwens, Bart; Schaye, Joop; Franx, Marijn
2015-06-01
The ATLAS3D Survey has reported evidence for a non-universal stellar initial mass function (IMF) for early-type galaxies (ETGs). The IMF was constrained by comparing stellar mass measurements from kinematic data with those from spectral energy distribution (SED) fitting. Here, we investigate possible effects of scatter in the reported stellar mass measurements and their potential impact on the IMF determination. We find that a trend of the IMF mismatch parameter with the kinematic mass-to-light ratio, comparable to the trend observed by Cappellari et al., could arise if the Gaussian errors of the kinematic mass determination are typically 30 per cent. Without additional data, it is hard to separate between the option that the IMF has a true large intrinsic variation or the option that the errors in the determination are larger than anticipated. A correlation of the IMF with other properties would help to make this distinction, but no strong correlation has been found yet. The strongest correlation is with velocity dispersion. However, it has a large scatter and the correlation depends on sample selection and distance measurements. The correlation with velocity dispersion could be partly caused by the colour-dependent calibration of the surface brightness fluctuation distances of Tonry et al. We find that the K-band luminosity-limited ATLAS3D Survey is incomplete for the highest M/L galaxies below 1010.3 M⊙. There is a significant IMF-velocity dispersion trend for galaxies with SED masses above this limit, but no trend for galaxies with kinematic masses above this limit. We also find an IMF trend with distance, but no correlation between nearest neighbour ETGs, which excludes a large environmental dependence. Our findings do not rule out the reported IMF variations, but they suggest that further study is needed.
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
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.
Properties of resonance wave functions.
NASA Technical Reports Server (NTRS)
More, R. M.; Gerjuoy, E.
1973-01-01
Construction and study of resonance wave functions corresponding to poles of the Green's function for several illustrative models of theoretical interest. Resonance wave functions obtained from the Siegert and Kapur-Peierls definitions of the resonance energies are compared. The comparison especially clarifies the meaning of the normalization constant of the resonance wave functions. It is shown that the wave functions may be considered renormalized in a sense analogous to that of quantum field theory. However, this renormalization is entirely automatic, and the theory has neither ad hoc procedures nor infinite quantities.
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.
Guided-wave-based damage detection in a composite T-joint using 3D scanning laser Doppler vibrometer
NASA Astrophysics Data System (ADS)
Kolappan Geetha, Ganesh; Roy Mahapatra, D.; Srinivasan, Gopalakrishnan
2012-04-01
Composite T-joints are commonly used in modern composite airframe, pressure vessels and piping structures, mainly to increase the bending strength of the joint and prevents buckling of plates and shells, and in multi-cell thin-walled structures. Here we report a detailed study on the propagation of guided ultrasonic wave modes in a composite T-joint and their interactions with delamination in the co-cured co-bonded flange. A well designed guiding path is employed wherein the waves undergo a two step mode conversion process, one is due to the web and joint filler on the back face of the flange and the other is due to the delamination edges close to underneath the accessible surface of the flange. A 3D Laser Doppler Vibrometer is used to obtain the three components of surface displacements/velocities of the accessible face of the flange of the T-joint. The waves are launched by a piezo ceramic wafer bonded on to the back surface of the flange. What is novel in the proposed method is that the location of any change in material/geometric properties can be traced by computing a frequency domain power flow along a scan line. The scan line can be chosen over a grid either during scan or during post-processing of the scan data off-line. The proposed technique eliminates the necessity of baseline data and disassembly of structure for structural interrogation.
Low-power 20-meter 3D ranging SPAD camera based on continuous-wave indirect time-of-flight
NASA Astrophysics Data System (ADS)
Bellisai, S.; Ferretti, L.; Villa, F.; Ruggeri, A.; Tisa, S.; Tosi, A.; Zappa, F.
2012-06-01
Three dimensional (3D) image acquisitions is the enabling technology of a great number of applications; culture heritage morphology study, industrial robotics, automotive active safety and security access control are example of applications. The most important feature is the high frame-rate, to detect very fast events within the acquired scenes. In order to reduce the computational complexity, Time-of-Flight algorithms for single sensor cameras are used. To achieve high-frame rate and high distance measurement accuracy it is important to collect the most part of the reflected light using sensor with very high sensitivity, allowing the implementation of a low-power light source. We designed and developed a single-photon detection based 3D ranging camera, capable to acquire distance image up to 22.5 m, with a resolution down to one centimeter. The light source used in this prototype employs 8 laser diodes sinusoidally modulated. The imager used in the application is based on Single-Photon Avalanche Diodes (SPADs) fabricated in a standard CMOS 0.35 μm technology. The sensor has 1024 pixels arranged in a 32x32 squared layout, with overall dimensions of 3.5mm x 3.5mm. The camera acquires 3D images through the continuous-wave indirect Time of Flight (cw-iTOF) technique. The typical frame-rate is 20 fps while the theoretical maximum frame-rate is 5 kfps. The precision is better than 5 cm within 22.5 m range, and can be effectively used in indoor applications, e.g. in industrial environment.
3D Printing Multi-Functionality: Embedded RF Antennas and Components
NASA Technical Reports Server (NTRS)
Shemelya, C. M.; Zemba, M.; Liang, M.; Espalin, D.; Kief, C.; Xin, H.; Wicker, R. B.; MacDonald, E. W.
2015-01-01
Significant research and press has recently focused on the fabrication freedom of Additive Manufacturing (AM) to create both conceptual models and final end-use products. This flexibility allows design modifications to be immediately reflected in 3D printed structures, creating new paradigms within the manufacturing process. 3D printed products will inevitably be fabricated locally, with unit-level customization, optimized to unique mission requirements. However, for the technology to be universally adopted, the processes must be enhanced to incorporate additional technologies; such as electronics, actuation, and electromagnetics. Recently, a novel 3D printing platform, Multi3D manufacturing, was funded by the presidential initiative for revitalizing manufacturing in the USA using 3D printing (America Makes - also known as the National Additive Manufacturing Innovation Institute). The Multi3D system specifically targets 3D printed electronics in arbitrary form; and building upon the potential of this system, this paper describes RF antennas and components fabricated through the integration of material extrusion 3D printing with embedded wire, mesh, and RF elements.
van de Kamp, Thomas; dos Santos Rolo, Tomy; Vagovič, Patrik; Baumbach, Tilo; Riedel, Alexander
2014-01-01
Digital surface mesh models based on segmented datasets have become an integral part of studies on animal anatomy and functional morphology; usually, they are published as static images, movies or as interactive PDF files. We demonstrate the use of animated 3D models embedded in PDF documents, which combine the advantages of both movie and interactivity, based on the example of preserved Trigonopterus weevils. The method is particularly suitable to simulate joints with largely deterministic movements due to precise form closure. We illustrate the function of an individual screw-and-nut type hip joint and proceed to the complex movements of the entire insect attaining a defence position. This posture is achieved by a specific cascade of movements: Head and legs interlock mutually and with specific features of thorax and the first abdominal ventrite, presumably to increase the mechanical stability of the beetle and to maintain the defence position with minimal muscle activity. The deterministic interaction of accurately fitting body parts follows a defined sequence, which resembles a piece of engineering. PMID:25029366
3D printing enables separation of orthogonal functions within a hydrogel particle.
Raman, Ritu; Clay, Nicholas E; Sen, Sanjeet; Melhem, Molly; Qin, Ellen; Kong, Hyunjoon; Bashir, Rashid
2016-06-01
Multifunctional particles with distinct physiochemical phases are required by a variety of applications in biomedical engineering, such as diagnostic imaging and targeted drug delivery. This motivates the development of a repeatable, efficient, and customizable approach to manufacturing particles with spatially segregated bioactive moieties. This study demonstrates a stereolithographic 3D printing approach for designing and fabricating large arrays of biphasic poly (ethylene glycol) diacrylate (PEGDA) gel particles. The fabrication parameters governing the physical and biochemical properties of multi-layered particles are thoroughly investigated, yielding a readily tunable approach to manufacturing customizable arrays of multifunctional particles. The advantage in spatially organizing functional epitopes is examined by loading superparamagnetic iron oxide nanoparticles (SPIONs) and bovine serum albumin (BSA) in separate layers of biphasic PEGDA gel particles and examining SPION-induced magnetic resonance (MR) contrast and BSA-release kinetics. Particles with spatial segregation of functional moieties have demonstrably higher MR contrast and BSA release. Overall, this study will contribute significant knowledge to the preparation of multifunctional particles for use as biomedical tools. PMID:27215416
van de Kamp, Thomas; dos Santos Rolo, Tomy; Vagovič, Patrik; Baumbach, Tilo; Riedel, Alexander
2014-01-01
Digital surface mesh models based on segmented datasets have become an integral part of studies on animal anatomy and functional morphology; usually, they are published as static images, movies or as interactive PDF files. We demonstrate the use of animated 3D models embedded in PDF documents, which combine the advantages of both movie and interactivity, based on the example of preserved Trigonopterus weevils. The method is particularly suitable to simulate joints with largely deterministic movements due to precise form closure. We illustrate the function of an individual screw-and-nut type hip joint and proceed to the complex movements of the entire insect attaining a defence position. This posture is achieved by a specific cascade of movements: Head and legs interlock mutually and with specific features of thorax and the first abdominal ventrite, presumably to increase the mechanical stability of the beetle and to maintain the defence position with minimal muscle activity. The deterministic interaction of accurately fitting body parts follows a defined sequence, which resembles a piece of engineering. PMID:25029366
A MATLAB function for 3-D and 4-D topographical visualization in geosciences
NASA Astrophysics Data System (ADS)
Zekollari, Harry
2016-04-01
Combining topographical information and spatially varying variables in visualizations is often crucial and inherent to geoscientific problems. Despite this, it is often an impossible or a very time-consuming and difficult task to create such figures by using classic software packages. This is also the case in the widely used numerical computing environment MATLAB. Here a MATLAB function is introduced for plotting a variety of natural environments with a pronounced topography, such as for instance glaciers, volcanoes and lakes in mountainous regions. Landscapes can be visualized in 3-D, with a single colour defining a featured surface type (e.g. ice, snow, water, lava), or with a colour scale defining the magnitude of a variable (e.g. ice thickness, snow depth, water depth, surface velocity, gradient, elevation). As an input only the elevation of the subsurface (typically the bedrock) and the surface are needed, which can be complemented by various input parameters in order to adapt the figure to specific needs. The figures are particularly suited to make time-evolving animations of natural processes, such as for instance a glacier retreat or a lake drainage event. Several visualization examples will be provided alongside with animations. The function, which is freely available for download, only requires the basic package of MATLAB and can be run on any standard stationary or portable personal computer.
NASA Astrophysics Data System (ADS)
Wang, S.; De Hoop, M. V.; Xia, J.; Li, X.
2011-12-01
We consider the modeling of elastic seismic wave propagation on a rectangular domain via the discretization and solution of the inhomogeneous coupled Helmholtz equation in 3D, by exploiting a parallel multifrontal sparse direct solver equipped with Hierarchically Semi-Separable (HSS) structure to reduce the computational complexity and storage. In particular, we are concerned with solving this equation on a large domain, for a large number of different forcing terms in the context of seismic problems in general, and modeling in particular. We resort to a parsimonious mixed grid finite differences scheme for discretizing the Helmholtz operator and Perfect Matched Layer boundaries, resulting in a non-Hermitian matrix. We make use of a nested dissection based domain decomposition, and introduce an approximate direct solver by developing a parallel HSS matrix compression, factorization, and solution approach. We cast our massive parallelization in the framework of the multifrontal method. The assembly tree is partitioned into local trees and a global tree. The local trees are eliminated independently in each processor, while the global tree is eliminated through massive communication. The solver for the inhomogeneous equation is a parallel hybrid between multifrontal and HSS structure. The computational complexity associated with the factorization is almost linear with the size of the Helmholtz matrix. Our numerical approach can be compared with the spectral element method in 3D seismic applications.
Hartzell, S.; Harmsen, S.; Williams, R.A.; Carver, D.; Frankel, A.; Choy, G.; Liu, P.-C.; Jachens, R.C.; Brocher, T.M.; Wentworth, C.M.
2006-01-01
A 3D seismic velocity and attenuation model is developed for Santa Clara Valley, California, and its surrounding uplands to predict ground motions from scenario earthquakes. The model is developed using a variety of geologic and geophysical data. Our starting point is a 3D geologic model developed primarily from geologic mapping and gravity and magnetic surveys. An initial velocity model is constructed by using seismic velocities from boreholes, reflection/refraction lines, and spatial autocorrelation microtremor surveys. This model is further refined and the seismic attenuation is estimated through waveform modeling of weak motions from small local events and strong-ground motion from the 1989 Loma Prieta earthquake. Waveforms are calculated to an upper frequency of 1 Hz using a parallelized finite-difference code that utilizes two regions with a factor of 3 difference in grid spacing to reduce memory requirements. Cenozoic basins trap and strongly amplify ground motions. This effect is particularly strong in the Evergreen Basin on the northeastern side of the Santa Clara Valley, where the steeply dipping Silver Creek fault forms the southwestern boundary of the basin. In comparison, the Cupertino Basin on the southwestern side of the valley has a more moderate response, which is attributed to a greater age and velocity of the Cenozoic fill. Surface waves play a major role in the ground motion of sedimentary basins, and they are seen to strongly develop along the western margins of the Santa Clara Valley for our simulation of the Loma Prieta earthquake.
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.
NASA Astrophysics Data System (ADS)
Trovato, Claudio; Aochi, Hideo; De Martin, Florent
2014-05-01
Understanding the source mechanism of long-period (LP) seismic signals on volcanoes is an important key point in volcanology and for the hazard forecasting. In the last decades, moment tensor inversions have led to various descriptions of the kinematic source mechanism. These inversions suppose a relatively simple structure of the medium. However, the seismic wave propagation in a realistic 3-D volcano model should be taken into account for understanding the complicated physical processes of magma and gas behaviors at depth. We are studying Etna volcano, Italy, to understand the volcanic processes during different stages of activity. We adopt a spectral element method (SEM), a code EFISPEC3D (De Martin, BSSA, 2011), which shows a good accuracy and numerical stability in the simulations of seismic wave propagation. First we construct the geometrical model. We use a digital elevation model (DEM) to generate finite element meshes with a spacing of 50 m on the ground surface. We aim to calculate the ground motions until 3 Hz for the shallowest layer with Vs = ~500 m/s. The minimal size of the hexahedral elements is required to be around 100 m, with a total number of elements n = ~2 10 ^ 6 for the whole model. We compare different velocity structure configurations. We start with a homogeneous medium and add complexities taking in account the shallow low velocity structure. We also introduce a velocity gradient towards depth. Simulations performed in the homogeneous medium turn in approximately 20 hours for calculations parallelized on 16 CPUs. Complex velocity models should take approximately the same time of computation. We then try to simulate the ground motion from the LP sources (0.1-1.5 Hz) obtained by the inversion for the Etna volcano in 2008 (De Barros, GRL, 2009 and De Barros, JGR, 2011). Some vertical and horizontal structures can be added to reproduce injected dikes or sills respectively.
Function of pretribosphenic and tribosphenic mammalian molars inferred from 3D animation.
Schultz, Julia A; Martin, Thomas
2014-10-01
Appearance of the tribosphenic molar in the Late Jurassic (160 Ma) is a crucial innovation for food processing in mammalian evolution. This molar type is characterized by a protocone, a talonid basin and a two-phased chewing cycle, all of which are apomorphic. In this functional study on the teeth of Late Jurassic Dryolestes leiriensis and the living marsupial Monodelphis domestica, we demonstrate that pretribosphenic and tribosphenic molars show fundamental differences of food reduction strategies, representing a shift in dental function during the transition of tribosphenic mammals. By using the Occlusal Fingerprint Analyser (OFA), we simulated the chewing motions of the pretribosphenic Dryolestes that represents an evolutionary precursor condition to such tribosphenic mammals as Monodelphis. Animation of chewing path and detection of collisional contacts between virtual models of teeth suggests that Dryolestes differs from the classical two-phased chewing movement of tribosphenidans, due to the narrowing of the interdental space in cervical (crown-root transition) direction, the inclination angle of the hypoflexid groove, and the unicuspid talonid. The pretribosphenic chewing cycle is equivalent to phase I of the tribosphenic chewing cycle, but the former lacks phase II of the tribosphenic chewing. The new approach can analyze the chewing cycle of the jaw by using polygonal 3D models of tooth surfaces, in a way that is complementary to the electromyography and strain gauge studies of muscle function of living animals. The technique allows alignment and scaling of isolated fossil teeth and utilizes the wear facet orientation and striation of the teeth to reconstruct the chewing path of extinct mammals. PMID:25091547
Function of pretribosphenic and tribosphenic mammalian molars inferred from 3D animation
NASA Astrophysics Data System (ADS)
Schultz, Julia A.; Martin, Thomas
2014-10-01
Appearance of the tribosphenic molar in the Late Jurassic (160 Ma) is a crucial innovation for food processing in mammalian evolution. This molar type is characterized by a protocone, a talonid basin and a two-phased chewing cycle, all of which are apomorphic. In this functional study on the teeth of Late Jurassic Dryolestes leiriensis and the living marsupial Monodelphis domestica, we demonstrate that pretribosphenic and tribosphenic molars show fundamental differences of food reduction strategies, representing a shift in dental function during the transition of tribosphenic mammals. By using the Occlusal Fingerprint Analyser (OFA), we simulated the chewing motions of the pretribosphenic Dryolestes that represents an evolutionary precursor condition to such tribosphenic mammals as Monodelphis. Animation of chewing path and detection of collisional contacts between virtual models of teeth suggests that Dryolestes differs from the classical two-phased chewing movement of tribosphenidans, due to the narrowing of the interdental space in cervical (crown-root transition) direction, the inclination angle of the hypoflexid groove, and the unicuspid talonid. The pretribosphenic chewing cycle is equivalent to phase I of the tribosphenic chewing cycle, but the former lacks phase II of the tribosphenic chewing. The new approach can analyze the chewing cycle of the jaw by using polygonal 3D models of tooth surfaces, in a way that is complementary to the electromyography and strain gauge studies of muscle function of living animals. The technique allows alignment and scaling of isolated fossil teeth and utilizes the wear facet orientation and striation of the teeth to reconstruct the chewing path of extinct mammals.
3D-Ultrastructure, Functions and Stress Responses of Gastropod (Biomphalaria glabrata) Rhogocytes
Kokkinopoulou, Maria; Güler, M. Alptekin; Lieb, Bernhard; Barbeck, Mike; Ghanaati, Shahram; Markl, Jürgen
2014-01-01
Rhogocytes are pore cells scattered among the connective tissue of different body parts of gastropods and other molluscs, with great variation in their number, shape and size. They are enveloped by a lamina of extracellular matrix. Their most characteristic feature is the “slit apparatus”, local invaginations of the plasma membrane bridged by cytoplasmic bars, forming slits of ca. 20 nm width. A slit diaphragm creates a molecular sieve with permeation holes of 20×20 nm. In blue-blooded gastropods, rhogocytes synthesize and secrete the respiratory protein hemocyanin, and it has been proposed–though not proven–that in the rare red-blooded snail species they might synthesize and secrete the hemoglobin. However, the cellular secretion pathway for respiratory proteins, and the functional role(s) of the enigmatic rhogocyte slit apparatus are still unclear. Additional functions for rhogocytes have been proposed, notably a role in protein uptake and degradation, and in heavy metal detoxification. Here we provide new structural and functional information on the rhogocytes of the red-blooded freshwater snail Biomphalaria glabrata. By in situ hybridization of mantle tissues, we prove that rhogocytes indeed synthesize hemoglobin. By electron tomography, the first three dimensional (3D) reconstructions of the slit apparatus are provided, showing detail of highly dense material in the cytoplasmic bars close to the slits. By immunogold labelling, we collected evidence that a major component of this material is actin. By genome databank mining, the complete sequence of a B. glabrata nephrin was obtained, and localized to the rhogocytes by immunofluorescence microscopy. The presence of both proteins fit the ultrastructure-based hypothesis that rhogocytes are related to mammalian podocytes and insect nephrocytes. Reactions of the rhogocytes to deprivation of food and cadmium toxification are also documented, and a possible secretion pathway of newly synthesized respiratory
3D-ultrastructure, functions and stress responses of gastropod (Biomphalaria glabrata) rhogocytes.
Kokkinopoulou, Maria; Güler, M Alptekin; Lieb, Bernhard; Barbeck, Mike; Ghanaati, Shahram; Markl, Jürgen
2014-01-01
Rhogocytes are pore cells scattered among the connective tissue of different body parts of gastropods and other molluscs, with great variation in their number, shape and size. They are enveloped by a lamina of extracellular matrix. Their most characteristic feature is the "slit apparatus", local invaginations of the plasma membrane bridged by cytoplasmic bars, forming slits of ca. 20 nm width. A slit diaphragm creates a molecular sieve with permeation holes of 20×20 nm. In blue-blooded gastropods, rhogocytes synthesize and secrete the respiratory protein hemocyanin, and it has been proposed-though not proven-that in the rare red-blooded snail species they might synthesize and secrete the hemoglobin. However, the cellular secretion pathway for respiratory proteins, and the functional role(s) of the enigmatic rhogocyte slit apparatus are still unclear. Additional functions for rhogocytes have been proposed, notably a role in protein uptake and degradation, and in heavy metal detoxification. Here we provide new structural and functional information on the rhogocytes of the red-blooded freshwater snail Biomphalaria glabrata. By in situ hybridization of mantle tissues, we prove that rhogocytes indeed synthesize hemoglobin. By electron tomography, the first three dimensional (3D) reconstructions of the slit apparatus are provided, showing detail of highly dense material in the cytoplasmic bars close to the slits. By immunogold labelling, we collected evidence that a major component of this material is actin. By genome databank mining, the complete sequence of a B. glabrata nephrin was obtained, and localized to the rhogocytes by immunofluorescence microscopy. The presence of both proteins fit the ultrastructure-based hypothesis that rhogocytes are related to mammalian podocytes and insect nephrocytes. Reactions of the rhogocytes to deprivation of food and cadmium toxification are also documented, and a possible secretion pathway of newly synthesized respiratory proteins
Anderson localization of matter waves in 3D anisotropic disordered potentials
NASA Astrophysics Data System (ADS)
Piraud, Marie
2016-05-01
We study quantum transport and Anderson localization of matterwaves in 3 dimensional correlated disorder, focusing on the effects of the anisotropy. Indeed, understanding the anisotropy effects is fundamental for experiments with ultracold atoms as well as for several other systems, such as electrons in MOSFETs, light in biological medium, liquid crystals. A major challenge is to understand whether the anisotropy of the diffusion tensor is altered by the interference terms at the origin of Anderson localization. In particular, its anisotropy at the mobility edge remains to be investigated. So far, all theoretical analysis have assumed - more or less implicitly - that the anisotropy of the diffusion tensor is preserved by interference effects, and have focussed on the vanishing of diffusion as a whole. In this talk, I will start by presenting the usual description of matterwave transport in disordered medium. I will then present our method to go beyond the standard self-consistent theory, which includes in particular the full anisotropic structure of the spectral function. It thus avoids the infrared divergence of the usual self-consistent theory and, most importantly, does not make any assumption on the anisotropy of the renormalized diffusion tensor when including quantum interference terms. Using a generic model of disorder with elongated correlations, we find that the diffusion tensor is strongly affected by the quantum interference terms and that the anisotropy strongly diminishes in the vicinity of the mobility edge. Our work paves the way to further investigation with speckle potentials, which are directly relevant to ultracold-atom experiments. It will permit comparison with previous predictions for the mobility edge and shed new light on ongoing experiments in the field of ultracold atoms.
High Rayleigh Number 3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Yuen (3), G. Wright, D.
2009-04-01
In the last quarter of a century many numerical methods, such as finite-differences, finite-volume, their yin-yang variants, finite-elements and pseudo-spectral methods have been used to study the problem of 3-D spherical convection. All have their respective strengths but also serious weaknesses, such as low-order and can involve high algorithmic complexity, as in triangular elements. Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. We have recently introduced a new grid/mesh-free approach, using radial basis functions ( RBFs) . It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and allows naturally node-refinement. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer( planetary surface ) and inner shell ( core-mantle boundary ). The velocity and stress components are expressed in terms of the scalar potential approach and the other remaining variable is the perturbed temperature field. We have studied the problem from the weakly nonlinear to a moderately nonlinear regime involving a Rayleigh number, about 1000 times super-critical. Both purely basal and partially internal -heating cases have been considered
High Rayleigh Number 3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Wright, G.; Yuen, D. A.
2009-04-01
In the last quarter of a century many numerical methods, such as finite-differences, finite-volume, their yin-yang variants, finite-elements and pseudo-spectral methods have been used to study the problem of 3-D spherical convection. All have their respective strengths but also serious weaknesses, such as low-order and can involve high algorithmic complexity, as in triangular elements. Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. We have recently introduced a new grid/mesh-free approach, using radial basis functions (RBFs). It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and allows naturally node-refinement. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer(planetary surface) and inner shell (core-mantle boundary). The velocity and stress components are expressed in terms of the scalar potential approach and the other remaining variable is the perturbed temperature field. We have studied the problem from the weakly onlinear to a moderately nonlinear regime involving a Rayleigh number, about 1000 times super-critical. Both purely basal and partially internal-heating cases have been considered.
(Bio)Chemical Tailoring of Biogenic 3-D Nanopatterned Templates with Energy-Relevant Functionalities
Sandhage, Kenneth H; Kroger, Nils
2014-09-08
The overall aim of this research has been to obtain fundamental understanding of (bio)chemical methodologies that will enable utilization of the unique 3-D nanopatterned architectures naturally produced by diatoms for the syntheses of advanced functional materials attractive for applications in energy harvesting/conversion and storage. This research has been conducted in three thrusts: Thrust 1 (In vivo immobilization of proteins in diatom biosilica) is directed towards elucidating the fundamental mechanism(s) underlying the cellular processes of in vivo immobilization of proteins in diatom silica. Thrust 2 (Shape-preserving reactive conversion of diatom biosilica into porous, high-surface area inorganic replicas) is aimed at understanding the fundamental mechanisms of shape preservation and nanostructural evolution associated with the reactive conversion and/or coating-based conversion of diatom biosilica templates into porous inorganic replicas. Thrust 3 (Immobilization of energy-relevant enzymes in diatom biosilica and onto diatom biosilica-derived inorganic replicas) involves use of the results from both Thrust 1 and 2 to develop strategies for in vivo and in vitro immobilization of enzymes in/on diatom biosilica and diatom biosilica-derived inorganic replicas, respectively. This Final Report describes progress achieved in all 3 of these thrusts.
3D-FE analysis of functionally graded structured dental posts.
Abu Kasim, Noor H; Madfa, Ahmed A; Hamdi, Mohd; Rahbari, Ghahnavyeh R
2011-01-01
This study aimed to compare the biomechanical behaviour of functionally graded structured posts (FGSPs) and homogenous-type posts in simulated models of a maxillary central incisor. Two models of FGSPs consisting of a multilayer xTi-yHA composite design, where zirconia and alumina was added as the first layer for models A and B respectively were compared to homogenous zirconia post (model C) and a titanium post (model D). The amount of Ti and HA in the FGSP models was varied in gradations. 3D-FEA was performed on all models and stress distributions were investigated along the dental post. In addition, interface stresses between the posts and their surrounding structures were investigated under vertical, oblique, and horizontal loadings. Strain distribution along the post-dentine interface was also investigated. The results showed that FGSPs models, A and B demonstrated better stress distribution than models C and D, indicating that dental posts with multilayered structure dissipate localized and interfacial stress and strain more efficiently than homogenous-type posts. PMID:22123011
Function and 3D Structure of the N-Glycans on Glycoproteins
Nagae, Masamichi; Yamaguchi, Yoshiki
2012-01-01
Glycosylation is one of the most common post-translational modifications in eukaryotic cells and plays important roles in many biological processes, such as the immune response and protein quality control systems. It has been notoriously difficult to study glycoproteins by X-ray crystallography since the glycan moieties usually have a heterogeneous chemical structure and conformation, and are often mobile. Nonetheless, recent technical advances in glycoprotein crystallography have accelerated the accumulation of 3D structural information. Statistical analysis of “snapshots” of glycoproteins can provide clues to understanding their structural and dynamic aspects. In this review, we provide an overview of crystallographic analyses of glycoproteins, in which electron density of the glycan moiety is clearly observed. These well-defined N-glycan structures are in most cases attributed to carbohydrate-protein and/or carbohydrate-carbohydrate interactions and may function as “molecular glue” to help stabilize inter- and intra-molecular interactions. However, the more mobile N-glycans on cell surface receptors, the electron density of which is usually missing on X-ray crystallography, seem to guide the partner ligand to its binding site and prevent irregular protein aggregation by covering oligomerization sites away from the ligand-binding site. PMID:22942711
Zhang, Yali; Xia, Lunguo; Zhai, Dong; Shi, Mengchao; Luo, Yongxiang; Feng, Chun; Fang, Bing; Yin, Jingbo; Chang, Jiang; Wu, Chengtie
2015-12-01
The hierarchical microstructure, surface and interface of biomaterials are important factors influencing their bioactivity. Porous bioceramic scaffolds have been widely used for bone tissue engineering by optimizing their chemical composition and large-pore structure. However, the surface and interface of struts in bioceramic scaffolds are often ignored. The aim of this study is to incorporate hierarchical pores and bioactive components into the bioceramic scaffolds by constructing nanopores and bioactive elements on the struts of scaffolds and further improve their bone-forming activity. Mesoporous bioactive glass (MBG) modified β-tricalcium phosphate (MBG-β-TCP) scaffolds with a hierarchical pore structure and a functional strut surface (∼100 nm of MBG nanolayer) were successfully prepared via 3D printing and spin coating. The compressive strength and apatite-mineralization ability of MBG-β-TCP scaffolds were significantly enhanced as compared to β-TCP scaffolds without the MBG nanolayer. The attachment, viability, alkaline phosphatase (ALP) activity, osteogenic gene expression (Runx2, BMP2, OPN and Col I) and protein expression (OPN, Col I, VEGF, HIF-1α) of rabbit bone marrow stromal cells (rBMSCs) as well as the attachment, viability and angiogenic gene expression (VEGF and HIF-1α) of human umbilical vein endothelial cells (HUVECs) in MBG-β-TCP scaffolds were significantly upregulated compared with conventional bioactive glass (BG)-modified β-TCP (BG-β-TCP) and pure β-TCP scaffolds. Furthermore, MBG-β-TCP scaffolds significantly enhanced the formation of new bone in vivo as compared to BG-β-TCP and β-TCP scaffolds. The results suggest that application of the MBG nanolayer to modify 3D-printed bioceramic scaffolds offers a new strategy to construct hierarchically porous scaffolds with significantly improved physicochemical and biological properties, such as mechanical properties, osteogenesis, angiogenesis and protein expression for bone tissue
NASA Astrophysics Data System (ADS)
Zhang, Yali; Xia, Lunguo; Zhai, Dong; Shi, Mengchao; Luo, Yongxiang; Feng, Chun; Fang, Bing; Yin, Jingbo; Chang, Jiang; Wu, Chengtie
2015-11-01
The hierarchical microstructure, surface and interface of biomaterials are important factors influencing their bioactivity. Porous bioceramic scaffolds have been widely used for bone tissue engineering by optimizing their chemical composition and large-pore structure. However, the surface and interface of struts in bioceramic scaffolds are often ignored. The aim of this study is to incorporate hierarchical pores and bioactive components into the bioceramic scaffolds by constructing nanopores and bioactive elements on the struts of scaffolds and further improve their bone-forming activity. Mesoporous bioactive glass (MBG) modified β-tricalcium phosphate (MBG-β-TCP) scaffolds with a hierarchical pore structure and a functional strut surface (~100 nm of MBG nanolayer) were successfully prepared via 3D printing and spin coating. The compressive strength and apatite-mineralization ability of MBG-β-TCP scaffolds were significantly enhanced as compared to β-TCP scaffolds without the MBG nanolayer. The attachment, viability, alkaline phosphatase (ALP) activity, osteogenic gene expression (Runx2, BMP2, OPN and Col I) and protein expression (OPN, Col I, VEGF, HIF-1α) of rabbit bone marrow stromal cells (rBMSCs) as well as the attachment, viability and angiogenic gene expression (VEGF and HIF-1α) of human umbilical vein endothelial cells (HUVECs) in MBG-β-TCP scaffolds were significantly upregulated compared with conventional bioactive glass (BG)-modified β-TCP (BG-β-TCP) and pure β-TCP scaffolds. Furthermore, MBG-β-TCP scaffolds significantly enhanced the formation of new bone in vivo as compared to BG-β-TCP and β-TCP scaffolds. The results suggest that application of the MBG nanolayer to modify 3D-printed bioceramic scaffolds offers a new strategy to construct hierarchically porous scaffolds with significantly improved physicochemical and biological properties, such as mechanical properties, osteogenesis, angiogenesis and protein expression for bone tissue
The performance of semilocal and hybrid density functionals in 3d transition-metal chemistry
NASA Astrophysics Data System (ADS)
Furche, Filipp; Perdew, John P.
2006-01-01
We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic π ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca2, Mn2, and Zn2 are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground
3D meshes of carbon nanotubes guide functional reconnection of segregated spinal explants.
Usmani, Sadaf; Aurand, Emily Rose; Medelin, Manuela; Fabbro, Alessandra; Scaini, Denis; Laishram, Jummi; Rosselli, Federica B; Ansuini, Alessio; Zoccolan, Davide; Scarselli, Manuela; De Crescenzi, Maurizio; Bosi, Susanna; Prato, Maurizio; Ballerini, Laura
2016-07-01
In modern neuroscience, significant progress in developing structural scaffolds integrated with the brain is provided by the increasing use of nanomaterials. We show that a multiwalled carbon nanotube self-standing framework, consisting of a three-dimensional (3D) mesh of interconnected, conductive, pure carbon nanotubes, can guide the formation of neural webs in vitro where the spontaneous regrowth of neurite bundles is molded into a dense random net. This morphology of the fiber regrowth shaped by the 3D structure supports the successful reconnection of segregated spinal cord segments. We further observed in vivo the adaptability of these 3D devices in a healthy physiological environment. Our study shows that 3D artificial scaffolds may drive local rewiring in vitro and hold great potential for the development of future in vivo interfaces. PMID:27453939
3D meshes of carbon nanotubes guide functional reconnection of segregated spinal explants
Usmani, Sadaf; Aurand, Emily Rose; Medelin, Manuela; Fabbro, Alessandra; Scaini, Denis; Laishram, Jummi; Rosselli, Federica B.; Ansuini, Alessio; Zoccolan, Davide; Scarselli, Manuela; De Crescenzi, Maurizio; Bosi, Susanna; Prato, Maurizio; Ballerini, Laura
2016-01-01
In modern neuroscience, significant progress in developing structural scaffolds integrated with the brain is provided by the increasing use of nanomaterials. We show that a multiwalled carbon nanotube self-standing framework, consisting of a three-dimensional (3D) mesh of interconnected, conductive, pure carbon nanotubes, can guide the formation of neural webs in vitro where the spontaneous regrowth of neurite bundles is molded into a dense random net. This morphology of the fiber regrowth shaped by the 3D structure supports the successful reconnection of segregated spinal cord segments. We further observed in vivo the adaptability of these 3D devices in a healthy physiological environment. Our study shows that 3D artificial scaffolds may drive local rewiring in vitro and hold great potential for the development of future in vivo interfaces. PMID:27453939
Template-free synthesis of functional 3D BN architecture for removal of dyes from water.
Liu, Dan; Lei, Weiwei; Qin, Si; Chen, Ying
2014-01-01
Three-dimensional (3D) architectures are of interest in applications in electronics, catalysis devices, sensors and adsorption materials. However, it is still a challenge to fabricate 3D BN architectures by a simple method. Here, we report the direct synthesis of 3D BN architectures by a simple thermal treatment process. A 3D BN architecture consists of an interconnected flexible network of nanosheets. The typical nitrogen adsorption/desorption results demonstrate that the specific surface area for the as-prepared samples is up to 1156 m(2) g(-1), and the total pore volume is about 1.17 cm(3) g(-1). The 3D BN architecture displays very high adsorption rates and large capacities for organic dyes in water without any other additives due to its low densities, high resistance to oxidation, good chemical inertness and high surface area. Importantly, 88% of the starting adsorption capacity is maintained after 15 cycles. These results indicate that the 3D BN architecture is potential environmental materials for water purification and treatment. PMID:24663292
Marsano, Anna; Conficconi, Chiara; Lemme, Marta; Occhetta, Paola; Gaudiello, Emanuele; Votta, Emiliano; Cerino, Giulia; Redaelli, Alberto; Rasponi, Marco
2016-02-01
In the past few years, microfluidic-based technology has developed microscale models recapitulating key physical and biological cues typical of the native myocardium. However, the application of controlled physiological uniaxial cyclic strains on a defined three-dimension cellular environment is not yet possible. Two-dimension mechanical stimulation was particularly investigated, neglecting the complex three-dimensional cell-cell and cell-matrix interactions. For this purpose, we developed a heart-on-a-chip platform, which recapitulates the physiologic mechanical environment experienced by cells in the native myocardium. The device includes an array of hanging posts to confine cell-laden gels, and a pneumatic actuation system to induce homogeneous uniaxial cyclic strains to the 3D cell constructs during culture. The device was used to generate mature and highly functional micro-engineered cardiac tissues (μECTs), from both neonatal rat and human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM), strongly suggesting the robustness of our engineered cardiac micro-niche. Our results demonstrated that the cyclic strain was effectively highly uniaxial and uniformly transferred to cells in culture. As compared to control, stimulated μECTs showed superior cardiac differentiation, as well as electrical and mechanical coupling, owing to a remarkable increase in junction complexes. Mechanical stimulation also promoted early spontaneous synchronous beating and better contractile capability in response to electric pacing. Pacing analyses of hiPSC-CM constructs upon controlled administration of isoprenaline showed further promising applications of our platform in drug discovery, delivery and toxicology fields. The proposed heart-on-a-chip device represents a relevant step forward in the field, providing a standard functional three-dimensional cardiac model to possibly predict signs of hypertrophic changes in cardiac phenotype by mechanical and biochemical co
Comparison and functional implications of the 3D architectures of viral tRNA-like structures.
Hammond, John A; Rambo, Robert P; Filbin, Megan E; Kieft, Jeffrey S
2009-02-01
RNA viruses co-opt the host cell's biological machinery, and their infection strategies often depend on specific structures in the viral genomic RNA. Examples are tRNA-like structures (TLSs), found at the 3' end of certain plant viral RNAs, which can use the cell's aminoacyl tRNA-synthetases (AARSs) to drive addition of an amino acid to the 3' end of the viral RNA. TLSs are multifunctional RNAs involved in processes such as viral replication, translation, and viral RNA stability; these functions depend on their fold. Experimental result-based structural models of TLSs have been published. In this study, we further examine these structures using a combination of biophysical and biochemical approaches to explore the three-dimensional (3D) architectures of TLSs from the turnip yellow mosaic virus (TYMV), tobacco mosaic virus (TMV), and brome mosaic virus (BMV). We find that despite similar function, these RNAs are biophysically diverse: the TYMV TLS adopts a characteristic tRNA-like L shape, the BMV TLS has a large compact globular domain with several helical extensions, and the TMV TLS aggregates in solution. Both the TYMV and BMV TLS RNAs adopt structures with tight backbone packing and also with dynamic structural elements, suggesting complexities and subtleties that cannot be explained by simple tRNA mimicry. These results confirm some aspects of existing models and also indicate how these models can be improved. The biophysical characteristics of these TLSs show how these multifunctional RNAs might regulate various viral processes, including negative strand synthesis, and also allow comparison with other structured RNAs. PMID:19144910
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITH TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITHOUT TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
Photoelectron wave function in photoionization: plane wave or Coulomb wave?
Gozem, Samer; Gunina, Anastasia O; Ichino, Takatoshi; Osborn, David L; Stanton, John F; Krylov, Anna I
2015-11-19
The calculation of absolute total cross sections requires accurate wave functions of the photoelectron and of the initial and final states of the system. The essential information contained in the latter two can be condensed into a Dyson orbital. We employ correlated Dyson orbitals and test approximate treatments of the photoelectron wave function, that is, plane and Coulomb waves, by comparing computed and experimental photoionization and photodetachment spectra. We find that in anions, a plane wave treatment of the photoelectron provides a good description of photodetachment spectra. For photoionization of neutral atoms or molecules with one heavy atom, the photoelectron wave function must be treated as a Coulomb wave to account for the interaction of the photoelectron with the +1 charge of the ionized core. For larger molecules, the best agreement with experiment is often achieved by using a Coulomb wave with a partial (effective) charge smaller than unity. This likely derives from the fact that the effective charge at the centroid of the Dyson orbital, which serves as the origin of the spherical wave expansion, is smaller than the total charge of a polyatomic cation. The results suggest that accurate molecular photoionization cross sections can be computed with a modified central potential model that accounts for the nonspherical charge distribution of the core by adjusting the charge in the center of the expansion. PMID:26509428
Poly(acrylic acid) brushes pattern as a 3D functional biosensor surface for microchips
NASA Astrophysics Data System (ADS)
Wang, Yan-Mei; Cui, Yi; Cheng, Zhi-Qiang; Song, Lu-Sheng; Wang, Zhi-You; Han, Bao-Hang; Zhu, Jin-Song
2013-02-01
Poly(acrylic acid) (PAA) brushes, a novel three dimensional (3D) precursor layer of biosensor or protein microarrays, possess high protein loading level and low non-specific protein adsorption. In this article, we describe a simple and convenient way to fabricate 3D PAA brushes pattern by microcontact printing (μCP) and characterize it with FT-IR and optical microscopy. The carboxyl groups of PAA brushes can be applied to covalently immobilize protein for immunoassay. Thriving 3D space made by patterning PAA brushes thin film is available to enhance protein immobilization, which is confirmed by measuring model protein interaction between human immunoglobulin G (H-IgG) and goat anti-H-IgG (G-H-IgG) with fluorescence microscopy and surface plasmon resonance imaging (SPRi). As expected, the SPRi signals of H-IgG coating on 3D PAA brushes pattern and further measuring specific binding with G-H-IgG are all larger than that of 3D PAA brushes without pattern and 2D bare gold surface. We further revealed that this surface can be used for high-throughput screening and clinical diagnosis by label-free assaying of Hepatitis-B-Virus surface antibody (HBsAb) with Hepatitis-B-Virus surface antigen (HBsAg) concentration array chip. The linearity range for HBsAb assay is wider than that of conventional ELISA method.
Neuronal-glial populations form functional networks in a biocompatible 3D scaffold.
Smith, Imogen; Haag, Marcus; Ugbode, Christopher; Tams, Daniel; Rattray, Marcus; Przyborski, Stefan; Bithell, Angela; Whalley, Benjamin J
2015-11-16
Monolayers of neurons and glia have been employed for decades as tools for the study of cellular physiology and as the basis for a variety of standard toxicological assays. A variety of three dimensional (3D) culture techniques have been developed with the aim to produce cultures that recapitulate desirable features of intact. In this study, we investigated the effect of preparing primary mouse mixed neuron and glial cultures in the inert 3D scaffold, Alvetex. Using planar multielectrode arrays, we compared the spontaneous bioelectrical activity exhibited by neuroglial networks grown in the scaffold with that seen in the same cells prepared as conventional monolayer cultures. Two dimensional (monolayer; 2D) cultures exhibited a significantly higher spike firing rate than that seen in 3D cultures although no difference was seen in total signal power (<50Hz) while pharmacological responsiveness of each culture type to antagonism of GABAAR, NMDAR and AMPAR was highly comparable. Interestingly, correlation of burst events, spike firing and total signal power (<50Hz) revealed that local field potential events were associated with action potential driven bursts as was the case for 2D cultures. Moreover, glial morphology was more physiologically normal in 3D cultures. These results show that 3D culture in inert scaffolds represents a more physiologically normal preparation which has advantages for physiological, pharmacological, toxicological and drug development studies, particularly given the extensive use of such preparations in high throughput and high content systems. PMID:26493605
Cloud 3D Effects Evidenced in Landsat Power Spectra and Autocorrelation Functions
NASA Technical Reports Server (NTRS)
Oreopoulos, Lazaros; Marshak, Alexander; Cahalan, Robert F.; Wen, Guoyong
1999-01-01
the spectral signatures of decorrelation between reflectance and optical depth at large scales becoming stronger as the magnitude of cloud top variations increase. Finally, the usefulness of power spectral analysis in evaluating the skill of novel optical depth retrieval techniques in removing 3D radiative effects is demonstrated. New techniques using inverse Non-local Independent Pixel Approximation (NIPA) and Normalized Difference of Nadir Reflectivity (NDNR) yield optical depth fields which better match the scale-by-scale variability of the true optical depth field.